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

1、电磁场实验报告实验一 模拟电偶极子的电场和等位线学院：电气工程及其自动化班级：学号：姓名:实验目的：1、了解并掌握 MATLAB 软件，熟练运用 MATLAB 语言进行数值运算2、熟练掌握电偶极子所激发出的静电场的基本性质3、掌握等位线与电力线的绘制方法实验要求：1、 通过编程，完成练习中的每个问题，熟练掌握MATLAB 的基本操作2、请将原程序以及运行结果写成 word 文档以方便检查实验内容: 一、相关概念回顾对于下图两个点电荷形成的电场其中距离分别为 *f- qx)2(y - qy)2, r2- q2x)2(y - q2y)2电场强度与电位的关系是E p 等位线函数为：(x, y, z)

2、 =C电力线函数为：ExEydx dy二、实验步骤1、打开 MATLAB 软件，新建命令文档并保存，并在文档中输入程序。两个电荷共同产生的电位为：厂亡1 1( )1 $2、 输入点电荷 q1 的坐标(qlx，q1y),以及 q1 所带的电量。调用 in put 函数。 如果不知道该函数的使用方法可在 MATLAB 命令行处键入 doc in put。3、输入点电荷 q1 的坐标(qlx, q1y),以及 q1 所带的电量。14、定义比例常系数=9e9，命令为 k=9e9。4n05、 定义研究的坐标系范围为5,5y一 5,5】,步长值为 0.1。6 将 x,y 两组向量转化为二维坐标的网点结构，

3、函数为meshgrid。命令为X,Y=meshgrid(x,y)，如果不知道该函数的使用方法可在 MATLAB 命令行处键入 docmeshgricL7、 计算任意一点与点电荷之间的距离r，公式为* = (x-q/)2 (y-qy)2，D =(x - q2X)2(y -q2y)2V=(丄一丄)8、计算由 q1,q2 两个点电荷共同产生的电势4n0129、注意，由于在 q1 和 q2 位置处计算电势函数为无穷大或者无穷小，因此要把 这两点去掉掉，以方便下面绘制等势线。具体命令可参考Vi nf1=fin d(V=i nf);V(Vi nf1)=NaN;Vi nf2=fi nd(V=-i nf);V

