时域有限差分法仿真一维TE波在分裂场完全匹配层【含源码】.docx

时域有限差分法仿真一维TE波在分裂场完全匹配层吸收边界条件下的传输一、时域有限差分法 (FDTD, Finite-Difference Time-Domain)FDTD是1966年K.S.Yee发表在AP上的一篇论文建立起来的，后被称为Yee网格空间离散方式核心思想是把带时间变量的Maxwell旋度方程转化为差分形式，模拟出电子脉冲和理想导体作用的时域响应号称目前计算电磁学界最受关注，最时髦的算法，但还在发展完善之中国外已有多种基于FDTD算法的电磁场计算的软件:XFDTD等。二、差分算法 三、时域有限差分算法步骤(1)采用一定的网格划分方式离散化场域；(2)对场内的偏微分方程及各种边界条件进行差分离散化处理，建立差分格式，得到差分方程组；(3)结合选定的代数方程组的解法，编制程序，求边值问题的数值解。四、吸收边界条件1、问题的提出在电磁场的辐射和散射问题中，边界总是开放的，电磁场占据无限大空间，而计算机内存是有限的，所以只能模拟有限空间。即：时域有限差分网格将在某处被截断。这要求在网格截断处不能引起波的明显反射，因而对向外传播的波而言，就像在无限大的空间传播一样，一种行之有效的方法是在截断处设置一种吸收边界条件。使传播到截断出的波被边界吸收而不产生反射。2、良匹配层基本原理1994年Beernger提出一种新颖的由吸收媒质构成的良匹配层(PML)的概念,这种人工设计的良匹配层由有耗,导电,导磁媒质组成,可吸收任意入射角,任意频率,任意偏振态的入射电磁波,且透射波以原来波速传播,且有相同的特征阻抗,在垂直入射方向衰减。根据推导，只要这种媒质的结构参数满足一定条件,这种媒质就可以起到完全吸收入射波的作用。五、计算所用到的电磁场公式根据麦克斯韦方程组和时域差分法，得到以下公式：六、仿真结果七、源代码%***********************************************************************% 时域有限差分法仿真一维TE波在分裂场完全匹配层吸收边界条件下的传输%***********************************************************************%%%Eref Ein Erans% PEC PML介质 PML PEC% +-----+---+------+---------------------+-----+------+--% z=0 z=Lz% r=1 Nref NinNtrans r=Nz+1%%***********************************************************************% 物理常量mu0 = 4e-7 * pi; % 磁导率c0 = 299792458; % 光速eps0 = 1/c0^2/mu0; % 电导率eta0 = sqrt(mu0/eps0); % 波导率 GHz = 10^9;mm = 10^(-3);%***********************************************************************% 参数初始化%***********************************************************************Lz = 1; % 长度单位：米Tmax = 4*Lz/c0; % 时间最大值Dz = 1*mm; % 空间网格尺寸R = 1/sqrt(3); %每一单位时间计算的网格单元数freq = 5*GHz; % 激励源的中心频率fmax = 9*GHz; % 画图的最大频率fmin = 1*GHz;Nz = round(Lz/Dz); % z轴上的单元格数量（round 四舍五入）n_pict = round(Nz/20); % 每n_pict步画一次图linew = 2; % 画图的线宽Dt = Dz*R/c0; % 时间步长Nt = round(Tmax/Dt); %