5用matlab解决电磁学中的电场问题

/用MAＴＬAB解决电磁学中电场问题摘要:二十一世纪以来，随着计算机技术的进一步提高以及电脑在国内的普遍使用，在很多教学领域运用计算机软件来辅助教学已经非常普遍，然而用电脑软件处理仿真电学及电磁学中电场及电势的研究却并不常见。在进行电场及电势问题的研究时，电场虽然是确实存在的,但是由于其抽象而不易被理解,而且即使在实验室中,我们也很难实现其理想化模型。这样就使得教师在教学过程中不能生动且直观的描述出来,学生在学习过程中也将很难理解和接受。本文利用计算机数学软件MAＴLＡＢ模拟仿真真空中自由电荷,电偶极子，带电细棒，尖端导体等附近静电场电场线及等势线分布图型，从而使其更加形象，生动，直观，更加便于学生学习,理解和接受，同时使的教师教学更加方便，快捷。关键字:MAＴLAB软件，静电场,尖端导体,尖端效应，电场强度,电势。目录1引言·····························································12用ＭAＴLAB处理静电场中的电场线和等势线的问题············1２.1用MATLAB模拟仿真真空中点电荷的电场线和等势线分布········1２．2用ＭATLＡＢ模拟仿真真空两个点电荷的电场线及等势线分布·······32.2.1用ＭATLAB模拟真空中两个同种点电荷的电场线及等势线分布·····３2。2.2用MATLＡＢ模拟真空中两个同号但不等量点电荷的电场线和电势分布··52.2.3用ＭATＬAＢ模拟电偶极子的电场线及等势线分布··································································63用ＭＡＴLAＢ模拟均匀带电细棒的电和分布情况···························8４用ＭATLAB模拟两个无限长导线的电位和电场分布······················９５尖端导体附近的电场及电势特点及其应用·····························１1５．1电场函数····················································11６总结·····························································137参考文献·························································1４１引言:二十一世纪以来，随着计算机技术的进一步提高以及电脑在国内的普遍使用,在很多教学领域运用计算机软件来辅助教学已经非常普遍，然而用电脑软件处理仿真电学及电磁学中电场及电势的研究却并不常见。我们在研究电磁学中的电场问题时，通常将自由电荷看作一个理想化的模型。电场虽然存在，但是人们是不能直接观察到的,并且这些理想化模型很难在实验室建立起来，这样就会显得空洞乏味，使得在教学的过程老师授课难以到位,学生在学习过程中不容易理解和接受。本文主要运用计算机数学软件ＭATLAＢ来模拟真空中自由电荷、带电细棒，电偶极子及尖端导体等产生的静电场在空间中的分布，并通过二维和三维图形模拟电场和电势的分布图像。我们把这些理想化的模型和抽象的概念通过计算机软件实现，既能使学生容易理解,也可以启发他们在学习过程中充分运用计算机软件资源,使得学习效率倍速提高。MATLAB作为一种科学高效的计算软件，其主要作用是进行矩阵的运算,并且还可以进行信息领域的处理。它比起其他计算机软件使用更加方便，输入更加简洁，运算更加高效,内容更加丰富，并且很容易由用户自行扩展.通过对MATLAB软件的运用，可以在学习过程中具有更高效率、更加富有创造性的进行科学计算２用MAＴLＡB处理静电场中的电场线和等势线的问题在电磁学中,点电荷，电场强度，以及点电荷在真空中电场的分布等,都是电磁学中重要的内容。在处理这些问题时我们一定要做到形象直观，使学生更易理解,进入二十一世纪以来,随着计算机技术的提高，作为二十一世纪的新型人才，我们必不可少的要用到一些计算机软件。下面我们就用MATLAB程序将这些抽象的概念模拟成形象的图像，清晰地展现在大家面前，以方便学习和交流。２。1用MATLAB模拟仿真真空中点电荷的电场线和等势线分布为了研究以及学习方便，我们假设真空中存在一点电荷，其所带电荷量为Ｑ，并规定其初始位置为坐标原点，从而建立平面直角坐标系,用数学软件ＭAＴLＡB模拟仿真出点电荷产生的电场线和等势线在真空中的分布图像。图1:真空中点电荷的电场线和等势线分布图2。2用MATＬＡB模拟仿真真空两个点电荷的电场线及等势线分布研究完真空中的点电荷后，我们不禁会想到如果是两个点电荷呢？他们的电场分布又会是何种情景，对于两个带等量同号的点电荷，由于它们的相对位置，电荷量的分布等因素的不同,用几何作图的方法描述出它们的电场线和等势线分布的立体图，就显得有一定困难。而通过运用数学软件MATＬAB,只要我们知道这两个点电荷的距离，明确它们所带的电荷量，就可以模拟出它们的电场线和等势线分布图像。使得我们更加形象直观的认识到电场及电势现象，从而在学习中更易理解和接受。2.２．１用MATLAB模拟真空中两个等量同号点电荷的电场线和电势分布根据我们高中及大学期间所学电磁学知识，如果我们知道点电荷在真空中所产生的电势分布为Ｖ,那么电场强度的大小为电势梯度的相反数，即：（1）为了研究的方便，我们往往规定真空中无穷远的电势为零,则它们所发出的的电场中，空间的电势分布为：(2）通过计算电势梯度，便可以得到两个点电荷的电场分布.在运用数学软件模拟两个点电荷电场线和等势线时，先运用指令mesｈｇrid（x，y)，它的作用是把分布在xy平面内的,一个区域内的所有场点都一一表示出来。设两个点电荷分别位于距离原点1ｍ处,为了有较好的直观效果，这里取库伦.图2：等量同号点电荷的电场线和等势线分布图2．2。2用MATLAB模拟真空中两个同号但不等量点电荷的电场线和电势分布由我们所学过的知识可知，同号但不等量电荷的电场和电势公式与等量同号电荷相同，下面我们就模拟同号但不等量点电荷的电场线和等势线分布。图3:等量同号点电荷的电场线和等势线分布图2.２.３用MATLAB模拟仿真电偶极子的电场线及等势线分布既然我们在前文已经做了同种电荷的电场及电势分布，那么对与异种电荷呢?在这里我们不妨再做一下它们的电场及电势分布图，首先我们假设这对异种电荷所带电荷量相等，并且处于真空中,此时我们习惯上叫这对异种电荷为电偶极子，下面我们就用数学软件MATLAB模拟这对电偶极子的电场线和等势线在真空中的分布情况。电偶极子的带电量为,正电荷的位置为（2,0)，负电荷的位置为(—2，0)。通过函数指令就能实现。图4：电偶极子的电场线和等势线分布图3用MＡTＬAB模拟均匀带电细棒的电场线和等势线分布情况假设电荷Q均匀分布在长为L的细棒上,求其在真空中的电场分布情况。根据我们学过的有关电场的知识,在真空中，对于均匀带电的细棒，根据电场的叠加原理,我们可以在一些特定位置通过数学积分求出电场强度的解析解。在细棒延长线上的电场强度为：（3)在细棒的垂直平分线上，与棒的距离为ｒ处的场强为:(4）然而对于其它位置的电场强度的求解，由于不具有对称性，很难得到解析解,这时我们则必须借助于计算机,将所得数字通过图形直观表达出来。图５：带电细棒的电场线和等势线分布图4用MATＬAB模拟两个无限长导线的电位及电力分布假设两根电荷线密度为,与ｚ轴平行，分别置于,,的平行无限长带电导线,那么其电位及电力分布可用一下程序实现:图6:无限长平行导线的电位分布图图7：无限长平行带电体电力线分布5尖端导体附近的电场特征及其应用５．1电场函数如下图所示:这是一尖端导体的示意图图８:尖端导体它的电场用函数可表示为如下三式：（７）（8）（9)其中,底面和顶点的距离为r，为半顶角,，。为了方便计算，我们把此尖端导体的电势假设为零，从而得出上式，式中Ａ为待定常量。对于尖端导体而言，我们认为,即。通过分析推到，得到上面三式运用MATLAＢ绘图结果，如图所示。图９：尖端导体的电场线和等势线分布图6总结:运用MATLＡB语言形象直观的模拟电磁学中的静电场问题，不但给学生的学习,老师的教学都带来很多的方便,而且将物理学中的抽象问题数字化、形象化，直观化,使人们更容易接受。通过对静电场问题的研究,用ＭATLAB解决物理中力、热、光、电，磁以及近代物理中的一些问题，对于学习效率的提高，科学研究方法的探寻，以及科学工作都有相当大的益处.进入新世纪以来,随着计算机技术的飞速发展，不再是古人所说隔行如隔山,而是各种学科之间的互相交融，因此对大学生的培养也提出了更加苛刻但却积极的要求。学生不仅要掌握本学科的知识,还需要有广博的知识视野和先进的学习工具载体，而其中计算机技术在理论课学习中的运用,在素质教育中起着举足轻重的作用。因此,计算机软件模拟技术将会很快运用于现代化教学,从而使教学过程更加方便快捷，教学效果更加形象直观，进一步提高我国现代化教育进程，做到与国际接轨。致谢：在这里，我首先向各位领导及老师致以崇高的敬意和深深的感谢。大学四年,我明白了许多道理，学会了许多知识,最重要的是我不再是一个只知道学习的呆子,我知道了如何学以致用。作为一名普通大学生，我完成了我应该完成的学业,我想说我爱你，大同大学.这篇论文是我的处女作，肯定会有许多不足之处,希望各位老师给予批评指正，敬礼!参考文献梁灿彬。《电磁学》，北京高等教育出版社，2０04马涛，《数字化大学物理》浙江，浙江大学出版社,2０08张德丰,《ＭATLAＢ数值计算方法》北京,机械工业出版社。2010刘定兴，胡先权，尖端导体表面附近的静电场研究【N】，重庆三峡学院学报，22（2006)10４-106张志涌,ｍatlａb教程，北京航空航天大学出版社，20１0，97—98（６）Ｃ.Zhｏng，Ｗ.B。Hu，Ｙ。F.Cheng，Ｏnｔhｅessentiａｌrｏleoｆcuｒreｎtdｅnsityinｅｌｅcｔrｏcataｌｙｔicａcｔivityoftheelｅctrｏdeposｉtedｐlatinumforoxidａtiｏnofａmmoｎia，Ｊ.PowerＳouｒceｓ,19６（２0１1）8064－80７２．陈怀深.《MATＬAB及其在理工课程中的应用指南》,西安电子科技大学出版社程守株。《普通物理学》,高等教育出版社苏金明,张莲花,刘波.MATLAB工具箱应用，北京,电子工业出版社，2002MＡTLＡB。help.美国。Ｍaｔh。work公司,2004Abｓtｒact:inthｅ21stcentury,alongwiththeｃoｍputerteｃｈnｏlｏgytoｆurthｅrimｐｒoｖeandthewidｅspreadｕseofｃompuｔeｒsathome,inaloｔofｔeａｃｈｉnguｓingｃomｐuｔersｏftwaretoａssｉstｉｎteaｃhinｇisveｒyｃomｍon,butｗithaｃｏmputｅrsoftwareinｔhｅｓｉmuｌaｔｉoｎelecｔricitｙandelecｔroｍaｇnetiｓmｅｌｅctrｉcｆieldandeｌeｃtricpoｔｅntiａｌoｆｔhereｓｅarcharｅnotcommon。Ｉｎthｅsｔudyｏfｅlectriｃｆieldanｄｐotentialprobleｍs,ａｌthoｕghｅlectricfielｄisreal，ｂｕtbecausｅofiｔsabｓtrａctandｄｉｆficｕlｔtｏｂeｕｎderstood，aｎｄeｖenｉntｈelab,weａlsodifficｕｌttｏrｅａｌｉzeitｓｉdealmodel.Ｔhisaｌlowsteaｃherｓiｎtheｔeａchingprocesstobevｉvidanｄintuｉｔivedescriptｉon，studｅntｓｉntheｌearｎｉｎgprocesswilｌbediｆficulttoundｅｒstaｎdａｎdａccepｔ.Thispapｅr，byusｉnｇmathemａticaｌｓofｔwareMＡＴLAＢcｏｍpuｔｅrsimulatｉoninavaｃuｕmfreｅcｈaｒge,electｒiｃdipole，ａｃhaｒgeｄthiｎrｏds,cuttｉng—edgeconductorsｓuchasｅlｅｃtrｉcfieldlinesanｄequipｏtｅntｉallinesnearｔｈｅelｅctrｏｓtaticfiｅlｄdistributionpaｔtｅrn，mａｋｉｎgitmｏreimage,vｉvid,iｎtｕitive，mｏreｆaciliｔatestudｅntｌeａｒniｎg，ｕndeｒstandaｎdacｃｅpt,atthesameｔｉｍｅ,ｍakethetｅachingmｏrecoｎvenｉｅnt,qｕicｋ．Ｋｅywords:ＭＡTLABsoｆｔware,tｈｅｅlｅctｒｏsｔaticfiｅlｄ，cuｔtinｇ-ｅdgｅｃoｎdｕctｏｒ，edｇｅeffecｔ，ｔheｅlectricfieｌdiｎｔeｎsity，eleｃtricpｏｔｅｎtial附录:图一单个点电荷程序如下：clｅａr;Ｅ0=8。85e－1２；c0=1/４/pi/E0；q=1.6＊10^(－1９）；xm=2.５；ym=2；x=linspaｃｅ(-xｍ，ｘm);ｙ=linspａcｅ（-ym，ym);[Ｘ，Y］=meshgrid(ｘ，y）；R=sqrt(Ｘ.^2+Y．^２）;U＝c0＊q．／R;ｕ＝１e—9：0.5e—9:5e—9；ｆｉguｒe（１)coｎtoｕr(X，Y,Ｕ,u）holｄonplot(０，0，’ｏ'，'mａｒkｅｒsize＇，１2)axisequalaｘistｉghtｔｉtlｅ（'单个点电荷的平面电场线与等势线',’fontsiｚe＇,12）；ｘlabｅｌ（’r'，’fontｓize'，12）;yｌabeｌ（'Ｅ(U)’,＇fontｓizｅ＇,１2）;图二等量同号点电荷程序如下：ｃｌear；clc;cｌoseall；E0=8.８5e—1２；c０=1／４／pi/E0；q=1。６*10^(-1９）;a=1；xm=２．5;ym=２;ｘ=liｎspaｃｅ(—xm,xm）；y＝linｓpacｅ(－yｍ,yｍ)；[X,Y]=mｅsｈｇriｄ(x，y）；Ｒ１=sqrt（（Ｘ+1）.^2+Ｙ.^2);R2=sqrｔ（(X－1）。^2+Y．^2）;U＝c0＊q./R1+c０＊q./R２；ｕ=1ｅ—9：0．5ｅ-９:5e－9;ｆｉgｕre(1)ｃｏｎtour(Ｘ，Ｙ,Ｕ,u）ｇridonlｅgeｎｄ（num2str(u’））ｈｏldonplｏt（—1,0,'o','ｍaｒkersize’，12）plot（1，０，＇o’，＇maｒｋｅrsｉzｅ',1２）axisｅqｕａｌｔightｔitｌe（＇等量同号点电荷的电场线与等势线＇,＇foｎtsｉｚe',1２）xlabｅｌ（'r’，'ｆontsize'，1２)ylaｂel（’E(u)’，’fｏntｓｉｚｅ’,１２)txt=[ｎuｍ2str（a)］；texｔ（－ｘｍ，—ｙm—０。3，txt，’fontsiｚｅ’，6）;图三同号但不等量点电荷程序如下clｅａｒ;clｃ；closeaｌl；Ｅ0=8。85ｅ－12;c0=1／4/pi/E0；q1＝1.6＊10^（—１9)；q2=2*ｑ1；a=1;xm=２．5；ym=2;x=linspacｅ(-xm，xm);ｙ=ｌinsｐace（-yｍ,ym)；［X，Y］=ｍｅsｈgriｄ（x,y）;R1=ｓqｒt(（X+1）.^2＋Ｙ．^２）;R2=sqrt(（X—１)．^２+Ｙ.^２)；Ｕ=c0＊q1．/R1+c0＊q2。／R2；u=１e—９:0．５ｅ－9：5e—９;figuｒe（1)contouｒ(Ｘ，Y，Ｕ,u）gridoｎlegeｎd(num2str（u'）)ｈoldｏnplｏｔ（—1,0，＇ｏ'，’mａrkeｒsｉｚe’，12)plot（1，0，＇o',’marｋｅrsizｅ’，12)axisｅｑuａlｔightｔitle('同号但不等量点电荷的电场线与等势线’,'foｎtsｉzｅ’,１２）ｘlａbel('ｒ'，'ｆｏntｓｉze',12）ylabel('E（u）’,＇fontsｉze＇，1２)tｘt=［num2sｔr(a）］；teｘｔ（—ｘm，－ym－０.3，txt,'ｆｏｎtsiｚe’,6);图四电偶极子的电场线及等势线程序如下cｌeａr;clf；ｑ=２ｅ－6；k=9e9；a=2。0；b＝0；x=—６：0.6:６；y＝x；［X，Ｙ]＝meshgrid(x,ｙ)；rｐ＝sqrt(（X—ａ）。^２＋（Y-b）．^2）；rm=sqrt（（X+a）。＾2+（Y+ｂ)。^２)；Ｖ=q＊k*（１./ｒp－１。/rｍ)；［Ex,Ｅy］=ｇradient(—V）;ＡE＝sqｒt（Ex。＾2+Ey.＾2)；Ex=Ex。/ＡE;Ey＝Ey.／AＥ；ｃv=ｌｉnsｐａｃe（miｎ（miｎ(V)）,max(ｍax（V)），５1）;contｏuｒ（Ｘ,Y，V,ｃv,＇r-')axis(’sｑｕａrｅ＇)title(’电偶极子的电场线与等势线＇，'fontsize'，12),hｏldｏnquiver(X，Y，Ｅｘ，Ey,0。６,’g’）ｐlot(a,b,＇bｏ'，a，b，＇ｇ+＇）plot（—a，—ｂ，'bo',－a,—b,’w-’)xlabel（'x')；ｙlａbｅｌ（'y’）;holdoff图5：带电细棒的电场线和等势线分布图程序如下ｃleaｒalllam＝le-9; ﻩep0=8．85*ｌｅ—1２；ﻩ c0＝lam/（4*pi*ep0)；ﻩﻩ％归并常数Lh=3； ﻩｘ＝-6。5：0.11：６.5;y=—5.5:0.１1：5．5；l=-Lｈ:０.1：Lh；［X，Y，Ｌ］＝ｍeshgrid(x，y,ｌ）;r=sｑrt（(Ｙ-l).^2+x。^2);dv＝ｃ0．/r；v=ｐi/４0*trａpz(dｖ,3); %求电势[Ex，Ey］=gradieｎt（-v,0。2）；ﻩ %求电场figureaxis(［-6，6，—5，5]）；L=lｉne（[0，0］,［-３,3］，'coｌor’，'ｒ’，’lineｓtylｅ’,'—'，＇ｌｉneｗｉｄth’，5。５）；ﻩﻩ％画带电棒holdoncoｎtour（X(：,:，1）,Y(：，:，１),v,［６，8，10,１2,１4，16，１8，20,２２,２４，２６,2８，３0,32］，’g’)hoｌdonsｘ=0．2；ｓy=[-３。２:，0．４：3。２];［Sx，Sｙ]＝meshgｒid（sx，sy）； ﻩ％计算电场线起点streａmｌｉne（X(:，：，1），