[求助]Sullivan Dennis书上1d fdtd请教

electromagnetic simulation using fdtd method (Sullivan Dennis) 书上的1.2题matlab代码如下，不知道仿真对没？请大佬指点 /~fu,2=7  
5xNOIOpDB  
TM_bu  
clear; 2.v`J=R  
KE=200; f.oY:3h:  
ex=1:KE+1; )-}<}< oO  
ex(1:KE+1)=0;
$-!7<a-  
hy=1:KE+1; iXXaB+w  
hy(1:KE+1)=0; oblw!)  
a=1:KE+1; e},:QL0X  
a(1:KE+1)=0; U-f8D  
kc = (KE / 2)+1;     Qx!Bf_,J  
t0 = 40.0;     1Ak0A6E  
spread = 12;     wq1s#ag<  
T = 0; _z.CV<  
Time1=40; 6I
z!_  
NSTEPS=Time1; Fy$C._C$  
M=moviein(Time1); j>v8i bS(  
for b = 1:Time1 Tkp"mT v?<  
     !/|B4Yv  
for n=1:NSTEPS 7X
>@r"9<  
     wnLpf  
    T=T+1; |Dq?<Ha  
for k=2:KE+1 (R RRG;*n#  
             dpt P(H  
  % ex(k) = ex(k) +0.5*(hy(k - 1) - hy(k)); 0%IZ -])  
     %ex(k) = ex(k) +1*(hy(k - 1) - hy(k)); \RFA?PuY  
   % ex(k) = ex(k) +0.25*(hy(k - 1) - hy(k)); /;21?o  
    % ex(k) = ex(k) +1.1*(hy(k - 1) - hy(k)); )fS6H<*  
    ex(k) = ex(k) +0.5*(hy(k - 1) - hy(k)); QAZs1;lU  
   ex(1)=0; <+b~E,  
   ex(201)=0; By}ZHK94I  
end 1L\r:mx3  
$`0,N_C<}  
if T<=100 %.\+j,G7  
pulse = exp(-0.5*(power((t0 - T) / spread, 2))); L6!Hv{ijn  
            ex(kc) = pulse; {cdrMP@""  
else Ip8:~Fl]  
    pulse = exp(-0.5*(power((t0 - T) / spread, 2))); @j%@Z  
    ex(kc) = ex(kc)+pulse; $~ d6KFT  
end +x4*T  
            
zqI|VH  
               pXh`o20I  
for k = 1:KE AAB_Ytf  
   JEZ0O&_R  
            %hy(k) = hy(k) + 0.5*(ex(k) - ex(k + 1)); J^~J&  
          %hy(k) = hy(k) + 1*(ex(k) - ex(k + 1)); uz=9L<$  
          % hy(k) = hy(k) + 0.25*(ex(k) - ex(k + 1)); [E2".F3  
            %  hy(k) = hy(k) + 1.1*(ex(k) - ex(k + 1)); w&]$!g4  
& .. P(C5@x(Z  
Jk{v(W#  

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

感觉仿真出来有数值色散，不知道哪里问题