利用streamline函数绘制电力线

利用streamline函数绘制电力线

​ 利用streamline函数绘制电场强度的分布情况，首先需要得到分布在二维和三维空间中的电场强度。比如，我们在这里以一对点电荷为例来说明情况。

​ 对于点电荷系统而言，在场中任意一点 P(x,y) 的电势可以表示为(在这里我们仅仅是为了展示如何绘制电力线，所以我们对系数$q_1/4\pi ϵ$ 不进行考虑)
$$\varphi =\frac{1}{\sqrt{(x+a{)}^2+y^2}}+\frac{Q^{\ast }}{\sqrt{(x-a{)}^2+y^2}}$$
其中，两个点电荷分别位于$(-a,0),(a,0)$ ,$Q^{\ast }=q_2/q_1$ ，可以计算得到在$xoy$平面内只有两个分量，且两个分量的表达式为
$$\begin{gathered} E_x=\frac{x+a}{[(x+a{)}^2+y^2{]}^{3/2}}+Q^{\ast }\frac{x-a}{[(x-a{)}^2+y^2{]}^{3/2}} \\ {E}_y=\frac{y}{[(x+a{)}^2+y^2{]}^{3/2}}+Q^{\ast }\frac{y}{[(x-a{)}^2+y^2{]}^{3/2}} \end{gathered}$$
当然，直接采用点电荷所产生的电场强度的公式，并利用叠加原理，同样可以得到上面的公式。

我们分别考虑两个点电荷系统异量同号和异量异号时候的电力线分布情况。不妨假设$Q^{\ast }=\pm 2$并且假定两个点电荷的位置位于$(\pm 1,0)$。对于初始电力线位置的考虑，因为电力线式起始于正电荷而终止于负电荷，当起始点距离电荷很近的时候，我们可以认为各个起始点绕电荷进行均匀分布。需要指出的是，若电力线分布在正电荷的周围，从正电荷出发，绘制流线的时候取场分布为$[E_x,E_y]$;若电力线分布在负电荷的周围，终止于负电荷，则绘制流线的时候取$[-E_x,-E_y]$,否则在使用streamline函数的时候无法正常的绘制流线。

首先绘制两个不等量同号电荷对应的电力线分布

程序的代码如下：

clear all
close all
x=-4:0.02:4;y=x;         %生成坐标
[X,Y]=meshgrid(x,y);     %形成网格坐标
R1=sqrt((X+1).^2+Y.^2);
R2=sqrt((X-1).^2+Y.^2);
phi=1./R1+2./R2;
[Ex,Ey]=gradient(-phi);  %利用梯度计算电场强度
hold on
r0=0.1;
th=20:20:360-20;
th=th*pi/180;            %角度转化为弧度
x1=r0*cos(th)-1;
y1=r0*sin(th); 
h=streamline(X,Y,Ex,Ey,x1,y1); %开始绘制流线
xr=-x1;yr=y1;
h=streamline(X,Y,Ex,Ey,xr,yr);
axis image;

程序完成运行之后可以得如下的电力线图像：

绘制两个不等量异号电荷对应的电力线分布

程序代码如下：

clear all
close all
x=-4:0.02:4;y=x;         %生成坐标
[X,Y]=meshgrid(x,y);     %形成网格坐标
R1=sqrt((X+1).^2+Y.^2);
R2=sqrt((X-1).^2+Y.^2);
phi=1./R1-2./R2;
[Ex,Ey]=gradient(-phi);  %利用梯度计算电场强度
hold on
r0=0.1;
th=20:20:360-20;
th=th*pi/180;            %角度转化为弧度
x1=r0*cos(th)-1;
y1=r0*sin(th); 
h=streamline(X,Y,Ex,Ey,x1,y1); %开始绘制流线
xr=-x1;yr=y1;
h=streamline(X,Y,-Ex,-Ey,xr,yr);
axis image;

程序完成运行之后得到的图像如下所示：

正如前面所示的那样，当我们使用streamline函数的时候，对于放置在右侧的负电荷，绘制流线图的时候，需要将计算所得到的电场强度分别取反，才能够绘制电力线出来。这是因为当我们采用异号的电荷的时候，电力线是需要进入到负电荷从而终止。