有限差分法分析电磁场边值问题定稿

西南科技大学本科生毕业论文 1第 1 章 绪 论1.1 电磁场理论产生的背景及其意义电磁场理论是人类探索自然活动的结晶和宝贵财富。人类认识电磁场运动规律的道路是漫长而曲折的。早在两千多年前，人类就有了关于磁石和摩擦起电的知识，我们祖先发明的指南针，为人类文明作出了不朽的贡献。但是，将电磁场现象系统地上升为理论的研究并加以应用则是 18 世纪中叶，特别是 19 世纪中叶以后的事情。17711773 年，卡文迪许（Henry Cavendish;1731_1810）进行了著名的静电实验，库伦（Chareles-Augustinde Coulomb，17361806）于 1785 年建立了关于静电和静磁的平方反比定律，这标志着电学和磁学定量研究的开始。此后，人们对电和磁现象进行了大量的观察和实验研究，其中，最著名的是伽伐尼（L.Calvani，17371798）在解剖青蛙是注意到青蛙腿的痉挛现象，从而发现电流；伏特（Alessandro Volt，17451827）用电化学方法产生了稳定的电流（即伏特电池） 。随后，欧姆（georg Simon Ohm，17891854）和基尔霍夫（Gustav Robert Kirchhoff，1824 1887）分别建立了后来用他们名字命名的电路定律。在很长的时期内，人们把电和磁看成是相互独立的现象，并不知道他们之间有什么联系。直到 1820 年奥斯特（Hans Christian Oersted，17771851）发现电流可使磁针偏转，级电流可产生磁力，才开始了将点与磁联系起来的研究。1825 年，安培（Andrc Maric Ampere，17751836）提出了确定两电流之间相互作用及载流导体能受到磁力作用的定律，即安培定律，毕奥（Biot）和萨法尔（Savart）确定了磁场和电流之间的定量关系，即毕奥-萨法尔定律。到此为止，人们一直都还是在静止的或恒定的状态下研究电磁现象。电磁学研究的一个重大发展史，1831 年法拉第（Michael Faraday，1791 1867）发现电磁感应现象，这是人们第一次对随时间变化的电磁场进行研究。电磁感应定律一方面推动了电磁在工程中的应用，另一方面它是电磁理论的一块基石。1864 年，麦克斯韦（James Clerk Maxwell，18311879）在总结前人发现的实验定律的基础上，进行了创造性的理西南科技大学本科生毕业论文 2论研究工作，建立了后来以他名字命名的麦克斯韦方程组，从而创立了完整的电磁理论体系 1。麦克斯韦电磁理论体系的建立，是 19 世纪人类文明史上的重大事件，它标志着人类文明迈进了电的时代。紧随其后，1866 年，西门子（William Siemens，18231883）发明了发电机；1876 年，贝尔（Alexander Graham Bell，18471922）发明了电话；1879 年，爱迪生（Thomas Alva Edison，18471931）成功地做了电磁波实验，对麦克斯韦方程组的正确性提供两块实验依据。赫兹实验后不到 6 年，意大利工程师马可尼（G.Marconi，18741937）和俄国的波波夫（A.S.Popov，18591906）分别实现了无线电远距离传播，并很快投入实际应用。其后，无线电报（1894 年） 、无线电广播（1906 年） 、导航（1911 年） 、微波通信（1933 年） 、雷达（1935 年）以及近代的无线电遥测、遥控、卫星通信、光纤通信等如雨后春笋般涌现出来。一个多世纪以来，由电磁学发展起来的现代电子技术已应用在电力工程、电子工程、通信工程、计算机技术等多学科领域。电磁理论已广泛应用于国防、工业、农业、医疗、卫生等领域，并深入到人们的日常生活中。今天，电磁场问题的研究及其成果的广泛运用，已成为人类社会现代化的标志之一。1.2 电磁场边值问题计算方法的重要性在一个电磁系统中，电场和磁场的计算对于完成该系统的有效设计师极端重要的。例如，在系统中，用一种绝缘材料是导体相互隔离是，就要保证电场强度低于绝缘介质的击穿强度。在磁力开关中，所要求的磁场强弱，应能产生足够大的力来驱动开关。在发射系统中进行天线的有效设计时，关于天线周围介质中电磁场分布的知识显然有实质性的意义。为了分析电磁场，我们可以从问题所涉及的数学公式入手。依据电磁系统的特性，拉普拉斯方程和泊松方程只能适合于描述静态和准静态（低频）运行条件下的情况。但是，在高频应用中，则必须在时域或频域中求解波动方程，以做到准确地西南科技大学本科生毕业论文 3预测电场和磁场，在任何情况下，满足边界条件的一个或多个偏微分方程的解，因此，计算电池系统内部和周围的电场和磁场都是必要的。电磁场理论早期主要运用在军事领域，其发展和无线电通信、雷达的发展史分不开的。现在，电磁场理论的应用已经遍及得学、生命科学和医学、材料科学和信息科学等几乎所有的科学技术领域。计算电磁场研究的内容涉及面很广，与电磁场工程、电磁场理论互相联系，互相依赖。对电磁场工程而言，计算电磁场要解决的是实际电磁场工程中越来越复杂的建模与仿真、优化设计等问题；而电磁场工程也为之提供实验结果，以验证其计算结果的正确性。对电磁场理论而言，计算电磁场可以为其研究提供进行复杂的数值及解析运算的方法、手段和计算结果；而电磁场理论这为计算电磁场问题提供了电磁规律、数学方程，进而验证其计算结果，计算电磁场对电磁场理论发展的影响不仅仅是提供一个计算工具而已 2，而是使整个电磁场理论发生了革命性的变化。毫不夸张地说，近二三十年来，电磁场理论的发展，无一不是与计算电磁场的发展相联系的。目前，计算电磁场已成为对复杂体系的电磁规律、电磁性质进行研究的重要手段，为电磁场理论的深入研究开辟了新的途径，并极大地推动了电磁场工程的发展。西南科技大学本科生毕业论文 4第 2 章 电磁场边值问题计算方法常用的计算电磁场边值问题的方法主要有两大类，其每一类又包含若干种方法，第一类是解析法；第二类是数值法。对于那些具有最简单的边界条件和几何形状规则的（如矩形、圆形等）问题，可用分离变量法和镜像法求电磁场边值问题的解析解（精确解） ，但是在许多实际问题中往往由于边界条件过于复杂而无法求得解析解。在这种情况下，一般借助于数值法求解电磁场的数值解。2.1 解析法电磁学是一门古老而又不断发展的学科。经典的数学分析方法是近百年来电磁学学科发展中的一个极为重要的方法。解析法包括建立和求解偏微分方程或积分方程。严格求解偏微分方程的经典方法是分离变量法 ,即在可分离变量的坐标系中求解 Maxwell 方程组或其退化形式 ,最后得到解析解。严格的求解积分方程的方法主要是变换数学法.例 1 一 无 限 长 直 接 地 金 属 槽 ,其 三 壁 电 势 为 零 ,顶盖与三壁绝缘且电势为V0sin x，其中 V0=100 V ,截面长宽分别为 a = 10 cm 和 b = 5 cm ,如图 2-1 所示 . a求金属槽内的电势 分布。图 2-1 金属槽截面西南科技大学本科生毕业论文 5分析 金属槽无限长 , 故槽内电势与坐标 z无关. 由于槽内各点上电荷密度 = 0 , 槽内电势满足二维直角坐标系中的拉普拉斯方程及其边界条件 : xaVyxyxbasin,0,0 2（2.1.1）（2.1.2）（2.1.3）（2.1.4）（2.1.5）应 用分离 变 量 法 , 得 到 满 足 方 程 ( 1 ) 和 边 界 条 件 式 (2)式 (4) 的解的形式为（2.1.6）aynshxAyxn1i,带入边界条件（5）得V0sin x= （2.1.7）aabnshxn1i比较系数得：， （2.1.8）abshA011n槽内电势的解析解为 10sin210sin,0 yhxayhxbVyx 解析法的优点是：可将解答表示为已知函数的显式，从而计算出精确的数值结果；可以作为近似解和数值解的检验标准；西南科技大学本科生毕业论文 6在解析过程中和在解的显式中可以观察到问题的内在联系和各个参数对数值结果所起的作用。但解析法也存在严重缺点，主要是，它仅能解决很少量的问题。事实上，只有在为数不多的坐标系中才能分离变量，而用积分方程法是往往求不出结果，致使分析过程既困难又复杂。例如，对标量赫姆霍兹方程，只有在 11 种坐标系下才能用分离变量法求解。如果边界面不是在 11 种坐标系中 1 个坐标系的 1 个坐标面或该坐标系的几个坐标面的组合，或者边界条件不是第一类边界条件（该标量在边界上的值为已知）或不是第二类边界条件（该标量在边界上沿法线方向的空间导数为已知） ，则分立变量就不能用。又如，只有当积分方程中的核是某些形式时，才能用变换数学法来严格求解。2.2 数值法数值法用高性能的计算机就可直接以数值的、程序的形式代替解析形式来描述电磁场的问题。在数值法中，通常以差分代替微分，用有限求和代替积分，这样，就将问题化为求解差分方程或代数方程问题。这方面的例子如有限差分法（FDM） 。数值法与解析法比较，在许多方面具有独特的优点。普适性强，用户拥有的弹性大。一个特定问题的边界条件、电气结构、激励等特性可以不编入基本程序，而由用户输入，更好的情况是通过图形界面输入。用户不必具备高度专业化的电磁场理论、数学及数值技术方面的知识就能用提供的程序解决实际问题。数值法的出现，使电磁场边值问题的分析研究从解析的经典方法进入到离散系统的数值分析方法，从而使许多解析法很难解决的复杂电磁场边值问题，有可能通过电磁场的计算机辅助分析获得很高精度的离散解（数值解） ，同时可极大地促进各种电磁场数值计算方法的发展。数值法的缺点是数据输入量大、计算量大、受硬件条件的限制。原则上，数值法可以求解具有任何复杂几何形状、复杂材料的电磁场工程问题。但是，在工程应西南科技大学本科生毕业论文 7用中，由于受计算机存储容量、执行时间以及解的数值误差等方面的限制，数值法在解大型复杂的电磁场工程问题时也难以完成任务。可以说数值法的发展大致分为两个阶段。起发展初期，是研究“解决得了”的问题，也就是研究该数值法能否应用于各个学科分支领域；而其发展后期，是研究“解决得好”的阶段，即探讨解决工程实际问题的各种改进方法、手段及相应的计算技术。近期的数值法研究中的大量工作都是为了实现这一目标。有的研究在小机器上计算大问题；有的研究减少内存占用，加快计算速度；还有的研究在一定程度上减少自由度和计算工作量；而最新的发展动向是研究高效的并行数值算法。西南科技大学本科生毕业论文 8第 3 章 有限差分法微分方程和积分微分方程数值解的方法。基本思想是把连续的定解区域用有限个离散点构成的网格来代替， 这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似，积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组 ，解此方程组就可以得到原问题在离散点上的近似解 3。然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。有限差分法的主要内容包括 4：如何根据问题的特点将定解区域作网格剖分；如何把原微分方程离散化为差分方程组以及如何解此代数方程组。此外为了保证计算过程的可行和计算结果的正确，还需从理论上分析差分方程组的性态，包括解的唯一性、存在性和差分格式的相容性、收敛性和稳定性。对于一个微分方程建立的各种差分格式，为了有实用意义，一个基本要求是它们能够任意逼近微分方程，这就是相容性要求。另外，一个差分格式是否有用，最终要看差分方程的精确解能否任意逼近微分方程的解，这就是收敛性的概念。此外，还有一个重要的概念必须考虑，即差分格式的稳定性。因为差分格式的计算过程是逐层推进的，在计算第 n1层的近似值时要用到第 n 层的近似值，直到与初始值有关。前面各层若有舍入误差，必然影响到后面各层的值，如果误差的影响越来越大，以致差分格式的精确解的面貌完全被掩盖，这种格式是不稳定的，相反如果误差的传播是可以控制的，就认为格式是稳定的。只有在这种情形，差分格式在实际计算中的近似解才可能任意逼近差分方程的精确解。关于差分格式的构造一般有以下 3 种方法。最常用的方法是数值微分法，比如用差商代替微商等。另一方法叫积分插值法，因为在实际问题中得西南科技大学本科生毕业论文 9出的微分方程常常反映物理上的某种守恒原理，一般可以通过积分形式来表示。此外还可以用待定系数法构造一些精度较高的差分格式。3.1 差分运算的基本概念有限差分法是指用差分来近似取代微分，从而将微分方程离散成为差分方程组。于是求解边值问题即转换成为求解矩阵方程 5。对单元函数 而言，取变量 的一个增量 = ，则函数 的增量可以表xfxxhxf示为= - (3.1.1)xff称为函数 的差分或一阶差分。函数增量还经常表示为xf= - (3.1.2)xf2hxf称为函数 的中心差分或一阶中心差分。xf函数一阶差分 与自变量增量 h 的比值 / 称为一阶差商。在一阶差分xfh运算中，它常用来近似函数 的一阶导数 。xf d函数 的二阶差商定义为xf2xfhxff(3.1.3) 2ff它常被用来近似函数 的二阶导数 。xf2/dxf我们还可以采用类似方法给出二阶以上差分的定义，并用它们来近似函数二阶以上的导数。但由于二阶以上的倒数在本次研究中没有用到，因此就不在赘述了。西南科技大学本科生毕业论文 103.2 有限差分法应用有限差分法基本思想是将场域离散为许多小网格，应用差分原理，将求解连续函数 的泊松方程的问题转换为求解网格节点上 的差分方程组的问题 67。现在，以静电场边值问题Fyx22)(sfL （3.2.1）（3.2.2）为例，说明有限差分法的应用。f(s)为边界点 s 的点函数，二位场域 D 和边界L 示于图 3.2-1 中。0yx02413DLhh图 3.2-1 有限差分的网格分割通常将场域分成足够小的正方形网格,网格线之间的距离为 ,节点 上的h4,3210电位分别用 和 表示。3210,4设函数 在 处可微，则沿 方向在 处的泰勒公式展开为xx0（3.2.3）nK nK0 00)(!将 和 分别代入式（3.2.3）,得13西南科技大学本科生毕业论文 11（3.2.4）0302001 )(!1)(!1)( xhxhx（3.2.5）0302003 )(!)(!)(由（3.2.4）-（3.2.5）得（3.2.6）hx2)(310（3.2.4）+（3.2.5）得（3.2.7）2301x2h0)(同理（3.2.8）hy2)(310（3.2.9）23012)(0y将式(3.2.7)、 （3.2.9）代入式（3.2.1） ，得到泊松方程的五点差分格式)(44 24321020321 Fh Fh当场域中 得到拉普拉斯方程的五点差分格式,（3.2.10）)(404 43210321 从这个公式我们可以看出，当我们将一个二维无源区场域剖分为一系列正方形网格时，场域内任何一个节点的电位都等于它周围四个节点电位的算术平均值。这就是规则正方形网格内某点的电位所满足的拉普拉斯方程的差分格式，或差分方程8。对于场域内的每一个结点，关系式（3.2.10）式都成立，都可以列出一个相同形式的差分方程。但是，对于近邻边界的结点，其边界不一定正好落在正方形网格的结点上，而西南科技大学本科生毕业论文 12可能如图 3.2-2 所示。其中 1、2 为边界线上的结点， p、 q 为小于 1 的正数。仿上所述，可推得对这些近邻边界结点的拉普拉斯方程的差分格式为 0)1(1)()1( 432 qpqp（3.2.11）式中： 1 和 2分别是给定边界条件函数 f (s)在对应边界点处的值，是已知的。0 123 4hhqhph图 3.2-2 近邻边界的结点3.3 边界条件的离散化处理4321h1h20图 3.3-1 边界条件的离散化处理若场域离散为矩形网格(如图 3.3-1 示)，差分格式为 9：(3.3.1)F2h1h1h1 04222 )(（1）第一类边界条件:给边界离散节点直接赋已知电位值（2）对称边界条件:合理减小计算场域，差分格式为：西南科技大学本科生毕业论文 13(3.3.2)(Fh241240图3.3-2 边界条件的离散化处理（3）第二类边界条件:边界线与网格线相重合的差分格式：(3.3.3)hffhn2102010,)(4)介质分界面衔接条件的差分格式(3.3.4) )1(44320 K其中 ba1444 A 4西南科技大学本科生毕业论文 14第 4 章 差分方程组的求解由上一章分析可以看出，对于场域内 D 的每一个节点，就有一个差分方程。场域内部节点的个数就等于差分方程的个数。若节点位于场域的边界 10，则这些边界节点的电位值由边界条件式给出。在对场域 D 内各个结点（包括所有场域内点和有关的边界结点）逐一列出对应的差分方程，组成差分方程组后，就可选择一定的代数解法，以算出各离散结点上待求的电位值。计算差分方程可以直接求解，也可以采用迭代法，相对而言，采用迭代法求解差分方程更受人们重视。这是因为差分方程组的系数一般是有规律的，且各个方程都很简单，包含的项数不多（最多不超过 5 项） ，因此，对于有限差分法，通常都采用逐次近似的迭代方法求解。4.1 简单迭代法例如如图 4.1-1 为一个正方形截面的无限长金属盒。盒子的两侧及底的电位为零，顶部电位为 1000V，求盒内的电位分布 11。000V10西南科技大学本科生毕业论文 15（1）在盒内取 33 个离散的电位节点由于场域内不存在电荷，其电位分布必满足拉普拉斯方程。在均匀剖分的条件下，其差分格式如式（3.2.10）所示。用简单迭代法求解 12，其步骤如下：第一步，在场域内部节点上选定电位初始值，为简单起见，可将它们都取为零，记为 = = =0，常称为零次解。01209第二步，将零次解代入差分方程式（3.2.10） ，得出诸内部节点电位值的一次解，它们为：= = =25014040201= = =25012051034= = =2501340620= = =01407015 4同理可得= = = = =015617819在求出一次解的 9 个节点电位值以后，原来零次解中的 9 个节点电位值将被一次解中的相应电位值所取代，在计算机的内存中不予保留，从而达到了节省存储空间的目的。第三步，重复上述步骤，令每一个内部节点上的二次解电位值等于该节点周围四个相邻节点（或边界点）一次解电位值的算术平均值，并用二次解电位值冲去内存中的原一次解电位值这样迭代一次又一次的继续下去，可望诸节点的电位值变化越来越小，这时可图 4.1-1 用有限差分法求金属盒内电位（）西南科技大学本科生毕业论文 16取这些节点上的电位值为该边值问题的数值解。迭代 39 次后的结果为=428V， =526V， =428V39139239=187V， =250V， =187V456=71V， =98V， =71V39739839利用 MATLAB 绘制等电位图 13，程序见附录，画出等电位图如下：图 4.1-2 简单迭代法（33）等电位分布图由得出的数值解可以看出，金属盒内点电位分布是越靠中间电位越高，越靠近金属盒顶部电位越高，这是由于金属盒底部和两边的电位都为零，而顶部最高。由此表明此方法计算出的电位值，符合金属盒内的电位分布情况。（2）在盒内取 55 个离散的电位节点上述例子是在金属盒内取 33 个离散的电位节点，若是取 55 个节点，如图4.1-3 所示：西南科技大学本科生毕业论文 17图 4.1-3 用有限差分法求金属盒内电位（55）计算步骤同（1） ，经过 91 次迭代之后，得到电位分布如下：0 100.0 1000.0 10000 10000 10000 00 468.7 629.2 669.4 629.2 468.7 00 245.5 378.8 419.3 378.8 245.5 00 134.6 221.2 250.0 221.2 134.6 00 71.8 121.2 138.4 121.2 71.8 00 31.3 53.5 61.3 53.5 31.3 00 0 0 0 0 0 0由以上矩阵得知，在盒内取 55 个节点，求出结果比取 33 个节点更具有连续性，精确度也越高，但是迭代次数也越多，计算所花的时间久越长。4.2 超松弛迭代法在迭代法的应用中，为加速迭代解收敛速度，一般采用的是超松弛迭代法 14。4.2.1 超松弛迭代法理论由于编写计算机程序的需要，每一网格结点的位置由双下标（ i， j）予以识别，如图 4.2-1 所示。对于差分方程（3.2.10）式，采用超松弛迭代法（规定迭代的运算顺序是：从左下角开始做起，即 i 小的先做；对固定的 i， j 小的先做。 ） ，则关于结点 0 迭代到第（ n+1）次时的近似值，应由如下迭代公式算得西南科技大学本科生毕业论文 18)(,)1(,)1(,)(1,)(,1)(,)1(, 44 njinjinjinjinjinjinji ) （4.2.1） ij2(i,j+1)0(i,j)3(i-1,j) 1(i+1,j)4(i,j-1)图 4.2-1 结点的双下标（ i， j）标号式中： 称为加速收敛因子，其取值范围是 1 2，当 2 时，迭代过程将不收敛。加速收敛因子 有一个最佳取值问题，但随具体问题而异。对于第一类边值问题，若一正方形场域由正方形网格分割（每边结点数为 m+1） ，则最佳收敛因子 0可按下式计算 sin120（4.2.2）在更一般的情况下， 只能凭借经验取值。0在超松弛迭代法的应用中，还必须涉及迭代解收敛程度的检验问题。对此，通常的处理方法是：迭代一直进行到所有内结点上相邻两次迭代解的近似值满足修正条件 Wnjiji),()1(,（4.2.3）时，终止迭代。即将式(12)作为检查迭代解收敛程度的依据。其中： W 是指定的最大允许误差。4.2.2 超松弛迭代法的 MATLAB 实现（1）程序框图西南科技大学本科生毕业论文 19启 动 给 定 边 值填 写 场 域 内 的 初 始 值叠 代 次 数 计 数 n:=0n:=n+1按 超 松 弛 迭 代 法 进 行一 次 迭 代 ， 求 )1(,nji所 有 内 点 相 邻 二 次 迭 代 值的 相 对 误 差 是 否 小 于 W？ 打 印 迭 代 次 数 n， 待 求数 值 解 )j,i(停 机 图 4.2-2 程序框图（2）MATLAB 程序代码如下 15：clearhx=5;hy=5; %设置网格节点数（33）v1=ones(hy,hx); %设置行列二维数组v1(1,:)=ones(1,hx)*1000; %上下两行的 Dirichlt 边界值v1(hy,:)=zeros(1,hx); %左右两列的 Dirichlt 边界值for i=1:hyv1(i,1)=0;v1(i,hx)=0;endv2=v1;maxw=1;w=0;n=0；a=1西南科技大学本科生毕业论文 20while(maxw1e-9) %由 V1 迭代，算出 V2，跌代精度为0.000000001n=n+1; %计算迭代次数maxw=0;for i=2:hy-1for j=2:hx-12(i,j)=v1(i,j)+(v1(i,j+1)+v1(i+1,j)+v2(i-1,j)+v2(i,j-1)-4*(v1(i,j)*a/4;%超松弛迭代公式w=abs(v2(i,j)-v1(i,j);if(wmaxw)maxw=w;endv(i,j)=v2(i,j)endendv1=v2;endn %输出迭代次数contour(v2,20) %画等电位线(3)计算结果如下：n=19v = 0 1000.0 1000.0 1000.0 00 428.6 526.8 428.6 00 187.5 250.0 187.5 00 71.4 98.2 71.4 00 0 0 0 0画出等电位分布图：西南科技大学本科生毕业论文 21图 4.2-3 超松弛迭代法（33）等电位分布图当取 55 个离散的电位节点时，只需将程序 hx=5;hy=5 改为 hx=7;hy=7，经过 54 次迭代之后计算出结果为：0 1000.0 1000.0 1000.0 1000.0 1000.0 00 468.7 629.2 669.4 629.2 468.7 00 245.5 378.8 419.3 378.8 245.5 00 134.6 221.2 250.0 221.2 134.6 00 71.8 121.2 138.4 121.2 71.8 00 31.3 53.5 61.3 53.5 31.3 00 0 0 0 0 0 0画出等电位图：西南科技大学本科生毕业论文 22图 4.2-4 超松弛迭代法（55）等电位分布图4.3 讨论与分析由（33）和（55）等电位分布图对比可知，两种取点方法求得等电位分布规律一致，但取的离散的点越多，求出的电位连续性越强。迭代收敛的速度与 有明显关系表1-1 迭代收敛的速度与 的关系收敛因子（ ） 1.0 1.1 1.15 1.2 1.3 1.5 1.8 2.0（33）迭代次数（ ）N39 29 23 19 24 41 122 发散当收敛因子 =1.2时迭代次数最少，计算所用时间最短。通过简单迭代法和超松弛迭代法对比发现，两种方法求出的数值解相同，作出的等电位线分布一样，但超松弛迭代法比简单迭代法收敛速度更快，迭代次数更少，计算时间更短。西南科技大学本科生毕业论文 23结 论电磁场边值问题的数值求解有很多种方法，有限差分法是数值计算中应用得最早而又相当简单，直观的一种方法 16。本文在讨论电磁场边值问题求解时，将电磁场进行了理想化处理，以一简单边界条件的电磁场边值问题为例，建立相应的数学模型，将场域离散为一些网格，运用差分原理，对场域内偏微分方程及边界上的边界条件进行差分离散化处理，在通过差分运算求出场域内电位值。可以通过上述分析得到这样一些有意义的结论：（1）使用有限差分法求解电磁场边值问题是可行的，只要将网格取得足够小，是可以将离散的点看成连续的。求出离散的数值解，更符合实际应用的需要，而且随着计算机技术的发展，求解差分方程的过程变得简单，使得有限差分法在电磁场问题计算中更具有优势。（2）超松弛迭代法的收敛因子 =1 时，相当于简单迭代法。收敛因子 的导入，使得超松弛迭代法比简单迭代法能更快的收敛，若 的值选取恰当，收敛速度还将加快。（3）场域内取的节点越多，计算结果就越精确，当节点划分足够多的时候，离散的点可以看作连续的。但节点划分越多，迭代次数就越多，计算量就越大，所以计算时应根据实际需要，划分合适的节点数。在电磁场计算中如何保证计算精度，减少计算工作量，提高计算速度，减少计算时间，是我们长期努力的方向。西南科技大学本科生毕业论文 24参 考 文 献1 冯慈璋 , 马西奎. 工程电磁场导论 M . 北京 : 高等教育出版社 , 2000 :32240 .2 徐立勤 曹伟电磁场与电磁波 M 北京 : 科学出版社,2006:102-107.3 何红雨. 电磁场数值计算法与 MATLAB 实现 M . 武汉 : 华中科技大学出版社 , 2004 : 4210 .4 余定峰 李 超. 有限差分法在静态电磁场计算中的应用 J .电子测试 2009(4),23-26.5 赵得奎 刘勇.MATLAB 在有限差分数值计算中的应用J.四川理工学院报,2005,18(4):61-64.6 Dlala , Emadl emad. Inverted and Forward Preisach Models for Numerical Analysis of Electromagnetic Field Problems J . IEEE Transactions on Magnetics , 2006 , 42 (8) :1 963 21 973.7 Kanai , Yasushi. Automatic mesh generation for 3D electro2 magnetic field analysis by FD - TD met hod J . IEEE Trans2 actions on Magnetics ,1998 Part 1 of 2 ,34 (5) :3 383.8 谢处方,饶克谨电磁场与电磁波 M 北京 : 高等教育出版社, 1999:96-100.9 钱焕延,计算方法 M 西安电子科技大学出版社,2007: 100-110.10 雷亚平,肖洪祥 匡晚成基于 MATLAB 的电磁场数值分析J .电子测试 2007,10