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

god

这是我为我女朋友的课程设计写的,如有不足,大家指点啊 1 v) F' ]3 E' e0 I1 e

t=input('请输入时间步长：'); # r) L2 I! E/ V. t3 Uh=input('请输入x轴步长：'); $ ~# \; y/ L; g- T" `r=t/h^2;a=1-2*r; A=[]; u=[];%第0时间层的数值 6 ]/ w; z- E' k" F% s: cm=1/h; 3 @. k f7 B% }j=0:h:1; 8 g, X( f, L1 I. Y4 hI=ones(size(j)); u=4*j.*(I-j); ; f% v% r) h r. o( R% t %第0时间层的值的计算 disp('5秒之后将会显示第0时间层的数值：') . x+ k+ s @- |+ `1 v/ w' E% ipause u ) m& z1 O3 J2 w1 @, pw=[];%第1时间层数值的计算 - _% w6 q, J+ u0 ? W7 |for j=2:m 6 L: N/ |3 i- H) Z+ [ w(j)=[r*u(j)+a*u(j-1)+r*u(j+1)]; 2 ^" k4 P6 W; r6 l" q7 L. M w(1)=0;w(m+1)=0; end %通过第0层数值来求第1层的计算 disp('5秒之后将会显示第1时间层的数值：') 2 v# Z, b1 X$ p" p0 Bpause 2 V6 x: N, T: T+ P6 l, rw q=[];%第2时间层的数值 for j=2:m q(j)=[r*w(j)+a*w(j-1)+r*w(j+1)]; / {" |5 Q8 e) U, X7 v6 N9 e: {4 K q(1)=0;q(m+1)=0; ) Y+ S4 M# l$ X0 b8 s; Jend %通过第1层数值来求第2层的计算 disp('5秒之后将会显示第2时间层的数值：') ! R" c4 w3 q# n4 @, { ^! ?pause q / J. G, v# E, _ ^, l# b: A, {f=[];%第3时间层的数值 ' O" n$ Q! S7 p* W* a1 ufor j=2:m - V) \7 ^6 m9 U& | f(j)=[r*q(j)+a*q(j-1)+r*q(j+1)]; 6 m; k+ b. H% { A% ~ f(1)=0;f(m+1)=0; end %通过第2层数值来求第3层的计算 disp('5秒之后将会显示第3时间层的数值：') 8 q8 T6 Z) k- t& Qpause f e=[];%第4时间层的数值 7 H W; M L) y$ a, nfor j=2:m 5 D+ ]( ?, Y0 L$ w" C/ q: D 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时间层的数值：') pause 1 @ Y/ Y8 R7 u" t# T5 E4 ke g=[];%第5时间层的数值 for j=2:m g(j)=[r*e(j)+a*e(j-1)+r*e(j+1)]; ' S7 T# L8 Z, W; S+ F1 f7 E g(1)=0;g(m+1)=0; , P4 I) w, R3 _* s `end %通过第4层数值来求第5层的计算 disp('5秒之后将会显示第5时间层的数值：') . i' k, }7 X6 h% L( r0 z' @pause g % C8 E0 u8 P4 o- T1 m' Bif r<=0.5 , o* A: g& R0 A" O- r disp('以上为上五层数据，因为r小于0.5，在稳定区域，所以方法稳定！') else disp('以上为上五层数据，因为r大于0.5，不在稳定区域，所以方法不稳定！') end 有用的,大家可以直接用啊

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

[em05][em05]

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