电偶极子的场

电偶极子的场

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); %正电荷在x=1.5;y=-1.5;负电荷在x=-1.5;y=1.5。

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');

clear;clf;q=2e-6;k=9e9;a=-2.2;b=2.2;x=-4:0.6:4;y=x;

>> [X,Y]=meshgrid(x,y); % 设置坐标网格

rp=sqrt((X-a).^2+(Y-b).^2);rm=sqrt((X+a).^2+(Y+b).^2); %正电荷在x=1.5;y=-1.5;负电荷在x=-1.5;y=1.5。

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