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

god

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

t=input('请输入时间步长：'); h=input('请输入x轴步长：'); r=t/h^2;a=1-2*r; A=[]; u=[];%第0时间层的数值 - m/ C8 k7 N& ]3 `m=1/h; 3 @2 q8 A; e+ Y, zj=0:h:1; I=ones(size(j)); u=4*j.*(I-j); 5 a3 [, |$ \1 p1 O %第0时间层的值的计算 disp('5秒之后将会显示第0时间层的数值：') pause * t& [3 ]' p9 h+ d( i0 qu w=[];%第1时间层数值的计算 for j=2:m w(j)=[r*u(j)+a*u(j-1)+r*u(j+1)]; & u; @/ h, w7 t9 G8 h0 @ w(1)=0;w(m+1)=0; end %通过第0层数值来求第1层的计算 disp('5秒之后将会显示第1时间层的数值：') . N! g9 n6 B5 |, [pause + D# p% B$ R8 a' k' ^6 Q, Cw $ _1 z" _& l3 Uq=[];%第2时间层的数值 for j=2:m - q: z4 w5 O, ^ T& X, Q q(j)=[r*w(j)+a*w(j-1)+r*w(j+1)]; ! l/ B! z! x! E# Z# N G! X# ] q(1)=0;q(m+1)=0; 3 y7 s9 v6 c" f' Rend %通过第1层数值来求第2层的计算 disp('5秒之后将会显示第2时间层的数值：') pause 3 I% K9 u% C1 f1 z) jq f=[];%第3时间层的数值 $ J) g2 X% }9 P! \- u% [8 Wfor j=2:m f(j)=[r*q(j)+a*q(j-1)+r*q(j+1)]; n6 E4 d! j/ O. e2 K: X f(1)=0;f(m+1)=0; ) k/ }5 n5 W: zend %通过第2层数值来求第3层的计算 3 e4 `7 ]# F! m! }4 E7 Y& ^disp('5秒之后将会显示第3时间层的数值：') pause / j3 R/ ~: b2 l# Z* M. r+ Y4 x$ ?f e=[];%第4时间层的数值 / b3 ?0 T2 G/ x+ h' K9 v. ofor 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层的计算 disp('5秒之后将会显示第4时间层的数值：') ' t1 R. b8 {( e) E; n4 j6 F' qpause e o; _$ M, E' C1 c3 O0 p% Z7 I# X$ k2 bg=[];%第5时间层的数值 % }7 f L& i$ `$ `2 {8 d8 s. C: Hfor j=2:m 6 U8 T* B. l- D' l$ \. z1 B* b g(j)=[r*e(j)+a*e(j-1)+r*e(j+1)]; g(1)=0;g(m+1)=0; $ x* [# K: a) Y5 f% ?" u6 |- Mend %通过第4层数值来求第5层的计算 disp('5秒之后将会显示第5时间层的数值：') % i( e" `+ F: C" J9 r2 L8 fpause g if r<=0.5 6 L, K9 z# w4 E( a' U. ?: N disp('以上为上五层数据，因为r小于0.5，在稳定区域，所以方法稳定！') else disp('以上为上五层数据，因为r大于0.5，不在稳定区域，所以方法不稳定！') end 有用的,大家可以直接用啊

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

[em05][em05] . ^# }( R: z1 Y8 d5 H* G7 {; c- o

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