[讨论]我这个二维FDTD的matlab程序为何会出现这个问题？

我这个二维FDTD的matlab程序为何会出现这个问题？中心放置一个高斯信号源Ez，长度为100dt，随着波的传播，外围是一个圆，这没错，可是为什么中间还有这个矩形？ f)N67z6  
)e]:T4*vo  
function fdtdex1 :n>:*e@w%  
c=3.0e8; 9Q].cDe[  
sigmac=0; YQe @C  
sigmam=0; )pJ}o&J  
epslon0=8.854e-12; b@5&<V;r2  
epslonr=1; d|DIqT~{W  
epslon=epslon0*epslonr; "eGS~-DVK  
mu0=1.2566306e-6; b\H(Lq17  
mur=1; %3"3OOT7  
mu=mu0*mur; QT^( oog=  
f=1e9; -hhE`Y  
lambda=c/f;
;41s&~eR  
Nlambda=5; ]:]2f9y  
dx=lambda/Nlambda; pmHd1 Wub  
S=0.7; %
l3RM*zb  
dt=dx/c/sqrt(3); p|Z"< I7p(  
CA=(1-sigmac*dt/2/epslon)/(1+sigmac*dt/2/epslon); 3/]J i^+  
CB=(dt/epslon/dx)/(1+sigmac*dt/2/epslon); !YAkHrF`[0  
DA=(1-sigmam*dt/2/mu)/(1+sigmam*dt/2/mu); 2c9@n9Vx3a  
DB=(dt/mu/dx)/(1+sigmam*dt/2/mu); +Y 3_)  
Ez=zeros(500,500);
uyAhN  
Hy=zeros(500,500); r|0C G^:C  
Hx=zeros(500,500); EL~s90C  
Centerx=ceil(500/2); c|ZZ+2IYd  
Centery=ceil(500/2); I%{U~  
for t=1:1000 6o6!Ol  
    Ez(2:499,2:499)=CA*Ez(2:499,2:499)+CB*(Hy(3:500,2:499)-Hy(2:499,2:499)+Hx(2:499,2:499)-Hx(2:499,3:500));  '6 w|z^  
    Ez(1,1:500)=0; N9{ivq|fO  
    Ez(1:500,1)=0; K\n %&w  
    Ez(500,1:500)=0; }L:LcM  
    Ez(1:500,500)=0; $0cMrf@  
    if (t<100) zz02F+H$Y  
        Ez(Centerx,Centery)=exp(-(t*dt-50*dt)^2*f^2); 2NE/ZqREg  
    end j>0S3P,  
    Hx(2:500,2:500)=DA*Hx(2:500,2:500)+DB*(Ez(2:500,1:499)-Ez(2:500,2:500)); +$8hTi,  
    Hy(2:500,2:500)=DA*Hy(2:500,2:500)+DB*(Ez(2:500,2:500)-Ez(1:499,2:500)); Na3tK}x  
    Emax=max(max(Ez)); '_(oa<g  
    Emin=min(min(Ez));
C|z`hNp  
    imagesc(Ez,[0 1]); -c|O!Lc-  
    title('Ez comp .. |%zhwDQ.  
AGrGZ7p]  

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

if (t<100) D0"+E*   
        Ez(Centerx,Centery)=exp(-(t*dt-50*dt)^2*f^2); pRt )B`#  
    end :()4eK/\  
>1joCG~  
有问题。 时间不够常。高斯信号还没有完全衰减。

向高手学习 71)HxC[6vA  
请问：CB=(dt/epslon/dx)/(1+sigmac*dt/2/epslon); 中的dt/epslon/dx为什么要除以dx？？ MhJ`>.z1  
请不吝赐教，谢谢！

准确的说应该是 dy. 楼主用的是均匀网哥。所以没有关系。

知道了，谢谢！