FDTD方法中的完美匹配层（PML）研究（Matlab代码实现）

 👨‍🎓个人主页：研学社的博客   

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

---

💥1 概述

在FDTD方法中，我们截断了问题空间。但是在截断时，我们面临着在其边界上的反射问题。吸收边界条件（ABC）但很难在FDTD方法中制备二维ABC和使用。
这个问题由PML解决。假设一个波向外传播，它最终会到达允许空间的边缘，这是由矩阵在程序中的尺寸决定的。此概率由此 PML 解决。

📚2 运行结果

运行视频：

FDTD方法中的完美匹配层（PML）(Matlab代码实现）

https://www.bilibili.com/video/BV1sh41157rf/

 FDTD方法中的完美匹配层（PML）(Matlab代码实现）_哔哩哔哩_bilibili

部分代码：

%%%%%_______________PML(Perfectly Matched Layer)_using_FDTD__________%%%%%%%
%_________________________SOUGATA_CHATTERJEE________________________________
%%%%%%%%%_________________SAMEER_KOLKATA_CENTER___________________%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%________02/03/2011_________%%%%%%%%%%%%%
ie=60;
je=60;
ic=fix(ie/2);
jc=fix(je/2);
epsz=8.8e-12;
ddx=0.1;
dt=ddx/6e8;
%%%%%%-----Time_interval-----------%%%%%%%%%%%%
nsteps=100;
%%%%%%-----Time_interval-----------%%%%%%%%%%%%
ez=zeros(ie,je);
dz=zeros(ie,je);
hy=zeros(ie,je);
ihx=zeros(ie,je);
ihy=zeros(ie,je);
hx=zeros(ie,je);
ga=ones(ie,je);
%%%%%%-----wave_specification-----------%%%%%%%%%%%%
t0=40.0;
spread=15.0;
%------------calculated the PML parameter--------------- 

🎉3 参考文献

部分理论来源于网络，如有侵权请联系删除。

[1]尤双双,李孝坤.基于PML的FDTD法正演模拟[J].中国锰业,2018,36(06):191-194.DOI:10.14101/j.cnki.issn.1002-4336.2018.06.045.

[2]Sougata Chatterjee (2023). Perfectly Matched Layer(PML) in FDTD method. 

[3]李志. 基于FDTD的高阶PML算法及其在集成电路仿真中的应用[D].天津工业大学,2021.DOI:10.27357/d.cnki.gtgyu.2021.000394.

🌈4 Matlab代码实现