实验八点电荷电场分布模拟MATLAB程序.docx

电力线图Q1=1,Q2=-4,d=2;x=-5:0.1:5;y=-5:0.1:5;[X,Y]=meshgrid(x,y);R1=sqrt(X.^2+(Y-d./2).^2);R2=sqrt(X.^2+(Y+d./2).^2);Z=Q1*(Y-d./2)./R1+Q2*(Y+d./2)./R2;C=-5:0.22:5;holdon[S,h]=contour(X,Y,Z,C);clabel(S);plot(0,d/2,marker,o,markersize,10,markerfacecolor,g,markeredgecolor,g);text(0,d/2,+,fontsize,20);plot(0,-d/2,marker,o,markersize,10,markerfacecolor,g,markeredgecolor,g);text(0,-d/2,-,fontsize,20);gridon,title(电力线图);电势分布图Q1=1,Q2=-4,d=2;x=-5:0.1:5;y=-5:0.1:5;[X,Y]=meshgrid(x,y);R1=sqrt(X.^2+(Y-d./2).^2);R2=sqrt(X.^2+(Y+d./2).^2);Z=(Q1./R1)+(Q2./R2);C=-4:0.55:1;[S,h]=contour(X,Y,Z,C);clabel(S);holdonplot(0,d/2,marker,o,markersize,10,markerfacecolor,g,markeredgecolor,g);text(0,d/2,+,fontsize,20);plot(0,-d/2,marker,o,markersize,10,markerfacecolor,g,markeredgecolor,g);text(0,-d/2,-,fontsize,20);gridon,title(电势分布图);电场矢量图Q1=1,Q2=-4,d=2;x=-5:0.1:5;y=-5:0.1:5;[X,Y]=meshgrid(x,y);R1=sqrt(X.^2+(Y-d./2).^2);R2=sqrt(X.^2+(Y+d./2).^2);Z=(Q1./R1)+(Q2./R2);C=-4:0.55:1;[Ex,Ey]=gradient(-Z);E=sqrt(Ex.^2+Ey.^2);Ex=Ex./E,Ey=Ey./E;holdon;plot(0,d/2,marker,o,markersize,10,markerfacecolor,g,markeredgecolor,g);text(0,d/2,+,fontsize,20);plot(0,-d/2,marker,o,markersize,10,markerfacecolor,g,markeredgecolor,g);text(0,-d/2,-,fontsize,20);holdon;H=quiver(X,Y,Ex,Ey,0.5)gridon,title(电场矢量图);