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

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
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);
contourf(X,Y,V,cv,'k-') % 用黑实线画填色等位线图
axis('square') % 在 Notebook 中,此指令不用 title('电偶极子的场 ');
hold on quiver(X,Y,Ex,Ey,0.7) % 第五输入宗量 0.7 使场强箭头长短适中。
电偶极子的等势线和电场线程序:
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);
q=1.6*10^(-19); %元电荷电量
r0=0.1; %电场线起点半径
theta=linspace(0,2*pi,13);
[x,y]=pol2cart(theta,a);
x=ver(x,y,0.5*x,0.5*y);
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 个电位值
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; %比例常数
title('单个点电荷的平面电场线与等势线','fontsize',20);%显示标题
xlabel('r','fontsize',16); %显示横坐标
ylabel('E(U)','fontsize',16); %显示纵坐标
相关文档
最新文档