为什么我的fdtd--mur边界条件没有用

clc ,AEaW  
clear (hFyp}jkk  
5tQZf'pHfd  
%定义参数 VFLW@  
                                                 %1.光及平板的相关性质 {'$+?V"&  
c=3e8;                                           %光速 w-wap  
mu=(4e-7)*pi;                                    %磁导率 n,8bQP=&  
epslon=8.84e-12;                                 %介电常数 )H|cri~D  
d=0.18;                                          %平板的x方向厚度   Gx/kel[Y}  
                                                 %2.传播过程的相关数据 @z1pE@7jK  
dz=0.015;                                        %Z轴方向最小间距 o#=@!m  
dt=dz/(2*c);                                     %变换时间的最小间隔 nX|]JW  
nz=400;                                          %Z轴方向上的中步数 $v0beN6MG  
nt=1500;                                         %T轴方向上的步数 _fE$KaP  
                                                 %3.输入函数的性质 =@y ?Np^A  
t=0.5e-9;                                        %高斯函数的分割周期 X&s"}Hf  
t0=3*t;                                          %3S置信度的3倍分割周期 {?hpW+1,#  
\zM3{{mV/  
cm=dt/(mu*dz);                                   %4.列式中的相关数据 $s"{C"4q  
ce=dt/(epslon*dz); 8c^Hfjr0  
                                     AOscewQ  
%初始化 9U]j@*QN  
ex(1:nz+1)=0;                                    %构建电场强度初始矩阵 x#c%+  
hy(1:nz)=0;                                      %构建磁强度初始矩阵 L wP  
exearly=zeros(1,nz+1); ;?.w!|6  
hyearly=zeros(1,nz); qEajT"?  
%循环 ~x6<A\  
for i=1:nt                                       %时间循环为nt /;clxtus  
   if i<=1000                                    %时间步长为1000 pUD(5v*0R  
       ex(1)=exp(-(i*dt-t0)^2/t^2);              %输入信号 ibn(eu<uW  
   end M" R=;n  
   for j=1:nz                                    %磁场及电场的变换关系（麦克斯韦变换关系） J`xCd/G  
      hy(j)=hy(j)+cm*(ex(j)-ex(j+1)); \k>1q/T0V  
   end ViYfK7Z  
   for k=2:nz Ze"m;T  
      ex(k)=ex(k)+ce*(hy(k-1)-hy(k));
Z+g9!
@'a  
+\)a p  

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