一个泊松方程的程序,大家一块研究啊
这是我为我女朋友的课程设计写的,如有不足,大家指点啊<P>t=input('请输入时间步长：');
h=input('请输入x轴步长：');
r=t/h^2;a=1-2*r;
A=[];
u=[];%第0时间层的数值
m=1/h;
j=0:h:1;
I=ones(size(j));
u=4*j.*(I-j);
%第0时间层的值的计算
disp('5秒之后将会显示第0时间层的数值：')
pause
u
w=[];%第1时间层数值的计算
for j=2:m
w(j)=;
w(1)=0;w(m+1)=0;
end %通过第0层数值来求第1层的计算
disp('5秒之后将会显示第1时间层的数值：')
pause
w
q=[];%第2时间层的数值
for j=2:m
q(j)=;
q(1)=0;q(m+1)=0;
end %通过第1层数值来求第2层的计算
disp('5秒之后将会显示第2时间层的数值：')
pause
q
f=[];%第3时间层的数值
for j=2:m
f(j)=;
f(1)=0;f(m+1)=0;
end %通过第2层数值来求第3层的计算
disp('5秒之后将会显示第3时间层的数值：')
pause
f
e=[];%第4时间层的数值
for j=2:m
e(j)=;
e(1)=0;e(m+1)=0;
end %通过第3层数值来求第4层的计算
disp('5秒之后将会显示第4时间层的数值：')
pause
e
g=[];%第5时间层的数值
for j=2:m
g(j)=;
g(1)=0;g(m+1)=0;
end %通过第4层数值来求第5层的计算
disp('5秒之后将会显示第5时间层的数值：')
pause
g
if r<=0.5
disp('以上为上五层数据，因为r小于0.5，在稳定区域，所以方法稳定！')
else
disp('以上为上五层数据，因为r大于0.5，不在稳定区域，所以方法不稳定！')
end
有用的,大家可以直接用啊</P>
<P>时间有限,其中的函数大家可以从程序中获知!!! </P>
[此贴子已经被作者于2004-12-27 21:52:53编辑过]
页:
[1]