实验八 点电荷电场分布模拟MATLAB程序

matlab点电荷的电势和电场解析式【原创实用版】目录一、引言二、点电荷的电势和电场解析式1.电势解析式2.电场强度解析式三、MATLAB 绘制点电荷的电势和电场分布图1.建立数据网格2.计算电势和电场强度3.绘制等势线和电场线四、结论正文一、引言在物理学中，点电荷是一种理想化的模型，用于研究静电场和电荷分布等问题。

对于点电荷，我们可以通过解析式来描述其电势和电场分布。

在本文中，我们将讨论如何使用 MATLAB 绘制点电荷的电势和电场分布图。

二、点电荷的电势和电场解析式1.电势解析式点电荷的电势解析式是根据库仑定律和电势定义推导得到的。

对于一个位于空间中任意位置的点电荷，其电势可以表示为：V(r) = k * q / r其中，V(r) 表示电势，k 为库仑常数，q 为电荷量，r 为距离点电荷的距离。

2.电场强度解析式根据电势的定义，电场强度与电势的关系为：E = -dV/dr对上式进行求导，可得电场强度的解析式为：E(r) = k * q / r^2其中，E(r) 表示电场强度。

三、MATLAB 绘制点电荷的电势和电场分布图1.建立数据网格为了绘制点电荷的电势和电场分布图，首先需要建立一个数据网格，用于表示空间中各点的坐标。

在 MATLAB 中，可以使用 meshgrid 函数建立数据网格。

例如，对于一个位于 x 轴上范围为 -5 到 5，y 轴上范围为 -5 到 5 的数据网格，可以编写如下代码：x = -5:0.1:5;y = -5:0.1:5;[x, y] = meshgrid(x, y);2.计算电势和电场强度利用电势和电场强度的解析式，可以计算数据网格中各点的电势和电场强度。

在 MATLAB 中，可以使用矩阵运算完成这一任务。

例如，对于一个位于 (x, y) 位置的点电荷，可以计算其电势和电场强度如下：z = k * q * (x^2 + y^2)^(-1/2);E = k * q / (x^2 + y^2)^(3/2);3.绘制等势线和电场线在计算出各点的电势和电场强度后，可以使用 MATLAB 的 plot3 函数绘制等势线和电场线。

龙源期刊网
基于MATLAB的点电荷的静电场模拟
作者：王文涛梅宇航
来源：《中学物理·高中》2013年第02期
1前言
静电场是继引力场之后，教学中又一次面临的有关“场”的问题，处于承上（力学）启下（电磁学）的重要位置，但由于静电场概念理解起来比较抽象，加之描述电场特性的电场强度，电势等物理量需要较强的空间想象能力和抽象思维能力，所以本节一直是教学的难点，MATLAB作为一种数学软件，计算命令丰富，作图功能强大，通过编程对静电场的模拟，可以提供给学生一个直观的物理图景，帮助学生正确理解静电场的概念和特性，还可以逐步引导学生加深对场的物质性认识.
2数学模型的建立
3.2等量异种点电荷的电场线与等势线
对程序适当修改，可得到等量异种点电荷的电场线和等势线分布图，见图3.
由图3可以看出，等量异种点电荷连线的中垂线是一等势线，所以沿着中垂线移动电荷到无穷远处，电场力不做功，因此若取无穷远处电势为零，则中垂线上各点的电势也为零.
3.3三维电势图
由图2、3可以很好的理解平面上点电荷的电场线和等势线的分布规律，但是要想形象的理解点电荷电势变化规律就需要借助于MATLAB的三维仿真功能了.
4结束语
本文首先建立一对电荷产生的电势和电场强度的数学模型，接着运用MATLAB绘图指令，画出了电场线和等势线的平面分布图以及电势的三维分布图，由以上分析不难看出，通过MATLAB编程模拟静电场，不受实验条件的限制，得到的图像细致直观，而且实验参数可以根据实验结果的需要随时调整，这在实际教学中对于帮助学生快速掌握静电场规律和提高学习兴趣有积极作用.。

利用MATLAB 模拟点电荷电势的分布一、目的1.熟悉单个点电荷及一对点电荷的电势分布情况；2.学会使用MATLAB 进行数值计算，并绘出相应的图形；二、原理根据库仑定律：在真空中，两个静止点电荷之间的作用力与这两个电荷的电量乘积成正比，与它们之间距离的平方成反比，作用力的方向在两个电荷的连线上，两电荷同号为斥力，异号为吸引力，它们之间的力F 满足：R RQ Q k F 221=（式1） 由电场强度E 的定义可知：R RkQ E 2=（式2） 对于点电荷，根据场论基础中的定义，有势场E 的势函数为R kQ U =（式3） 在MATLAB 中，由以上公式算出各点的电势U ，可以用MATLAB 自带的库函数绘出相应的电势分布情况。

