exm040522_2.m

来自「这本书是matlab经典书籍」· M 代码 · 共 15 行

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('\fontname{隶书}\fontsize{22}偶极子的场'),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