作者uniqueco (宝)
看板Grad-ProbAsk
标题[理工] [工数]-pde
时间Tue Feb 9 18:42:05 2010
d^2 u d^2u
---- = 9 ---- - e^-x 0<x<4 t>0
dt^2 dx^2
BC: u(0,t)=u(4,t)=0 t≧0
du(x,0)
IC: u(x,0)=sin(πx) , --------=0 0<x<4
dt
这题我用sin展开
算到
2 2 -4 n
d^2 nπ (4-4e (-1) )*64
------bn = -9----bn - -----------------------------------------
dt^2 16 nπ*(16-(n^2)(π^2))*3(n^2)(π^2)
答案超丑的..
请问有人能算算看吗?
没算过过程那麽丑的pde..
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 114.41.75.32
1F：→ honestonly:我算过..真的那麽丑.. 02/09 18:44
2F：→ honestonly:因为 解特解的关系 囧rz.. 02/09 18:44
3F：→ honestonly:不过 好像怪怪的耶.. 02/09 18:45
4F：→ honestonly:刚刚翻了一下题库班那本里面有 周易的答案 有个地方 02/09 18:48
5F：→ honestonly:负号 没加 所以後面都差一个负号orz.. 02/09 18:48
6F：→ honestonly:我觉得原po应该有算错.. 我记得很丑 但是没这麽丑XDDD 02/09 18:48
7F：→ honestonly:原PO可以考虑把e^x展开那项 写成f(x)之类的 02/09 18:50
8F：→ honestonly:之後计算对t微分的特解的时候 对x没有影响orz 02/09 18:50
9F：→ uniqueco:谢谢 我在重算看看=.= 02/09 19:08
10F：→ CRAZYAWIND:这题我是用分态下去算的@@~ 我印象中我没算出那麽丑= = 02/09 19:25
11F：→ uniqueco:分态@@? 那是什麽呀... 02/09 19:32
12F：推 CRAZYAWIND:令u(x.t) = ψ(x.t) + φ(x) 02/09 19:34
13F：→ CRAZYAWIND:这样就可以得到跟t有关的暂态解 跟t无关的稳态解 02/09 19:34
14F：→ uniqueco:嗯嗯@@ 谢谢 02/09 19:35
15F：推 hihaka2001:请问喔 好像会做到n=1的时候 那n不为1要讨论吗 02/09 22:29
16F：→ hihaka2001:sorry是n=4 02/09 22:29
17F：推 muwalker:周易答案背错了 他讲这题的时候α正负号不太对XDD 02/11 14:22
18F：→ honestonly:原来他有讲这题..我都在放空XDDD 02/11 18:49
19F：→ honestonly:上次还被不认识的说 所以你都在发呆罗...? XDDD 02/11 18:49