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

god

这是我为我女朋友的课程设计写的,如有不足,大家指点啊 6 U) T" Y' H* O- H

t=input('请输入时间步长：'); ) w* O, L; D5 a( gh=input('请输入x轴步长：'); & x4 y$ r8 N6 F2 q& W9 j, Z0 L/ or=t/h^2;a=1-2*r; A=[]; " |) c; A$ P2 A+ v$ Tu=[];%第0时间层的数值 8 E- c% s& t9 ~) W0 Bm=1/h; : N. h! e8 g, y* C0 mj=0:h:1; I=ones(size(j)); 4 s) y$ }. k0 |( u) s u=4*j.*(I-j); %第0时间层的值的计算 disp('5秒之后将会显示第0时间层的数值：') 1 W+ t' L& A0 h* t1 Gpause , R( {5 n; N6 E. Vu / j$ {. D! L, \# m7 v5 o' N R6 |w=[];%第1时间层数值的计算 for j=2:m w(j)=[r*u(j)+a*u(j-1)+r*u(j+1)]; w(1)=0;w(m+1)=0; * P; B _% `" v" c4 t" kend %通过第0层数值来求第1层的计算 disp('5秒之后将会显示第1时间层的数值：') / Q: ~; O1 M! m/ Dpause w ) ]4 L. s5 L' Fq=[];%第2时间层的数值 for j=2:m - C `# t& ^1 J+ f# g q(j)=[r*w(j)+a*w(j-1)+r*w(j+1)]; q(1)=0;q(m+1)=0; end %通过第1层数值来求第2层的计算 disp('5秒之后将会显示第2时间层的数值：') pause q f=[];%第3时间层的数值 + T- I* r3 R; O, p# Ffor j=2:m 9 p) w0 b3 T; G- z8 s* n f(j)=[r*q(j)+a*q(j-1)+r*q(j+1)]; f(1)=0;f(m+1)=0; end %通过第2层数值来求第3层的计算 6 S6 R2 L# l4 p' Y- Odisp('5秒之后将会显示第3时间层的数值：') pause f e=[];%第4时间层的数值 - B7 r6 p! U. @for j=2:m e(j)=[r*f(j)+a*f(j-1)+r*f(j+1)]; , c x2 s& e& J0 f$ v3 w. \2 ` e(1)=0;e(m+1)=0; 9 M( c' x7 Y4 N- P: ]end %通过第3层数值来求第4层的计算 / a# I5 m- \) R& v8 Fdisp('5秒之后将会显示第4时间层的数值：') ( u+ y; X3 p: q) Z5 K, \1 N* w4 ppause 8 I) \* j% ~3 e$ J: O1 re g=[];%第5时间层的数值 for j=2:m 6 r6 X, g3 b# O g(j)=[r*e(j)+a*e(j-1)+r*e(j+1)]; g(1)=0;g(m+1)=0; 4 H3 T# X/ G+ b! D# o% e$ Hend %通过第4层数值来求第5层的计算 disp('5秒之后将会显示第5时间层的数值：') pause g if r<=0.5 disp('以上为上五层数据，因为r小于0.5，在稳定区域，所以方法稳定！') ' |2 B; y! F v Q% E/ Aelse . j4 z2 ~7 |% n$ N disp('以上为上五层数据，因为r大于0.5，不在稳定区域，所以方法不稳定！') - h1 [" D7 V! F/ k7 c; {end & z: S* p) o5 m8 T6 m 9 s! q/ P: c$ K$ V& G% m有用的,大家可以直接用啊

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

[em05][em05]

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