电磁场 计算机辅助分析 有限差分.

电磁场的数学建模与解答技巧电磁场是电荷和电流所产生的相互作用效应，它在工程学、物理学以及计算机模拟中都扮演着重要角色。

为了更好地理解和分析电磁场，数学建模和解答技巧是必不可少的。

本文将从电磁场的数学建模入手，介绍几种常用的数学建模方法，并给出解答技巧的实例。

一、电磁场的数学建模方法之一：微分方程微分方程是描述电磁场的一种常用数学工具。

通常，通过麦克斯韦方程组可以得到电磁场满足的偏微分方程。

对于静电场，可以使用拉普拉斯方程描述，表示为：∇²ϕ = -ρ/ε₀其中ϕ是电势，ρ是电荷密度，ε₀是真空介电常数。

对于静磁场，则可以使用斯托克斯方程描述，表示为：∇×B = μ₀J其中B是磁感应强度，J是电流密度，μ₀是真空磁导率。

通过求解这些微分方程，可以得到电磁场的分布情况。

二、电磁场的数学建模方法之二：有限元法有限元法是一种常用的数值解法，可用于求解任意形状的电磁场问题。

该方法将电磁场区域划分为有限个小单元，并在每个小单元内以多项式函数逼近电磁场的分布。

通过建立离散的代数方程组，并求解该方程组，可以得到电磁场的近似解。

三、电磁场的数学建模方法之三：有限差分法有限差分法是一种离散方法，通过将连续的电磁场问题转化为离散的代数问题进行求解。

该方法将连续的电磁场区域划分为网格，并在每个网格节点上进行逼近。

通过近似微分算子，将偏微分方程转化为差分方程，并通过迭代求解差分方程得到电磁场的解。

四、电磁场解答技巧实例为了更好地展示电磁场解答技巧，以下给出一个实例。

考虑一个带有一根无限长直导线的无限大平面问题。

已知导线的电流密度为I，求解该情况下的磁场分布。

根据安培环路定理，可以得到这个问题的微分方程为：∇×B = μ₀Iδ(x)δ(y)ez其中δ表示狄拉克δ函数，ez表示z轴方向上的单位向量。

通过对微分方程进行求解，可以得到在导线周围的磁场强度为：B = μ₀I/2πr其中r表示距导线的径向距离。

电磁场与电磁波实验有限差分法作者: 日期:电磁场与电磁波实验报告实验项目：有限差分法一、实验目的及要求1学习有限差分法的原理与计算步骤；2、学习用有限差分法解静电场中简单的二维静电场边值问题；3、学习用Matlab语言描述电磁场与电磁波中内容，用matlab求解问题并用图形表示出了，学习matlab语言在电磁波与电磁场中的编程思路。

二、实验内容理论学习：学习静电场中边值问题的数值法中的优先差分法的求解知识；实践学习：学习用matlab语言编写有限差分法计算二维静电场边值问题；三、实验仪器或软件Matlab7.0电脑四、实验原理有限差分法的基本思想将计算场域划分成网格，把求解场域内连续的场分布用求解网格节点上的离散数值解来代替；即用网格节点的差分方程近似代替场域内的偏微分方程来求解。

简单迭代法小(°)先对场域内的节点赋予初始值㈡，这里上标(0)表示第°次近似值，即初始值。

然后再按照：VUi]进行反复迭代。

若当第N次迭代结束后，所有内节点相邻两次迭代值之间的绝对误差小于事先给定的精度，则迭代停止。

MAX①:N)- ①:N‘)W初始值的赋予是任意的；赋予初始值后，请按“从左到右、从下到上”的固定顺序依次计算各节点值; 当所有节点都算完一遍后，再用它们的新值代替旧值，即完成一次迭代。

