利用MATLAB计算电磁场有关分布概要

电磁场实验报告实验一模拟电偶极子的电场和等位线学院：电气工程及其自动化班级：学号：姓名:实验目的：1、了解并掌握MATLAB 软件，熟练运用MATLAB 语言进行数值运算。

2、熟练掌握电偶极子所激发出的静电场的基本性质3、掌握等位线与电力线的绘制方法实验要求：1、通过编程，完成练习中的每个问题，熟练掌握MATLAB 的基本操作。

2、请将原程序以及运行结果写成word 文档以方便检查实验内容：一、相关概念回顾对于下图两个点电荷形成的电场两个电荷共同产生的电位为：2101201211()4π4πp r r q q r r r r φεε-=-= 其中距离分别为22111()()r x q x y q y =-+-，22222()()r x q x y q y =-+- 电场强度与电位的关系是p φ=-∇E 等位线函数为: (,,)x y z C φ=电力线函数为：d d yx E E x y=二、实验步骤1、打开MATLAB 软件，新建命令文档并保存，并在文档中输入程序。

2、输入点电荷q1的坐标（q1x ，q1y ）, 以及q1所带的电量。

调用input 函数。

如果不知道该函数的使用方法可在MATLAB 命令行处键入 doc input 。

3、输入点电荷q1的坐标（q1x ，q1y ）, 以及q1所带的电量。

4、定义比例常系数90194πe ε=， 命令为 k=9e9。

5、定义研究的坐标系范围为[][]5,5,5,5x y ∈-∈-,步长值为0.1。

6、将x,y 两组向量转化为二维坐标的网点结构，函数为meshgrid 。

命令为[X,Y]=meshgrid(x,y)，如果不知道该函数的使用方法可在MATLAB 命令行处键入 doc meshgrid 。

7、计算任意一点与点电荷之间的距离r，公式为1r =，2r =8、计算由q1,q2两个点电荷共同产生的电势01211()4πq V r r ε=-9、注意，由于在q1和q2位置处计算电势函数为无穷大或者无穷小，因此要把这两点去掉掉，以方便下面绘制等势线。

具体命令可参考Vinf1=find(V==inf); V(Vinf1)=NaN;Vinf2=find(V==-inf); V(Vinf2)=NaN;如果是可以解释这四句话的原理，可以有加分！10、根据天长强度与电位函数的关系φ=-∇E ，可直接计算E ，调用gradient 函数。

如果不知道该函数的使用方法可在MATLAB 命令行处键入 doc gradient 。

参考命令为[Ex,Ey]=gradient(-V)11、计算E的模值q =E Ex.^212、计算电场强度的单位矢量，x x e E =E ，y y e E =E，注意在计算时运算要加点，Ey=Ey./ Eq13、生成你要绘制的等位线的数量与每条等位线上的电位值 cv=linspace(min(min(V)),max(max(V)),49)该命令表示在最大电位与最小电位之间插入49个点，形成一个向量cv 14、绘制等位线contourf (X,Y,V,cv,'k-')如果不知道该函数的使用方法可在MATLAB命令行处键入doc contourf。

15、进行一些修饰axis('square')title('\fontname{Impact}\fontsize{16}³¡ÓëµÈλÏß');hold on16、绘制电场线quiver(X,Y,Ex,Ey,0.5)如果不知道该函数的使用方法可在MATLAB命令行处键入doc quiver。

17进行一些修饰plot(q1x,q1y,'wo')plot(q2x,q2y,'ws')xlabel('x')ylabel('y')hold off18、结果验证（1）q1x=1,q1y=0,q1=4e-9; q1x=-1,q1y=0,q2=-4e-9（2）q1x=1,q1y=1,q1=10e-9; q1x=-1,q1y=-1,q2=-4e-9（3）q1x=1,q1y=1,q1=100e-9; q1x=-1,q1y=-1,q2=100e-9三、开放性试验画出电偶极子的等位线和电力线 ( r >>d )在球坐标系中,通过用二项式展开，又有r >>d ，得2101201211()4π4πp r r q q r r r r φεε-=-= 22111()()r x q x y q y =-+-12222(cos )4d r r rd θ=++ 12221(cos )4d r r rd θ=+-用二项式展开，又有r >>d ，得 1cos 2d r r θ=- 2cos 2dr r θ=+所以 2200cos 4π4πr p qd r r θφεε⋅==p ep=qd,表示电偶极矩（dipole moment)，方向由-q 指向 +q 。

等位线方程 ( 球坐标系 ) ：cos r C θ=d d r r r E E θθ= 30(2cos sin )4πp r q rθφθθε=-∇=+E e e将E θ和E r 代入E 线方程有2sin r D θ=Q1x=1 Q1y=2 Q1=10 Q2x=1 Q2y=-2 Q2=10Q1x=1 Q1y=2 Q1=10 Q2x=1 Q2y=-2 Q2=-10Q1x=1 Q2y=2 Q1=10 Q2x=-1 Q2y=-2 Q2=10实验二 MATLAB电磁场有限元计算实验目的：4、了解有限元算法的原理，熟练运用MATLAB环境的PDE工具。

