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

god

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

t=input('请输入时间步长：'); 2 Z$ S0 g s8 d9 z9 [h=input('请输入x轴步长：'); r=t/h^2;a=1-2*r; 1 m, w& M8 Y0 c3 f0 QA=[]; & {' x1 c o, w+ c$ ^u=[];%第0时间层的数值 ' E& Y& r0 d) R c1 i* r. B* h; X1 ~m=1/h; : I' l- r' i! d/ {& T' \( I7 Vj=0:h:1; I=ones(size(j)); ( D4 I9 |+ [& A) c* L0 z u=4*j.*(I-j); 7 |5 O9 B) w8 H8 X! e8 K %第0时间层的值的计算 ) R: v* W" R, c4 |) @disp('5秒之后将会显示第0时间层的数值：') pause K: r5 M! W, _u w=[];%第1时间层数值的计算 + Q% C, d0 H4 }7 }- d1 zfor j=2:m w(j)=[r*u(j)+a*u(j-1)+r*u(j+1)]; w(1)=0;w(m+1)=0; end %通过第0层数值来求第1层的计算 disp('5秒之后将会显示第1时间层的数值：') pause ' D) E6 I5 u; E7 aw q=[];%第2时间层的数值 for j=2:m q(j)=[r*w(j)+a*w(j-1)+r*w(j+1)]; ! l9 H. d( u, }; `" ` q(1)=0;q(m+1)=0; end %通过第1层数值来求第2层的计算 disp('5秒之后将会显示第2时间层的数值：') pause " t* U; k1 [0 d. \& k' `5 ?) J) |q . z) ?1 i% Y/ r% j- w7 if=[];%第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; 1 Q2 d- u+ @* P6 g: U% r8 nend %通过第2层数值来求第3层的计算 3 V7 G4 |1 o2 s" Xdisp('5秒之后将会显示第3时间层的数值：') ; }. W! V' d. q! Y# ipause f ; ~( t- b! L9 p1 i, r* [% u- V8 `e=[];%第4时间层的数值 ( _: y% s N$ M' X6 d$ d+ qfor 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层的计算 ; S/ q! C/ i% g8 edisp('5秒之后将会显示第4时间层的数值：') pause , {3 U3 y: T! ae % `& d: z- q) H) }* Rg=[];%第5时间层的数值 for j=2:m * D4 J9 j: r4 ~$ ~ i0 ^+ I g(j)=[r*e(j)+a*e(j-1)+r*e(j+1)]; g(1)=0;g(m+1)=0; end %通过第4层数值来求第5层的计算 * J9 l8 l# Y" g. F; l* d3 i& wdisp('5秒之后将会显示第5时间层的数值：') & z. Y8 L+ N( r/ hpause 7 `* W* j4 Q7 Og 1 Z X% o0 z; p( e0 ]& Iif r<=0.5 ) N% P( O, @- N2 j disp('以上为上五层数据，因为r小于0.5，在稳定区域，所以方法稳定！') + }! O& w0 u! @! Z4 z$ [4 v' H3 b# Velse disp('以上为上五层数据，因为r大于0.5，不在稳定区域，所以方法不稳定！') end 4 B X, ?! h3 ^9 f) n2 `' s4 H ! U8 Y! Z+ V6 ] H9 S有用的,大家可以直接用啊

( z2 V1 S' Y' x8 w% y% i; C; z7 z

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

[em05][em05]

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