电磁场大作业2.docVIP

电磁场大作业 一、题目 采用计算机编程绘制双电荷静电系统的电力线分布和等位面分布（二维xy平面分布），讨论不同的电荷值q1和q2，以及不同的间距d的特性。 基础知识及公式推导 由库伦定律：在真空中，两个静止点电荷之间的作用力与这两个电荷的电量乘积成正比，与它们之间距离的平方成反比，作用力的方向在两个电荷的连线上，两电荷同号为斥力，异号为吸力，它们之间的力F 满足： 有电场强度E的定义可知： 对于点电荷,根据场论基础中的定义,有势场E的势函数为: 并且 三、matlab编程实现 1、两个同号等量电荷对的电场线和等势线 q=1; %电量比 xm=2.5; ym=2; x=linspace(-xm,xm); y=linspace(-ym,ym); [X,Y]=meshgrid(x,y); R1=sqrt((X+1).^2+Y.^2); %第一个正电荷到场点的距离 R2=sqrt((X-1).^2+Y.^2); %第二个正电荷到场点的距离 U=1./R1+q./R2; u=1:0.5:4; figure %创建图形窗口 contour(X,Y,U,u) %画等势线 grid on %加网格 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,Ey]=gradient(-U,x(2)-x(1),y(2)-y(1)); dth1=20; th1=(dth1:dth1:180-dth1)*pi/180; r0=0.1; x1=r0*cos(th1)-1; y1=r0*sin(th1); streamline(X,Y,Ex,Ey,x1,y1) streamline(X,-Y,Ex,-Ey,x1,-y1) dth2=dth1/q; th2=(180-dth2:-dth2:dth2)*pi/180; x2=r0*cos(th2)+1; 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,16) %显示横坐标 ylabel(E(U),fontsize,16) %显示纵坐标 txt=[电荷比:\itQ\rm_2/\itQ\rm_1= num2str(q)]; text(-xm,-ym-0.3,txt,fontsize,16) 图像： 2、两个不等量同号电荷对的电场线和等势线（以电荷比为3为例） q=3; %电量比 xm=2.5; ym=2; x=linspace(-xm,xm); y=linspace(-ym,ym); [X,Y]=meshgrid(x,y); R1=sqrt((X+1).^2+Y.^2); %第一个正电荷到场点的距离 R2=sqrt((X-1).^2+Y.^2); %第二个正电荷到场点的距离 U=1./R1+q./R2; u=1:0.5:4; figure %创建图形窗口 contour(X,Y,U,u) %画等势线 grid on %加网格 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,Ey]=gradient(-U,x(2)-x(1),y(2)-y(1)); dth1=20; th1=(dth1:dth1:180-dth1)*pi/180; r0=0.1; x1=r0*cos(th1)-1; y1=r0*sin(th1); streamline(X,Y,Ex,Ey,x1,y1) streamline(X,-Y,Ex,-Ey,x1,-y1) dth2=dth1/q; th2=(180-dth2:-dth2:dth2)*pi/180; x2=r0*cos(th2)+1; 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,16) %显示横坐标 ylabel(E(U),fontsize,16) %显示纵坐标 tx