实验1(静电磁场).doc

点电荷系的电场分布的可视化 实验目的：MATLAB程序演示了2个点电荷组成的点电荷系（含电偶极子）激发的电势分布立体图和场强分布图，结合图像探讨了电荷对称分布情况下，点电荷系的等势面及电场分布的特点以及场强与电势之间的关系，运用MATLAB对其进行可视化。 引言：电场强度、电势是描述静电场属性的重要物理量，利用等势面和电场线可以形象地描述静电场。但是有些电荷体系其电势的分布和等势面是相当复杂抽象的，人们通常只能靠大脑凭空想象。通过运用MATLAB强大的数学运算能力和绘图功能，可以对2个点电荷系统、多个点电荷系统的电势分布和矢量场图进行可视化，使得电场中电势的分布直观明了。 1.材料与方法： 1.1 基本原理 考虑真空中的一点电荷系，如果电场是由n个点电荷、、……所激发，电场空间中的某点的电势由电势叠加原理公式求得:1 式中是场点到点电荷的距离. 在空间直角坐标系下，设点坐标为，所在位置的坐标为，则在点处的电势为 2 其中为空间中的任意一点。观察上式显然是一个三元函数，再加上电势 ，共4个变量，则不能在三维空间中描述它们的关系。我们可以通过以下的简单方法处理： 考虑平面内的电势分布，则点在平面内的电势为......... 3 由在静电场中电场强度可以表示成电势的负梯度，即或写成。则在平面上的电场强度等于电势的负梯度可以写成为如下形式：................4 根据式，我们对2个点电荷编写相应的MATLAB程序不难得到电势分布的三维立体图和电场矢量分布图. 当时，我们来研究其电场强度和电势的分布情况，通过MATLAB程序来实现这种情况下电场强度和电势分布的可视化。 1.2运用MATLAB程序设计与实现 1．2.1两个等量同种电荷 我们先来研究由2个点电荷组成的点电荷系针对带等量同种电荷，为了计算的方便，在作图时，我们可以取这两个点电荷的电量大小为，编写MATLAB程序（见附录MATLAB程序一）得到如下图形： 图1.两个等量同种电荷电势分布和矢量场图 1.2.2两个等量异种电荷 .我们再来研究由2个点电荷组成的点电荷系针对带等量异种电荷，为了计算的方便，在作图时，我们可以取这两个点电荷的电量大小为，编写MATLAB程序（见附录MATLAB程序二）得到如下图形： 图2.两个等量异种电荷电势分布和矢量场图 1.2.3电偶极子当这两个点电荷带等量异种电荷并且间距很小时所组成的系统就成为了电偶极子。电偶极矩表示电偶极子的大小和空间取向。 设电偶极子的两电荷和间的距离为，求距离电偶极子远处的电势和场强的分布。 1.电势分布 设场点P到的距离分别为和，则单独存在时P点的电势为............................5 根据电势叠加原理，有..................6 电偶极子的中点O到场点P的距离为r, 按题意，于是有： , ; 则在场点的电势可以写成.........7 其中、 又 、将它们代入U的表达式，可得................8 电偶极子在远处的性质是由它的电偶极矩决定的。 2 场强分布 由于轴对称性，U与方位角无关，采用平面极坐标系，其极轴沿电矩p，原点O位于电偶极子的中心。场强E的两个分量分别为：............... 在电偶极子的延长线上，有： 在电偶极子的中垂面上有MATLAB编程思路： 第一步：设定相关物理量的符号表示； 第二步：编好点电荷系电场强度、电势的公式； 第三步：根据点电荷系电场强度、电势的公式运用MATLAB的相关函数plot,mesh等模拟出点电荷系电场的立体图样； 通过编写MATLAB程序三（见附录程序三），我们可以对电偶极子电场的电势分布进行可视化，得到如下图形： 1.2.4多点电荷系 通过MATLBA程序我们还可以实现多个点电荷系的电场和电势分布的可视化，只需要对程序一中的“v1./x-3.^2+y.^2.^0.5+1./x+3.^2+y.^2.^0.5”语句按照各个点电荷的位置代入式（3）稍加修改即可得到不同点电荷系的电场和电势的分布图。如以下是3个、4个、5个、8个、9个点电荷系的电场与电势分布图形，并且改变点电荷的正负和带电量，其电场、电势分布也要发生变化。参考文献： [1] 钟季康 鲍鸿吉 . 大学物理习题计算机解法——MATLAB编程应用[M]. 机械工业出版社, 2008.1 [2] 胡盘新，钟季康．在大学物理教材中引入计算机数值解的尝试[J]．物理与工程，2006，162：47-50． [3 ] 　张志涌. 精通MATLAB615 版[M] . 北京:北京航空航天大学出版社,2003. [4] 路宏敏，赵永久，朱满座.电磁场与电磁波基础[M].科学出版社 2006.9 [5] 郭硕鸿.电动力学[M]