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

god

这是我为我女朋友的课程设计写的,如有不足,大家指点啊 ) z/ j- D! |8 S, `3 Q% I* p% ?, [

t=input('请输入时间步长：'); h=input('请输入x轴步长：'); r=t/h^2;a=1-2*r; A=[]; u=[];%第0时间层的数值 , q3 q' m1 Z* q D! ?0 E. s: Bm=1/h; 9 |6 ~. ^* @4 d; J- Pj=0:h:1; $ q9 L1 d* H$ n2 |8 i- Y) @0 Q/ O; DI=ones(size(j)); u=4*j.*(I-j); + J/ M7 C" Q1 {6 f %第0时间层的值的计算 8 N: G+ `- _" ]& e+ a' g7 M3 Pdisp('5秒之后将会显示第0时间层的数值：') pause * _) g& w7 F. n: w# ju 3 A1 p8 Y, d6 t! o7 t1 Sw=[];%第1时间层数值的计算 . S6 Z$ d9 T D* ~& Yfor j=2:m w(j)=[r*u(j)+a*u(j-1)+r*u(j+1)]; ) _) I4 a& w) m( B5 ]8 h w(1)=0;w(m+1)=0; end %通过第0层数值来求第1层的计算 8 F. [' e0 |" C4 ^disp('5秒之后将会显示第1时间层的数值：') 5 p! u4 V V7 M6 qpause w q=[];%第2时间层的数值 for j=2:m * E4 h0 j. p( C1 I- n q(j)=[r*w(j)+a*w(j-1)+r*w(j+1)]; q(1)=0;q(m+1)=0; * r/ q0 `& d; X* N+ |2 ? B8 Pend %通过第1层数值来求第2层的计算 : ]0 t t/ o; r0 O. T8 Pdisp('5秒之后将会显示第2时间层的数值：') 7 p$ l" F: N$ A: V' T& G r# Spause q f=[];%第3时间层的数值 % t6 Z2 j+ E1 `' d' p* x' m) {2 Ifor j=2:m f(j)=[r*q(j)+a*q(j-1)+r*q(j+1)]; , l9 d% F% \4 l, F f(1)=0;f(m+1)=0; end %通过第2层数值来求第3层的计算 ! w3 C7 Z: f/ Mdisp('5秒之后将会显示第3时间层的数值：') pause : q1 e p9 z, j) ?# Rf ) X2 Q8 w8 ?9 D0 o; s. o$ A3 l# Ye=[];%第4时间层的数值 & v1 J7 S. w* t3 efor j=2:m e(j)=[r*f(j)+a*f(j-1)+r*f(j+1)]; e(1)=0;e(m+1)=0; 1 P5 E4 f1 i9 b$ @& ^end %通过第3层数值来求第4层的计算 disp('5秒之后将会显示第4时间层的数值：') pause / o. O- k) s: D/ ?3 i7 Re g=[];%第5时间层的数值 for j=2:m 5 K. }: A% K; q6 ~: I g(j)=[r*e(j)+a*e(j-1)+r*e(j+1)]; g(1)=0;g(m+1)=0; ( ^2 {1 L# x+ N3 i& }8 j t$ Bend %通过第4层数值来求第5层的计算 disp('5秒之后将会显示第5时间层的数值：') pause g if r<=0.5 9 b* p4 k6 ^: M disp('以上为上五层数据，因为r小于0.5，在稳定区域，所以方法稳定！') $ J8 L2 g( W. q* I' m3 Oelse disp('以上为上五层数据，因为r大于0.5，不在稳定区域，所以方法不稳定！') end 2 f; A6 _% [2 e5 V4 a1 N. P有用的,大家可以直接用啊

; t* |7 N, W/ ~ e- U

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

[em05][em05]

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