作者lhz135 (PPP)
看板Grad-ProbAsk
标题[理工]请帮我看看这pde这样令对不对?
时间Thu Mar 19 23:25:22 2009
ρ为偏微
ρW ρ^2W W(x,0)=T(x) , W(0,t)=0 ,W(L,t)=0
── =a ───+C
ρt ρx^2 a,C 是常数
C
解:令W=u(x,t)+──x(x-L) W(0,t)=u(0,t)=0
2a W(L,t)=u(L,t)=0 C
W(x,0)=u(x,0) + ──x(x-L)
2a
C
u(x,0)=T(x)-──x(x-L)
2a
令u=XT
nπx nπ
得Xn=cn sin── Tn=bn exp-(──)^2t u=ΣBn X T
L L
2 L C nπx
Bn=─∫ T(x)-──x(x-L) sin── dx
L 0 2a L
请教 这样令 这样解答 对吗?? 谢谢指教~~~
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 203.67.141.162
1F：推 jsdre:不是要令暂态稳态解吗!? 03/19 23:29
2F：推 gcaa:边界条件都为0 所以可以不用另 直接复利叶sine展开即可 03/19 23:32
3F：→ gcaa:or 直接用拉式转换解 03/19 23:33
4F：推 shinyhaung:我记得是令W=u(x,t)+v(x,t)对吧 同一楼 03/19 23:38
5F：推 atlaskuma:他令了阿 只是他先把他解出来了 03/19 23:39
6F：→ shinyhaung:是喔 没看清楚XD 话说楼上正取了没? 03/19 23:42
7F：推 atlaskuma:没>"< 现在在拼中央 03/19 23:45
8F：→ atlaskuma:都只有备取~"~ 03/19 23:46
9F：→ Zrst:这样令没错，解应该也没错 不过我是用看的 错了不负责 03/19 23:49
10F：推 jsdre:18298 有没有人会阿 > < 03/19 23:52
11F：推 atlaskuma:那题你的题目应该有打错 因为(a)(b)小题是不同的矩阵 03/19 23:55
12F：→ jsdre:矩阵!? 03/20 00:00
13F：→ atlaskuma:sorry 我看错篇了XD 03/20 00:01
14F：→ lhz135:阿~~那麽到底对不对阿@@?? 03/20 00:04
15F：→ shinyhaung:就是对了吧XDDDD 03/20 00:09
16F：→ lhz135:感恩~~~ 03/20 00:10
17F：→ ohya74921:因为是均质且边界符合Sturm Lioville边界，所以系统一定 03/20 10:18
18F：→ ohya74921:由特徵状态构成，令 w=Σcn(t) sin (nπx/L) 带入原方程 03/20 10:20
19F：推 ohya74921:得到Σ[cn'(t)+ (n^2π^2a/L^2) cn(t)]sin(nπx/L) = C 03/20 10:24
20F：→ ohya74921:等号两边相等有 c0'(t)=C & cn'(t)+n^2π^2a/L^2cn(t)=0 03/20 10:26
21F：→ ohya74921:sorry 没有C0(t),因为特徵含数不包括n=0 03/20 10:34
22F：→ ohya74921:所以是 cn'(t)+n^2π^2a/L^2cn(t) = C 03/20 10:35