电磁波简单一维FDTD的MATLAB实现

第一次写博客，发表自己的学习成果，如有不对之处，希望各位评论指出。
内容如题，公式为直角坐标系中FDTD一维公式。

代码

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%标题：一维FDTD MATLAB模拟
%版本：MATLAB R2018b Update4 
%时间：2019.7.3 13:47:30
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clc;clear all;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Z向TEM波
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
kzmax=200;                                    %模拟空间长度
tmax=200;                                     %模拟时间长度
ex=zeros(200,200);                            %声明电场，电场矢量为x向
hy=zeros(200,200);                            %声明磁场，磁场矢量为y向
epsz=8.854187817e-12;                         %真空介电常数
muz=(4*pi)*10^(-7);                           %真空磁导率
for t=2:tmax-1                                %时间步进
    ex(1,t)=exp(-(t-1-30)*(t-1-30)/(100));    %高斯脉冲函数在z为1处