oujizidianchang.m

此段程序描述了电偶极子在静态的情况下的电场、电势分布图。 比较简单注重在一些基本语句的使用

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