电磁场理论试验一.DOC

CENTRAL SOUTH UNIVERSITY 题 目 利用Matlab模拟点电荷电场的分布 姓 名 刘翔 学 号 0909081215 班 级 电气0802 任课老师 李义勇 实验日期 2010-10-25 电磁场理论 实验一 　　　　　——利用Matlab模拟点电荷电场的分布 实验目的： 熟悉单个点电荷及一对点电荷的电场分布情况； 学会使用Matlab进行数值计算，并绘出相应的图形； 实验原理： 根据库伦定律：在真空中，两个静止点电荷之间的作用力与这两个电荷的电量乘积成正比，与它们之间距离的平方成反比，作用力的方向在两个电荷的连线上，两电荷同号为斥力，异号为吸力，它们之间的力F满足： (式1) 由电场强度E的定义可知： (式2) 对于点电荷，根据场论基础中的定义，有势场E的势函数为 (式3) 而 (式4) 在Matlab中，由以上公式算出各点的电势U，电场强度E后，可以用Matlab自带的库函数绘出相应电荷的电场分布情况。 实验内容： 画单个点电荷的平面电场线与等势线，正点电荷与负点电荷任选一个作图； 易证： （式5） 即电势在空间某点的梯度方向与该点场强方向相反。 根据（式5），可画出正点电荷的等势线如下 其中箭头方向为电场方向。 源程序如下： claer [x,y]=meshgrid(-4:0.1:4,-4:0.1:4); z=1./sqrt(x.^2+y.^2+0.01); [px,py]=gradient(-z); contour(x,y,z,1:0.5:4); hold on quiver(x,y,px,py,b) axis equal tight 那么，在平面上确定的一个点，可得到该点的梯度方向，即电场方向的反方向。那么根据平面上点的坐标和梯度方向则可以画出点电荷的电场线。如下： 源程序如下： clear xm=2; ym=2; x=linspace(-xm,xm); y=linspace(-ym,ym); [X,Y]=meshgrid(x,y); R1=sqrt(X.^2+Y.^2); U=1./R1; %设kQ=1 figure contour(X,Y,U,1:0.5:4) grid on hold on [Ex,Ey]=gradient(-U,x,y); %求电势的梯度 dth1=20; th1=(dth1:dth1:180-dth1)*pi/180; r0=0.1; x1=r0*cos(th1); y1=r0*sin(th1); streamline(X,Y,Ex,Ey,x1,y1) streamline(X,-Y,Ex,-Ey,x1,-y1) axis equal tight 画一对点电荷的平面电场线与等势线，可以在一正一负，两个负电荷和两个正电荷之中任选一组； 选取空间两个正点电荷。 根据（式5），则对于空间一点的电势，有 （式六） 其中U为点电荷u在点产生的电势，V为点电荷v在点产生的电势。显然空间中点电势的梯度方向为该点场强的反方向。 当两正点电荷点所带电量相等时，等势线与电场线如下图所示， 当两正点电荷所带电量不等时，等势线与电场线如下图所示， 由图中所示，可见带点量大的电荷（图中右侧电荷）的电场线较带电量小的电荷密，电场强度为零的点偏离中心靠近带点量小的电荷。 源程序如下： clear q=2; xm=2; ym=2; x=linspace(-xm,xm); y=linspace(-ym,ym); [X,Y]=meshgrid(x,y); R1=sqrt((X+1).^2+Y.^2+0.01); R2=sqrt((X-1).^2+Y.^2+0.01); U=1./R