4、(Vi nf2)=NaN;如果是可以解释这四句话的原理，可以有加分！10、 根据天长强度与电位函数的关系 E，可直接计算 E,调用 gradient 函数。如果不知道该函数的使用方法可在MATLAB 命令行处键入 doc gradient。参考命令为Ex,Ey=gradie nt(-V)E |=E2+E2、亠,11、 计算 E 的模值Eqx Ey，注意在计算时运算要加点，Ex.A212、 计算电场强度的单位矢量，$ = 丿|E，勺=| E，注意在计算时运算要 力卩点，Ey=Ey./ Eq13、生成你要绘制的等位线的数量与每条等位线上的电位值cv=li nspace (min(min( V),m

5、ax(max(V),49)该命令表示在最大电位与最小电位之间插入49 个点，形成一个向量 cv14、绘制等位线con tourf (X,Y,V,cv,k-)如果不知道该函数的使用方法可在MATLAB 命令行处键入 doc contourf。15、进行一些修饰axis(square)title(fontnamelmpactfontsize163 ?o?(ie?);hold on16、绘制电场线quiver(X,Y,Ex,Ey,0.5)如果不知道该函数的使用方法可在MATLAB 命令行处键入 doc quiver。17 进行一些修饰plot(q1x,q1y,wo)plot(q2x,q2y,ws)x

6、label(x)ylabel(y)hold off18、结果验证(1) q1x=1,q1y=0,q1=4e-9; q1x=-1,q1y=0,q2=-4e-9场与竽位线I I L j I L J I i I I J J i* I * * 8 I* d ! I I I-1 I H M 4 I昇卜和I tifjfjJjfj I I h 4 I L d I iI I h I Li J I I I I il I I I J Z * i* J i IIMEII汕I.川汕加；加: 出卅|：常1忖I陶川齢加;I h、 L LI I F I I b J j. i!d I I . M k S I:HIEilJif

7、iiSJ；:二亠：二”：; ;- - s-fe-s-fe- : : n:n:kuku nunu. . ; ;二；：二:：;:k:%v戶 x忙 r. i h h J ? h ! riI j I：卜f| ! !屮f八F 11!. i . i I i I 1 i J /i r r j j rH Jzz iP i i-WS.4i、lTt i aI I i IFUi J申八2 * * iim * *壮讥計:冷I$;弭$门弭： $爲弗：髭書醫鲁力*1 I hLhll沾h、L、 、 、,ZZX/jFJ JJ ? J .J |b| |J |物由仃曰川加出曲帝眶岳mHS1!Z Frqrjjrirji. iiJi

8、aiibii.|ih I f r F.%JLXli I * F * 2 Pi* 2 r m承承j,j_i. r .2 r、I i 1 I 4 i r r f i* 1 j * H I .KS-S-K* * *f / t r r J f ii ! i . . i r i川ii铝ii报讯ry x r-= -riUWFmmmiLiLUTiL L I I.t i i i i i J Ji、.m m % i, i i 5 h b 1 Ii丨ii V %u 1 fc. rrf jmmFFkrr”i*J;F. . .B.s r .- / IZ. r r . .I-/ p；i- , , j r I I J r

9、 I I , h ! *1 5、.、.、 r.i r r s* ? I I* 1 I r ! I I I 1 *1、 .“UP.FFUrj-.i-. |耳1|1厂|3、&frpJ r i1I I I S * ”J*.-J*J . J,”阳ym T f H I f M H卄* F-r.r.hn I*4 i i* I t I 1 1 *rrFPui / i i r J I I . 1 i 丿t*i* .m m i qi”Fi i ii i 11u it“ r 2 -. 1.1. _ _ F ffi壮f Ef *4 It I I訂II、)m,r I .1 .1 r i 3 N -I I I

10、f丨iI勺i勺 IIi, L * Ff zz-.U r. *、.、w l-T- 1s. - x:;：；：汁：讥： ： ；壮讥讥：松：汀；汀汀鳥;2衣】：&乂苫严“it tssiu1 LI4r A1 . V R；I I I I ! I . J F J f f T f f f fjf fJ*J* f r J*J* f1| Isr . /. f r r | , | | 1 L I h hJ t r f J d I I J I. J I LILI I him h、IS L h X % h V h F* u、i ! E百i j衬I i F *、L h Ji、wnyr- | . f打|11 i k

11、 i *.i b-!-! ! - - yi*:!|-|J|Jlnil|faLLILfa-L-|LlX.h*i lUUl 4 XXX . /；11 j i J j 11 , j , i r, i. . i, , i,sk . . n I. L 、 a11f f f J! f J 1I d I i I L I I H I fa * LILI、*%*、m y . j |- |- J I : I ,1 4 ( , I I, I I -1 I L -I ! L L * L Lb 1. X k V/ H Jri rJ-d J I .1UI 1kI n I*iL Lk X h * zp* mrrir i i

12、 b d i. i h h i h i卄.i /t dfwjr.1 . -i L I. -i kO i i -iHir：-(2) q1x=1,q1y=1,q 1=10e-9; q1x=-1,q1y=-1,q2=-4e-9(3) q1x=1,q1y=1,q1=100e-9; q1x=-1,q1y=-1,q2=100e-9场与等楼线5-1023三、开放性试验p=汁_(1_丄)4n0r1r24n0n r2riKx-qx)2(y -qy)2用二项式展开，又有 rd，得A 二 rdcosra 二 rdco s2 2p=qd,表示电偶极矩（dipole moment），方向由-q 指向 +q。 等位”2丄

13、(2匚rdcoEri所以厂普2一4n0r4n0r画出电偶极子的等位线和电力线(rd )线方程（球坐标系）：场与等位线二-崖：I:1si-505xr =C、cosdr _ rdE?Epq3(2cos versin v ej4氓r将 Ee和 Er 代入 E 线方程有r = D sin2Q1x=1Q1y=2 Q 仁 10Q2x=1Q2y=-2Q2=10场与等位线Q1x=1 Q1y=2 Q 仁 10Q2x=1 Q2y=-2 Q2=-10r ff r jrj* - J-ll ll-mr- L_rrL_r_T-ra* l!”/-AfcF* 打 a*J尸|SV9rarFVMWWYHHHM7,TJ-r FrF

14、* .-.-_-.-i_-.-.-_-.-_-.-_-n_r-j-.-5!昶WJ.- 总 4 宀71.亠AJK*AJ.7 _ i:EEI_ ”“T .7仝j* . l mrLf g_rTL-Fnrr -,亠w ja-_XV- - _ rhWWAWVAl h-J f-ww孑X.W 5-5531 _sl.xc-l-诧讥铁彊岭枣蔡-11tcjtr*董K釁姿釁崔*盐恥鳌董石X令i沁奈至說5g 谯濺-:.vj-;-n-rc - w总M枣弐时J5 Ke龍- 4A %Lm% Lcl说診5翘洪捡買席期的秒”+毬轻“ 匸-S?;.mc.h.-.hnk-.*- n c ” ，很幣isg洛r溟虫覚宴戻壮；：翌卫-_

15、.* 2*2 -,-=2嚴煞扯敢尋芒一卞瓷養区:览佯昭毬迟迟 2-一、irLJ-fL 5S-iL-X誥逼盖裟-、， ! i I ! .1 dL L住蹩矍養養;: I : I : ; r I :L L L321” 0寸JJm-if* SSS? JJ-S5 J? J-.-lx.-ry .J*X_ J- -,a-uXF -. Jx J)-B#y一亦娄隽抵訂垃恣铸十密裟I场与等位线5Q1x=1 Q2y=2 Q 仁 10 Q2x=-1 Q2y=-2 Q2=1053212实验二 MATLAB 电磁场有限元计算实验目的：4、了解有限元算法的原理，熟练运用 MATLAB 环境的 PDE 工具5、熟练运用 PD

16、E 工具分析简单的电磁场边值问题。实验内容:在电磁场的计算中，仅对那些具有最简单边界条件和场域几何形状规则的 问题才有解析解，多数问题的求解必须用数值计算的方法，其场域分布的数值计 算内容是学习难点。本实验将有限元法和Matlab 结合起来对电磁场教学中的电位分布问题进行计算。结果表明使用Matlab 对有限元分析编程中的矩阵进行处有限元法是以变分原理和剖分插值为基础的一种数值计算方法，其基本思想有限元简介(Nonodeelement（单元）理，程序设计清晰简便，易于理解和实是将场域方程等价为一个条件变分问题，然后由条件变分问题对场域进行剖分离 散为方程组进行求解。对于一个电场来说，其储能总是

17、趋于最小，这样变分法的泛 函和电场的储能就联系起来了。对于边界为 L 的无源空气介质二维静电场中 厂 个圭寸闭场域 S 内的等价能量泛函可以写为：在有限元分析中，将所研究的区域 S 划分成有限的 n 个三角形网格单元。 对应m 个节点，ds 为单元 e 的面积。对任意三角形单元 e 中任一点的电位可以认 为由该三角形的三个节点（分别设为 i、j、k）上的电位 u 随该点坐标 x、y 变 化而线性确定。因此，对于单元 e 构造插值函数：Ue = CtU i + Ct U + OkEZk =力Oh ZZh ft= Lih其中 ah 称为形状函数。那么有插值函数的一阶偏导数为：V*矗加L亍站hTjr

18、k出H Obh2 ij,i从而得到能量函数 We:+（爲訓）M则将单元 e 中的能量函数 We 对每一个节点电位 ul （ I = i, j , k） 求一阶偏导数 dw.（u）_ r（hii.irixi mi 1 z RT_表示为矩阵形式有：然后进行总体合成，将各单元的能量函数对同一节点的电位一阶偏导数相加得所要求解的线性方程组。由以上分析，可知在该场域内电场有限元数学模型为：IKJ二（）式中 U 为 n 个节点处的待求电位，K 为 n 阶矩阵。最后进行强加边界条件处理 消去已知电位节点在系数矩阵中所在的行和列，得到简化后的方程，继而可以对 电位进行求解。流程框图如下图所示。LkllAln(

19、) *V4=0-A ii 1/Limi-三IIr4，得:()fl 01+it=Tj,k巾7*D= pVx/= J7X =().VZ?= 0静态场满足上方基本方程，式中D 为电位移，为电荷密度，H 为磁场强度，J为电流密度，E 为电场强度，B 为磁感应强度.对于恒定的电场：式中电位满足泊松(Poisson 方程：V 二利用上述方程，再加上边界条件，利用 Matlab 中的偏微分工具箱，即可求解带 电体周围空间的电场分布.输入 pdetool 可进入软件环境。、静电场仿真对于不存在电荷的空间部分有电荷体密度为零方程：,上式退化为拉普拉斯(Lap lace)幵始两点电荷的电场：两等值异号点电荷单位，

20、两者间距为 1,求其电势分布整个求解域取中心为原点，半径为 2 的圆,两空间电荷点位置为（-0.5,0）和（0.5,0）作为一种近似，画一个尽量小的圆,取半径为 0.05.大圆的边界条件是 Di richlet 边界条 件，取 h= 1, r= 0,这种做法是模拟远处的电势为零.由于大圆与小圆之间的区域没 有电荷，满足 Laplace 方程，因此在选择方程时选取 Elliptic（椭圆）方程，其方程类 型为：一 V （ c V I ） +au = / *取系数为 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,表示求解区域为二者之差.矩形所 有的