单个点电荷的平面电场线与等势线

1、单个点电荷的平面电场线与等势线a=0.2;k=9*10人9;%比例常数q=1.6*10人(-19); %元电荷电量r0=0.1;%电场线起点半径theta=linspace(0,2*pi,13);x,y=pol2cart(theta,a);x=x;0.05*x;y=y;0.05*y;quiver(x,y,0.5*x,0.5*y)plot(x,y)hold onu=k*q/r0;u1=linspace(1,3,7)*u;x=linspace(-0.1,0.1,100);X,Y=meshgrid(x);r1=sqrt(X.A2+Y.A2);U=k.*q./r1;contour(X,Y,U,u1)t

2、itle(单个点电荷的平面电场与电场线,fontsize,20)%显示标题 xlabel(r,fontsize,16) %显示横坐标 ylabel(E(U),fontsize,16) %显示纵坐标2、画一对点电荷的平面电场线与等势线程序代码如下：%同号点电荷对的电场线和等势线clear%清除变量%横坐标范围%横坐标范围%横坐标向量%横坐标范围%横坐标范围%横坐标向量%纵坐标向量%设置坐标网点%第一个正电荷到场点的距离%第二个正电荷到场点的距离%计算电势%等势线的电势向量ym=2;x=linspace(-xm,xm); y=linspace(-ym,ym); X,Y=meshgrid(x,y);

3、R1=sqrt(X+1).A2+Y.A2);R2=sqrt(X-1).A2+Y.A2);U=1./R1+q./R2;u=1:0.5:4;figure contour(X,Y,U,u) grid on%创建图形窗口%画等势线%加网格%图例%保持图像%图例%保持图像%画水平线%画竖直线dthl=20;%左边电场线角度间隔legend(num2str(u)hold on plot(-xm;xm,0;0)plot(0;0,-ym;ym)plot(-1,0,o,MarkerSize,12)%画第一个正电荷plot(1,0,o,MarkerSize,12)%画第二个正电荷Ex,E y=gradient(-

4、U,x(2)-x(l),y(2)-y(l);%用电势梯度求场强的两个分量thl=(dthl:dthl:l80-dthl)*pi/l80;%电场线的起始角度r0=0.l;%电场线起点半径xl=r0*cos(thl)-l;%电场线的起点横坐标yl=r0*sin(thl);%电场线的起点纵坐标streamline(X,Y,Ex,Ey,xl,yl)%画左上电场线streamline(X,-Y,Ex,-Ey,xl,-yl) %画左下电场线 dth2=dthl/q;%右边电场线角度间隔th2=(l80-dth2:-dth2:dth2)*pi/l80;%电场线的起始角度x2=r0*cos(th2)+l;%电场线的起点横坐标y2=r0*sin(th2);%电场线的起点纵坐标streamline(X,Y,Ex,Ey,x2,y2)%画右上电场线streamline(X,-Y,Ex,-Ey,x2,-y2)%画右下电场线axis equal tight%使坐标刻度相等title(等量同号点电荷的电场线和等势线,fontsize,20)%显示标题 xlabel(r,fontsize,l6)%显示横坐标ylabel(E(U),fontsize,l6) %显示纵坐标txt=电荷比: