资源描述：
《标量场和静电场的仿真分析报告.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、实验一标量场和静电场的仿真分析实验目的：1、掌握标量场的梯度计算方法，理解梯度的物理意义；2、掌握点电荷产生的电场特性；3、熟悉多个极性不同的点电荷形成的点电荷系的电场特性；4、熟悉电偶极子的电场特性以及电力线和等位线的关系。实验原理：1，标量场：标量场u沿指定方向的变化率就是标量场在该方向的方向导数就是标量场的梯度，计算公式为由此就可以计算出标量场的梯度。2．静电场：电位表达式这就是等位面方程．在Matlab中可求解该方程并用极坐标作图，即可得到电场的等位线图(也可画出三维立体等位面图，如图4)．电位求
2、出后就可得到球坐标系下电场E的表达式实验步骤：1、标量场的梯度和等位线仿真（1）建立梯度的数学模型（2）利用matlab软件进行仿真（3）观察并分析仿真图中梯度和等位线之间的相互关系2、点电荷系的电位分布（1）建立点电荷系电位的数学模型（2）利用matlab软件进行仿真（3）观察并分析仿真图中电位分布的特点3、电偶极子的场（1）建立电偶极子的电位和电场的数学模型（2）利用matlab软件进行仿真（3）观察并分析仿真图中电位线和电力线的特点和关系实验报告要求：（1）抓仿真程序结果图（2）理论分析与讨论[例1
3、] 求二维标量场u(r)=y2-x的梯度。理解梯度的关键是：梯度是与等值面垂直的一个矢量。Matlab程序如下：[x,y]=meshgrid(-2:.2:2,-2:.2:2);%设定坐标x，y的围z=y.^2-x;%计算公式[px,py]=gradient(z,.2,.2);contour(z)%求梯度holdonquiver(px,py)%绘制梯度图像holdofftitle('等值线与梯度');%图像标题理论分析：由实验原理中梯度的概念和计算公式，按照题目的要求就可以得到上图所示的图形。根据图形分析可
4、得到，场中每一点处的梯度垂直于过该点的等值面，且指向函数增大的方向，也就是说，梯度就是该等值面的法相矢量。[例2](1)2个等量同号点电荷组成的点电荷系的电势分布图为了方便求解，令则：clearv='1./((x-3).^2+y.^2).^0.5+1./((x+3).^2+y.^2).^0.5';%读入电势计算方程xmax=10;%x轴的坐标最大值ymax=10;%y轴的坐标最大值ngrid=30;xplot=linspace(-xmax,xmax,ngrid);%绘图区域、网格线设定[x,y]=mesh
5、grid(xplot);%生成二维网格vplot=eval(v);%执行输入的电势计算方程[explot,eyplot]=gradient(-vplot);%计算电场强度clf;subplot(1,2,1),meshc(vplot);%画含等势线的三维曲面xlabel('x');ylabel('y');zlabel('电位');subplot(1,2,2),axis([-xmaxxmax-ymaxymax])cs=contour(x,y,vplot);%画等势线clabel(cs);holdon;%在等势
6、线上编号quiver(x,y,explot,eyplot)%用箭头描述矢量场xlabel('x');ylabel('y');holdoff;分析：(1)如右图所示，电场线从正电荷出发，终止在无穷远处。电场线与等势线垂直，任何两条电场线都不相交。(2)电势较高的等势线分别包围着电荷，电势较低等势线包围着两个电荷。电场强度大的地方，电场线较密，等势线也较密。(3)当两个电荷的电量相等时，电场线和等势线对中垂线是对称的。(4)根据左图所示，点电荷场强的分量Ex在电荷附近特别大。在点电荷附近的右侧，Ex的方向沿x
7、轴正向，在点电荷附近的左侧，Ex的方向沿x轴负向，因此在点电荷的右侧形成高峰，左侧形成深谷。(5)根据左图所示，点电荷场强的分量Ey在电荷附近也特别大。在点电荷附近的前方，Ey的方向沿y轴正向，在点电荷附近的后方，Ey的方向沿y轴负向，因此在点电荷的前方形成高峰，后方形成深谷。(2)2个等量异号点电荷组成的点电荷系的电势分布图clearv='1./((x-3).^2+y.^2).^0.5-1./((x+3).^2+y.^2).^0.5';%读入电势计算方程xmax=10;%x轴的坐标最大值ymax=10;
8、%y轴的坐标最大值ngrid=30;xplot=linspace(-xmax,xmax,ngrid);%绘图区域、网格线设定[x,y]=meshgrid(xplot);%生成二维网格vplot=eval(v);%执行输入的电势计算方程[explot,eyplot]=gradient(-vplot);%计算电场强度clf;subplot(1,2,1),meshc(vplot);%画含等势线的三维曲面xlabel('x');yla
