作者iyenn (晓风)
看板Grad-ProbAsk
标题Re: [理工] [工数]-拉普拉斯解ODE
时间Tue Sep 22 10:24:29 2009
※ 引述《abcxyz123 (阿里八八与四十大盗)》之铭言：
: 85年交大土木 喻老三本中的51页例13
: 1-t^2/π^2 0<t<π
: y"+ y = f(t) f(t)= 〈 y(0)=0 y'(0)=0
: 0 t>π
: -2cost/π^2 - t^2/π^2 + 2/π^2 0<t<π
: 解答是 y(t) = 〈
: (1-4/π^2)cost + 2sint/π^2 t>π
: ------------------------------------------------------
: 没算前看解答，看到 t>π 0的地方有解答觉得怪怪@@ 不知题目错还是我错
: 我有先解
: y"+ y = (1-t^2/π^2)[u(t)-u(t-π)]
: 两端取L-T 整理 :
: y(s) =[1/s(s^2+1)-2/(π^2)s^3(s^2+1)]
: +e^-πs[2/(π^2)s^3(s^2+1)+(2/πs(s^2+1)]
: 然後做到这我就停了...
: 因为s^2+1除下去要做部分分式 次方好高 感觉像当兵@@
: 不知道是我解错还是有更快的方法
: 谢谢~
看解答的表示方法应该是用摺积摺出来的
(s^2+1)Y=F(S)
Y=F(S)H(S)
h(t)=sint
y=f(t)*sint
t φ^2
y= ∫(1 - -----)sin(t-φ)dφ
0 π^2
t 1 t
=cos(t-φ)| - ----{[φ^2-2]cos(t-φ)+[2φ]sin(t-φ)}|
0 π^2 0
1
=1-cost - ---- {[t^2-2]+2cost}
π^2
-2cost t^2 2
= ------ - ----- + ----- +1-cost (0<t<π)
π^2 π^2 π^2
~~~~~~~~~~~~~~~~~~~~~~~~ ~~~~~~
跟你给的一样 这项应该是解答漏了.(也可能是我错-,- )
信号与系统还老师很久了,摺积不太熟
可能摺的过程去闪到腰吧XD
不过对照你的Y(S)项逆回去应该是解答有错
=>
π 1 π
y=cos(t-φ)| - ----{[φ^2-2]cos(t-φ)+[2φ]sin(t-φ)}|
0 π^2 0
1
=-2cost - -----{-[π^2-2]cost+2cost-2πsin(t)}
π^2
-4 2
={-1 + ----}cost + ----sin(t) ,t>π
π^2 π
就当我骗骗p币吧XD
--
路遥知马力,脸书见人心
---> FB Online XD
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 123.193.214.165
1F：推 youmehim:t>π的解法好新奇 09/22 10:56
2F：→ youmehim:就是非齐次项只作用t=0~π的意思嘛? 09/22 10:58
3F：→ iyenn:新奇?= =? 09/22 11:08
4F：→ iyenn:因为函数相乘只有在0~pi有值 09/22 11:09
5F：推 youmehim:..看来我折积要重念了 09/22 11:25
6F：推 abcxyz123:喔~原来是这样阿 感谢大大 09/22 22:10