电磁场数值分析
- 格式:docx
- 大小:1.31 MB
- 文档页数:20
电磁场数值分析方法及其应用电磁场是无处不在的,它在我们的日常生活中也发挥着极其重要的作用,比如说电视、手机、电脑和家用电器等等。
由于电磁现象的特殊性质,使得电磁场的理论计算非常困难,因此需要引入数值计算方法,对电磁场进行模拟分析,这就是电磁场数值分析方法的基本概念。
一、电磁场数值分析方法简介1. 经典电磁场理论在介绍电磁场数值分析方法之前,我们需要先了解一下经典电磁场理论,也即麦克斯韦方程组。
麦克斯韦方程组描述了电磁场的本质规律,包括电场E、磁场B、电荷密度ρ和电流密度J等四个基本物理量。
这些物理量之间的关系是非常复杂的,因此对于麦克斯韦方程组的求解,需要引入数值计算方法。
2. 电磁场数值计算方法电磁场数值计算方法是指采用离散化方法,将复杂的连续介质分割成有限的、简单的小单元,通过在每个小单元内求解基本电磁场变量的数值解,再通过数值方法进行拼合,最终得到求解区域内的电磁场分布特征。
3. 数值计算方法分类目前常用的电磁场数值计算方法主要包括有限元法、时域有限差分法、频域有限差分法、矩量法等等。
这些方法各有特点,适用于不同的电磁问题求解。
二、电磁场数值分析方法应用1. 微波器件设计微波器件中电磁场的分布特征是十分重要的,它决定了微波器件的性能。
采用电磁场数值分析方法可以清晰地描述微波场的分布特征,从而进行优化和改进设计,提高微波器件的性能。
2. 汽车电磁兼容性分析汽车中各类电子设备的数量越来越多,它们之间的干扰和互相影响也越来越严重。
采用电磁场数值分析方法可以对汽车中的电磁问题进行深入分析,确定干扰成因,从而提出解决方案。
3. 太阳能电池板设计太阳能电池板在光电转化过程中,需要考虑光的反射、折射和吸收等问题。
而这些问题都涉及到电磁场的分布特征。
因此,采用电磁场数值分析方法可以对太阳能电池板的设计进行优化,并提高其能量转换效率。
三、结论电磁场数值分析方法是一种强大的工具,它可以帮助我们深入了解电磁场的本质规律,并对各类电磁问题进行分析和优化设计。
数值分析方法在电磁场计算中的应用电磁场是物理学中最重要的一部分之一,它广泛应用于现代工业、交通、通信、能源和医疗设备等领域。
因此,研究电磁场的行为对于建立新技术和改进现有技术非常重要。
不过由于电磁场是一个非线性的动态系统,因此分析它的行为非常困难。
为了解决这个问题,我们需要数值分析方法来帮助我们更好地理解电磁场的行为。
电磁场的计算方法有很多种,常见的有有限元法、有限差分法等等。
本文将着重介绍有限差分法在电磁场计算中的应用。
有限差分法是经典的数值计算方法,它是一种数值求解偏微分方程的方法。
它的基本原理是将要求解的偏微分方程转化为差分方程,然后利用计算机来求解这个差分方程。
有限差分法的求解过程是离散化的,因此它更便于计算机的处理,同时它的数值误差也比较小。
有限差分法在电磁场计算中的应用非常广泛。
我们可以利用有限差分法来计算电磁场的强度、分布、辐射等参数。
下面我们将介绍一些在电磁场计算中使用有限差分法的实例。
首先,我们来看一个简单的电磁场问题:平面电容器之间的电场强度。
在这个问题中,我们需要求解电场的分布情况。
我们可以利用有限差分法来求解这个问题。
将计算区域离散化成若干个网格点,然后利用电场的高斯定理,将它的积分式子转化为差分式子,最后用差分方程来求解电场值。
在电磁场计算中,还有一些需要注意的问题。
首先是边界条件的处理。
由于有限差分法是一种离散的方法,因此我们需要在计算区域的外部放置边界条件。
这些边界条件包括电场的值、电势的值、电荷密度等等。
其次是计算精度的问题。
由于有限差分法是一种数值方法,因此它的计算精度有时会受到误差的影响。
我们可以通过适当地选择网格点的数量和大小来提高计算精度。
总体来说,有限差分法在电磁场计算中的应用非常广泛,并且具有很好的计算效果。
在实际应用中,我们需要根据具体问题选择合适的数值计算方法,并且在计算时注意处理边界条件和计算精度的问题。
电磁场数值计算与分析技术研究1. 研究背景电磁场是物理学中重要的研究领域,涉及到电磁波传播、电磁辐射、电磁场对物质的影响等多个方面。
在现代科学技术中,电磁场的应用十分广泛,如无线通信、电子设备、雷达测量等。
而电磁场数值计算与分析技术则是电磁场研究中的基础工具,它能够通过计算机模拟的方式帮助我们快速地了解电磁场的特性,分析电磁场对物体的影响。
2. 电磁场数值计算的方法电磁场数值计算的方法主要分为两类,即有限元法和有限差分法。
这两种方法在具体应用中各有优缺点。
有限元法是一种适用于复杂结构的数值计算方法,它将电磁场模型划分为有限个小的单元,然后在每个单元内进行计算,最后整合得到整个模型的计算结果。
有限元法的优点在于它能够处理各种复杂结构,如非线性材料、异形结构等,并且具有精度高、计算速度快等特点。
但是,有限元法的计算成本比较高,需要大量的计算资源,并且需要较高的计算技术水平。
有限差分法是一种比较简单的数值计算方法,它将空间分为一个个离散的网格,然后通过在不同的网格点上进行计算,得到整个空间内的电磁场分布。
有限差分法的优点在于它很容易实现且计算速度快,但是对于复杂的结构和材料效应处理能力较弱,并且需要网格的密度比较高才能够得到比较精确的结果。
3. 电磁场数值计算技术的应用电磁场数值计算技术的应用非常广泛,其中包括电磁波传播、电磁场对物体的影响、电磁设备设计等。
在电磁波传播方面,电磁场数值计算技术可以通过计算电磁波在空间中的传播路径、干扰区域等,来帮助无线通信等领域的设计和优化。
在电磁场对物体的影响方面,电磁场数值计算技术可以帮助我们计算电磁场对物体的激发情况,例如电磁波照射在人体上的吸收情况等,这对于电磁辐射防护等领域非常重要。
在电磁设备设计方面,电磁场数值计算技术可以帮助我们了解电磁场在设备内的分布情况,优化电磁场对设备的影响,提高设备的性能和可靠性。
4. 电磁场数值计算技术的未来发展随着计算机技术的不断进步,电磁场数值计算技术也在不断发展。
电磁场数值分析电和磁现象在自然界普遍存在,两者相互依存形成一个不看分割的整体。
电能产生磁,磁能生电。
很早以前人们就注意到电现象和磁现象,但是两者之间的这种相互联系在很长的一段时间内都没有被人们认识。
直到奥斯特首先发现了通电直导线周围存在磁场这一现象人们才开始把电和磁放在一起来研究。
然而这个时候人们依然没有办法揭示电和磁中间的秘密,只是停留在实验研究阶段,没有形成科学的理论。
1831年法拉第发现了电磁感应定律,从此电和磁的计算可以量化了,人类历史也开启了一个新的时代—电气时代。
由于法拉第的杰出工作,电和磁不再是不可触摸的了,人们已经掌握了运用它的钥匙。
在法拉第之后,另一位杰出的科学家麦克斯韦则更进一步,建立了麦克斯韦方程组,电和磁的理论已经到了相当完美的程度。
现代电机,不管结构多么复杂,都是基于法拉第电磁感应定律和麦克斯韦方程组的原理来运行的,其电和磁的相关量都可以利用这两个定律来进行精确地分析,在设计电机时,我们也是基于这两个定律对电机的电磁过程来进行精确的设计,从而设计出理想的电机。
学会电磁场分析,主要是基于麦克斯韦方程组的相关计算,对电机的学习非常重要。
它为我们今后的学习打下基础。
在学习过程中,主要要把握以下几个度之间的关系:梯度、旋度、散度,这三者的变换正体现了电和磁之间的转换。
一基本原理电磁场的内在规律由电磁场基本方程组—麦克斯韦(Maxwell )方程组表达。
这些方程是由麦克斯韦对大量实验结果及基本概念进行了数学加工和推广归纳而成的。
麦克斯韦方程组是分析和计算电磁场问题的出发点,它既可写成微分形式,又可写成积分形式。
微分形式的麦克斯韦方程组为 t DJ H ∂∂+=⨯∇(1) t BE ∂∂-=⨯∇(2) 0=⋅∇B(3) ρ=⋅∇D (4)式中,E 为电场强度(V/m );B 为磁感应强度(T );D 为电位移矢量(C/m 2);H 为磁场强度(A/m );J 为电流密度(A/m 2);ρ为电荷密度(C/m 2)。
高压输电线路的电磁场数值分析与优化设计1、背景随着经济社会的快速发展,能源的需求日益增长,高压输电线路作为能源输送的主要方式,也在不断地发展。
然而,高压输电线路会产生强烈的电磁场,可能对人体和环境造成潜在的危害,因而有必要对其电磁场进行数值分析与优化设计。
2、电磁场的概念电磁场是指电场和磁场在空间中的分布情况。
无论是直流系统还是交流系统,都会产生电磁场。
电磁场是可以测量的物理量,可以用电场强度和磁场强度来描述。
3、高压输电线路的电磁场高压输电线路的电磁场主要来自于通过导线的电流和被感应的感应电流引起的磁场。
电磁场的强度与电流的大小、线路的设计参数、地形地貌、气象因素等因素有关。
一旦高压输电线路产生了电磁场,就会对环境造成潜在的影响,例如会干扰物体或设备的正常运行,甚至对人员、动植物产生潜在的影响。
4、电磁场的数值分析方法在高压输电线路的设计和规划中,需要对其电磁场进行数值分析。
目前基于计算机的电磁场数值模拟方法可以预测电磁场的强度和时空分布特性。
其中,有限元法和有限差分法是比较常用的数值方法。
在进行数值模拟时,需要准确获取线路的参数,例如线路的频率、形状、材料等信息。
同时,也需要合理的设置模拟的网格和起始条件等信息,以及进行合理的模拟和分析。
5、高压输电线路的电磁场优化设计在高压输电线路设计中,也需要考虑电磁场的影响,进行电磁场优化设计。
通过改变线路的参数,例如线路高度、跨距距离、导线间距、杆塔间距、地面导电率等设计参数,可以减少电磁辐射,降低电磁场的强度。
此外,线路的布局、转角及地形等也需要合理考虑。
通过合理的电磁场设计,可以减少对土地、水资源、环境和人类的影响,最大程度地保障公共利益和民众健康。
6、结论高压输电线路的电磁场是一项需要注意和研究的问题。
通过对电磁场的数值分析和优化设计,可以减少对环境和人身体带来的影响。
此外,随着能源的需求不断增加,高压输电线路的优化设计也需要不断完善,以满足能源输送的需求。
湖北工业大学研究生考试答题纸考试科目工程电磁场数值计算研究生姓名陈天丽学号120130104任课教师邹玲教授学院、专业电气与电子工程学院成绩二0一四年6 月19日《工程电磁场数值计算》课程学习总结这一学期的工程电磁场数值计算学完了,在老师的教导下以及与同学的课堂交流中我学习了很多很多东西,接下来我将从以下七个方面来总结以下这一学期我们学习的东西。
1.高斯消元法 1.1高斯消元法概念高斯消除法是求解线性代数方程组最古老的方法之一。
它不仅容易在计算机上实现,同时,又是构造其他方法的基础。
基本思想:按序逐次消去未知量,把原来的方程化为等价的三角形方程组,或者说,用矩阵行的初等变换将系数矩阵A 约化为简单三角形矩阵;然后按相反方向顺序向上回代求解方程组。
一.下面以一个例子来说明高斯消除法的计算过程。
123123123234 6 (1)352 5 (2)433032 (3)x x x x x x x x x ⎧++=⎪++=⎨⎪++=⎩ 将上述方程写成矩阵形式23463525433032⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦(1)以第一行为基底,消元:12121132*==k k k 131311422*===k k k (2)第二行减去第一行乘以12*k21211112332()02**=+∙=+⨯-=k k k k222212123153()22**=+∙=+⨯-=k k k k23231312324()42**=+∙=+⨯-=-k k k k221312356()42**=+∙=+⨯-=-p p k k(3)同理,第三行减去第一行乘以13*k31311113442()02**=+∙=+⨯-=k k k k32321213433()32**=+∙=+⨯-=-k k k k333313134304()222**=+∙=+⨯-=k k k k331334326()202**=+∙=+⨯-=p p k p变形后矩阵变为234600.544032220⎡⎤⎢⎥--⎢⎥⎢⎥-⎣⎦(4)同理,以第二行为基地,消元:232322360.5*-===-k k k 323212233(3)0.5()00.5**=+∙=--⨯-=k k k k 33331313322(4)()20.5**=+∙=--⨯-=-k k k k331323320(4)()40.5**=+∙=--⨯-=-p p k k再次变形后的矩阵为234600.544004﹣﹣﹣2﹣⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦对应的方程为1232340 (1)++=x x x 230.54 4 (2)-=-x x 32 4 (3)-=-x解得3212813x x x ⎧=⎪=⎨⎪=-⎩二.有限元的方程组的求解方法归纳:13121110112223202122001020300n n n n n n n n n k k p k k p k k k k p k k k k ϕϕϕ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦高斯法如下:以第一行为基底消元:11ij ijp p k *=1111j jk k k *=第二 行减去第一行乘12k *第n 0行减去第一行乘01n k *同理有如下通式111111ii i i i p p p k p p k k **=-∙=-∙111111j ij ij j jij i k k k k kk k k **=-∙=-∙1.2列主元消除法一.基本实例 二.基本思想 给出增广矩阵111211,1212222,112,1a ,b =n n n n n n nnn n a a a a a a a A aa a a +++⎡⎤⎢⎥⎢⎥⎡⎤⎢⎥⎣⎦⎢⎥⎢⎥⎣⎦用增广矩阵表示方程组,在增广矩阵上进行计算,其计算步骤是: (1) 选1,111a max i i i na ≤≤=,交换第1行和第1i ,然后进行消元得,()()()()()()()()()()()()()()111111121n 1,1111111212222,11111n12,1a ,b =n n n n nn n n a a a a a a a A a a a a +++⎡⎤⎢⎥⎢⎥⎡⎤⎢⎥⎣⎦⎢⎥⎢⎥⎢⎥⎣⎦(2) 选()21,1i22a max i i n a ≤≤=,交换第2行和第2i ,然后进行消元,得()()22,b A ⎡⎤⎣⎦依次类推,每次消元前都要换行取最大的列元素为主元 三.列主元消去法技巧和注意在消元过程中适当选取主元素是十分必要的。
工程电磁场数值分析解读工程电磁场数值分析是一种应用有限元法来计算和解决电磁场问题的方法。
该方法通过将电磁场的连续性方程离散化为有限个小单元,再通过求解矩阵方程组来获取数值解。
这种分析方法能够定量计算电磁场的分布和特性,并为工程设计和优化提供重要的参考依据。
对于电磁场数值分析的解读,可以从以下几个方面进行讨论:首先,可以对电磁场的分布进行解读。
通过数值计算,可以得到电磁场在不同位置的数值结果,可以用来表示电磁场的强弱、方向和空间分布特性。
可以对电磁场的分布情况进行比较和分析,以评估电磁场的均匀性和一致性,为设计提供优化方案。
其次,可以对电磁场的特性进行解读。
通过数值计算,可以计算并分析电磁场的一些重要参数,如电场强度、磁场强度、电势、电感、电容等。
这些参数能够揭示电磁场的基本特性,并对电磁设备和系统的工作性能进行评估和优化。
另外,可以对电磁场的影响进行解读。
电磁场数值分析能够计算出电磁场对物体的作用效果,如力、热、电磁感应等。
通过对电磁场的影响进行解读,可以预测电磁场对设备、器件和系统的影响,并为电磁兼容性设计提供技术支持。
此外,还可以对电磁场数值分析方法和结果的准确性进行解读。
电磁场数值分析是一种近似求解的方法,所得数值结果可能与实际情况存在一定差异。
因此,在解读时需要对数值结果进行验证和确认,通过模型实验或其他可靠手段来验证模型的准确性和可靠性。
总之,工程电磁场数值分析是一种重要的工程设计方法,能够定量计算和解决电磁场问题。
通过对电磁场分布、特性和影响等方面进行解读,可以为工程设计和优化提供重要的参考依据。
同时也需要关注分析方法的准确性和结果的可靠性,以确保分析结果的准确性。
电磁场数值模拟技术的研究与应用电磁场是物理学中重要的概念之一,主要描述有电荷的物质间存在的相互作用。
电磁场的研究在生产生活中具有广泛的应用。
例如,电磁波可以作为通讯媒介,还能够在医学领域中被应用于成像等方面。
因此,电磁场数值模拟技术的研究和应用具有非常重要的意义。
一、电磁场数值模拟技术电磁场数值模拟技术主要是指对电磁场的分析和计算。
目前主要采用有限元方法、有限差分法、矢量势法、边界元法等来进行电磁场数值模拟。
1、有限元方法有限元方法是利用算法将大的物理问题分解成一些小的简单问题,然后进行求解的数值分析方法。
它的优点在于可以分析复杂的结构和三维几何形态,也可以用复杂的边界条件解决问题,因此能够在避免实验中的高成本和高风险的情况下,更好地掌握电磁场的性质。
2、有限差分法有限差分法是采用差分方程模拟物理问题的一种方法,具有简单易学,适用性强等优点。
它可以通过计算网格内某一点的导数或者数量值,来计算整个网格的电磁特性,所有计算数据都是通过互相关的差分公式来预测。
在建立数值模拟模型时,有限差分法非常适用于模型边界不规则,模型尺寸比较大的情况。
3、矢量势法矢量势法是一种用于求解电磁场分布的迭代计算方法。
该方法基于矢量势函数,其优点在于可以加速计算速度。
在进行矢量势方法计算时,需要对场区域进行离散化处理,然后以矢量势为主要电磁学量,通过数值迭代方法反复计算直到最终的模型求解出结论。
4、边界元法边界元法是基于经典场论,利用边界元法与物理问题的相关性,将物理问题进行离散化,通过求解离散化后的物理问题模型而达到计算结果的数值方法。
边界元法包括时间域和频域两种形式,它克服了有限元方法和有限差分法等方法的缺点,也能解决更加复杂的边界条件问题。
但它的计算较为复杂,其计算量也会随着模型的复杂程度增加而增加。
二、电磁场数值模拟的应用1、电磁场模拟在电磁波通信中的应用电磁波通信是一个广泛的领域,电磁场模拟可以用来研究和优化通信信道的性能。
电磁场数值分析引言电磁场是物理学中一个重要的研究领域,涉及到各种现实世界中的物理现象,如电磁感应、电磁波传播等。
为了更好地理解和研究电磁场,数值分析成为一种重要的工具。
本文将介绍电磁场数值分析的基本概念、方法和应用。
电磁场基本概念电磁场指的是由电荷和电流引起的电场和磁场的组合。
电场是由电荷引起的一种物理场,其描述了电荷间的相互作用。
磁场则是由电流引起的一种物理场,其描述了电流的磁性效应。
电磁场的数值分析主要涉及以下概念:1.电场强度:指在某一点产生的电场的强度,通常用矢量表示。
2.磁场强度:指在某一点产生的磁场的强度,也通常用矢量表示。
3.电势:指在某一点产生的电场对单位正电荷所做的功。
4.磁感应强度:指在某一点产生的磁场对单位正电荷所做的功。
电磁场数值分析方法电磁场数值分析基于数值计算方法,通过离散化的方式将连续的电磁场问题转化为离散的数值问题。
常用的电磁场数值分析方法包括有限差分法(Finite Difference Method, FDM)、边界元法(Boundary Element Method, BEM)、有限元法(Finite Element Method, FEM)等。
有限差分法有限差分法是一种基于差分近似的数值计算方法,将连续的变量离散化为有限个节点上的变量。
在电磁场数值分析中,有限差分法通常用于解决电场或磁场的分布问题。
该方法将空间离散化为网格,通过差分近似计算相邻节点间的电势或磁感应强度。
边界元法边界元法是一种基于积分方程的数值计算方法,将连续的物理场问题转化为边界上的积分方程。
在电磁场数值分析中,边界元法通常用于解决边界值问题,如电势或磁场在给定边界上的分布。
该方法通过将边界上的物理量表示为边界上的基本解的线性组合,通过求解线性方程组得到物理量的数值解。
有限元法有限元法是一种基于变分原理的数值计算方法,将连续的问题离散化为有限个元素上的问题。
在电磁场数值分析中,有限元法通常用于解决较为复杂的问题,如非线性材料的电磁场问题。
无刷直流电机电磁场数值分析及参数优化设计无刷直流电机是一种应用广泛的电机类型,它具有高效、低噪音、轻巧等优点,在工业制造、家用电器、医疗设备等领域中有广泛的应用。
其中,电磁场是无刷直流电机最核心的部分,其性能直接影响着无刷直流电机的工作效率、功率、寿命等方面,因此对其进行电磁场数值分析及参数优化设计便显得尤为重要。
一、无刷直流电机电磁场模型在进行无刷直流电机电磁场数值分析及参数优化设计前,首先需要了解其电磁场模型。
无刷直流电机的电磁场模型可以理解为一个二维、稳态、线性、磁滞、非线性耗散的电磁场问题,其中会涉及到各种不同的物理量。
1、电机几何模型首先,我们需要将无刷直流电机的几何模型转化为计算机可计算的三维有限元模型。
一般来讲,采用商业软件如ANSYS、COMSOL等进行建模。
其次,在进行建模时需要考虑电机叶片、铁芯、永磁体、电线等各种结构对磁场分布的影响。
2、动、静子系统的等效电路模型其次,应该对无刷直流电机的动子和静子系统的等效电路进行建模,以求解其特定电路方程。
动子系统一般由电阻、电感、电源电压、反电动势等元件组成;静子系统一般由电阻、电感、电源电压、永磁体等元件组成。
3、麦克斯韦方程最后,针对定子和转子的磁场建立麦克斯韦方程,对电机内的电磁场分布进行描述。
二、无刷直流电机电磁场数值分析1、网格剖分进行无刷直流电机电磁场数值分析时,需要将其几何模型转化为计算机可计算的三维有限元模型,并将其划分为网格。
网格剖分的质量将直接影响计算结果的准确度和效率。
因此,在进行剖分时需要兼顾划分精度和计算复杂度,提高计算速度与精度的平衡。
2、边界条件的设置进行无刷直流电机电磁场数值分析时,设置正确的边界条件也十分重要。
边界条件需要考虑电机结构、电路结构及其材料等多方面因素,不同条件的设置会直接影响计算结果的准确度。
3、电机工作状态下的电磁场分析在进行电磁场数值分析时,需要考虑无刷直流电机不同工作状态下的电磁场分析,如静态状态、动态状态以及电机启动和停止过程中的电磁场变化等。
电磁场数值分析《电磁场数值分析》(作业)--- 2016学年 ---学院:学号:姓名:联系方式:任课教师:2016年6月6日作业1一个二维正方形(边长a=10mm)的静电场区域,电位边界条件如图所示(单位:V),求区域内的电位分布。
要求用超松弛迭代法求解差分方程组进行计算。
代码:hx=11;hy=11;v1=zeros(hy,hx);v1(hy,:)=ones(1,hx)*100;v1(1,:)=ones(1,hx)*50;for i=1:hy;v1(i,1)=0;v1(i,hx)=100;endw=2/(1+sin(pi/(hx-1)));maxt=1;t=0;v2=v1;n=0;while(maxt>1e-6)n=n+1;maxt=0;for i=2:hy-1;for j=2:hx-1;v2(i,j)=(1-w)*v1(i,j)+w*(v1(i+1,j)+v1(i,j+1)+v2(i-1,j )+v2(i,j-1))/4;t=abs(v2(i,j)-v1(i,j));if (t>maxt)maxt=t;endendendv1=v2;endsubplot(1,2,1)mesh(v2)axis([0,11,0,11,0,100])subplot(1,2,2)contour(v2,20)结果:coef1=dt/(mu0*dx);coef2=dt/(eps0*dx);coef3=(c*dt-dx)/(c*dt+dx);ezold=ez;for now=1:timestep;hx=hx-coef1*(ez(:,2:ymesh+1)-ez(:,1:ymesh));hy=hy+coef1*(ez(2:xmesh+1,:)-ez(1:xmesh,:));ez(2:xmesh,2:ymesh)=ez(2:xmesh,2:ymesh)-...coef2*(hx(2:xmesh,2:ymesh)-hx(2:xmesh,1:ymesh-1))-...coef2*(hy(2:xmesh,2:ymesh)-hy(1:xmesh-1,2:ymesh));ez(1,:)=ezold(2,:)+coef3*(ez(2,:)-ezold(1,:));ez(xmesh+1,:)=ezold(xmesh,:)+coef3*(ez(xmesh,:)-ezold (xmesh+1,:));ez(:,1)=ezold(:,2)+coef3*(ez(:,2)-ezold(:,1));ez(:,ymesh+1)=ezold(:,ymesh)+coef3*(ez(:,ymesh)-ezold (:,ymesh+1));ez(xmesh/2+1,ymesh/2+1)=sin(now*dt*2*pi*c/25.0); mesh(ez)pause(0.01)ezold=ez;end结果:作业3基于Pocklington方程用MoM分析半波对称振子天线:观察天线线径和分段数目分别取不同值对天线阻抗和辐射特性的影响(半径分别取0.001λ,0.0001λ,0.00001λ,分段数取11,21,31)代码:%%初始化参数c=3e-8;r=1;f=c/r;w=2*pi*f;e0=8.85e-12;u0=4*pi*1e-7;a=0.0001*r;L=0.5*r;k=2*pi/r;N=31;dl=L/(N+1);l=L/2-dl/2;lz=-l:dl:1;lzs=lz(1:N);lzm=lz(1:N)+dl/2;lze=lz(2:N+1);%%阻抗矩阵元素求解fi=log(dl/a)/(2*pi*dl)-k/(4*pi)*1j;fi_1=exp(-k*dl*1j)/(4*pi*dl);fi_2=exp(-k*2*dl*1j)/(8*pi*dl);z=ones(N,N);for m=1:Nfor n=1:Nif m==nfi1=fi;fi2=fi_1;fi3=fi_2;z(m,n)=((k^2*dl^2-2)*fi1+fi2+fi3);elseif abs(m-n)==1fi1=fi_1;fi2=fi;fi3=fi_2;z(m,n)=((k^2*dl^2-2)*fi+fi2+fi3);elsefi1=exp(-k*abs(m-n)*dl*1j)/(4*pi*abs(m-n)*dl);fi2=exp(-k*abs(m+1-n)*dl*1j)/(4*pi*abs(m+1-n)*dl);fi3=exp(-k*abs(n+1-m)*dl*1j)/(4*pi*abs(n+1-m)*dl); z(m,n)=((k^2*dl^2-2)*fi+fi2+fi3);endendend%%电压矩阵求解V=zeros(N,1);V((N+1)/2)=-1*(1j*w*e0);I=z\V;Z_in=1/I((N+1)/2);disp(['输入阻抗=',num2str(Z_in)])I_amp=abs(I);Max=max(I_amp);Iunit2=[0;I_amp/Max(1);0];figure(1)h=0:dl/r:L/r;Ithe=sin(pi*h*r/L);plot(h,Iunit2,'b',h,Ithe,'r','linewidth',2)legend('pocklinton','解析值')grid onxlabel('电长度')ylabel('归一化电流')%%方向图theta=0:0.01:2*pi;abs_f=zeros(1,length(theta));for n=1:1:Nabs_f=abs_f+I(n)*exp(k*(n*dl-L/2)*cos(theta)*1j);endabs_f=abs(sin(theta)*dl.*abs_f);Max_f=abs(sum(I)*dl);Far_patten2=abs_f/Max_f(1);theta_2=0:0.1:2*pi;Far_theory=abs((cos(k*(L/2)*cos(theta_2))-cos(k*L/2)) ./sin(theta_2));figure(2)polar(theta,Far_patten2,'-b')hold onpolar(theta_2,Far_theory,'or')hold offlegend('pocklinton','解析值')title('半波阵子天线E面方向图')figure(3)polar(theta,ones(1,length(theta)),'-b')title('半波阵子天线H面方向图')%%半波阵子增益I_in=I((N+1)/2);A=(w*u0)^2/(4*pi*sqrt(u0/e0)*real(Z_in)*(abs(I_in))^2 );G_theta=A*abs_f.^2;Max_gain=max(G_theta);Max_gain_dB=10*log10(Max_gain);disp(['半波阵子增益=',sprintf('%.4fdB',Max_gain_dB)])结果:作业4基于电场积分方程用MoM分析对称振子天线:计算振子总长度分别为0.25λ ,0.5λ,λ,1.5λ时,振子的输入阻抗和E面方向图。
代码:lamda=1;a=0.0001;me=8.85e-12;mu=4*pi*(1e-7);arg=2*pi*(3e-8)/lamda;L=0.2*lamda;k=2*pi/lamda;N=21;dL=L/(N+1);l=L/2-dL/2;lz=-l:dL:1;lzs=lz(1:N);lzm=lz(1:N)+dL/2;lze=lz(2:N+1);for m=1:Nfor n=1:Nif n==mFmnmm=(1/(2*pi*dL))*log(dL/a)-1j*k/(4*pi); Fmnee=(1/(2*pi*dL))*log(dL/a)-1j*k/(4*pi); Fmnss=(1/(2*pi*dL))*log(dL/a)-1j*k/(4*pi); Fmnse=exp(-1j*k*dL)/(4*pi*dL);Fmnes=exp(-1j*k*dL)/(4*pi*dL);elseif abs(n-m)==1Fmnmm=exp(-1j*k*dL)/(4*pi*dL);Fmnee=exp(-1j*k*dL)/(4*pi*dL);Fmnss=exp(-1j*k*dL)/(4*pi*dL);if n>mFmnse=exp(-1j*k*2*dL)/(4*pi*2*dL);Fmnes=exp(1/(2*pi*dL))*log(dL/a)-1j*k/(4*pi);elseFmnes=exp(-1j*k*2*dL)/(4*pi*2*dL);Fmnse=exp(1/(2*pi*dL))*log(dL/a)-1j*k/(4*pi);endelsenum=abs(n-m);Fmnmm=exp(-1j*k*num*dL)/(4*pi*num*dL);Fmnee=exp(-1j*k*num*dL)/(4*pi*num*dL);Fmnss=exp(-1j*k*num*dL)/(4*pi*num*dL);if n>mFmnse=exp(-1j*k*(num+1)*dL)/(4*pi*(num+1)*dL);Fmnes=exp(-1j*k*(num-1)*dL)/(4*pi*(num-1)*dL);elseFmnes=exp(-1j*k*(num+1)*dL)/(4*pi*(num+1)*dL);Fmnse=exp(-1j*k*(num-1)*dL)/(4*pi*(num-1)*dL);endendz(m,n)=1j*arg*mu*dL*dL*Fmnmm+(1/(1j*arg*me))*(Fmnee-F mnes-Fmnse-Fmnss);endendV=zeros(N,1);fedp=(N+1)/2;V(fedp)=1;I=linsolve(z,V);Z=V(fedp)/I(fedp);theta=0:pi/100:2*pi;ftheta=0;for m=1:length(theta)F(m)=0;for n=1:NF(m)=F(m)+I(n)*exp(1j*k*(n-fedp)*dL*cos(theta(m)));endendF=abs(F);F1=F/max(F);polar(theta,F1,'b.-')结果:作业5对课程的建议、自己的收获等。