[已解决]matlab编的FDTD不收敛,求助

刚编好了计算区域,还没有编边界条件,试了一下,但是总是不收敛,望大家赐教.源程序如下,能运行,但不收敛(FDTD-2D TE) `y*o-St3  
clear;clc; 7b hJt_`Q  
$@eFSA5k,7  
e$4$G<8;y  
% 迭代求解电场和磁场 G0Wv=tX|  
c.Do b?5  
for t=1:H K)nn;j=  
    pulse=exp(-(t*dt-t0)^2/T^2);%高斯脉冲激励 QvDD   
    Hz(IC,JC)=pulse+Hz(IC,JC);%软源 e"]DIy4s  
- %`iLu  
for i=2:IE-1 e)kVS}e?  
    for j=2:JE-1 Pi'[d7o  
        Hz(i,j)=Hz(i,j)-h0*(Ey(i,j)-Ey(i-1,j)-Ex(i,j)+Ex(i,j-1));%z方向磁场 2@?\"kR"!  
    end QmY1Bn?s  
end Qe4"a*l-r  
X@^"@  
for i=1:IE o0FVVSl  
    for j=1:JE-1 mAuN* (  
        Ex(i,j)=Ex(i,j)-e0*(Hz(i,j+1)-Hz(i,j));%x方向电场 }!1pA5x$  
    end ;2Za]%'  
end <lf6gb  
?T7ndXX  
for i=1:IE-1 >OW>^%\!1  
    for j=1:JE .*c%A^>  
      &nb .. i4\m/&of3y  
r=:o$e  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册

递推公式错了。 z{ydP Ra  
Ex(i,j)=Ex(i,j)-e0*(Hz(i,j+1)-Hz(i,j));%x方向电场 Th\t6K~  
U8s&5~IPn  
====>
&W:R#/|  
Tv`-h  
Ex(i,j)=Ex(i,j)+e0*(Hz(i,j+1)-Hz(i,j));%x方向电场

成功了,万分感谢,不过相当惭愧,犯了个低级错误.

写code都这样的，很多笔误什么的。把x写成y，在比较长的code里面很常见的。我曾经调过一个星期的程序，结果发现在一个地方笔误了。   

个人建议你，下一步，把单向差分改成中间差分吧，这样误差小一点，和葛书上的公式保持一致，这样看资料习惯点。

     仔细点、认真点，就行了