电偶极子的等时线和电场线

电偶极子的等势线和电场线程序：
clear;clf;
q=2e-6;
k=9e9;
a=1.5;
b=-1.5;
x=-6:0.6:6;
y=x; [X,Y]=meshgrid(x,y); % 设置坐标网点
rp=sqrt((X-a).^2+(Y-b).^2);
rm=sqrt((X+a).^2+(Y+b).^2);
V=q*k*(1./rp-1./rm); % 计算电势
[Ex,Ey]=gradient(-V); % 计算场强
AE=sqrt(Ex.^2+Ey.^2);
Ex=Ex./AE;Ey=Ey./AE;% 场强归一化，使箭头等长
cv=linspace(min(min(V)),max(max(V)),49);% 产生 49 个电位值
contourf(X,Y,V,cv,'k-') % 用黑实线画填色等位线图
axis('square') % 在 Notebook 中，此指令不用 title('电偶极子的场 ');
hold on quiver(X,Y,Ex,Ey,0.7) % 第五输入宗量 0.7 使场强箭头长短适中。
plot(a,b,'wo',a,b,'w+') % 用白线画正电荷位置
plot(-a,-b,'wo',-a,-b,'w-') % 用白线画负电荷位置
xlabel('x');ylabel('y');
hold off;

%单个点电荷的平面电场线与等势线
a=0.2;
k=9*10^9; %比例常数
q=1.6*10^(-19); %元电荷电量
r0=0.1; %电场线起点半径
theta=linspace(0,2*pi,13);
[x,y]=pol2cart(theta,a);
x=[x;0.05*x];
y=[y;0.05*y];
quiver(x,y,0.5*x,0.5*y);
plot(x,y);
hold on;
u=k*q/r0;
u1=linspace(1,3,7)*u;
x=linspace(-0.1,0.1,100);
[X,Y]=meshgrid(x);
r1=sqrt(X.^2+Y.^2);
U=k.*q./r1;
contour(X,Y,U,u1);
title('单个点电荷的平面电场线与等势线','fontsize',20);%显示标题
xlabel('r','fontsize',16); %显示横坐标
ylabel('E(U)','fontsize',16); %显示纵坐标