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

.电磁场实验报告实验一模拟电偶极子的电场和等位线学院：电气工程及其自动化班级：学号： 姓名:;.实验目的：1、了解并掌握 matlab软件，熟练运用 matlab语言进行数值运算。2、熟练掌握电偶极子所激发出的静电场的基本性质3、掌握等位线与电力线的绘制方法实验要求：1、通过编程，完成练习中的每个问题，熟练掌握matlab的基本操作。2、请将原程序以及运行结果写成word 文档以方便检查实验内容：一、相关概念回顾对于下图两个点电荷形成的电场两个电荷共同产生的电位为：q11 ()p4 rrqr2r14r r01201 2其中距离分别为 r(xq x)2( yq y)2 ， r(xq x)2( yq y)2111222电场强度与电位的关系是e p等位线函数为 :( x,y, z)c电力线函数为：exeydxdy二、实验步骤1、打开 matlab软件，新建命令文档并保存，并在文档中输入程序。2、输入点电荷 q1 的坐标（ q1x，q1y） , 以及 q1 所带的电量。调用input 函数。如果不知道该函数的使用方法可在matlab命令行处键入doc input。3、输入点电荷 q1 的坐标（ q1x，q1y） , 以及 q1 所带的电量。4、定义比例常系数19e9 ， 命令为 k=9e9。4 05、定义研究的坐标系范围为x5,5 , y5,5,步长值为 0.1。6、将x,y两组向量转化为二维坐标的网点结构，函数为meshgrid。命令为x,y=meshgrid(x,y) ，如果不知道该函数的使用方法可在matlab 命令行处键入doc meshgrid。7、计算任意一点与点电荷之间的距离r，公式为 r(xq x)2( yq y)2 ，111222r(xq x)2( yq y) 2vq( 11 )8、计算由 q1,q2 两个点电荷共同产生的电势4 0r1r29、注意，由于在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)e q11、计算 e 的模值22eexy，注意在计算时运算要加点，ex.212、计算电场强度的单位矢量，ex加点， ey=ey./ eqexe ， eyeye ，注意在计算时运算要13、生成你要绘制的等位线的数量与每条等位线上的电位值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(fontnameimpactfontsize163 ?); 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三、开放性试验画出电偶极子的等位线和电力线( rd )在球坐标系中 ,通过用二项式展开，又有rd，得pq11 ()4 rrqr2r14 r r01201 2r( xq x)2( yq y)21112r(r 2drd241cos) 2d 22r1( r41r dc o s 2 )用二项式展开，又有rd，得r1r所以d cos 2qdcosr2rp erd c o s 2p4 r 24 r 200p=qd, 表示电偶极矩（ dipole moment)，方向由 -q 指向 +q。等位线方程( 球坐标系) ：rccosd rr de ree pq 4 0 r3(2cosersine )将 e和 er 代入 e 线方程有rd sin2.q1x=1q1y=2q1=10q2x=1q2y=-2q2=10;.q1x=1 q1y=2 q1=10q2x=1 q2y=-2 q2=-10q1x=1 q2y=2 q1=10q2x=-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 构造插值函数 :其中 ah 称为形状函数。那么有插值函数的一阶偏导数为:从而得到能量函数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 有如下关系 :ba故可简化为椭圆方程 :2aj画出大小为 2* 2 的矩形 r1,两导线用直径为0. 2、相距 0. 8 的两个圆表示 . 矩形的边界条件是 di richlet 边界条件 ,取 h= 1, r= 0。这种做法是模拟远处的磁场势为零.在设置方程类型时,选取应用模式为elliptic( 椭圆)方程, 对于矩形其它部分系数取mangetostatics.故在选择方程时选取=1、j=0.在表示导线的圆内