一个泊松方程的程序,大家一块研究啊

god

这是我为我女朋友的课程设计写的,如有不足,大家指点啊

t=input('请输入时间步长：'); h=input('请输入x轴步长：'); $ t+ n6 D& q- Z& q+ tr=t/h^2;a=1-2*r; : r1 A2 ~+ M- V U9 L2 rA=[]; / a( W$ \) V, `) b8 gu=[];%第0时间层的数值 8 D! Y( C) k' @3 o7 G& Z& @m=1/h; j=0:h:1; I=ones(size(j)); u=4*j.*(I-j); 8 s) F0 m- q( m' G %第0时间层的值的计算 ! M( _/ r/ ]. t1 s+ G8 Y7 j2 `disp('5秒之后将会显示第0时间层的数值：') pause u w=[];%第1时间层数值的计算 4 C7 C/ R2 Z0 L3 _: K/ zfor j=2:m w(j)=[r*u(j)+a*u(j-1)+r*u(j+1)]; : F! B& t6 x: B$ G& q* y. @9 | w(1)=0;w(m+1)=0; : J' Y4 w! N# ~4 t, w. rend %通过第0层数值来求第1层的计算 ) \! J' o# p8 f" m; Odisp('5秒之后将会显示第1时间层的数值：') pause w q=[];%第2时间层的数值 4 w# G8 W- v- V0 w& V% Z0 T3 ofor j=2:m q(j)=[r*w(j)+a*w(j-1)+r*w(j+1)]; ^* T/ L4 [" ~5 o q(1)=0;q(m+1)=0; + |2 T/ f, n( X2 O+ [0 k+ V6 Z- bend %通过第1层数值来求第2层的计算 1 \/ T3 X6 i! M' A/ E& ^) Tdisp('5秒之后将会显示第2时间层的数值：') ' Y7 V9 b2 k) }pause q - v& ]! I1 P& e# q0 f% r( D) Jf=[];%第3时间层的数值 for j=2:m $ v3 M7 Y( C. ~8 i& I f(j)=[r*q(j)+a*q(j-1)+r*q(j+1)]; 9 \8 b& a8 U" ` y' K6 v# \7 F f(1)=0;f(m+1)=0; 5 {6 {% i( \, c# U) [end %通过第2层数值来求第3层的计算 , Y" m. z7 L0 N" H; `disp('5秒之后将会显示第3时间层的数值：') pause f $ C% k" G! j0 v9 V! x# ~e=[];%第4时间层的数值 for j=2:m . s; G! u: ~9 l' [1 m3 i1 X4 t e(j)=[r*f(j)+a*f(j-1)+r*f(j+1)]; ( B1 E, Q. K2 @: \) } e(1)=0;e(m+1)=0; ( |( w1 M" D/ }- qend %通过第3层数值来求第4层的计算 ' }2 F& p x- R8 r5 W. [+ ?disp('5秒之后将会显示第4时间层的数值：') $ g7 s- B& M8 |6 ?) Bpause " b1 g" o* d: o2 S+ Ce 0 Q% y. n7 N- t ` F& @g=[];%第5时间层的数值 for j=2:m / Y5 {3 M5 |. K; W+ T g(j)=[r*e(j)+a*e(j-1)+r*e(j+1)]; % f' h7 ]/ Y1 U3 y# P: \; J g(1)=0;g(m+1)=0; * p% V2 v' P7 qend %通过第4层数值来求第5层的计算 disp('5秒之后将会显示第5时间层的数值：') pause : M0 J2 n! G9 Og 1 t& e1 _3 |. ~if r<=0.5 " W0 z% z* m5 M5 a disp('以上为上五层数据，因为r小于0.5，在稳定区域，所以方法稳定！') else * z+ u8 Q% J+ j# B disp('以上为上五层数据，因为r大于0.5，不在稳定区域，所以方法不稳定！') / L" Q" t Z( b `) }$ k0 e N8 Vend 有用的,大家可以直接用啊

5 G7 Y) }/ V p, H

时间有限,其中的函数大家可以从程序中获知!!!

[em05][em05]

[此贴子已经被作者于2004-12-27 21:52:53编辑过]