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

god

这是我为我女朋友的课程设计写的,如有不足,大家指点啊 / D# G! u- F& J9 d) S 5 k, B3 S }3 q

t=input('请输入时间步长：'); 1 I! t G( ]( @, |( th=input('请输入x轴步长：'); r=t/h^2;a=1-2*r; A=[]; 2 z8 l1 d6 T, |4 ]2 ^3 {' Z1 Mu=[];%第0时间层的数值 m=1/h; j=0:h:1; 9 I; Z+ e: L7 o+ m7 OI=ones(size(j)); u=4*j.*(I-j); %第0时间层的值的计算 disp('5秒之后将会显示第0时间层的数值：') pause . N, O! A' L) p5 N! @$ Bu 9 A) P, { e- I, I* ~! M+ }# ww=[];%第1时间层数值的计算 & l1 p( [3 \/ Y! nfor j=2:m " m, x/ ]' t5 I0 e w(j)=[r*u(j)+a*u(j-1)+r*u(j+1)]; w(1)=0;w(m+1)=0; # \0 t# @* o! \, jend %通过第0层数值来求第1层的计算 disp('5秒之后将会显示第1时间层的数值：') % d* [" t! X% Ppause & p' D7 E) g3 B/ w& c( Z( W) T7 kw , c5 C1 n9 z# X# m: k$ X& Mq=[];%第2时间层的数值 - k+ A1 v7 x d( F: K7 e: C8 `for j=2:m 5 d0 p |2 p, h( O( H7 ~5 l% A q(j)=[r*w(j)+a*w(j-1)+r*w(j+1)]; q(1)=0;q(m+1)=0; 6 N( Y) d, y1 \( z( Rend %通过第1层数值来求第2层的计算 disp('5秒之后将会显示第2时间层的数值：') pause q : Y+ g E6 h3 V+ m5 Tf=[];%第3时间层的数值 for j=2:m f(j)=[r*q(j)+a*q(j-1)+r*q(j+1)]; 0 V! {- T& x/ D5 N+ k( f/ e f(1)=0;f(m+1)=0; end %通过第2层数值来求第3层的计算 0 \( T- N' A0 M8 a3 n1 l' h2 jdisp('5秒之后将会显示第3时间层的数值：') 7 O J- S0 T2 k8 {0 o+ Zpause f e=[];%第4时间层的数值 for j=2:m ! Q$ m M2 |% @ e(j)=[r*f(j)+a*f(j-1)+r*f(j+1)]; 4 S( Y8 ~" n3 B8 Z$ k- \1 y: Q e(1)=0;e(m+1)=0; end %通过第3层数值来求第4层的计算 disp('5秒之后将会显示第4时间层的数值：') 4 `! y3 e7 }, H5 @& ]pause " ?: l, r) |& Z7 \e g=[];%第5时间层的数值 / i, R* x4 }) O+ e. l6 C. d- ^for j=2:m g(j)=[r*e(j)+a*e(j-1)+r*e(j+1)]; ( B5 L/ t8 w2 w& C g(1)=0;g(m+1)=0; end %通过第4层数值来求第5层的计算 disp('5秒之后将会显示第5时间层的数值：') pause g if r<=0.5 9 U! g; u8 t1 { disp('以上为上五层数据，因为r小于0.5，在稳定区域，所以方法稳定！') else 8 G5 b1 E# a& [; J0 H% T1 \# a8 g# I disp('以上为上五层数据，因为r大于0.5，不在稳定区域，所以方法不稳定！') end ! D" g8 b' B. R$ `- C有用的,大家可以直接用啊

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

[em05][em05]

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