中南大学电磁场理论实验2.docVIP

题 目 姓 名 学 号 班 级 任课老师 实验日期 电磁场理论 实验二 　　　　　——利用Matlab模拟点电荷电场的分布 实验目的： 熟悉单个点电荷及一对点电荷的电场分布情况； 学会使用Matlab进行数值计算，并绘出相应的图形； 实验原理： 根据库伦定律：在真空中，两个静止点电荷之间的作用力与这两个电荷的电量乘积成正比，与它们之间距离的平方成反比，作用力的方向在两个电荷的连线上，两电荷同号为斥力，异号为吸力，它们之间的力F满足： (式1) 由电场强度E的定义可知： (式2) 对于点电荷，根据场论基础中的定义，有势场E的势函数为 (式3) 而 (式4) 在Matlab中，由以上公式算出各点的电势U，电场强度E后，可以用Matlab自带的库函数绘出相应电荷的电场分布情况。 三、实验内容： 根据库伦定律，利用Matlab强大的绘图功能画出单个点电荷的电场分布情况，包括电力线和等势面。 根据库伦定律,利用Matlab强大的绘图功能画出一对点电荷的电场分布情况,包括电力线的分布和等势面。 实验代码1： r0=0.12; % 射线的半径 th=linspace(0,2*pi,13); % 电力线的角度 [x,y]=pol2cart(th,r0); % 将极坐标转化为直角坐标 x=[x;0.1*x]; % 插入x的起始坐标 y=[y;0.1*y]; % 插入y的起始坐标 plot(x,y,b) % 用蓝色画出所有电力线 grid on % 加网格 Hold on % 保持图像 plot(0,0,o,MarkerSize,12) % 画电荷 xlabel(x,fontsize,16) % 用16号字体标出X轴 ylabel(y,fontsize,16) % 用16号字体标出Y轴 k=9e9; % 设定k值 q=1e-9; % 设定电荷电量 r0=0.1; % 设定最大等势线的半径 u0=k*q/r0; % 算出最小的电势 u=linspace(1,3,7)*u0; % 求出各条等势线的电势大小 x=linspace(-r0,r0,100); % 将X坐标分成100等份 [X,Y]=meshgrid(x); % 在直角坐标中形成网格坐标 r=sqrt(X.^2+Y.^2); % 各个网格点到电荷点的距离 U=k*q./r; % 各点的电势 contour(X,Y,U,u) % 画出点电荷的电势面 title(正电荷的电场线和等势线,fontsize,20) %显示标题 实验代码2： k=9e9 q1=1e-9; % 设置正电荷电量 q2=-1e-9; % 设置负电荷电量 a=0.02; % 设置电荷到原点的距离 xx0=0.05; % 设置X轴的范围 yy0=0.04; % 设置Y轴的范围 x=linspace(-xx