五、实验步骤复习理论知识；编写matlab程序；六、结果分析与问题讨论1、程序:clearX=[0,0,0,0,0;0,25,25,25,0;0,50,50,50,0;0,75,75,75,0;100,100,100,100,100]Pot=[0,0];for i=2:4for j=2:4(i ,Pptx(1 ；j2,=(X(!-.1)j)+xe k1)+X3+1)2X0+1))4'Pot(1)=abs(PotX(i-1,j-1)-X(i,j));'''Pot(2)=max(Pot)endendX(2:4,2:4)=PotXnum=1;while(max(1000.*Pot)>1) Pot(2)=0;for i=2:4for j=2:4声PotX(i-1,j-1)=(X(i-1,j)+X(i,j-1)+X(i+1,j)+X(i,j+1))/4Pot(1)=abs(PotX(i-1,j-1)-X(i,j));Pot(2)=max(Pot)endendX(2:4,2:4)=PotXnum=nu m+1endsurf([0:4],[0:4],X);shadi ng in terpcolorbar('horiz')title(' 有限差分法计算电位图');2、运行结果X =0 0 0 0 00 25 25 25 00 50 50 50 00 75 75 75 0100 100 100 10C 1 100%第一次迭代PotX =18.7500Pot =6.2500 6.2500PotX =7.1440 9.8230 7.144018.7515 25.0023 18.751542.8583 52.6801 42.8583Pot =0.3815 0.7629%第28次迭代X =0 0 0 0 00 7.1440 9.8230 7.14400 100.0000 100.0000 100.0000 100.0000 100.0000 num =283、波形图matlab 软件在使用有限差分法研究静电场边值问题中有着重要的作用，它能够快捷有效 并且准确的解决边值问题，是解决计算相对复杂问题的有效工具。

电磁波时域有限差分方法电磁波时域有限差分方法是一种在计算电磁波传播过程中广泛使用的数值模拟方法。

它通过将电磁场的时域偏导数转化为差分形式进行离散计算，从而得到电磁场的时域响应。

这种方法在电磁波仿真、电磁辐射、雷达散射以及通信系统设计等领域具有重要的应用价值。

时域有限差分方法的理论基础是电磁波的麦克斯韦方程组。

通过将麦克斯韦方程组进行离散化，将时域偏导数转化为差分形式，并使用合适的差分格式来近似电场和磁场的时域分布。

通过迭代计算离散化后的麦克斯韦方程组，可以得到电磁场在时域上的演化过程。

具体来说，时域有限差分方法的基本步骤如下：1. 网格划分：首先对仿真区域进行网格划分，将空间离散为有限的小单元。

典型的网格划分包括一维、二维和三维的情况。

2. 差分格式选择：根据实际问题选择合适的差分格式，如中心差分格式、向前差分格式或向后差分格式等。

差分格式的选择会直接影响计算结果的准确性和稳定性。

3. 时间步长确定：为了保证计算结果的稳定性，需要根据空间离散步长和电磁波传播速度来确定合适的时间步长。

时间步长的选择需要满足稳定性条件。

4. 初始条件和边界条件设定：在仿真开始前，需要设定初始条件和边界条件。

初始条件指定电磁场在仿真区域内的初始分布，而边界条件则决定了电磁场与仿真区域边界的相互作用关系。

5. 迭代求解：通过迭代计算离散化的麦克斯韦方程组，可以得到电场和磁场在时域上的演化过程。

每一次迭代都涉及更新电场和磁场的数值。

时域有限差分方法相比其他电磁波计算方法具有一定的优势。

首先，它能够模拟电磁场的时域响应，对于短脉冲信号或非稳态过程的仿真非常有用。

其次，它在空域和频域上的计算误差相对较小，并且可以处理各种不规则形状的仿真区域。

此外，时域有限差分方法还可以结合其他方法，如有限元方法和边界元方法，进行更精确的仿真计算。

虽然时域有限差分方法在电磁波仿真中取得了显著的成果，但它也存在一些局限性。

首先，它的计算速度相对较慢，特别是在三维仿真中。

计算电磁场理论中的有限差分法与有限元法电磁场理论是电磁学的重要组成部分，研究电磁场的分布和变化规律对于解决实际问题具有重要意义。

