作者abcxyz123 (阿里八八与四十大盗)
看板Grad-ProbAsk
标题[理工] [工数]-拉普拉斯解ODE
时间Mon Sep 21 23:45:25 2009
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除下去要做部分分式 次方好高 感觉像当兵@@
不知道是我解错还是有更快的方法
谢谢~
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 122.121.207.170
※ 编辑: abcxyz123 来自: 122.121.207.170 (09/21 23:47)
1F：→ hihaka2001:我的感觉直接对这个曲L(1-t^2/π^2)[u(t)-u(t-π)] 09/22 00:25
2F：→ abcxyz123:对阿= ="我是这样做 不过有四项没错吧? 09/22 00:31
3F：→ abcxyz123:後面算到次方数这麽高 解答感觉不多 所以觉得怪 09/22 00:31
4F：推 fonlintw0621:不要算了 放到最後算 有够难算@@ 09/22 00:33
5F：推 squallting:四楼好建议 09/22 00:46
6F：→ abcxyz123:噗.交大土木真的很喜欢出四次方以上的部分分式... 09/22 00:48
7F：→ iyenn:照解的形式,可用摺叠去试试. 09/22 01:21