投稿 电偶极子和电四极子电场的计算机模拟.docVIP

电偶极子和电四极子电场的计算机模拟 王明美 （合肥师范学院电子信息工程学院，安徽，合肥，230601） 摘要：使用MATLAB对由电偶极子和电四极子的电场强度和电势进行数值计算，给出其分布，并绘制出相应电场线和等势线的图形。 关键词：电学系统，电场线，等势面，MATLAB绘图 MATLAB在大学物理中具有广泛的应用[1]-[7]，本文使用MATLAB对电偶极子和电四极子电场进行模拟,图像比较清晰，程序相对简单且适于推广。 1.电偶极子的电场和电势 对于两个点电荷系统，设两个点电荷相距为，所带电量分别为和。 由电场强度的叠加原理可知： （1） 有势场的势函数为 (2) 而 (3) 为了便于数值计算，电势可化为 ， (4) 其中，x* = x/a，y* = y/a，U0 = kQ1/a，q* = q2/q1。U0是q1在原点产生的电势，是电荷比。 场强可表示为 ， 即 ， (5) 其中，E0 = U0/a，U* = U/U0。同理可得 ， (6) 电场强度可通过约化电势计算。电场线可根据电场强度绘制。电场线从正电荷发出，当电场线的起点离电荷很近时，可认为起点绕电荷是均匀分布的。 在MATLAB中，仍然可以由以上公式算出各点的电势U，电场强度E后，再用MATLAB的库函数绘制出相应的电场线和等势线的分布情况。 电偶极子的电场线和等势线的画法（等量异号点电荷对q2:q1=1）程序如下： clear %清除变量 q=1; %电量比 xm=2.5; %横坐标范围 ym=2; %横坐标范围 x=linspace(-xm,xm); %横坐标向量 y=linspace(-ym,ym); %纵坐标向量 [X,Y]=meshgrid(x,y); %设置坐标网点 R1=sqrt((X+1).^2+Y.^2); %第一个正电荷到场点的距离 R2=sqrt((X-1).^2+Y.^2); %第二个正电荷到场点的距离 U=1./R1-q./R2; %计算电势 u=-4:0.5:4; %等势线的电势向量 figure %创建图形窗口 contour(X,Y,U,u,--) %画等势线 hold on %保持图像 plot(-1,0,o,MarkerSize,12) %画正电荷 plot(1,0,o,MarkerSize,12) %画负电荷 [Ex,Ey]=gradient(-U,x(2)-x(1),y(2)-y(1));%用电势梯度求场强的两个分量 r0=0.1; %电场线起点半径 a=20; %电场线角度间隔 b=(a:a:360-a)*pi/180; %正电荷电场线的起始和终止角度及其步长 x2=r0*cos(b)-1; %电场线的起点横坐标 y2=r0*sin(b); %电场线的起点纵坐标 streamline(X,Y,Ex,Ey,x2,y2) %画正电荷电场线 c=(-90-a:a:90+a)*pi/180; %负电荷电场线的起始和终止角度及其步长 x3=r0*cos(c)+1; %电场线的起点横坐标 y3=r0*sin(c); %电场线的起点纵坐标 streamline(X,Y,-Ex,-Ey,x3,y3) %画负电荷电场线 axis equal tight