电偶极子周围电场分布

k=9e9;

q1=-1e-9;

q2=-q1;

a=0.2;

xx0=0.5;

yy0=0.5;

x=linspace(-xx0,xx0,20);

y=linspace(-yy0,yy0,20);

[X,Y]=meshgrid(x,y);

r1=sqrt((X-a).^2+Y.^2);

r2=sqrt((X+a).^2+Y.^2);

U=k*q1./r1+k*q2./r2;

u0=50;

u=linspace(u0,-u0,11);

contour(X,Y,U,u,'k-');

Grid on

Hold on

plot(a,0,'o','MarkerSize',12);

plot(-a,0,'o','Markersize',12);

xlabel('x','fontsize',16);

ylabel('y','fontsize',16);

[Ex,Ey]=gradient(-U);

E=sqrt(Ex.^2+Ey.^2);

Ex=Ex./E;

Ey=Ey./E;

quiver(X,Y,Ex,Ey);

tetle('一对电荷的等势线图和电场线图','fontize',20);