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

刚编好了计算区域,还没有编边界条件,试了一下,但是总是不收敛,望大家赐教.源程序如下,能运行,但不收敛(FDTD-2D TE) 3T?f5+@I  
clear;clc; )d:K:YXt  
F;-90w  
6e<^oH  
% 迭代求解电场和磁场 ><K!~pst}  
1|]xo3j"'  
for t=1:H -,|ha>r  
    pulse=exp(-(t*dt-t0)^2/T^2);%高斯脉冲激励 VVFV8T4  
    Hz(IC,JC)=pulse+Hz(IC,JC);%软源 %E aE,  
yjFe'  
for i=2:IE-1 r^Y~mq  
    for j=2:JE-1 Xo8DEr  
        Hz(i,j)=Hz(i,j)-h0*(Ey(i,j)-Ey(i-1,j)-Ex(i,j)+Ex(i,j-1));%z方向磁场 VXwPdMy*L  
    end A4 5m)wQ  
end GkciA{  
P:^=m*d  
for i=1:IE IkU|W3Vo  
    for j=1:JE-1 St2Q7K5s{  
        Ex(i,j)=Ex(i,j)-e0*(Hz(i,j+1)-Hz(i,j));%x方向电场  $WR?  
    end ~{P:sjsU  
end a*s\Em7f  
I-bF{  
for i=1:IE-1 L1{T ?aII  
    for j=1:JE 9ad6uTc  
      &nb .. @/ k x er  
FQ!Oxlq,Q  

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

递推公式错了。 ?bI?GvSh  
Ex(i,j)=Ex(i,j)-e0*(Hz(i,j+1)-Hz(i,j));%x方向电场 mp:m`sh*i  
'HB~Dbq`V  
====> O]ZC+]}/  
h'!V8'}O?  
Ex(i,j)=Ex(i,j)+e0*(Hz(i,j+1)-Hz(i,j));%x方向电场

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

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

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

     仔细点、认真点，就行了