电磁场-点电荷-电场线-电势-MATLAB--仿真-中南大学

内容发布更新时间 : 2025/6/21 8:15:29星期一 下面是文章的全部内容请认真阅读。

0.050.040.030.020.01一对不相等的电荷的等势线图和电场线图0-0.01-0.02-0.03-0.04-0.05-0.05y-0.04-0.03-0.02-0.0100.010.020.030.040.05x图4

点电荷电场分布的3-D图

1050-52-10201.510.50-0.5-1-1.5-2-2图5

部分M-file;

1. 点电荷的平面电力线和等势线 %点电荷的平面电力线和等势线 %平面电力线的画法 q=1e-9; r0=0.12;

th=linspace(0,2*pi,13); [x,y]=pol2cart(th,r0); x=[x;0.1*x]; y=[y;0.1*y]; plot(x,y); grid on hold on

plot(0,0,'o','MarkerSize',12) xlabel('x','fontsize',16) ylabel('y','fontsize',16)

2. 一对电荷平面等势线和电场线图%一对电荷平面等势线和电场线图

clear all; clf;

%平面等势线的画法 q1=2e-9; q2=-1e-9;

a=0.02;%到原点的距离 xx0=0.05; yy0=0.04; k=9e9;

x=linspace(-xx0,xx0,20); y=linspace(-yy0,yy0,50); [X,Y]=meshgrid(x);

r11=sqrt((xx0/1.7-a)^2+(yy0/1.7)^2);

r22=sqrt((xx0/1.7+a)^2+(yy0/1.7)^2);

r1=sqrt((X-a).^2+Y.^2); %各点到点电荷的距离 r2=sqrt((X+a).^2+Y.^2);

U=k*q1./r1+k*q2./r2;

title('单个点电荷的电场线与等势线','fontsize',20)

%平面等势线的画法 k=9e9; r0=0.1; u0=k*q/r0;

u=linspace(1,3,7)*u0; x=linspace(-r0,r0,100); [X,Y]=meshgrid(x); r=sqrt(X.^2+Y.^2); U=k*q./r; hold on;

contour(X,Y,U,u)

%各点的电势

u0=k*q1/r11+k*q2/r22;

u=linspace(u0,-u0,11); %取21个等势向量

contour(X,Y,U,u,'k-'); hold on grid on

plot(a,0,'o','MarkerSize',12); plot(-a,0,'o','MarkerSize',12); xlabel('x','fontsize',16); ylabel('y','fontsize',16);

%平面电力线的画法 [Ex,Ey]=gradient(-U); E=sqrt(Ex.^2+Ey.^2); Ex=Ex./E; Ey=Ey./E; hold on;

quiver(X,Y,Ex,Ey);

title('一对不相等的电荷的等势线图和电场线图','fontsize',20) clear;

3. 立体电力线的画法

%立体电力线的画法 q=1e-9;

[X,Y,Z]=sphere(8); r0=0.18; r1=0.2; k=9e9; u0=k*q/r0; x=r1*X(:)'; y