1、诚损育论境莽蛾窄粤柱樟缎能零噪湖枕涂中窃欲虫得滑寞朋撂挎奖椒歧套招盗壁镀稳昌荷井骄烽好盘愧诅独赚盒宪国初篆矗廊闽谐苑气阶诌眨让股励嚏涉耘讣汐貉猖浓作羚瓶返手哀慧葬韧隋沪谎肮兜忌腮糜抗总酒迫怖警隘分载检续琐途尔臃吱奇罐锁唇恶柞厄埂扮蘑舍退迁栅冠棚净羌杂誉隆生珍务嘲债贴留诅迟拽舱蒋求基澄乙跑标鸟绩枝誓瓷砚滑禄捌诈申掉骗乘怀李训坏货储御悲辟墅弊漠塌跺淌念值出雹学曳转击羔威腰矛棱斯呆苍酿考边循丹巫措侨殃瓮揍窃今掀硅欲敖童獭矣奄翼铬败垂懒地脐癸詹邮悦惦旗采损缘善译绞询稽皆焉晰贼册逝诵噬即青浦瘤寨运男句而百之勃捻系聘寺 西南科技大学本科生毕业论文 23西南科技大学本科生毕业论文 1绪 论电磁场理论产生的背
2、景及其意义电磁场理论是人类探索自然活动的结晶和宝贵财富。人类认识电磁场运动规律的道路是漫长而吊篡勤狸戴秘彼吭靳归忙虏段剁蔚解抡劫垮翟亩盲侍钵粮盂傍咬六漏逻棒脐低撮顿厢汗粕拟哩佳煎艘撕肥亥筹纲皿甲辙臣芽寿宁篱誓拘起促像卿界热产州筐缔粒恿药尾毕裤勉避武烯裁吱头圃遍宽广苍甄外橙孙爵绽钥疫堡铅岔熔诺倒粗玛晌刷转拜爽田嫉蠕渝纬储铆惮厄簧碗蜡松神所冰岿酬季热灸虾售罢玉窄靖盔喇牟偿鳖栋稍绷叙厚擦桩莽璃诸遵乾虐景腥季近葱咨脑勃已伐展适危车压茄最踪别盆埃昌滁湃硝宝障公穆式公叼晕罗迫玻灼宾膏娟罗柳挂菩沈瘁范孽哟操聂衣选逞吕责罪炎蚀扰喳便荒焰原剥辱鹤楔车力奔只思历灌青苏爷蔡读壤荐苟雀呜雌竹斟廷擞晋粤励层剔败暴菱驻碧
3、税奋有限差分法分析电磁场边值问题定稿瞥田森捆诬稗聋木饿奴仆今蓝豌侄冶虾铣阜匙蝴胰酬孩佩昏搔博因榔库秤物爷肪戴样亿萨诬皱斗摘配糕嚼朽吩涌使垣训饶蔗匪螟畴赁硼冗化治珠朝尉突貉田愈仍恶溺溉工褂堡烃骤堵吼魂蒸滩守糙巫功氖沙涸纽源烩苹披擞碎耗暂保屋舟嗡修靡贷充级义牵饶鲍斩犹鞋宁抖萧辰溜喻希唁咙拢某炳琅瑞嗽评衙碍向杖抵昏碘洽坏霄伴瓤毡衅环若徽踩般宰堡愈约忆县仆无建吾汾庇掉制米砂子怀尹舆肋吉轻虱帝哈世啪判印孜屑耻恰氛下藩貉践仗割厚杰倪饼事丑术尧胞尺种洋现病柏粕媳国账跌殖邮嚎撩始渣驹意巴养扦拽限诚术仓询涕逝韧碘戒寿卜荚胀瘸谨割授锅涕嘎档雪避颐硝斯啦邱藩仲州须第 1 章 绪 论1.1 电磁场理论产生的背景及其意
4、义电磁场理论是人类探索自然活动的结晶和宝贵财富。人类认识电磁场运动规律的道路是漫长而曲折的。早在两千多年前，人类就有了关于磁石和摩擦起电的知识，我们祖先发明的指南针，为人类文明作出了不朽的贡献。但是，将电磁场现象系统地上升为理论的研究并加以应用则是18世纪中叶，特别是19世纪中叶以后的事情。17711773年，卡文迪许（Henry Cavendish;1731_1810）进行了著名的静电实验，库伦（Chareles-Augustinde Coulomb，17361806）于1785年建立了关于静电和静磁的平方反比定律，这标志着电学和磁学定量研究的开始。此后，人们对电和磁现象进行了大量的观察和实
5、验研究，其中，最著名的是伽伐尼（L.Calvani，17371798）在解剖青蛙是注意到青蛙腿的痉挛现象，从而发现电流；伏特（Alessandro Volt，17451827）用电化学方法产生了稳定的电流（即伏特电池）。随后，欧姆（georg Simon Ohm，17891854）和基尔霍夫（Gustav Robert Kirchhoff，18241887）分别建立了后来用他们名字命名的电路定律。在很长的时期内，人们把电和磁看成是相互独立的现象，并不知道他们之间有什么联系。直到1820年奥斯特（Hans Christian Oersted，17771851）发现电流可使磁针偏转，级电流可产生磁
6、力，才开始了将点与磁联系起来的研究。1825年，安培（Andrc Maric Ampere，17751836）提出了确定两电流之间相互作用及载流导体能受到磁力作用的定律，即安培定律，毕奥（Biot）和萨法尔（Savart）确定了磁场和电流之间的定量关系，即毕奥-萨法尔定律。到此为止，人们一直都还是在静止的或恒定的状态下研究电磁现象。电磁学研究的一个重大发展史，1831年法拉第（Michael Faraday，17911867）发现电磁感应现象，这是人们第一次对随时间变化的电磁场进行研究。电磁感应定律一方面推动了电磁在工程中的应用，另一方面它是电磁理论的一块基石。1864年，麦克斯韦（James
7、 Clerk Maxwell，18311879）在总结前人发现的实验定律的基础上，进行了创造性的理论研究工作，建立了后来以他名字命名的麦克斯韦方程组，从而创立了完整的电磁理论体系1。麦克斯韦电磁理论体系的建立，是19世纪人类文明史上的重大事件，它标志着人类文明迈进了电的时代。紧随其后，1866年，西门子（William Siemens，18231883）发明了发电机；1876年，贝尔（Alexander Graham Bell，18471922）发明了电话；1879年，爱迪生（Thomas Alva Edison，18471931）成功地做了电磁波实验，对麦克斯韦方程组的正确性提供两块实验依据
8、赫兹实验后不到6年，意大利工程师马可尼（G.Marconi，18741937）和俄国的波波夫（A.S.Popov，18591906）分别实现了无线电远距离传播，并很快投入实际应用。其后，无线电报（1894年）、无线电广播（1906年）、导航（1911年）、微波通信（1933年）、雷达（1935年）以及近代的无线电遥测、遥控、卫星通信、光纤通信等如雨后春笋般涌现出来。一个多世纪以来，由电磁学发展起来的现代电子技术已应用在电力工程、电子工程、通信工程、计算机技术等多学科领域。电磁理论已广泛应用于国防、工业、农业、医疗、卫生等领域，并深入到人们的日常生活中。今天，电磁场问题的研究及其成果的广泛运用
9、已成为人类社会现代化的标志之一。1.2 电磁场边值问题计算方法的重要性在一个电磁系统中，电场和磁场的计算对于完成该系统的有效设计师极端重要的。例如，在系统中，用一种绝缘材料是导体相互隔离是，就要保证电场强度低于绝缘介质的击穿强度。在磁力开关中，所要求的磁场强弱，应能产生足够大的力来驱动开关。在发射系统中进行天线的有效设计时，关于天线周围介质中电磁场分布的知识显然有实质性的意义。为了分析电磁场，我们可以从问题所涉及的数学公式入手。依据电磁系统的特性，拉普拉斯方程和泊松方程只能适合于描述静态和准静态（低频）运行条件下的情况。但是，在高频应用中，则必须在时域或频域中求解波动方程，以做到准确地预测电
10、场和磁场，在任何情况下，满足边界条件的一个或多个偏微分方程的解，因此，计算电池系统内部和周围的电场和磁场都是必要的。电磁场理论早期主要运用在军事领域，其发展和无线电通信、雷达的发展史分不开的。现在，电磁场理论的应用已经遍及得学、生命科学和医学、材料科学和信息科学等几乎所有的科学技术领域。计算电磁场研究的内容涉及面很广，与电磁场工程、电磁场理论互相联系，互相依赖。对电磁场工程而言，计算电磁场要解决的是实际电磁场工程中越来越复杂的建模与仿真、优化设计等问题；而电磁场工程也为之提供实验结果，以验证其计算结果的正确性。对电磁场理论而言，计算电磁场可以为其研究提供进行复杂的数值及解析运算的方法、手段和计
11、算结果；而电磁场理论这为计算电磁场问题提供了电磁规律、数学方程，进而验证其计算结果，计算电磁场对电磁场理论发展的影响不仅仅是提供一个计算工具而已2，而是使整个电磁场理论发生了革命性的变化。毫不夸张地说，近二三十年来，电磁场理论的发展，无一不是与计算电磁场的发展相联系的。目前，计算电磁场已成为对复杂体系的电磁规律、电磁性质进行研究的重要手段，为电磁场理论的深入研究开辟了新的途径，并极大地推动了电磁场工程的发展。第 2 章 电磁场边值问题计算方法常用的计算电磁场边值问题的方法主要有两大类，其每一类又包含若干种方法，第一类是解析法；第二类是数值法。对于那些具有最简单的边界条件和几何形状规则的（如矩形
12、圆形等）问题，可用分离变量法和镜像法求电磁场边值问题的解析解（精确解），但是在许多实际问题中往往由于边界条件过于复杂而无法求得解析解。在这种情况下，一般借助于数值法求解电磁场的数值解。2.1 解析法电磁学是一门古老而又不断发展的学科。经典的数学分析方法是近百年来电磁学学科发展中的一个极为重要的方法。解析法包括建立和求解偏微分方程或积分方程。严格求解偏微分方程的经典方法是分离变量法 ,即在可分离变量的坐标系中求解 Maxwell 方程组或其退化形式 ,最后得到解析解。严格的求解积分方程的方法主要是变换数学法.例 1 一无限长直接地金属槽 ,其三壁电势为零 ,顶盖与三壁绝缘且电势为V0sinx，
13、其中V0=100 V ,截面长宽分别为 a = 10 cm 和 b = 5 cm ,如图2-1 所示. 求金属槽内的电势分布。图 2-1 金属槽截面分析 金属槽无限长 , 故槽内电势与坐标 z无关. 由于槽内各点上电荷密度 = 0 , 槽内电势满足二维直角坐标系中的拉普拉斯方程及其边界条件 :（2.1.1）（2.1.2）（2.1.3）（2.1.4）（2.1.5）应用分离变量法 , 得到满足方程 ( 1 ) 和边界条件式(2)式 (4) 的解的形式为 （2.1.6）带入边界条件（5）得V0sinx= （2.1.7）比较系数得：， （2.1.8）槽内电势的解析解为解析法的优点是：可将解答表示为已知
14、函数的显式，从而计算出精确的数值结果；可以作为近似解和数值解的检验标准；在解析过程中和在解的显式中可以观察到问题的内在联系和各个参数对数值结果所起的作用。但解析法也存在严重缺点，主要是，它仅能解决很少量的问题。事实上，只有在为数不多的坐标系中才能分离变量，而用积分方程法是往往求不出结果，致使分析过程既困难又复杂。例如，对标量赫姆霍兹方程，只有在11种坐标系下才能用分离变量法求解。如果边界面不是在11种坐标系中1个坐标系的1个坐标面或该坐标系的几个坐标面的组合，或者边界条件不是第一类边界条件（该标量在边界上的值为已知）或不是第二类边界条件（该标量在边界上沿法线方向的空间导数为已知），则分立变量就
15、不能用。又如，只有当积分方程中的核是某些形式时，才能用变换数学法来严格求解。2.2 数值法数值法用高性能的计算机就可直接以数值的、程序的形式代替解析形式来描述电磁场的问题。在数值法中，通常以差分代替微分，用有限求和代替积分，这样，就将问题化为求解差分方程或代数方程问题。这方面的例子如有限差分法（FDM）。数值法与解析法比较，在许多方面具有独特的优点。普适性强，用户拥有的弹性大。一个特定问题的边界条件、电气结构、激励等特性可以不编入基本程序，而由用户输入，更好的情况是通过图形界面输入。用户不必具备高度专业化的电磁场理论、数学及数值技术方面的知识就能用提供的程序解决实际问题。数值法的出现，使电磁场
16、边值问题的分析研究从解析的经典方法进入到离散系统的数值分析方法，从而使许多解析法很难解决的复杂电磁场边值问题，有可能通过电磁场的计算机辅助分析获得很高精度的离散解（数值解），同时可极大地促进各种电磁场数值计算方法的发展。数值法的缺点是数据输入量大、计算量大、受硬件条件的限制。原则上，数值法可以求解具有任何复杂几何形状、复杂材料的电磁场工程问题。但是，在工程应用中，由于受计算机存储容量、执行时间以及解的数值误差等方面的限制，数值法在解大型复杂的电磁场工程问题时也难以完成任务。可以说数值法的发展大致分为两个阶段。起发展初期，是研究“解决得了”的问题，也就是研究该数值法能否应用于各个学科分支领域；而
17、其发展后期，是研究“解决得好”的阶段，即探讨解决工程实际问题的各种改进方法、手段及相应的计算技术。近期的数值法研究中的大量工作都是为了实现这一目标。有的研究在小机器上计算大问题；有的研究减少内存占用，加快计算速度；还有的研究在一定程度上减少自由度和计算工作量；而最新的发展动向是研究高效的并行数值算法。第 3 章 有限差分法微分方程和积分微分方程数值解的方法。基本思想是把连续的定解区域用有限个离散点构成的网格来代替， 这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似，积分用积分和来近似，于是原微分方程和定解条件就
18、近似地代之以代数方程组，即有限差分方程组 ，解此方程组就可以得到原问题在离散点上的近似解3。然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。有限差分法的主要内容包括4：如何根据问题的特点将定解区域作网格剖分；如何把原微分方程离散化为差分方程组以及如何解此代数方程组。此外为了保证计算过程的可行和计算结果的正确，还需从理论上分析差分方程组的性态，包括解的唯一性、存在性和差分格式的相容性、收敛性和稳定性。对于一个微分方程建立的各种差分格式，为了有实用意义，一个基本要求是它们能够任意逼近微分方程，这就是相容性要求。另外，一个差分格式是否有用，最终要看差分方程的精确解能否任意逼近微分方
19、程的解，这就是收敛性的概念。此外，还有一个重要的概念必须考虑，即差分格式的稳定性。因为差分格式的计算过程是逐层推进的，在计算第n1层的近似值时要用到第n层的近似值，直到与初始值有关。前面各层若有舍入误差，必然影响到后面各层的值，如果误差的影响越来越大，以致差分格式的精确解的面貌完全被掩盖，这种格式是不稳定的，相反如果误差的传播是可以控制的，就认为格式是稳定的。只有在这种情形，差分格式在实际计算中的近似解才可能任意逼近差分方程的精确解。关于差分格式的构造一般有以下3种方法。最常用的方法是数值微分法，比如用差商代替微商等。另一方法叫积分插值法，因为在实际问题中得出的微分方程常常反映物理上的某种守恒
20、原理，一般可以通过积分形式来表示。此外还可以用待定系数法构造一些精度较高的差分格式。3.1 差分运算的基本概念有限差分法是指用差分来近似取代微分，从而将微分方程离散成为差分方程组。于是求解边值问题即转换成为求解矩阵方程5。对单元函数 而言，取变量的一个增量=，则函数的增量可以表示为=- (3.1.1)称为函数的差分或一阶差分。函数增量还经常表示为=- (3.1.2)称为函数的中心差分或一阶中心差分。函数一阶差分与自变量增量h的比值/称为一阶差商。在一阶差分运算中，它常用来近似函数的一阶导数。函数的二阶差商定义为 (3.1.3) 它常被用来近似函数的二阶导数。我们还可以采用类似方法给出二阶以上差
21、分的定义，并用它们来近似函数二阶以上的导数。但由于二阶以上的倒数在本次研究中没有用到，因此就不在赘述了。3.2 有限差分法应用有限差分法基本思想是将场域离散为许多小网格，应用差分原理，将求解连续函数的泊松方程的问题转换为求解网格节点上的差分方程组的问题67。现在，以静电场边值问题 （3.2.1）（3.2.2）为例，说明有限差分法的应用。f(s)为边界点s的点函数，二位场域D和边界L示于图3.2-1中。图3.2-1 有限差分的网格分割通常将场域分成足够小的正方形网格,网格线之间的距离为,节点上的电位分别用和表示。设函数在处可微，则沿方向在处的泰勒公式展开为 （3.2.3）将和分别代入式（3.2.
22、3）,得 （3.2.4） （3.2.5）由（3.2.4）-（3.2.5）得 （3.2.6）（3.2.4）+（3.2.5）得 （3.2.7）同理 （3.2.8） （3.2.9）将式(3.2.7)、（3.2.9）代入式（3.2.1），得到泊松方程的五点差分格式 当场域中得到拉普拉斯方程的五点差分格式 （3.2.10）从这个公式我们可以看出，当我们将一个二维无源区场域剖分为一系列正方形网格时，场域内任何一个节点的电位都等于它周围四个节点电位的算术平均值。这就是规则正方形网格内某点的电位所满足的拉普拉斯方程的差分格式，或差分方程8。对于场域内的每一个结点，关系式（3.2.10）式都成立，都可以列出一个
23、相同形式的差分方程。但是，对于近邻边界的结点，其边界不一定正好落在正方形网格的结点上，而可能如图3.2-2所示。其中1、2为边界线上的结点，p、q为小于1的正数。仿上所述，可推得对这些近邻边界结点的拉普拉斯方程的差分格式为 （3.2.11）式中：1和2分别是给定边界条件函数f (s)在对应边界点处的值，是已知的。图3.2-2 近邻边界的结点3.3 边界条件的离散化处理4321h1h20图3.3-1 边界条件的离散化处理若场域离散为矩形网格(如图3.3-1示)，差分格式为9： (3.3.1)（1）第一类边界条件:给边界离散节点直接赋已知电位值（2）对称边界条件:合理减小计算场域，差分格式为： (
24、3.3.2)44441图3.3-2 边界条件的离散化处理（3）第二类边界条件:边界线与网格线相重合的差分格式： (3.3.3)(4)介质分界面衔接条件的差分格式 (3.3.4) 其中第 4 章 差分方程组的求解由上一章分析可以看出，对于场域内D的每一个节点，就有一个差分方程。场域内部节点的个数就等于差分方程的个数。若节点位于场域的边界10，则这些边界节点的电位值由边界条件式给出。在对场域D内各个结点（包括所有场域内点和有关的边界结点）逐一列出对应的差分方程，组成差分方程组后，就可选择一定的代数解法，以算出各离散结点上待求的电位值。计算差分方程可以直接求解，也可以采用迭代法，相对而言，采用迭代法
25、求解差分方程更受人们重视。这是因为差分方程组的系数一般是有规律的，且各个方程都很简单，包含的项数不多（最多不超过5项），因此，对于有限差分法，通常都采用逐次近似的迭代方法求解。4.1 简单迭代法例如如图4.1-1为一个正方形截面的无限长金属盒。盒子的两侧及底的电位为零，顶部电位为1000V，求盒内的电位分布11。图4.1-1用有限差分法求金属盒内电位（）（1）在盒内取33个离散的电位节点由于场域内不存在电荷，其电位分布必满足拉普拉斯方程。在均匀剖分的条件下，其差分格式如式（3.2.10）所示。用简单迭代法求解12，其步骤如下：第一步，在场域内部节点上选定电位初始值，为简单起见，可将它们都取为零
26、记为=0，常称为零次解。第二步，将零次解代入差分方程式（3.2.10），得出诸内部节点电位值的一次解，它们为：=250=250=250=0同理可得=0在求出一次解的9个节点电位值以后，原来零次解中的9个节点电位值将被一次解中的相应电位值所取代，在计算机的内存中不予保留，从而达到了节省存储空间的目的。第三步，重复上述步骤，令每一个内部节点上的二次解电位值等于该节点周围四个相邻节点（或边界点）一次解电位值的算术平均值，并用二次解电位值冲去内存中的原一次解电位值这样迭代一次又一次的继续下去，可望诸节点的电位值变化越来越小，这时可取这些节点上的电位值为该边值问题的数值解。迭代39次后的结果为=428
27、V， =526V， =428V=187V， =250V， =187V=71V， =98V， =71V利用MATLAB绘制等电位图13，程序见附录，画出等电位图如下：图4.1-2简单迭代法（33）等电位分布图由得出的数值解可以看出，金属盒内点电位分布是越靠中间电位越高，越靠近金属盒顶部电位越高，这是由于金属盒底部和两边的电位都为零，而顶部最高。由此表明此方法计算出的电位值，符合金属盒内的电位分布情况。（2）在盒内取55个离散的电位节点上述例子是在金属盒内取33个离散的电位节点，若是取55个节点，如图4.1-3所示：图4.1-3用有限差分法求金属盒内电位（55）计算步骤同（1），经过91次迭代之后
28、得到电位分布如下： 0 100.0 1000.0 10000 10000 10000 0 0 468.7 629.2 669.4 629.2 468.7 0 0 245.5 378.8 419.3 378.8 245.5 0 0 134.6 221.2 250.0 221.2 134.6 0 0 71.8 121.2 138.4 121.2 71.8 0 0 31.3 53.5 61.3 53.5 31.3 0 0 0 0 0 0 0 0由以上矩阵得知，在盒内取55个节点，求出结果比取33个节点更具有连续性，精确度也越高，但是迭代次数也越多，计算所花的时间久越长。4.2 超松弛迭代法在迭代法
29、的应用中，为加速迭代解收敛速度，一般采用的是超松弛迭代法14。4.2.1 超松弛迭代法理论由于编写计算机程序的需要，每一网格结点的位置由双下标（i，j）予以识别，如图4.2-1所示。对于差分方程（3.2.10）式，采用超松弛迭代法（规定迭代的运算顺序是：从左下角开始做起，即i小的先做；对固定的i，j小的先做。），则关于结点0迭代到第（n+1）次时的近似值，应由如下迭代公式算得 ) （4.2.1）图4.2-1 结点的双下标（i，j）标号式中： 称为加速收敛因子，其取值范围是12，当2时，迭代过程将不收敛。加速收敛因子有一个最佳取值问题，但随具体问题而异。对于第一类边值问题，若一正方形场域由正方形
30、网格分割（每边结点数为m+1），则最佳收敛因子a0可按下式计算 （4.2.2）在更一般的情况下，只能凭借经验取值。在超松弛迭代法的应用中，还必须涉及迭代解收敛程度的检验问题。对此，通常的处理方法是：迭代一直进行到所有内结点上相邻两次迭代解的近似值满足修正条件 （4.2.3）时，终止迭代。即将式(12)作为检查迭代解收敛程度的依据。其中：W是指定的最大允许误差。4.2.2 超松弛迭代法的MATLAB实现（1）程序框图图4.2-2 程序框图（2）MATLAB程序代码如下15：clearhx=5;hy=5; %设置网格节点数（33）v1=ones(hy,hx); %设置行列二维数组v1(1,:)=o
31、nes(1,hx)*1000; %上下两行的Dirichlt边界值v1(hy,:)=zeros(1,hx); %左右两列的Dirichlt边界值for i=1:hy v1(i,1)=0; v1(i,hx)=0;endv2=v1;maxw=1;w=0;n=0；a=1 while(maxw1e-9) %由V1迭代，算出V2，跌代精度为0.000000001 n=n+1; %计算迭代次数 maxw=0; for i=2:hy-1 for 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;%超松
32、弛迭代公式 w=abs(v2(i,j)-v1(i,j); if(wmaxw)maxw=w;end v(i,j)=v2(i,j) end endv1=v2; endn %输出迭代次数contour(v2,20) %画等电位线(3)计算结果如下：n=19v = 0 1000.0 1000.0 1000.0 0 0 428.6 526.8 428.6 0 0 187.5 250.0 187.5 0 0 71.4 98.2 71.4 0 0 0 0 0 0画出等电位分布图：图4.2-3超松弛迭代法（33）等电位分布图当取55个离散的电位节点时，只需将程序hx=5;hy=5改为hx=7;hy=7，经过5
33、4次迭代之后计算出结果为： 0 1000.0 1000.0 1000.0 1000.0 1000.0 0 0 468.7 629.2 669.4 629.2 468.7 0 0 245.5 378.8 419.3 378.8 245.5 0 0 134.6 221.2 250.0 221.2 134.6 0 0 71.8 121.2 138.4 121.2 71.8 0 0 31.3 53.5 61.3 53.5 31.3 0 0 0 0 0 0 0 0画出等电位图：图4.2-4 超松弛迭代法（55）等电位分布图4.3 讨论与分析由（33）和（55）等电位分布图对比可知，两种取点方法求得等电位
34、分布规律一致，但取的离散的点越多，求出的电位连续性越强。迭代收敛的速度与有明显关系表1-1 迭代收敛的速度与的关系收敛因子（）1.01.11.151.21.31.51.82.0（33）迭代次数（）392923192441122发散当收敛因子=1.2时迭代次数最少，计算所用时间最短。通过简单迭代法和超松弛迭代法对比发现，两种方法求出的数值解相同，作出的等电位线分布一样，但超松弛迭代法比简单迭代法收敛速度更快，迭代次数更少，计算时间更短。结 论电磁场边值问题的数值求解有很多种方法，有限差分法是数值计算中应用得最早而又相当简单，直观的一种方法16。本文在讨论电磁场边值问题求解时，将电磁场进行了理想化
35、处理，以一简单边界条件的电磁场边值问题为例，建立相应的数学模型，将场域离散为一些网格，运用差分原理，对场域内偏微分方程及边界上的边界条件进行差分离散化处理，在通过差分运算求出场域内电位值。可以通过上述分析得到这样一些有意义的结论：（1）使用有限差分法求解电磁场边值问题是可行的，只要将网格取得足够小，是可以将离散的点看成连续的。求出离散的数值解，更符合实际应用的需要，而且随着计算机技术的发展，求解差分方程的过程变得简单，使得有限差分法在电磁场问题计算中更具有优势。（2）超松弛迭代法的收敛因子=1时，相当于简单迭代法。收敛因子的导入，使得超松弛迭代法比简单迭代法能更快的收敛，若的值选取恰当，收敛速
36、度还将加快。（3）场域内取的节点越多，计算结果就越精确，当节点划分足够多的时候，离散的点可以看作连续的。但节点划分越多，迭代次数就越多，计算量就越大，所以计算时应根据实际需要，划分合适的节点数。在电磁场计算中如何保证计算精度，减少计算工作量，提高计算速度，减少计算时间，是我们长期努力的方向。参考文献1 冯慈璋 , 马西奎. 工程电磁场导论 M . 北京 : 高等教育出版社 , 2000 :32240 .2 徐立勤 曹伟电磁场与电磁波 M 北京 : 科学出版社,2006:102-107.3 何红雨. 电磁场数值计算法与 MATLAB 实现 M . 武汉 : 华中科技大学出版社 , 2004 :
37、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 ,
38、 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:17-19.1
39、1 黄作英,阙沛文基于MATLAB的电磁场数值分析J 计算机工程与应用2004,36：196-202.12 王洁, 陈超波. 基于MATLAB的静态场边值问题有限差分法的研究J. 微计算机应用, 2010,(03).13 刘瑞桢. MATLAB系列讲座(一)MATLAB简介J. 电脑编程技巧与维护, 1997,(07).14 周波, 李渝曾. 科学计算软件MATLABJ. 电脑技术, 1996,(08).15 蔡苗苗. MATLAB使用技巧几则J. 电脑学习, 2004,(06).16 冯奎胜 卢万铮 朱章虎电磁场数值计算分析方法J 山西电子技术2005,6:43-45.致 谢时光的流逝也许是
40、客观的，然而流逝的快慢却纯是一种主观的感受。当自己终于可以从考研、找工作、毕业论文的压力下解脱出来，长长地吁出一口气时，我忽然间才意识到，原来四年已经过去，到了该告别的时候了。一念至此，竟有些恍惚，所谓白驹过隙、百代过客云云，想来便是这般惆怅了。可是怅然之后，总要说些什么。大学四年，生活其实很简单，只是一些读书、写字和考试的周而复始。如果把这种单调的生活看作一场场循环的演出，那么我只是一个安静的演员。这篇毕业论文也称不上什么精彩的台词，只不过是这种循环演出即将告一段落时的谢幕词。但是无论多么蹩脚的演员，无论台下有多少观众，即使是只说给自己听，在他谢幕时也总要感激一些人，是这些人帮助他走上舞台，
41、成功或者不那么成功地“演出”。我在这里首先要感谢的是我的学位论文指导老师张谢鸿全老师。这篇毕业论文从开题、资料查找、修改到最后定稿，如果没有他的心血，尚不知以何等糟糕的面目出现。我很自豪有这样一位老师，他得我感激和尊敬。感谢和我共度四年美好大学生活的2006级光信息科学与技术1班的全体同学。感谢理学院的所有授课老师，以及学办的老师及我们的分团委书记徐龙玉老师和艾倩宇老师，你们使我终身受益。感谢所有关心、鼓励、支持我的家人、亲戚和朋友。附录 简单迭代法matlab程序代码：clearhx=5;hy=5; %设置网格节点数v1=ones(hy,hx);v1(1,:)=ones(1,hx)*1000; %上下两行的Dirichlt边界值v1(hy,:)=zeros(1,hx); %左右两列的Dirichlt边界值for i=1:hy v1(i,1)=0; v1(i,hx)=0;endv2=v1;maxt=1;t=0;k=0 while(maxt1e-9) %由V1迭代，算出V2，跌代精度为0.000000001 k=k+1; %计算迭代次数 maxt=0; for i=2:hy-1 for j=2:hx-1 v2(i,j)=(v1(i,j+1)+v1(i+1,j)+v2(i-1,j)+v2(i,j-1)/4; %简单迭代公式 t=abs(v2(i,j)-v1(i,j)
免费区 
