电磁场实验指导书及实验报告.docxVIP

电磁场理论实验一 利用Mat lab模拟点电荷电场的分布 实验目的： 熟悉单个点电荷及一对点电荷的电场分布情况； 学会使用Mat lab进行数值计算，并绘出相应的图形； 实验原理： 根据库伦定律：在真空中，两个静止点电荷之间的作用力与这两个电荷的电 量乘积成正比，与它们之间距离的平方成反比，作用力的方向在两个电荷的连线 上，两电荷同号为斥力，异号为吸力，它们之间的力F满足： Ff警R R~ （式1） 由电场强度E的定义可知： E烽 R1 （式2） 对于点电荷，根据场论基础中的定义，有势场E的势函数为 uQ （式 3） R 而 E=-VU （式 4） 在Mat lab中，由以上公式算出各点的电势U,电场强度E后，可以用Mat lab 自带的库函数绘出相应电荷的电场分布情况。 实验容： 1.单个点电荷 *点电荷的平面电力线和等势线 真空中点电荷的场强大小是E=kq /r2 ,其中k为静电力恒量，q为电量，r 为点电荷到场点P（x,y）的距离。电场呈球对称分布，取电量q 0,电力线是以 电荷为起点的射线簇。以无穷远处为零势点，点电荷的电势为U=kq /r,当U取 常数时，此式就是等势面方程.等势面是以电荷为中心以r为半径的球面。 ?平面电力线的画法 在平面上，电力线是等角分布的射线簇，用MATLAB画射线簇很简单。取射 线的半径为(都取国际制单位)rO=O. 12,不同的角度用向量表示(单位为弧度) th=linspace(0,2*pi, 13)。射线簇的终点的直角坐标为： [x,y]=pol2cart (th,rO) ?插入x的起始坐标x=[x； 0. l*x].同样插入y的起始 坐标，y二[y； 0. l*y], x和y都是二维数组，每一列是一条射线的起始和终止 坐标。用二维画线命令plot(x.y)就画出所有电力线。 ?平面等势线的画法 在过电荷的截面上，等势线就是以电荷为中心的圆簇，用MATLAB画等势 线更加简单。静电力常量为k=9e9,电量可取为q=le- 9；最大的等势线的半径 应该比射线的半径小一点? r0=0. Io其电势为u0=k8q /r0?如果从外到里取7条 等势线，最里面的等势线的电势是最外面的3倍，那么各条线的电势用向量表 示为：u=：linspace(l,3,7)*u0。从- rO到rO取偶数个点，例如100个点,使 最中心点的坐标绕过0,各点的坐标可用向量表示：x=linspace(- r0,r0,100), 在直角坐标系中可形成网格坐标：[X,Y]=meshgrid(x)°各点到原点的距离为： r=sqrt(X. 2+Y. 2),在乘方时，乘方号前面要加点，表示对变量中的元素进行 乘方计算。各点的电势为U=k8q. /r,在进行除法运算时，除号前面也要加点， 同样表示对变量中的元素进行除法运算。用等高线命令即可画出等势线 contour(X,Y,U,u),在画等势线后一般会把电力线擦除，在画等势线之前插入 如下命令hold on就行了。平面电力线和等势线如图1,其中插入了标题等等。 越靠近点电荷的中心，电势越髙，电场强度越大，电力线和等势线也越密。 单个点电荷的电场线与等势线 图1 源程序： %点电荷的平面电力线和等势线 %平面电力线的画法 q=le-9； rO=O. 12； th=linspace(0,2*pi 113)； [x,y]=pol2cart(th,rO)； x=[x；0. l*x]: y=[y；O. l*y]: plot(x,y)； grid on hold on plot(0,0,o,MarkerSize,12) xlabel(x,fontsize,16) ylabel(y,fontsize,16) title(r单个点电荷的电场线与等势线fontsize1,20) %平面等势线的画法 k=9e9； r0=0. 1； u0=k*q/r0； u=linspace(lt3t7) *u0； x=linspace(~rO,rO,100)； [X,Y]=meshgrid(x)； r=sqrt (X. 2+Y. 2)； U=k*q. /r； hold on； contour(X,Y,U,u) clear； *点电荷的立体电力线和等势面 ?立体电力线的画法 先形成三维单位球面坐标，绕z轴一周有8条电力线[X,Y,Z]=sphere(8), 每维都是9X9的网格矩阵，将X化为行向量，就形成各条电力线的终点x坐 标x=ro=X(：)z ,其他两个坐标也可同样形成终点坐标y=ro+Y(：)r , z=ro+Z(：) o对x坐标插入原点x=[x (zeros(size(x))],其他两个坐标如下形 成 y=[y(zeros(size(y))], z=[z(zeros(size(z))],用三维画线命令