2维FDTD+摩尔吸收边界，遇金属块反射

2维FDTD全波数值算法+摩尔吸收边界，基于matlab模拟二维平面波的空间传递，并且设置金属块模拟波的反射过程

clear all;

u=4.*pi.*10.^（-7）;
e=8.854187817.*10.^（-12）;
cig=0;
c=3*10^8;

Imax=200;
Jmax=200;
Kmax=200;
Nmax=800;
d=0.001;
dt=d/2/c;

s=（c*dt-d）/（c*dt+d）;                               %边界计算公式系数（相邻）
s1=（c*dt-sqrt（2）*d）/（c*dt+sqrt（2）*d）;              %边界计算公式系数（斜对角）

Ca=1;
Cb=dt/e/d;
Da=1;
Db=dt/u/d;

Ex=zeros（ImaxJmax）;
Ey=zeros（ImaxJmax）;
Ez=zeros（ImaxJmax）;
Hx=zeros（ImaxJmax）;
Hy=zeros（ImaxJmax）;
Hz=zeros（ImaxJmax）;

A=zeros（ImaxJmax）;                                 %三组记录上一时刻的数组
B=zeros（ImaxJmax）;
C=zeros（ImaxJmax）;
for n=1:Nmax
    
    
     A=Ex;                                          %记录上一时刻Ex    
    for x=1:60                                      %计算Ex（1-60） 
        for y=2:Jmax
         Ex（xy）=Ca*Ex（xy）+Cb*（Hz（xy）-Hz（xy-1））;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       3481  2017-04-26 14:57  FDTD2withM.m

----------- ---------  ---------- -----  ----

                 3481                    1