在计算电磁场中，有限差分法和有限元法是两种常用的数值计算方法。

本文将从理论原理、应用范围和优缺点等方面对这两种方法进行探讨。

有限差分法是一种将连续问题离散化的方法，通过将连续的电磁场分割成网格，然后在每个网格上进行离散计算。

这种方法的基本思想是将微分方程转化为差分方程，然后利用差分方程进行求解。

有限差分法的优点是简单易懂，计算过程直观，适用于各种电磁场问题的求解。

然而，由于差分法中的网格离散化会引入一定的误差，所以在计算精度上存在一定的限制。

与有限差分法相比，有限元法是一种更加精确的数值计算方法。

有限元法将电磁场问题的求解区域划分为有限个小单元，然后在每个小单元上建立适当的插值函数，通过求解代数方程组得到电磁场的近似解。

有限元法的优点是可以处理复杂的几何形状和材料特性，适用于各种边界条件和非线性问题。

然而，有限元法的计算过程相对较为复杂，需要对问题进行合理的离散化和网格划分，同时对于大规模问题，计算量也较大。

在实际应用中，根据具体问题的特点和求解要求，选择合适的数值计算方法是十分重要的。

对于简单的电磁场问题，如一维导线的电流分布，可以选择有限差分法进行求解。

而对于复杂的电磁场问题，如三维空间中的电磁波传播，有限元法更适合。

此外，有限差分法和有限元法还可以结合使用，通过将两种方法的优点相结合，提高计算精度和效率。

除了理论原理和应用范围，有限差分法和有限元法的优缺点也值得关注。

有限差分法的优点是简单易懂，计算过程直观，而且对于一些简单问题可以得到较为准确的结果。

然而，由于差分法中的网格离散化会引入一定的误差，对于复杂问题的求解精度有限。

相比之下，有限元法可以处理复杂的几何形状和材料特性，适用于各种边界条件和非线性问题，计算精度较高。

然而，有限元法的计算过程相对复杂，需要对问题进行合理的离散化和网格划分，同时对于大规模问题计算量较大。

西南科技大学本科生毕业论文第 1 章绪论1.1电磁场理论产生的背景及其意义电磁场理论是人类探索自然活动的结晶和宝贵财富。

人类认识电磁场运动规律的道路是漫长而曲折的。

早在两千多年前，人类就有了关于磁石和摩擦起电的知识，我们祖先发明的指南针，为人类文明作出了不朽的贡献。

但是，将电磁场现象系统地上升为理论的研究并加以应用则是18世纪中叶，特别是19世纪中叶以后的事情。

1771——1773年，卡文迪许（Henry Cavendish;1731_1810）进行了著名的静电实验，库伦（Chareles-Augustinde Coulomb，1736——1806）于1785年建立了关于静电和静磁的平方反比定律，这标志着电学和磁学定量研究的开始。

此后，人们对电和磁现象进行了大量的观察和实验研究，其中，最著名的是伽伐尼（L.Calvani，1737——1798）在解剖青蛙是注意到青蛙腿的痉挛现象，从而发现电流；伏特（Alessandro Volt，1745——1827）用电化学方法产生了稳定的电流（即伏特电池）。

随后，欧姆（georg Simon Ohm，1789——1854）和基尔霍夫（Gustav Robert Kirchhoff，1824——1887）分别建立了后来用他们名字命名的电路定律。

在很长的时期内，人们把电和磁看成是相互独立的现象，并不知道他们之间有什么联系。

直到1820年奥斯特（Hans Christian Oersted，1777——1851）发现电流可使磁针偏转，级电流可产生磁力，才开始了将点与磁联系起来的研究。

1825年，安培（Andrc Maric Ampere，1775——1836）提出了确定两电流之间相互作用及载流导体能受到磁力作用的定律，即安培定律，毕奥（Biot）和萨法尔（Savart）确定了磁场和电流之间的定量关系，即毕奥-萨法尔定律。

