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

god

这是我为我女朋友的课程设计写的,如有不足,大家指点啊 2 o5 n: Q$ n( p; u0 x

t=input('请输入时间步长：'); h=input('请输入x轴步长：'); 3 f) L/ L6 h' r$ r$ l3 {r=t/h^2;a=1-2*r; / T8 y, s3 T2 RA=[]; . a+ D; W1 W, Tu=[];%第0时间层的数值 ! \8 y+ R) k Am=1/h; ; o. G! S( H2 {3 Sj=0:h:1; + K, y# N/ i+ d( O, N9 h8 YI=ones(size(j)); ! f6 Q1 _ t5 \- f4 @1 ~* Z u=4*j.*(I-j); %第0时间层的值的计算 disp('5秒之后将会显示第0时间层的数值：') M9 G6 @$ y. Vpause u * H6 }# ^9 E5 I4 d8 N# T _5 O9 pw=[];%第1时间层数值的计算 . T$ p# ^8 P# L* x, a* W, \7 q8 yfor j=2:m w(j)=[r*u(j)+a*u(j-1)+r*u(j+1)]; 0 z5 m9 A; V: ?$ B" ^# S w(1)=0;w(m+1)=0; end %通过第0层数值来求第1层的计算 ! ?+ D" Z8 w8 c0 Z, G ?6 h6 \+ |disp('5秒之后将会显示第1时间层的数值：') + _& y3 `9 h. S+ Ppause 0 c) T1 g1 I+ b, \6 C1 }w q=[];%第2时间层的数值 for j=2:m q(j)=[r*w(j)+a*w(j-1)+r*w(j+1)]; q(1)=0;q(m+1)=0; end %通过第1层数值来求第2层的计算 3 r8 E D# e) Odisp('5秒之后将会显示第2时间层的数值：') 0 n/ J1 [ _6 }pause q 2 U( V: }# e' n% Cf=[];%第3时间层的数值 for j=2:m f(j)=[r*q(j)+a*q(j-1)+r*q(j+1)]; f(1)=0;f(m+1)=0; end %通过第2层数值来求第3层的计算 disp('5秒之后将会显示第3时间层的数值：') 6 i" \# X2 z! n" qpause 4 n7 `9 T2 f6 ~" X) V. Sf e=[];%第4时间层的数值 0 h4 f: K# \; W- {' H' k. }for j=2:m e(j)=[r*f(j)+a*f(j-1)+r*f(j+1)]; e(1)=0;e(m+1)=0; end %通过第3层数值来求第4层的计算 1 k! N+ x) S# e7 Adisp('5秒之后将会显示第4时间层的数值：') 4 k, B) ]2 R- S" jpause 4 a4 K5 A! E$ O3 z& C1 T' me g=[];%第5时间层的数值 + |% a4 B) t0 e6 g( H; mfor j=2:m g(j)=[r*e(j)+a*e(j-1)+r*e(j+1)]; g(1)=0;g(m+1)=0; 5 B5 d) I5 r, d+ iend %通过第4层数值来求第5层的计算 disp('5秒之后将会显示第5时间层的数值：') ; E( b4 ^; v- L' E$ t( y6 \, G6 o6 Ypause g 8 i4 _9 P1 j! k, T* L: F5 Z0 }if r<=0.5 / Q1 C& c$ q0 ], R" \' a- U! { disp('以上为上五层数据，因为r小于0.5，在稳定区域，所以方法稳定！') else disp('以上为上五层数据，因为r大于0.5，不在稳定区域，所以方法不稳定！') end \) f- R( Z- O: P 有用的,大家可以直接用啊

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

[em05][em05]

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