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

god

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

t=input('请输入时间步长：'); h=input('请输入x轴步长：'); ; r/ }! f! _5 |. D- ?r=t/h^2;a=1-2*r; 5 ~. ?' E+ G) A+ @" LA=[]; " J4 w4 X1 S% b; e! I1 ?. Gu=[];%第0时间层的数值 m=1/h; j=0:h:1; I=ones(size(j)); 6 O. z9 }5 x6 T. U4 n: k, H u=4*j.*(I-j); 2 z# t6 g# J7 z1 @5 R %第0时间层的值的计算 1 F1 M5 f7 v) C+ w, B+ g/ k4 n7 }) {disp('5秒之后将会显示第0时间层的数值：') : b$ f, Z) K. Ypause $ k, f; ?: H8 T7 Cu 8 L6 {: }$ _( y" I; m0 A1 t& Mw=[];%第1时间层数值的计算 for j=2:m # v) j4 Q7 b% w2 K, |1 | w(j)=[r*u(j)+a*u(j-1)+r*u(j+1)]; w(1)=0;w(m+1)=0; end %通过第0层数值来求第1层的计算 ) K; y9 E f j1 O1 wdisp('5秒之后将会显示第1时间层的数值：') - X( K, Q: z1 r$ Y5 x2 hpause w : f2 z8 [8 X9 o+ Oq=[];%第2时间层的数值 / Q* B% H! i9 z9 m2 Tfor j=2:m - z( T& Y! ^% ~ X# R A2 I q(j)=[r*w(j)+a*w(j-1)+r*w(j+1)]; q(1)=0;q(m+1)=0; 1 @7 R A' p4 ~( eend %通过第1层数值来求第2层的计算 disp('5秒之后将会显示第2时间层的数值：') pause % D8 s+ D/ ]6 O7 m9 \ N3 {q 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时间层的数值：') . P: n8 l8 R- I! ]pause - @1 `* G2 {& m1 x8 w( d# ~f : O% t) K5 E+ I9 ie=[];%第4时间层的数值 for j=2:m ! G5 Z' _: l' I, e# S, n; g e(j)=[r*f(j)+a*f(j-1)+r*f(j+1)]; e(1)=0;e(m+1)=0; : y' r a4 n, ~end %通过第3层数值来求第4层的计算 " v+ O; T$ @ K+ K: T) pdisp('5秒之后将会显示第4时间层的数值：') pause e ) F0 o# W+ L! w" }4 O: B1 J8 e8 ~3 {g=[];%第5时间层的数值 7 M1 h* g. f) [# Q) Vfor j=2:m g(j)=[r*e(j)+a*e(j-1)+r*e(j+1)]; ' h7 i9 [: D( M, [6 P g(1)=0;g(m+1)=0; T: h4 e/ W% M7 F4 Hend %通过第4层数值来求第5层的计算 - X F3 ?9 i4 x, ]. ndisp('5秒之后将会显示第5时间层的数值：') 5 E# G" j4 T; K" npause % @3 W2 G' N! H5 i, N' Sg 7 |1 u/ Z0 ^' p9 J3 H" Wif r<=0.5 disp('以上为上五层数据，因为r小于0.5，在稳定区域，所以方法稳定！') ; N+ E4 u! ~! s# a3 ^: Felse disp('以上为上五层数据，因为r大于0.5，不在稳定区域，所以方法不稳定！') # N8 U/ [& N n3 N2 Jend / ^3 W- f9 n% { 1 O/ r9 D/ |8 d$ q) x N( G有用的,大家可以直接用啊

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

[em05][em05] ; E! e# _- Z8 V: N& x2 t4 P( ?' Y

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