三、MATLAB 基本语法(一)标识符与数标识符是标志变量名、常量名、函数名和文件名的字符串的总称。

(二)矩阵及其元素的赋值赋值就是把数赋予代表常量或变量的标识符。

MATLAB 中的变量或常量都代表矩阵，标量应看作1×1价的矩阵。

赋值语句的一般形式为变量=表达式（或数）列如，输入语句a=[1 2 3;4 5 6;7 8 9]则显示结果为a=1 2 34 5 67 8 9输入 x=[1 2 3 4 5 6 7 8 9]结果为x= 1 2 3 4 5 6 7 8 9可以看出，矩阵的值放在方括号中，同一行中各元素之间以逗号或空格分开，不同行的元素以分号隔开。

语句的结尾可用回车或逗号“，”，此时会立即显示运算结果；如果不希望显示结果，就以分号“；”结尾再回车，此时运算仍然执行，只是不作显示。

变量的元素用圆括号“（）”中的数字（也称为下标）来注明，一维矩阵（也称数组）中的元素用一个下标表示，二维矩阵可有两个下标数，以逗号分开。

在MATLAB中可以单独给元素赋值，例如，a(2,3)=6,x(2)=2等。

(三)元素群运算把n×m矩阵中的每个元素当作对象，成群地执行某种运算，称为元素群运算。

基于Matlab模拟点电荷电场线和等势线1. 引言1.1 背景介绍电场理论是物理学中的重要概念，描述了在空间中存在的电荷所产生的相互作用力。

点电荷模型是电场研究中常用的简化模型，通过模拟点电荷的分布和运动，可以很好地描述电场的特性。

在现实生活中，我们经常会遇到点电荷电场的问题，比如电荷在空间中的分布及其对周围环境的影响。

基于Matlab的数值模拟方法可以帮助我们更好地理解电场的特性。

通过模拟点电荷的分布情况，我们可以绘制出电场线和等势线，从而直观地展示电场的分布情况和强度。

这不仅有助于理论研究，还可以在工程实践中提供重要参考。

通过基于Matlab的点电荷电场线和等势线模拟，我们可以更深入地探讨电场的性质，为相关领域的研究和应用提供支持和指导。

【字数：205】1.2 研究意义电场是物理学中非常重要的概念之一，它描述了空间中各点所受电荷作用力的性质。

而点电荷则是电荷密度在空间中极小的模型，通过研究点电荷的电场线和等势线的分布情况，可以帮助我们更好地理解电场的性质和规律。

基于Matlab进行点电荷电场线和等势线的模拟，不仅可以直观地展示电场和电势在空间中的分布情况，还可以通过调整参数来研究不同条件下电场和电势的变化规律。

研究点电荷电场线和等势线的分布对于学术研究和工程应用具有重要意义。

在学术研究方面，通过对电场线和等势线的模拟分析，可以深入探讨电场的特性和规律，进一步推动电磁学理论的发展。

在工程应用方面，电场线和等势线的模拟可以帮助工程师设计和优化电子元件、电路和传感器等设备，从而提高其性能和稳定性。

深入研究基于Matlab模拟点电荷电场线和等势线的方法和应用具有重要的理论和实际意义。

1.3 研究目的研究目的是在Matlab环境下通过模拟点电荷的电场线和等势线，深入探讨电荷在空间中产生的电场分布情况，以及不同点电荷配置对电场线和等势线的影响。

通过研究电场线和等势线的形态和分布规律，可以更好地理解电荷之间的作用关系，为进一步研究静电场提供依据。

matlab点电荷的电势和电场解析式
一、引言
在电磁学中，点电荷的电势和电场分析是基础内容。

Matlab作为一种数学计算工具，可以方便地帮助我们求解这些问题。

本文将介绍如何利用Matlab求解点电荷的电势和电场，以及它们之间的关系。

二、Matlab中点电荷电势的解析式
1.点电荷电势公式
点电荷电势公式为：
φ= k * Q / r
其中，φ表示电势，k为库仑常数，Q为电荷量，r为距离点电荷的距离。

2.Matlab代码实现
```Matlab
clear;
phi = 1/4 * pi * epsilon_0 * Q / r;
```
三、Matlab中点电荷电场的解析式
1.点电荷电场公式
点电荷电场公式为：
E = k * Q / r^2
其中，E表示电场强度，k为库仑常数，Q为电荷量，r为距离点电荷的距离。

2.Matlab代码实现
```Matlab
clear;
E = k * Q / (r.^2);
```
四、点电荷电势和电场的关系
电势是描述电场在空间中的分布情况，而电场强度则表示单位正电荷所受到的力。

点电荷的电势和电场强度之间的关系为：
φ= - gradient(E)
即电势是电场强度的负梯度。

