差分解偏微分方程问题出错求助

aiguo

用有限差分法,构造差分格式时，用的最简单的构造方法，将微分方程中的微分改为差分，即：
dx1/dt=[x1(j,n+1)-x(j,n)]/t
dx1/dz=[x1(j+1,n)-x(j,n)]/h;
h为时间步长和空间步长
程序中取时间步长为t=0.005; 空间步长为h=0.1;
全部程序如下：
%空间有J-1份和时间有N-1份，但是J和N个节点
%% 1.初始化
N=201; % time iteration number
J=51; % z iteration number
t=0.005;h=0.1;
x1(J,N)=0;x2(J,N)=0;x3(J,N)=0;%建立初始矩阵
%初始条件
for j=1:J
x1(j,1)=0;x2(j,1)=0;x3(j,1)=1.9e-19;
end
%边界条件
for n=1:N
x1(1,n)=0;x2(1,n)=0;x3(1,n)=0; %初始边界
x1(J,n)=0;x2(J,n)=0;x3(J,n)=0; %终止边界
end
%除了初始时刻和两边边界外，满足的差分方程
for n=2:N
for j=2:J-1
x1(j,n+1)=t*(x2(j,n)-2e3*x1(j,n)+1.2e-12*x3(j,n)*(x2(j,n)-x1(j,n))-5*(x2(j+1,n)-x2(j,n))/h)+x1(j,n);
x2(j,n+1)=t*(1-x1(j,n)-2*x2(j,n)-1.2e-12*x3(j,n)*(x2(j,n)-x1(j,n))-5*(x2(j+1,n)-x2(j,n))/h)+x2(j,n);
x3(j,n+1)=t*(1.2e8*x3(j,n)*(x2(j,n)-x1(j,n))-2e4*x3(j,n))+x1(j,n);
end;
end;
x1
结果错误！麻烦高手帮忙看一下！有什么问题！
另外我还有几个问题请教：1）用差分格式时，步长如何选择？
2）对于边界条件如何定，比较模糊，有高手能不能给点建议？
这些问题估计是我们这些初接触差分解偏微分方程和MATLAB的菜鸟特别想知道的，如果能够解答，应该比较有价值的
各位高手对于这方面的主题能否总结一下

[ 本帖最后由 eight 于 2007-10-4 18:52 编辑 ]

花如月

问题比较专业，建议找本书籍看看。原理懂了，就是编程的问题了。你的结果不对，说明程序有原理性错误而不是语法错误。外行是一点忙都帮不上的，既然知道结果错了，就找导致错误出现的原因。一步一步来吧，程序都是一条一条写出来的。

xjzuo

程序存在一些问题，边界和初值也古怪－－－请将原问题用word贴一下．
另：差分格式步长的选择需要一定的技巧，与边界条件一样，需要视具体问题而定．例如显式差分与隐式差分步长的选择要求就不一样，热传导情形与FDTD也不一样．．．

aiguo

楼上的xjzuo大哥,谢谢您的指点
能否告之您的QQ,我想向您更具体的请教
我的QQ:393326327

aiguo

楼上大哥要的WORD文档

花如月

讨论版不支持你这种做法。如果需要用qq讨论问题，可以去math群里讨论（加群的办法公告里有）。

xjzuo

相当奇怪的一些人------碰到过很多次了，叫他们用word贴出原问题及参数,目的是想看看原公式及其背景，以便提出建议，可他们都只是将贴出的代码copy一遍------伤脑筋．．．
例如本贴，连处理什么系统、对象都不提一下．．．

[ 本帖最后由 xjzuo 于 2007-10-8 23:43 编辑 ]