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

刚编好了计算区域,还没有编边界条件,试了一下,但是总是不收敛,望大家赐教.源程序如下,能运行,但不收敛(FDTD-2D TE) _p*9LsN$L  
clear;clc; `3F#k[IR  
BX?DI-o^h  
o^@"eG
$,  
% 迭代求解电场和磁场 45c?0tj  
j9NF|  
for t=1:H vN{@c(=g  
    pulse=exp(-(t*dt-t0)^2/T^2);%高斯脉冲激励 TN0KS]^A3  
    Hz(IC,JC)=pulse+Hz(IC,JC);%软源 Cq2Wpu-u  
E%`J=C}  
for i=2:IE-1 mU #F>  
    for j=2:JE-1 ]j3>=Jb;  
        Hz(i,j)=Hz(i,j)-h0*(Ey(i,j)-Ey(i-1,j)-Ex(i,j)+Ex(i,j-1));%z方向磁场 n4kq=Z%
 
    end bgorW"'  
end x]c8?H9,&  
H.ZIRt!RB  
for i=1:IE >R F|Q  
    for j=1:JE-1 M1kA-Xr  
        Ex(i,j)=Ex(i,j)-e0*(Hz(i,j+1)-Hz(i,j));%x方向电场 >$
ZG=&  
    end 5.6tVr  
end UGPDwgq\v  
Y*c]C;%=  
for i=1:IE-1 ~*Y#Y{  
    for j=1:JE -$Z1X_~;)<  
      &nb .. 5tlRrf  
=K`.$R  

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

递推公式错了。 jW:7PS  
Ex(i,j)=Ex(i,j)-e0*(Hz(i,j+1)-Hz(i,j));%x方向电场 *u$aItx  
` i[26Qb  
====> E'WXi!>7p  
)O[8 D  
Ex(i,j)=Ex(i,j)+e0*(Hz(i,j+1)-Hz(i,j));%x方向电场

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

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

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

     仔细点、认真点，就行了