五、结论
通过Matlab编程，我们可以方便地求解点电荷的电势和电场。

掌握这些解析式，有助于我们在实际问题中分析电场分布和电势变化，为后续的电磁学学习打下基础。

利用Mat lab損拟点电荷电场的分布一・实验目的：1. 烬思融个点电命及时点电&的电场分布愴况i2. 儒会便HI 计卸.并绘出Hl应的图移二・实验原理：眾厳冷伦;口人作何空中.曲个»itA电尙Z何的作用力与这构个电荷的电fit蔡枳成正It.弓它的平方谥反比.作用力的方向金电倚的连段1・曲电斥力.wy W力.它们2何的力$滑足*4式U山电场誉咬［的ill文顼知*(式2)<1 TA电荷.根卅场论垩的中的迄义.<1的场［的的晦数为（/•学R(A 3)向 E.-0U d（i M4lUt> P.由以上公式W Hl ft AM电钓U・电场新唱（右.可以用Malhb门谐的相应电荷的电场分衛情况.三.实匕内容1. ■草个点电背的平■电场线9等勞纽尊祈线就乂以电荷为中心・用MalUb価零铃歿电加曲札鼎电力用3 为k・9・t••电St可取为q・“g 般大的*勢銭的Y径凶逐比射线的丫栓小 A. r^Ql.H电势为屿二丄%・如果从外到中茶等野线.MVlfi的邯针找的电5迄*外面的护乩騒么缶*饯的电紡用向吊丧不切—亦刑“（1以7）・％・从"判巾丸偶数个点.RtaiooV点.传嵐中心点的生轿慢ilo・/点的坐杯町用向IB灰示I x./imparr(-j;.G.IOO).在血fl!樂标系中町形阪期悟世标：［儿町二林心皿(町・*点到廩点的为：F二儿八2・丫厂2・fiMaUA中进行喉方运"时・桑方号曲面更加点.戏示对交■中的元It透务彙方计算・备点的电势为(/“S "同什饱.住进h»iAizi»W.聲号前面也"加点.冋什住不时变鍛中的兀素进打除決运A用等矗线命令出帑勞线. 節图谕EKRWtaF：■■个迄电"0・2】■■••icr 肌■比■常■q・1.6・W“” Qit电•电■ rO-O.l;■电场纽g戊丫怜thota-llnspacetO^^^plUS)； [x9y]-pcl2cart(th«ta fl aU x>lxj0.05«x]j y-(y;0.05e y：； quiwr<Mry.O.S*x.0.5*yI plotlx«y) hold on u-k*q/rO|ul-lm&p»c4( X v3*7)*u;x-Lln5pAC«(-0.1>0.19100)| |X,YI “・*hgr idf M);rX-Bqrt(x.e2»Y.M2>;U-k.•q-/rl;contourfX^Y^U.ulI电背馆丫血电场Mft*, v fontBixeS20l>U^bS xl«fc*ll*r\*font*iie\lS>tU>b«*kyUbell •t<U)\t font»ite\16l2. Hi 一对走电債的平Ifc 电场嫂与羚毎绘 程序代刑如Fi电&林的电场絃和线■电■比〈焼•!小曲电■比点电價H 釣电址线和*铃銀只鬲占* qgtUM 》x-Ue>sp4C«(-x».xa); y-lin»p4c«( -yw.ya)： !X«Y)-TC9hgrld<x r y>2 Rl-3qrtHX«l>.*2<Y.*2)； R2-flqrtllX-l>.A2<Y.*2); U-l./RUq./R2; u ・l:0.5:4; figurecoAtourIX,Y r U e uigrid on l«q«nd(nuB^str1u*)> bold on plol<|-xjT>;xn}. *0;01» ploKIOrOUI ywuynH plot<-l«Q» *o*,^Kark^rStx*4 ・12) pl^Kl.O. e o*»<Nerk«rSia«* «12>tEx,IyJ^radl«nt(-U f x(2>-xm <y<2>-y<ll>MR1 电付 H 反欢第・的卿个分・ dehl-20| ・4垃电场纽角用・(■thl-<dthl ：dthl ：ie0-dthl)*pl/160; ♦电f 的 rO-O.U«l-rO-c© ・2bl >-l;Q 电场线的■堡标■电场4的q-1; xr>2«5; 眄2$■■帘体沟■电勢MHi«itra««u«BUM»ifUMIUfll i**ra：个壬电丄yl-rO a iln<thlMAtreABlXne(X.Y«Ex9Ey.x2.yl) ■•庄卜电初i&treanIlne(X.-Y«£x,-Ey,xl.-yl> ■・圧*电场红dth2^dthl/qi itiiH电你傀仪但*th2-<180-dtb2:-clth2:dth2rpX/ie0; ■电场n«lCteftrtx2«rO*coB<th2Hl; ■电场线钟V力■上”y2«rO a s:n<th2»; ♦电绻很的atr«aBllne(X.Y.b v Ey.x2r y2lstr«Mlina(X v-Y«Kx0-Ky r x2.-*y2) tH/iF电场幼«xl> eqS tl<3httitlec电场岐xlabcK a r\ e fMt91ze\14> QU联■住毎ylabclfl e E<U)\ e fontslzo\l()nct-l •卍Utt八仇Q\g2八让S«ul・・ n®2atr(ql IM«»tttAt»Nt(* m. /M-0.3r txt«*fonts&ae9«1€)' SI示电*比耿厂I靱厂"卜出点电的W的电场线和馬势统如图? ffi/ii：K2 - 电背的平面电场爼与粤竹怨“£・护三眄（1）甲个电備的；［M电场分布如闺3所不ffi 3 MX个电苗的立体电场分術畀汗代田如F：个电績“"电场仔令k«5•10*Sjq-10A|-^);r0-0.1;uO-k e q/rO|[X•丫“[•■phoir・ W e・rO・)U : I *iy-rO e Y( 11 ■匹•M・2( :l • jx«f X；v.»ro4(•&>•<«) Hiy-lyII«roI) J ；!•(*;x«roMl ・(*♦(*) I I; plot3(x.y«t); hold ©<iu・l"・pec・(1.3・5)・uOH)C・Y・Zl・ sph«rv;r-«e q./\>;ZIX<OAYcfll-nanjfor 1-1x5 surttrf ll A X*rlll<Y«rU)<ZI♦n<1shading int<»rpUtleC*个电紆訝代电场分命•••“"■“■••20八/乐标11 xlaMirxS e：onts：z«\X«) yla^X(e y^a：ontslzo\2«>zlab«：( *x\e:Gnt&：2«S16> 护警牛*（2）需■同号点电債时的电场理咬分It的占血设两个点电爸的电At为Q.场APd. r）的场色的舅分St为场強的y分■为g严咯mq♦聖■■&•HWHftiX 系M坨MS 磁*・<0•■[("釧7丁厂[(—盯・>丁(6b) 4%；匕足買的令确It融v的n^6t：匕是■的偶常放・足y的命的畝・匕和 &的空阿分布比牧乂余•需©通过■而相僅找乂不兴分布《1律・取匕・kQ/『为电场期电场強度釣分■町衣示为Z)尸〃九”♦/-yr八【注・<・广严）・(63)图点电荷时的电场侵度分■的曲囱axis tight%«KMi理庠代码如Fl电紳H 的电场無电分■的tlAiW 电场乞*分・的•如1cle«rrl3--(tx<D.-2<y.*2l.-<3/2r ； 左山喊点的护寓的 £»^字符席r23-•dx -X>.-2^y.-2H-(3/2r ；%«*M6边用內f)■禹的二次方字符“Ex-ir»Xln«<rix*l)./- «13 ・4<* 11./- r23|)；mam* By-lnXlM<(v y./v rl3 •*/./• r23)l; %«>»« y 5f ffl-16;■字It 大小 ■・ *kCHI眄2・5『 x«linspac<(-xn«xn 9501; ylin»p«c«(-ywi«yn tf 40li (X.Y1-Mah9rid|x 9y>i subplot 1123); surf(x 0y«Ex(X«YI) box on tltl«(• T ・HI 号炉KG 场・dtJt'E $t*AdD*• 'fontsixc 4.：aHxUbell «fa>*41 爪・卷你yUb«)r\Ky/a\a fMteU9\r«>tUb«ll •MtK.x/MQ^Xrtn - - *2*. •fonltU*' •"八41 示鼻维蒔 •Xia tiahttKIhMl subplot < 122)i tMtfnman 2 sutr (x.y«£y(X.Y|) ■•■Mbox onalatoell ^ltx/a*«^fontsixe* .fa) ■里示*■标 Qll 示 a*u四.实匕总结Ihr 电场不业.換不忆 它不ft 好通的“三物质雾謀由尿7\分子构 成.也没有可见的形态.fiKHW 可以護检測的运动速度.能■和动占有空 刚.M 斡真实的客或仔任・实lAVkAMimvhABiM*M«aai tta*絵中通过仿真软件MATIAB绘出的电场（或电势）的分布怕・讣我们对电场这艸桁喷右了屯律的峪斤认识.用MATIAB 101 HI的立体用也更冇利『对电场的nw.对丁对应如识的理解和吸ftwitt大的ffiitt.在以噸的学刃中•我仅只是佚用MATLAB的litfl计氛的功絶•通过这个实勉对于MATLAB强大的仿血功能有r出加渾対的r*i.为滋圧次的学列此软件开r -个很好的头.4il MAUAB ■出的电场线和聲勞找能U澤我们对电场的了酬. 任角闍的辻程中・个电術电■相等时•电场线和第的线对中*线业対称的.出芍个点背电■不H1尊时.电场线势找对中•役圧不对片的•但足电场找和等的线仍堆4111的.MU.咬心地鴉謝，老帅构朱帅兄在实购叩给卩的IB牙！。

基于Matlab模拟点电荷电场线和等势线张雷【期刊名称】《《科技视界》》【年(卷),期】2019(000)034【总页数】2页(P70,76)【关键词】点电荷; 电场线; 等势线; Matlab【作者】张雷【作者单位】西安文理学院陕西西安 710065【正文语种】中文【中图分类】O441.1-4; G6421 点电荷的电场点电荷是电磁学中最简单的一种理想化模型。

单个点电荷在空间某点所激发的电场强度大小E=q/4πεor2 和电势U=q/4πεor，用电场线和等势面直观形象地来表示这两个物理量。

电场线从正电荷 ( 无穷远)出发终止于负电荷（无穷远），电场线的切线方向代表该点的场强方向，电场线的疏密来表示电场的强弱。

电场线和等势面处处垂直。

然而在真空中某点所激发的电场通常是由点电荷组或者连续带电体所激发。

对于点电荷组在空间某点所激发的电场强度和电势，我们可以根据电场强度矢量叠加和电势标量叠加进行数值求解。

由于电场线和等势面的抽象性和复杂性，此时就很难想象出空间的分布。

Matlab 有强大的数学运算能力和绘图功能，能直接利用软件提供的函数，编写较简单的程序即可解决相应的物理问题。

已有文献[1 -3]对电偶极子、电四极子，多个点电荷系统的电势分布和矢量场图进行计算仿真。

根据库仑定律，一试验电荷qo 与点电荷Q 相距为r 时，受到的静电力为F，根据电强度的定义，点电荷Q在r 处产生的电场强度为其中k 是静电力常量，k = 9×109N·m2/C2。

2 利用maltab 编写程序利用matlab 编写程序，并运行程序得。

3 分析总结通过借助matlab 软件编写程序能够很好的模拟点电荷的电场线和等势面，运行程序得到的图像和我们所做的静电场模拟实验相吻合。

图1 点电荷的电场强度和电势曲线从图1 可以得到。

（1）ro 表示参考距离。

（2）与电场强度相比，在参考距离之内，电势随距离的减小而增加得较慢。

用MABLAB描述三维空间中点电荷的电势分布摘要：MATLAB语言是一种科学计算语言，所运用的范围非常广阔，利用MATLAB 语言可以解决诸多问题，如本文所要解决的点电荷的三维空间电势的分布，电势是点电荷电场的一个重要性质，了解电势的性质对理解电场很有帮助，所以，本论文重点讨论用MATLAB语言描绘电势在三维空间中的分布问题。

关键词：MATLAB 点电荷电势引言：MATLAB语言是一种科学计算语言，它功能强，效率高，而且简单易学，所以，所被运用的范围非常大，能够解决线性代数，自动控制理论，数字信号处理，时间序列分析，动态系统仿真，图像处理等诸多问题，MATLAB语言简单易懂，丰富多样，尤其能够数据可视化，以达到真正理解的目的。

用MATLAB语言来解决很多物理领域的问题，是一种很好很简捷的方法，而且丰富的二维，三维图像信息能够证明模糊的疑问，使问题清晰化，简明化。

接下来，我们将运用MATLAB语言通过编程来解决三维空间中点电荷的电势分布问题。

一程序构思在高中和大学，我们都大量地学习过关于电场这方面的知识，电场是一种特殊的物质，它看不见，摸不到，但是却存在于真实的空间中，为了研究它的性质，我们需要用一些实型来真实地表达出，于是，人们假象出了电场线这一物质，电场线不是随意地画出的，而是依据了一定的规律，从高中和大学的课本上都可知，电场线的方向是根据电势的某个趋势而改变的，电场线的方向也是电势逐渐降低的方向，而且正电荷和负电荷的电势分布又有所不同，正点电荷是随着距离的增大逐渐减小，负点电荷的电势是随着距离的增大逐渐增大，由此看来，电势最能提体现电场的特点，要想研究电场的性质，首先需要知道电势的分布情况，因此，在这篇文章中重点讨论不同情况下点电荷的电势分布。

而且运用MATLAB语言将电势的分布生动地体现在三维空间中，来证明电势的这个分布特点，二相关的物理知识及MATLAB知识由物理公式可知电场中电势的公式为U=q/4πεr，跟据此公式可以求出距点荷距离为r 处的电势；其中还须知道此点电荷的空间坐标，其中，r=2020)()(y y x x -+-。

利用MATLAB模拟点电荷电力线和等势面摘要：本文依据探索点电荷的电力线和等势线的实验理论，着重阐述利用Matlab这一计算机软件模拟和绘制点电荷的电力线和等势面的分布（包括单个正点电荷和一对等量相异点电荷），形象直观展现了点电荷在其周围场中产生的电场线以及等势面的分布状况，这为与电荷有关的理论与实验提供了帮助，并且可以通过Matlab编程描绘这样直观可视的图样总结出不同情况的点电荷的电力线和等势面的分布规律。

关键字：电场；Matlab；试探电荷；电力线；等势线；模拟图样0引言在真空中的点电荷，在周围产生了电场。

电场既看不见也摸不着，人们在探究电场的性质及分布状况的时候，用试探电荷在电场中的表现及分布状况来表征电场的状况，然后在纸张上画出电荷的电场线和等势面。

该过程不仅繁琐艰难，而且误差较大。

如果我们利用Matlab这一计算机软件进行模拟该过程，将会给我们带来极大的便利，所绘制得到的电场线和等势线面不仅生动直观，而且较接近实际情况。

正因为如此，随着当今科技的迅速发展，Matlab这一计算机软件也进行了持续的开发，并得到了广泛的应用。

1物理依据本次试验的目的，就在于熟悉单个点电荷及一对点电荷的电场分布情况，并且学会使用Matlab进行数值编程与计算，并绘出相应的图像。

根据库伦定律：在真空中，两个静止点电荷之间的作用力与这两个电荷的电量乘积成正比，与它们之间距离的平方成反比，作用力的方向在两个电荷的连线上，两电荷同号为排斥力，异号为吸引力，它们之间的作用力F满足[1]：F=k Q1Q2(1)R2，ε0称为介电常数，根据电场力的概念：其中k=14πε0F=EQ(2)跟库伦定律的表达式（1）对比，可以得到点电荷产生的电场强度E为[2](3)E=kQR2(3)式即为点电荷产生的电场的势函数。

我们以(3)式作为目标函数对电场利用Matlab 进行模拟。

在Matlab 中，由以上公式算出各点的电场强度E 后，可以用Matlab 自带的库函数绘出相应电荷的电场分布情况。

基于Matlab模拟点电荷电场线和等势线作者：张雷来源：《科技视界》2019年第34期【摘要】在静电场中引入电场强度和电势后，通过等电势线图和场强分布图可以具体的描述静电场这种抽象的物质场。

利用计算机技术来模拟静电场等物质场逐渐成为趋势，本文介绍了如何利用Matlab软件模拟点电荷的电场线和等势线。

【关键词】点电荷;电场线;等势线;Matlab中图分类号： O441.1-4;G642 文献标识码： A 文章编号： 2095-2457（2019）34-0070-001DOI：10.19694/ki.issn2095-2457.2019.34.028Simulation of Electric Field and Equipotential line of Point Charge with MatlabZHANG Lei（School of mechanical and material engineering，xi’an university of arts and s ciences，Xi’an Shaanxi 710065， China）【Abstract】The electrostatic field which is a matter field so the abstract to understand can be described in specific equipotential line and electric field line. The use of computer technology to simulate the electrostatic field gradually become a trend. This paper introduces how to use Matlab software to simulate point charge electric field lines and equipotential line.【Key words】Point charge; The electric field lines; Equipotential line; Matlab1 点电荷的电场点电荷是电磁学中最简单的一种理想化模型。

利用MATLAB模拟点电荷电力线和等势面摘要：本文依据探索点电荷的电力线和等势线的实验理论，着重阐述利用Matlab这一计算机软件模拟和绘制点电荷的电力线和等势面的分布（包括单个正点电荷和一对等量相异点电荷），形象直观展现了点电荷在其周围场中产生的电场线以及等势面的分布状况，这为与电荷有关的理论与实验提供了帮助，并且可以通过Matlab编程描绘这样直观可视的图样总结出不同情况的点电荷的电力线和等势面的分布规律。

关键字：电场；Matlab；试探电荷；电力线；等势线；模拟图样0引言在真空中的点电荷，在周围产生了电场。

电场既看不见也摸不着，人们在探究电场的性质及分布状况的时候，用试探电荷在电场中的表现及分布状况来表征电场的状况，然后在纸张上画出电荷的电场线和等势面。

该过程不仅繁琐艰难，而且误差较大。

如果我们利用Matlab这一计算机软件进行模拟该过程，将会给我们带来极大的便利，所绘制得到的电场线和等势线面不仅生动直观，而且较接近实际情况。

正因为如此，随着当今科技的迅速发展，Matlab这一计算机软件也进行了持续的开发，并得到了广泛的应用。

1物理依据本次试验的目的，就在于熟悉单个点电荷及一对点电荷的电场分布情况，并且学会使用Matlab进行数值编程与计算，并绘出相应的图像。

根据库伦定律：在真空中，两个静止点电荷之间的作用力与这两个电荷的电量乘积成正比，与它们之间距离的平方成反比，作用力的方向在两个电荷的连线上，两电荷同号为排斥力，异号为吸引力，它们之间的作用力F满足[1]：F=k Q1Q2(1)R2，ε0称为介电常数，根据电场力的概念：其中k=14πε0F=EQ(2)跟库伦定律的表达式（1）对比，可以得到点电荷产生的电场强度E为[2](3)E=kQR2(3)式即为点电荷产生的电场的势函数。

我们以(3)式作为目标函数对电场利用Matlab 进行模拟。

在Matlab 中，由以上公式算出各点的电场强度E 后，可以用Matlab 自带的库函数绘出相应电荷的电场分布情况。

电磁学1、点电荷的电场研究真空中，两个带正电的点电荷，在电量相同和电量不同情况下的电场分布。

V=V i+V2= q i+—^ ，E=-▽ V小％* 4“0「22、程序实现主程序文件名为point.mclear allep0=8.85*le-12; %真空中的电容率cO=F(4*pi*epO);e=1.6e-10;h=0.018;x=-0.5:h:0.5;y=-0.5:h:0.5;str{1}=两同号等量点电荷'str{2}=两同号不等量点电荷'[X,Y]=meshgrid(x,y);q=[e;1.9*e];for i=1:2V=c0*e./sqrt((X+0.2).A2+Y.A2)+c0.*q(i)./sqrt((X-0.2).A2+Y.A2); %求电势[Ex,Ey]=gradie nt(-V,h); % 求电场figure(i)cou nter(X(:,:,1),Y(:,:,1),V；・％等势面[20,-20,19,-19,18,-18,17,-17,16,-16,15,-15,14,-14,13,-13,12,-12,11,-11,10,-10"; Axis([-0.38,0.38,-0.28,0.28])hold onphi=0:pi/17:2*pi; %以下画电场线sx1=0.2+0.01*cos(phi);sy 1=0.01*si n( phi);streamli ne(X(:,:,1),Y(:,:,1),Ex,Ey,sx1,sy1);hold onsx2=-0.2+0.01*cos(phi);sy2=0.01*si n( phi);streamli ne(X(:,:,1),Y(:,:,1),Ex,Ey,sx2,sy2);title(str(i))text(-0.215,0,'+：'fo ntsize；20); %标示点电荷text(0.185,0,'+：'fo ntsize,20);end二、带电细棒的电场1、若电荷Q均匀分布在长为L的细棒上，求真空中，带电细棒的电场在xy平面内的分布情况。

用MATLAB解决电磁学中电场问题摘要：二十一世纪以来，随着计算机技术的进一步提高以及电脑在国内的普遍使用，在很多教学领域运用计算机软件来辅助教学已经非常普遍，然而用电脑软件处理仿真电学及电磁学中电场及电势的研究却并不常见。

在进行电场及电势问题的研究时，电场虽然是确实存在的，但是由于其抽象而不易被理解，而且即使在实验室中，我们也很难实现其理想化模型。

这样就使得教师在教学过程中不能生动且直观的描述出来，学生在学习过程中也将很难理解和接受。

本文利用计算机数学软件MATLAB模拟仿真真空中自由电荷，电偶极子，带电细棒，尖端导体等附近静电场电场线及等势线分布图型，从而使其更加形象，生动，直观，更加便于学生学习，理解和接受，同时使的教师教学更加方便，快捷。

关键字：MATLAB软件，静电场，尖端导体，尖端效应，电场强度，电势。

目录1 引言 (1)2用MATLAB处理静电场中的电场线和等势线的问题 (1)2.1用MATLAB模拟仿真真空中点电荷的电场线和等势线分布 (1)2.2用MATLAB模拟仿真真空两个点电荷的电场线及等势线分布 (3)2.2.1用MATLAB模拟真空中两个同种点电荷的电场线及等势线分布 (3)2.2.2用MATLAB模拟真空中两个同号但不等量点电荷的电场线和电势分布..5 2.2.3用MATLAB模拟电偶极子的电场线及等势线分布 (6)3用MATLAB模拟均匀带电细棒的电和分布情况 (8)4用MATLAB模拟两个无限长导线的电位和电场分布 (9)5尖端导体附近的电场及电势特点及其应用 (11)5.1电场函数 (11)6总结 (13)7参考文献 (14)1引言：二十一世纪以来，随着计算机技术的进一步提高以及电脑在国内的普遍使用，在很多教学领域运用计算机软件来辅助教学已经非常普遍，然而用电脑软件处理仿真电学及电磁学中电场及电势的研究却并不常见。

我们在研究电磁学中的电场问题时，通常将自由电荷看作一个理想化的模型。

ｌＳｓ～１００２—４９５６ＣＮｌｌ一２０３４／Ｔ实验技术与管理Ｅ邛耐ｍｅⅡ血Ｔ即ｈｎｏｌｏ酣帅ｄＭ跚ａｐ咖ｅｎｔ第２４卷第１０期２００７年１０月ｖ０１．２４ＮｏｌＯ０ｅｔ２００７基于ＭＡＴＬＡＢ构建点电荷系的电势与电场强度分布图孔祥鲲１”，原立格。

，杨宏伟’（１＿南京农业大学理学院。

江苏南京２１００９５；２．解放军镇江船艇学院基础部，江苏镇汪２１２００３）摘要：用ＭＡｒⅡ‘ＡＢ程序演示了多个点电荷组成的点电荷系激发的电势分布立体图和场强分布图，结合图像探讨了电荷对称分布情况下，点电荷系的等势面及电场分布的特点以及场强与电势之间的关系。

关键词：ＭＡｌＩＡＢ；点电荷系；等势面；电场强度；电势分布中图分类号：１Ｐ３９１．９文献标识码：Ａ文章编号：１００２—４９５６（２００７）１０．００７５．０４Ｔｈｅｄｉｓｔｒｉｂｕｔｉｏｎｄｉａｇｍｍｓｏｆｅｌｅｃｔｒｉｃｐｏｔｅｎｔｉａｌａｎｄｅｌｅｃｔｒｉｅｆｉｅｌｄｉｎｔｅｎｓｉｔｙｉｎｄｏｔ－ｃｈａ唔ｅｓｙｓｔｅｍｂａｓｅｄｏｎＭＡ’ｒＬＡＢＫＯＮＧＸｉａｎｇ．ｋｕｎｌ”，ＹＵＡＮ“一ｇｅｌ，ＹＡＮＧＨｏ“ｇ．ｗｅｉ‘（１．ｃｏｕｅｇｅｏｆｓｅｊ曲ｃｅ，Ｎ＆ｎｊｊｎｇＡ画ｃｌｌｌｔｕＩａｌＵ上ｌｉ坩鸬竹，ＮａｎＪｉｎｇ２１００９５，Ｃｌｌｉ吡；２．Ｄ印心址０ｆＦｏｍｌｄａｔｉ加ＺＩｌｅｎｊｉａＩｌｇｗａｔｅｒｃｒａｆｔｃｏｕ８９ｅ０ｆⅡ陀ＰＬＡ，ｚｈｅｎｊｉ卸ｇ２１２００３，ｃｈｉ矾）Ａｂｉ喇：Ｂｙ珊岫ＭＡＴｕＢ，“８ｐ啦ｒ画ｖｅｓｓｉｔｙｉｎｄｏｔ－ｃｈａ蜡ｅｔｌＩｅａｒｅｔｌｌｅｄｉ“ｂｕ６０ｎｄｉａ酽籼０ｆｅＪｅｃ咄ｐｏｌｅｎⅡａｌ柚ｄｅｌｅｃｍ丘ｅｌｄｉｎｔｅｎ．ｓｙ８ｔｅｍ，Ｗｉｔｌｌ∞：ａｍｐｌｅｓ０ｆⅡ”ｄ０Ｉ—ｃｈａ‘ｇｅｐｍｒｇｙ８ｔｅｍ，吐ｌｅｃｈａｒａｃｌｅ＾ＢⅡｃｓ０ｆｅ叩ｉｐｏｔｅｎｔｉａＩｓｕｄ配ｅ，６ｅｌｄｄｉ嘶ｂＩｌ６叩０ｆｔ｝ｌｅｄｅ嘶ｃｉｎｔ哪畸ａｒＩｄｔｈｅｒｅｌａｌｉｏｎＢｈｉｐ‰ｅｅｎ岫ｄｃｐ０忙Ⅱｔ柚ａｎｄｅｌｅｃ伍ｃ６ｅＩｄｉｎＩｅ嘣ｔｙｄｉ∞Ｌｌｓ埘．ｐｏｔｅｎｔｉａｌｄｉｓｔｒｉｂＩｌ—Ｋｅｙ帅ｒ出：ＭＡＩＬＡＢ；ｄ０卜ｃｈ嘲∞ｓ”ｋｍ；ｅｑｌｌｉｐ０Ｉｅ而ａｌ８毗ｋｅ；ｅｌ哪一ｃ矗ｅｌｄｉｎｔｅｎ嘶；ｅｌｅ砌ｃｔｉｏｎ电场强度、电势是描述静电场属性的重要物理量，剽用等势面和电场线可以形象地描述静电场。