时域有限差分法的Matlab仿真.pdf

时域有限差分法的Matlab仿真

关键词:Matlab矩形波导时域有限差分法

摘要：介绍了时域有限差分法的基本原理，并利用Matlab仿真，对矩形波导谐振
腔中的电磁场作了模拟和分析。

关键词：时域有限差分法；Matlab；矩形波导；谐振腔

目前，电磁场的时域计算方法越来越引人注目。时域有限差分（Finite
DifferenceTimeDomain，FDTD）法［1］作为一种主要的电磁场时域计算方法，
最早是在1966年由K.S.Yee提出的。这种方法通过将Maxwell旋度方程转化为
有限差分式而直接在时域求解，通过建立时间离散的递进序列，在相互交织的网
格空间中交替计算电场和磁场。经过三十多年的发展，这种方法已经广泛应用到
各种电磁问题的分析之中。

Matlab作为一种工程仿真工具得到了广泛应用［2］。用于时域有限差分法，可以
简化编程，使研究者的研究重心放在FDTD法本身上，而不必在编程上花费过多
的时间。
下面将采用FDTD法，利用Matlab仿真来分析矩形波导谐振腔的电磁场，说明
了将二者结合起来的优越性。

1FDTD法基本原理

时域有限差分法的主要思想是把Maxwell方程在空间、时间上离散化，用差分
方程代替一阶偏微分方程，求解差分方程组，从而得出各网格单元的场值。FDTD
空间网格单元上电场和磁场各分量的分布如图1所示。

电场和磁场被交叉放置，电场分量位于网格单元每条棱的中心，磁场分量位
于网格单元每个面的中心，每个磁场（电场）分量都有4个电场（磁场）分量环
绕。这样不仅保证了介质分界面上切向场分量的连续性条件得到自然满足，而且
还允许旋度方程在空间上进行中心差分运算，同时也满足了法拉第电磁感应定律
和安培环路积分定律，也可以很恰当地模拟电磁波的实际传播过程。
1.1Maxwell方程的差分形式

旋度方程为：

将其标量化，并将问题空间沿3个轴向分成若干网格单元，用Δx，Δy和Δz
分别表示每个网格单元沿3个轴向的长度，用Δt表示时间步长。网格单元顶点
的坐标（x，y，z）可记为：

其中：i，j，k和n为整数。
同时利用二阶精度的中心有限差分式来表示函数对空间和时间的偏导数，即
可得到如下FDTD基本差分式：

由于方程式里出现了半个网格和半个时间步，为了便于编程，将上面的差分
式改写成如下形式：

其中：ε和μ分别为介质的介电常数和磁导率。
1.2数值色散及稳定性条件

为了减小数值色散，在选取空间网格尺寸时，应满足λ≥10Δ，
min
Δ=min(Δx,Δy,Δz)，λ是被研究媒质空间的最小波长值。由此可以看出：
min
减小网格尺寸可以减小数值色散，但是会引起计算存储量的增大，因此需综合考
虑，权衡处理。

为了使数值计算稳定，时间步长的选择应满足：

2矩形波导谐振腔模型及仿真分析

矩形波导谐振腔是由两端短路的一段金属波导构成，如图2所示。本文选用的
是WJB100型矩形波导，即宽度a=22.86mm，高度b=10.16mm，选取谐振腔长度

l=50.80mm。由微波理论［3］可知，当b<a<l时，TE模的谐振波长最长，是矩形
101
波导谐振腔的主模。TE模在矩形腔的3个方向都不传输能量，呈驻波分布，并且
101
电场分量只有E分量，在腔体中央最强；磁场有H和H两个分量，在腔壁附近
yxz
最强，腔体中央为0。由于矩形波导谐振腔有着广泛的应用，并且大多工作在主模
状态，这里选择微分高斯脉冲以激励起TE模。
101

决定网格单元的尺寸和时间步长取Δx=Δy=Δz=Δ=1.27mm，则空间网格
数为18×8×40，取Δt=Δ/(2c)2.118ps。

设置激励源微分高斯脉冲的表达式为：

其优点是不含零频率分量，时域波形和频谱如图3所示。为了在谐振腔中激励
起TE模，并且抑制其他高次模，选择线源脉冲，使之在腔内xz平面中心处沿
101
y轴方向分布，并选择合适的t和τ值。经过反复试验，取t=95.316ps，
00
τ=285.948ps。

设置边界条件对于矩形波导谐振腔，腔体的6个面都是金属，为方便起见，
这里假设都是理想导体，即腔内导体边界上的所有切向电场分量为0，所有法向
磁场分量为0。

矩形谐振腔的谐振波长为：

则TE模的谐振波长为λ=41.69mm，对应的谐振频率f=7.19GHz。
10100
图4是Matlab仿真计算出的矩形波导谐振腔的谐振频率，很明显，与理论值非
常接近，由此可以确定确实激励起了TE模，并且较好地抑制了高次模。
101
图5是程序运行N=3000步时，腔体中央平面上的各分量瞬时分布图，可
以看出，符合电磁场理论。

3结语

以上结合FDTD和Matlab对矩形波导谐振