5、熟练运用PDE工具分析简单的电磁场边值问题。

实验内容：一、有限元简介在电磁场的计算中, 仅对那些具有最简单边界条件和场域几何形状规则的问题才有解析解, 多数问题的求解必须用数值计算的方法,其场域分布的数值计算内容是学习难点。

本实验将有限元法和Matlab 结合起来对电磁场教学中的电位分布问题进行计算。

结果表明使用Matlab对有限元分析编程中的矩阵进行处理,程序设计清晰简便,易于理解和实现。

node(节点)element(单元)有限元法是以变分原理和剖分插值为基础的一种数值计算方法,其基本思想是将场域方程等价为一个条件变分问题,然后由条件变分问题对场域进行剖分离散为方程组进行求解。

对于一个电场来说,其储能总是趋于最小,这样变分法的泛函和电场的储能就联系起来了。

对于边界为L 的无源空气介质二维静电场中,一个封闭场域S 内的等价能量泛函可以写为:在有限元分析中,将所研究的区域S划分成有限的n 个三角形网格单元。

对应m个节点, ds 为单元e的面积。

对任意三角形单元e 中任一点的电位可以认为由该三角形的三个节点(分别设为i、j、k) 上的电位u 随该点坐标x、y 变化而线性确定。

因此, 对于单元e 构造插值函数:其中a h称为形状函数。

那么有插值函数的一阶偏导数为:从而得到能量函数We:则将单元e中的能量函数We 对每一个节点电位ul ( l = i, j , k)求一阶偏导数, 得:表示为矩阵形式有：然后进行总体合成, 将各单元的能量函数对同一节点的电位一阶偏导数相加, 获得所要求解的线性方程组。

由以上分析,可知在该场域内电场有限元数学模型为:式中U 为n 个节点处的待求电位, K 为n 阶矩阵。

最后进行强加边界条件处理, 消去已知电位节点在系数矩阵中所在的行和列, 得到简化后的方程,继而可以对电位进行求解。

流程框图如下图所示。

二、静电场仿真静态场满足上方基本方程，式中 D 为电位移, 为电荷密度, H 为磁场强度, J 为电流密度, E为电场强度, B为磁感应强度.对于恒定的电场:式中电位满足泊松( Poisson)方程:对于不存在电荷的空间部分有电荷体密度为零,上式退化为拉普拉斯( Laplace) 方程:利用上述方程, 再加上边界条件, 利用Matlab 中的偏微分工具箱, 即可求解带电体周围空间的电场分布.输入pdetool可进入软件环境。

两点电荷的电场：两等值异号点电荷单位, 两者间距为1,求其电势分布.整个求解域取中心为原点,半径为 2 的圆,两空间电荷点位置为(-0.5,0)和( 0.5,0),作为一种近似,画一个尽量小的圆,取半径为0.05.大圆的边界条件是Di richlet边界条件,取h= 1, r= 0,这种做法是模拟远处的电势为零.由于大圆与小圆之间的区域没有电荷, 满足Laplace 方程, 因此在选择方程时选取Elliptic(椭圆)方程,其方程类型为:取系数为c= 1, a= 0, f= 0.在表示点电荷的小圆内, 我们认为电荷是均匀分布的, 满足Poisson 方程, 在选择方程时也取Elliptic方程, 取系数为c= 1, a= 0, f= 0. 2.其两点电荷电势分布上图所示,电力线用箭头表示.三、静电场中的导体问题描述: 在电场强度为E 的静电场中放置一根无限长的导体,研究截面上的电势分布。

首先画一个2*2的矩形R1,然后在中心原点画半径为0. 3的圆E1.然后将Set formula对话框中的公式改为R1-E1,表示求解区域为二者之差.矩形所有的边界条件是Dirichlet边界条件, 取h=1, r= y.而在圆的边界取h=1, r=0.由于求解域没有电荷,因此在选择方程时选取Elliptic(椭圆)方程,系数为c=1, a=0, f=0.其电势分布如下图所示,电力线用箭头表示.四、两根载流长直导线的磁场问题描述: 两根载流长直导线,相距为0.8,导线直径为0.2, 求电流引起的磁场. 从麦克斯韦(Maxwell)方程组出发,其磁场强度B和磁感应强度H的关系为:μ=B H磁场势A 与B 有如下关系: =∇⨯B A故可简化为椭圆方程:2μ∇=-A J画出大小为2* 2的矩形R1,两导线用直径为0. 2、相距0. 8 的两个圆表示.矩形的边界条件是Di richlet 边界条件,取h= 1, r= 0。

这种做法是模拟远处的磁场势为零.在设置方程类型时, 选取应用模式为Mangetostatics.故在选择方程时选取Elliptic(椭圆)方程, 对于矩形其它部分系数取μ=1、J=0.在表示导线的圆内,取μ= 1, J=1.两根载流长直导线的磁场势和磁力线如下图所示,磁力线用箭头表示.实验要求：a) 有限元法计算电磁场问题的基本思想是什么？b) 求解静电场两个等量异号带电导线电场位函数及电场强度。

c) 求解静电场中的导体等等势面，写出导体与电介质分界面上的衔接条件，并与实验结果对比。