实验1(静电磁场)点电荷的可视化处理资料

1、精选优质文档-倾情为你奉上电磁场与微波技术实验报告（I） 课程实验：点电荷系的电场分布的可视化 班 级： 姓 名： 指导老师： 实验日期： 2015.10.15 专心-专注-专业点电荷系的电场分布的可视化1、 实验目的：1、 用MATLAB程序演示了2个点电荷组成的点电荷系（含电偶极子）激发的电势分布立体图和场强分布图；2、 结合图像探讨了电荷对称分布情况下，点电荷系的等势面及电场分布的特点以及场强与电势之间的关系；3、 熟悉单个点电荷及一对点电荷的电场分布情况；4、 学会使用Matlab进行数值计算，并绘出相应的图形，运用MATLAB对其进行可视化处理。引言：电场强度、电势是描述静电场属性的

2、重要物理量，利用等势面和电场线可以形象地描述静电场。但是有些电荷体系其电势的分布和等势面是相当复杂抽象的，人们通常只能靠大脑凭空想象。通过运用MATLAB强大的数学运算能力和绘图功能，可以对2个点电荷系统、多个点电荷系统的电势分布和矢量场图进行可视化，使得电场中电势的分布直观明了。 二、实验原理考虑真空中的一点电荷系，如果电场是由n个点电荷、所激发，电场空间中的某点的电势由电势叠加原理公式求得: (1)式中是场点到点电荷的距离.在空间直角坐标系下，设点坐标为，所在位置的坐标为，则在点处的电势为 (2)其中为空间中的任意一点。观察上式显然是一个三元函数，再加上电势 ，共4个变量，则不能在三维空间

3、中描述它们的关系。我们可以通过以下的简单方法处理：考虑平面内的电势分布，则点在平面内的电势为 . (3)由在静电场中电场强度可以表示成电势的负梯度，即或写成。则在平面上的电场强度等于电势的负梯度可以写成为如下形式： .(4)根据式，我们对2个点电荷编写相应的MATLAB程序不难得到电势分布的三维立体图和电场矢量分布图.当时，我们来研究其电场强度和电势的分布情况，通过MATLAB程序来实现这种情况下电场强度和电势分布的可视化。三MATLAB编程步骤1、设定相关物理量的符号表示；2、编好点电荷系电场强度、电势的公式；3、根据点电荷系电场强度、电势的公式运用MATLAB的相关函数plot,mesh等

4、模拟出点电荷系电场的立体图样；通过编写MATLAB程序，我们可以生成相关的MATLAB图像进行可视化，得到图形。四、实验内容1、两个等量同种电荷我们先来研究由2个点电荷组成的点电荷系针对带等量同种电荷，为了计算的方便，在作图时，我们可以取这两个点电荷的电量大小为，编写MATLAB程序（见附录MATLAB程序一）得到如下图形：程序为：clearv=1./(x-3).2+y.2).0.5+1./(x+3).2+y.2).0.5; % 读入电势计算方程xmax=10; ymax=10; ngrid=30;xplot=linspace(-xmax,xmax,ngrid); %绘图区域、网格线设定x,y

5、=meshgrid(xplot); %生成二维网格vplot=eval(v); %执行输入的电势计算方程explot,eyplot =gradient(-vplot); %计算电场强度clf;subplot(1,2,1),meshc(vplot); % 画含等势线的三维曲面xlabel(x);xlabel(x);ylabel(y);zlabel(电势);subplot(1,2,2),axis(-xmax xmax -ymax ymax)cs=contour(x,y,vplot); %画等势线clabel (cs) %在等势线上编号hold on;quiver(x,y,explot,eyplot

6、) %用箭头描述矢量场xlabel(x);ylabel(y);hold off;MATLAB图形为：图1.两个等量同种电荷电势分布和矢量场图2、两个等量异种电荷.我们再来研究由2个点电荷组成的点电荷系针对带等量异种电荷，为了计算的方便，在作图时，我们可以取这两个点电荷的电量大小为，编写MATLAB程序（见附录MATLAB程序二）得到如下图形：MATLAB程序：cleard=120; %设定循环次数for i=1:d y(i)=(i-d/2); %建立y坐标的数组 yy=y(i); %将当前的位置的y数组值赋给yy变量 for n=1:d x(n)=(n-d/2); %建立x坐标的数组 xx=x

