fdtd 磁场 matlab,矩形谐振腔电磁场的FDTD分析和Matlab仿真

矩形谐振腔电磁场的FDTD分析和Matlab仿真

摘要：目前，电磁场的时域计算方法越来越引人注目。这种方法已经广泛应用到各种电磁问题的分析之中。而将Matlab作为一种仿真工具，用于时域有限差分法，可以简化编程，使研究者重心放在FDTD本身上，而不必在编程上花费过多的时间。本课题通过用FDTD方法计算矩形谐振腔电磁场分布，并用Matlab 进行仿真。

关键词：时域有限差分法，Matlab仿真，矩形谐振腔

1.引言

时域有限差分法(Finite-Dfference Time-Domain Method)是求解电磁问题的一种数值技术，是在1966年由K.S.Yee第一次提出的。FDTD法直接将有限差分式代替麦克斯韦时域场旋度方程中的微分式，得到关于场分量的有限差分式，用具有相同电参量的空间网格去模拟被研究体，选取合适的场始值和计算空间的边界条件，可以得到包括时间变量的麦克斯韦方程的四维数值解，通过傅里叶变换可求得三维空间的频域解。时域有限差分法突出的优点是所需的存储量及计算时间与N成正比，使得很多复杂的电磁场计算问题成为可能，用时域有限差分法容易模拟各种复杂的结构，使得用其他方法不能解决的问题有了新的处理方法。

本文主要讨论如何用Matlab语言来编写FDTD的吸收边界条件以及编程时应注意的问题。

2时域有限差分法的基本理论

2.1 时域有限差分法的简介

1966年K.S.Yee首次提出了一种电磁场数值计算的新方法——时域有限差分(Finite-Dfference Time-Domain Method)方法。对电磁场E、H分量在时间和空间上采取交替抽样的离散方式，每一个E(或H)场分量四周有四个H(或E)场分量环绕，应用这种离散方式将含时间变量的麦克斯韦旋度方程转化为一组差分方程，并在时间轴上逐步推进地求解空间电磁场。Yee提出的这种抽样方式后来被称为Yee元胞。

FDTD方法是求解麦克斯韦方程的直接时域方法。在计算中将空间某一样本点的电场(或磁场)与周围格点的磁场(或电场)直接相关联，且介质参数已赋值给空间每一个元胞，因此这一方法可以处理复杂形状目标和非均匀介质物体的电磁散射、辐射等问题。同时FDTD的随时间推进可以方便地给出电磁场的时间演化过程，在计算机上以伪彩色方式显示，这种电磁场可视化结果清楚的显示了物理过程，便于分析和设计。

2.2 FDTD数值计算的优势

FDTD算法，其空间节点采用Yee元胞的方法，电场和磁场节点空间与时间上都采用交错抽样，因而使得麦克斯韦旋度方程离散后构成显示差分方程，相比较宇前面的波动方程求解，计算等到大大简化。由于FDTD采用吸收边界条件的