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

god

这是我为我女朋友的课程设计写的,如有不足,大家指点啊 4 B; `0 @) W- \/ B& f

t=input('请输入时间步长：'); h=input('请输入x轴步长：'); r=t/h^2;a=1-2*r; A=[]; ) w$ J8 H. V/ U6 {u=[];%第0时间层的数值 m=1/h; ' k$ N1 c4 c' ?& z' Ej=0:h:1; I=ones(size(j)); u=4*j.*(I-j); %第0时间层的值的计算 disp('5秒之后将会显示第0时间层的数值：') 0 S0 d {" A0 z8 J+ H* Q* Kpause u - ^5 W2 i8 f: I8 e2 H4 m$ Bw=[];%第1时间层数值的计算 * s/ g9 X( s/ P6 }+ s& i9 k9 _1 T) k3 P1 Xfor j=2:m w(j)=[r*u(j)+a*u(j-1)+r*u(j+1)]; w(1)=0;w(m+1)=0; * u: \! K- [/ _' ^4 V/ r% send %通过第0层数值来求第1层的计算 disp('5秒之后将会显示第1时间层的数值：') ; W6 s: {) H% u: \' d/ g5 Z5 }pause 6 I0 y4 \6 A7 s" Y* }) o; rw 6 I' t1 U7 g2 r a+ o: Iq=[];%第2时间层的数值 for j=2:m , E1 r* X& J6 F# w q(j)=[r*w(j)+a*w(j-1)+r*w(j+1)]; S' E7 {4 m5 ], g! A9 j q(1)=0;q(m+1)=0; 7 C6 r2 `* Y' x4 p/ a6 [end %通过第1层数值来求第2层的计算 6 Z% d& V$ H3 ^$ {& w7 Gdisp('5秒之后将会显示第2时间层的数值：') pause q / @: R& U" k4 j# |0 `f=[];%第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时间层的数值：') 8 g* b/ |2 J2 v3 ppause + |! Y9 v. w1 D: \" J/ a. }f # ?0 ~3 Z: {; [e=[];%第4时间层的数值 4 T7 G m6 _' D7 h5 i+ Bfor j=2:m e(j)=[r*f(j)+a*f(j-1)+r*f(j+1)]; $ v. ?( V; J. z2 ~; Q e(1)=0;e(m+1)=0; end %通过第3层数值来求第4层的计算 disp('5秒之后将会显示第4时间层的数值：') pause 5 a/ e$ ?6 s2 oe 2 E {% }* f- i4 zg=[];%第5时间层的数值 ' T) Q& R7 g$ O" bfor j=2:m g(j)=[r*e(j)+a*e(j-1)+r*e(j+1)]; g(1)=0;g(m+1)=0; end %通过第4层数值来求第5层的计算 : x' ?$ [5 \6 C) n2 O! `: Q zdisp('5秒之后将会显示第5时间层的数值：') ( j" G: O! l: |) y# u3 @' _pause : @- j9 F; \8 J; a( w3 T" Qg if r<=0.5 disp('以上为上五层数据，因为r小于0.5，在稳定区域，所以方法稳定！') % v) K. n8 x4 J% v9 _& ~8 o, I' @( n& b" Uelse disp('以上为上五层数据，因为r大于0.5，不在稳定区域，所以方法不稳定！') & J k$ y, O: \( uend 2 ]" b8 m: U( C. N$ D0 F有用的,大家可以直接用啊

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

[em05][em05]

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