到此为止，人们一直都还是在静止的或恒定的状态下研究电磁现象。

Research Institute of RF & Wireless Techniques School of Electronic and Information EngineeringSouth China University of Technology褚庆昕华南理工大学电子与信息学院计算电磁场第7讲频域有限差分法(2)-分界面的差分网格Research Institute of RF & Wireless Techniques 分界面与网格线重合 分界面与网格线呈对角线 三角形分界面微带问题第7讲内容7分界面与网格线重合Research Institute of RF & Wireless Techniques先假设为Research Institute of RF & Wireless Techniques7-2 分界面对网格呈对角线Research Institute of RF & Wireless Techniques7Research Institute of RF & Wireless Techniques7Research Institute of RF & Wireless TechniquesResearch Institute of RF & Wireless TechniquesResearch Institute of RF & Wireless TechniquesResearch Institute of RF & Wireless TechniquesResearch Institute of RF & Wireless Techniques 在差分法中难度最大的就是边界（包括分界）的差分处理。

如果有边界条件依据，则可得到令人信服的格式。

有些情况下，边界条件不清楚（例如角域中尖点处的边界条件），便有较大的随意性。

实验一 用有限差分法解静电场边值问题一、目的1．掌握有限差分法的原理与计算步骤； 2．理解并掌握求解差分方程组的超松弛迭代法，分析加速收敛因子α的作用； 3．学会用有限差分法解简单的二维静电场边值问题，并编制计算程序。

二、方法原理有限差分法是数值计算中应用得最早而又相当简单、直观的一种方法。

应用有限差分法通常所采取的步骤是：⑴ 采用一定的网格分割方式离散化场域。

⑵ 进行差分离散化处理。

用离散的、只含有限个未知数的差分方程组，来近似代替场域内具有连续变量的偏微分方程以及边界上的边界条件（也包括场域内不同媒质分界面上的衔接条件）。

⑶ 结合选定的代数方程组的解法，编制计算机程序，求解由上面所得对应于待求边值问题的差分方程组，所得解答即为该边值问题的数值解。