7、(n); %将当前的位置的x数组值赋给xx变量 r=sqrt(xx2+yy2)+eps; %计算当前的径向位置 V(i,n)=10*(xx/r)/(r2-0.25*xx2/(r2); %计算电势 endendmesh(x(30:70),y(30:70),V(30:70,30:70) %画出普通的三维曲面网格图xlable(x);ylable(y);zlable(V);MATLAB图形为：图2.两个等量异种电荷电势分布和矢量场图3、 电偶极子 当这两个点电荷带等量异种电荷并且间距很小时所组成的系统就成为了电偶极子。电偶极矩表示电偶极子的大小和空间取向。设电偶极子的两电荷和间的距离为，求距离电偶极

8、子远处的电势和场强的分布。(1).电势分布设场点P到的距离分别为和，则单独存在时P点的电势为 .(5)根据电势叠加原理，有 .(6)电偶极子的中点O到场点P的距离为r, 按题意，于是有：,;则在场点的电势可以写成 .(7)其中、又 、将它们代入U的表达式，可得 .(8)电偶极子在远处的性质是由它的电偶极矩决定的。( 2 ) 场强分布(9)由于轴对称性，U与方位角j无关，采用平面极坐标系，其极轴沿电矩p，原点O位于电偶极子的中心。场强E的两个分量分别为：(10) .在电偶极子的延长线上，有：(11) (12) 在电偶极子的中垂面上有 (13) (14) 3.1、电偶极子电势在MATLAB中的模拟

9、：MATLAB程序：q=1;d=2;e0=8.*10.-12;x=-3:0.1:3;y=-3:0.1:3;x,y=meshgrid(x,y);z=q.*(1./sqrt(y-1).2+x.2)-1./sqrt(y+1).2+x.2)./(4*pi*e0);mesh(x,y,z);图像：图3 电偶极子电势在MATLAB中的模拟3.2、电偶极子电场MATLAB中的模拟：q=1;d=2;e0=8.*10.-12;x=-3:0.1:3;y=-3:0.1:3;x,y=meshgrid(x,y);z=q.*(1./sqrt(y-1).2+x.2+0.01)-1./sqrt(y+1).2+x.2+0.01)

10、./(4*pi*e0);contour(x,y,z);px,py=gradient(z);holdonstreamslice(x,y,px,py,k)图形： 图4 电偶极子电场MATLAB中的模拟：4、多点电荷系通过MATLBA程序我们还可以实现多个点电荷系的电场和电势分布的可视化，只需要对程序一中的“v=1./(x-3).2+y.2).0.5+1./(x+3).2+y.2).0.5”语句按照各个点电荷的位置代入式（3）稍加修改即可得到不同点电荷系的电场和电势的分布图。如以下是3个、4个、5个、8个、9个点电荷系的电场与电势分布图形，并且改变点电荷的正负和带电量，其电场、电势分布也要发生变化。

11、 图5.三个电荷电势分布和矢量场图图6.四个电荷电势分布和矢量场图图7.五个电荷电势分布和矢量场图图9.九个电荷电势分布和矢量场图图8.八个电荷电势分布和矢量场图5、 实验总结：由于电场看不见，摸不着，它不像普通的“三态”物质那样由原子、分子构成，也没有可见的形态，但其具有可以被检测的运动速度、能量和动量，占有空间，是一种真实的客观存在。电场强度、电势是描述静电场属性的重要物理量，利用等势面和电场线可以形象地描述静电场。但是有些电荷体系其电势的分布和等势面是相当复杂抽象的，人们通常只能靠大脑凭空想象。通过运用MATLAB强大的数学运算能力和绘图功能，可以对2个点电荷系统、多个点电荷系统的电势分布和矢量场图进行可视化，使得电场中电势的分布直观明了。实验中通过仿真软件MATLAB绘出的电场（或电势）的分布图，让我们对电场这种物质有了更深的感性认识，用MATLAB画出的立体图也更有利于对电场的理解，对于对应知识的理解和吸收有很大的帮助。在以前的学习中，我仅只是使