现在，以静电场边值问题⎪⎩⎪⎨⎧==∂∂+∂∂)2()()1(02222s f D y x Lϕϕϕ中在为例，说明有限差分法的应用。

f (s )为边界点s 的点函数，二位场域D 和边界L 示于图5.1-1中。

x图5.1-1 有限差分的网格分割1． 离散化场域应用有限差分法时，首先需从网格划分着手决定离散点的分布方式。

通常采用完全有规律的方式，这样在每个离散点上可得出相同形式的差分方程，有效地提高解题速度。

如图5.1-1所示，现采用分别与x ，y 轴平行的等距（步距为h ）网格线把场域D 分割成足够多的正方形网格。

各个正方形的顶点（也即网格线的交点）称为网格的结点。

这样，对于场域内典型的内结点0，它与周围相邻的结点1、2、3和4构成一个所谓对称的星形。

2．差分格式造好网格后，需把上述静电场边值问题中的拉普拉斯方程（1）式离散化。

设结点0上的电位值为ϕ0。

结点1、2、3和4上的电位值相应为ϕ1、ϕ2、ϕ3和ϕ4，则基于差分原理的应用，拉普拉斯方程（1）式在结点0处可近似表达为ϕ1+ϕ2+ϕ3+ϕ4-4ϕ1=0 （3）这就是规则正方形网格内某点的电位所满足的拉普拉斯方程的差分格式，或差分方程。

电磁场计算方法与技术的研究电磁场是物质界中的一个现象，在物理学的研究中占有极为重要的地位。

随着科技的不断发展，电磁场计算方法和技术也得到了越来越广泛的应用。

本文将探讨电磁场计算的方法和技术。

一、电磁场的基本性质电磁场是一种能量和动量的传递媒介，它是由电荷和电流产生的。

它包括电场和磁场，它们之间的关系描述了电磁场本身的性质。

电场和磁场的产生和变化都受到麦克斯韦方程组的描述。

在实际研究中，我们需要根据麦克斯韦方程组建立数学模型，通过数值计算的方法来得到电磁场的分布规律。

二、电磁场计算的方法1. 有限差分法有限差分法是一种常用的电磁场计算方法，它基于离散化的思想，将空间和时间离散化为等间距的点，通过有限差分近似代替微分运算。

这种方法可以用于求解各种边值问题，包括静电场、静磁场和电磁波等。

有限差分法具有高效、简单、精度较高等特点，但是其精确度受到空间剖分的影响，如果空间间距过大，会导致精度不足。

2. 有限元法有限元法也是一种常用的电磁场计算方法，它是一种离散化方法，将连续的物理场离散为有限个离散元素。

通过有限元法可以得到各个单元内部的电磁场，从而推导出连续介质内部电磁场的分布规律。

3. 边界元法边界元法是一种将边界上的物理量加以积分来求解内部物理场的方法。

它将物理场分为积分区域内、积分区域外两部分，利用格林公式和边界条件将积分区域外的物理量转化为边界上的物理量，从而求解内部物理场。

边界元法具有高效、高精度、适应性强等优点。

三、电磁场计算的技术1. 计算机辅助设计软件计算机辅助设计软件可以用来模拟、计算电磁场分布规律，可以方便地建立模型，进行计算和分析，从而进行电磁场的优化设计。

2. 多物理场仿真软件多物理场仿真软件可以模拟多种物理场之间的耦合和相互作用，包括电场、磁场、温度场等。

它们不仅可以对各种物理场进行计算和分析，还可以对各种现象进行多方位的仿真，提高了计算精度和准确性。

3. 电磁场测量和检测仪器电磁场测量和检测仪器可以用来检测和测量电磁场的分布规律。

用有限差分方法求解微波电磁场问题本章主要内容是说明用差分法求解在微波器件和微波技术中常常遇见的一些偏微分方程的边值问题。

我们知道，很多给定边界条件的偏微分方程的求解相当复杂。

除少数情况外，要求它的精确解是颇为困难的，一般采用近似方法。

有限差分法就是经常采用的一种近似方法，它是用离散的、含有有限个未知数的差分方程去替代连续变量的微分方程，并把相应的差分方程的解作为该边值问题数值形式的近似解。

1 用差分方程解拉普拉斯方程在微波系统中很多问题，例如同轴线的台阶电容、谐振腔隙缝处的漏散电容、微带线的特性阻抗等，要求出它们的值，首先就要找出这些线或谐振腔内静电电位分布，这些电位分布是满足拉普拉斯方程的。

用差分方法解拉普拉斯方程是很方便的，所以我们开始就讨论它。

将拉普拉斯方程化成差分方程的方法在很多书上都可找到[6, 7]，下面将列出公式而不作推导，仅对差分方程的求解过程作一些简单介绍。

一、基本差分公式我们要求的电位函数u ，它在区域D 内满足下面的拉普拉斯方程02222=∂∂+∂∂yux u (1-1) 在边界上S ，它服从以下条件：()p f u S = (1-2)式中()p f 为边界点p 的函数。

这类问题一般称为第一类边值问题或称狄里赫利问题。

为了用差分方法求解电位分布，先在y x -平面分别作两族平行于x 轴和y 轴的直线,线间的距离为h ，于是各直线的x 和y 坐标分别为：jh y ih x j i == ;式中j i ,为正整数，取值1、2、……。

这样区域D 就被许多边长为h 的正方形所覆盖，在图1-1中示出了这种情况。

各正方形的顶点被称为网格的节点，从图可以看到，各节点所处位置有所不同。

一些节点（例如a 节点）恰落在边界上S ，我们把它叫做边界节点。

有些节点到边界的距离不足h （例如节点b ）,这些节点叫做不规则节点。

但是大部分节点到边界的距离大于h ，例如图上的0点，它们属于规则节点。

差分法就是求这些离散节点处u 的近似值。

电磁场的数值计算方法与应用引言：电磁场是物理学中一个重要的研究领域，它涉及到电磁波、电磁感应等多个方面。

为了更好地理解和应用电磁场，科学家们开发了各种数值计算方法。

本文将介绍电磁场的数值计算方法及其应用。

一、有限差分法有限差分法是一种常用的数值计算方法，它将连续的电磁场问题离散化为离散的网格点问题。

通过在网格点上近似计算电场和磁场的导数，可以得到电场和磁场在空间中的分布情况。

有限差分法的优点是简单易懂，适用于各种电磁场问题的求解。

例如，可以利用有限差分法计算电磁波在介质中的传播，或者计算导体中的电磁感应现象。

二、有限元法有限元法是一种广泛应用于工程领域的数值计算方法，它可以用于求解各种复杂的电磁场问题。

有限元法将电磁场问题离散化为一系列的小区域，称为有限元。

通过在每个有限元上近似计算电场和磁场的分布，可以得到整个电磁场的数值解。

有限元法的优点是适用于各种不规则形状的区域，可以处理复杂的边界条件和材料特性。

例如，可以利用有限元法分析电磁场在电机中的分布，或者计算电磁屏蔽结构的性能。

三、边界元法边界元法是一种特殊的数值计算方法，它将电磁场问题转化为在边界上求解的问题。

边界元法通过在边界上近似计算电场和磁场的分布，可以得到整个电磁场的数值解。

边界元法的优点是可以减少计算的自由度，提高计算效率。

例如，可以利用边界元法计算电磁波在散射体上的散射现象，或者计算导体表面的电磁场分布。

四、数值计算方法在电磁场问题中的应用数值计算方法在电磁场问题中有着广泛的应用。

例如，在通信领域中，可以利用数值计算方法分析电磁波在天线和传输线中的传播特性，以及在无线通信系统中的传播损耗和干扰现象。

在电力系统中，可以利用数值计算方法分析电磁场对输电线路和变压器的影响，以及计算电力设备的电磁兼容性。

在电子设备设计中，可以利用数值计算方法分析电磁场对电路元件的耦合和干扰，以及计算电磁屏蔽结构的性能。

总之，数值计算方法在电磁场问题的研究和应用中发挥着重要的作用。

第一章 有限差分法一元函数泰勒公式：设函数()f x 在0x 处的某邻域内具有1n +阶导数，则对该邻域异于0x 的任意点x ，在0x 与x 之间至少存在一点ξ，使得()20000000()()()()()()()()()2!!n n n f x f x f x f x f x x x x x x x R x n '''=+-+-+⋅⋅⋅+-+其中，(1)10()()()(1)!n n n f R x x x n ξ++=-+ 二元函数的泰勒公式：设函数(,)z f x y =在点00(,)x y 的某一邻域内连续且有直到1n +阶连续偏导数，00(,)x h y k ++为此邻域内任意点，则有00000020000100(,)(,)()(,)1()(,)2!1()(,)!1()(,)(1)!n n f x h y k f x y h k f x y x yh k f x y x y h k f x y n x y h k f x h y k n x yθθ+∂∂++=++∂∂∂∂+++⋅⋅⋅∂∂∂∂++∂∂∂∂+++++∂∂式中01θ<<；0000(,)0(,)mm m p p m p m x y p m p p f h k f x y c h kxy x y --=⎛⎫∂∂∂+= ⎪∂∂∂∂⎝⎭∑1.利用泰勒展开求不等间距的差分格式。

（1）2x y ϕ∂∂∂ （2）33xϕ∂∂解：(1)2600104001040101040010400010041()()2!11()()(,)!(1)!n n h h h h x y x yh h h h x h y h n x y n x yϕϕϕϕϕϕθθ+∂∂∂∂=+-++-++⋅⋅⋅∂∂∂∂∂∂∂∂+-++-+-+∂∂+∂∂（1.1）250020120020120102012002012000200121()()2!11()()(,)!(1)!n n h h h h x y x yh h h h x h y h n x y n x y ϕϕϕϕϕϕθθ+∂∂∂∂=+++++⋅⋅⋅∂∂∂∂∂∂∂∂++++++∂∂+∂∂（1.2）280090110090110109011009011000900111()()2!11()()(,)!(1)!n n h h h h x y x yh h h h x h y h n x y n x y ϕϕϕϕϕϕθθ+∂∂∂∂=+--+--+⋅⋅⋅∂∂∂∂∂∂∂∂+--+----∂∂+∂∂（1.3）270010030010030101003001003000100031()()2!11()()(,)!(1)!n n h h h h x y x yh h h h x h y h n x y n x yϕϕϕϕϕϕθθ+∂∂∂∂=+-+-+⋅⋅⋅∂∂∂∂∂∂∂∂+-+-+-∂∂+∂∂（1.4）()()()()605080702010400104020201200201202090110090110010*********()()2!1()()2!1()()2!1()(2!A B C D A h h A h h x y x y B h h B h h x y x y C h h C h h x y x y D h h D h h x y x ϕϕϕϕϕϕϕϕϕϕϕϕϕϕϕ-+-+-+-∂∂∂∂=-++-+∂∂∂∂∂∂∂∂++++∂∂∂∂∂∂∂∂+--+--∂∂∂∂∂∂∂∂+-+-∂∂∂∂20)y ϕ （1.5）0102090100412011032222010209010222204012011030000Ah Bh Ch Dh Ah Bh Ch Dh Ah Bh Ch Dh Ah Bh Ch Dh -+-+=⎧⎪+--=⎪⎨+++=⎪⎪+++=⎩ （1.6） 有线性代数知识可知，有多种差分格式。

计算电磁学中的时域有限差分法的数值特性分析及应用摘要时域有限差分法(Finite Difference Time Domain，FDTD)是解决电磁问题非常有效的一种数值方法。

本文先介绍了FDTD的基本原理，分析了FDTD解的稳定性和数值色散分析，然后用FDTD求解电磁散射问题，吸收边界条件的设置起着关键性作用。

通过时间和空间上的递推算法对FDTD中的两种吸收边界条件：Mur吸收边界条件和完全匹配层（PML）的吸收效果进行了比较和分析。

同时，引入参数对PML 的差分方程进行了优化，避免了将电磁场分裂为两个分量进行计算，进而降低了计算内存开销。

实验结果证明PML具有更优越的吸收性能。

关键词：计算电磁学；时域有限差分法（FDTD）；吸收边界条件1.绪论1.1 电磁场数值计算方法概述自1873年麦克斯韦建立电磁场基本方程以来，电磁理论和应用的发展已经有一百多年的历史，Maxwell方程组的提出对于科学技术的发展具有重要的推动作用。

解析法、近似法、数值法共同构成求解Maxwell方程组的主要手段[1]。

在现代电磁场工程中，由于问题的复杂性，要求得到封闭形式的解已不可能，就是半解析的近似方法也只能在个别问题中得到有限的应用，能够较广泛发挥作用的，只有各种数值方法。

随着计算机技术的发展，诞生了一门解决复杂电磁理论和工程问题的应用科学——计算电磁学[2,3]。

最近几十年，各具优势和特色的新颖算法层出不穷相继提出。

在经历了理论和实践两方面检验的基础上，一些有生命力的数值计算方法取得长足进步，应用范围不断拓展。

关于电磁场数值计算方法如图1所示：图1 电磁场数值计算方法分类1.2 FDTD研究背景FDTD是电磁场数值计算中一种有效的方法。

在1966年K.S.Yee发表的著名论文“Numerical solution of initial boundary value problems involving Maxwell’s equation in isotropic Media”中，用后来被称为Yee氏网格的空间离散方式，把带有时间变量Maxwell方程转化为差分方程，诞生了后来被称作FDTD的一种新的电磁场数值解法[4]。