电磁场实验指导书及实验报告
CENTRAL SOUTH UNIVERSITY
题目利用Matlab模拟点电荷电场的分布姓名xxxx
学号xxxxxxxxxx
班级电气xxxx班
任课老师xxxx
实验日期2010-10
电磁场理论 实验一
——利用Matlab 模拟点电荷电场的分布
一．实验目的：
1．熟悉单个点电荷及一对点电荷的电场分布情况； 2．学会使用Matlab 进行数值计算，并绘出相应的图形；
二．实验原理：
根据库伦定律：在真空中，两个静止点电荷之间的作用力与这两个电荷的电量乘积成正比，与它们之间距离的平方成反比，作用力的方向在两个电荷的连线上，两电荷同号为斥力，异号为吸力，它们之间的力F 满足：
R R Q Q k F ?
212=
(式1)
由电场强度E 的定义可知：
R R kQ E ?
2
= (式2)
对于点电荷，根据场论基础中的定义，有势场E 的势函数为
R kQ
U = (式3)
而 U E -?= (式4)
在Matlab 中，由以上公式算出各点的电势U ，电场强度E 后，可以用Matlab 自带的库函数绘出相应电荷的电场分布情况。
三．实验容：
1. 单个点电荷
点电荷的平面电力线和等势线
真空中点电荷的场强大小是E=kq /r^2 ,其中k 为静电力恒量, q 为电量, r 为点电荷到场点P(x,y)的距离。电场呈球对称分布, 取电量q> 0, 电力线是以
电荷为起点的射线簇。以无穷远处为零势点, 点电荷的电势为U=kq /r,当U 取常数时, 此式就是等势面方程.等势面是以电荷为中心以r 为半径的球面。
◆平面电力线的画法
在平面上, 电力线是等角分布的射线簇, 用MATLAB 画射线簇很简单。取射线的半径为( 都取国际制单位) r0=0.12, 不同的角度用向量表示( 单位为弧度) th=linspace(0,2*pi,13)。射线簇的终点的直角坐标为:
[x,y]=pol2cart(th,r0)。插入x 的起始坐标x=[x; 0.1*x].同样插入y 的起始坐标, y=[y; 0.1*y], x 和y 都是二维数组, 每一列是一条射线的起始和终止
坐标。用二维画线命令plot(x,y)就画出所有电力线。
◆平面等势线的画法
在过电荷的截面上, 等势线就是以电荷为中心的圆簇, 用MATLAB 画等势
线更加简单。静电力常量为k=9e9, 电量可取为q=1e- 9; 最大的等势线的半径应该比射线的半径小一点? r0=0.1。其电势为u0=k8q /r0。如果从外到里取7 条等势线, 最里面的等势线的电势是最外面的3 倍, 那么各条线的电势用向量表
示为: u=linspace(1,3,7)*u0。从- r0 到r0 取偶数个点, 例如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, 其中插入了标题等等。越靠近点电荷的中心, 电势越高, 电场强度越大, 电力线和等势线也越密。
x
y
单个点电荷的电场线与等势线
图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)
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)
clear;
点电荷的立体电力线和等势面
◆立体电力线的画法
先形成三维单位球面坐标, 绕z 轴一周有8 条电力线[X,Y,Z]=sphere(8),每维都是9×9 的网格矩阵, 将X 化为行向量, 就形成各条电力线的终点x 坐
标x=r0=X(:)′, 其他两个坐标也可同样形成终点坐标y=r0+Y(:)' ,
z=r0+Z(:)' 。对x坐标插入原点x=[x(zeros(size(x))], 其他两个坐标如下形成y=[y(zeros(size(y))], z=[z(zeros(size(z))], 用三维画线命令
plot3(x,y,z), 就画出所有电力线。
◆立体等势面的画法
画5 条等势面时, 各面的电势为u=linspace(1,3,5)+u0, 各等势面的半径为r=k6q. /u, 其中第一个球面的半径为rr=r(1)。三维单位球面的坐标可由[X,Y,Z]=sphere 命令形成, 每维都是21×21 的网格矩阵, 由于外球会包围球, 因此把球面的四分之一设为非数, 表示割去该部分Z(X<0&Y<0)=nan. 用曲面命
令可画出第一个曲面surf(rr6X,rr6Y,rr6Z), 只要取不同的半径就能画出不同
的等势面.为了使等势面好看, 可设置一个颜色浓淡连续变化的命令shading interp。点电荷的立体电力线和等势面如图2, 旋转图片可从不同的角度观察。
x
正电荷电场线等势面的三维图形
y
z
图2
源程序：
%立体电力线的画法 q=1e-9;
[X,Y,Z]=sphere(8); r0=0.18; r1=0.2; k=9e9; u0=k*q/r0; x=r1*X(:)'; y=r1*Y(:)'; z=r1*Z(:)';
x=[x;zeros(size(x))]; y=[y;zeros(size(y))]; z=[z;zeros(size(z))]; plot3(x,y,z) hold on;
%立体等势线之画法
u=linspace(1,3,5)*u0; %画5 条等势面时, 各面的电势为
电磁场HFSS实验报告
实验一 T形波导的内场分析 实验目的 1、熟悉并掌握HFSS的工作界面、操作步骤及工作流程。 2、掌握T型波导功分器的设计方法、优化设计方法和工作原理。实验仪器 1、装有windows 系统的PC 一台 2、HFSS15.0 或更高版本软件 3、截图软件 实验原理 本实验所要分析的器件是下图所示的一个带有隔片的T形波导。其中，波导的端口1是信号输入端口，端口2和端口3是信号输出端口。正对着端口1一侧的波导壁凹进去一块，相当于在此处放置一个金属隔片。通过调节隔片的位置可以调节在端口1传输到端口2，从端口1传输到端口3的信号能量大小，以及反射回端口1的信号能量大小。 T形波导
实验步骤 1、新建工程设置： 运行HFSS并新建工程：打开HFSS 软件后，自动创建一个新工程：Project1，由主菜单选File\Save as ，保存在指定的文件夹内，命名为Ex1_Tee；由主菜单选Project\ Insert HFSS Design，在工程树中选择HFSSModel1，点右键，选择Rename项，将设计命名为TeeModel。 选择求解类型为模式驱动（Driven Model）：由主菜单选HFSS\Solution Type ，在弹出对话窗选择Driven Model 项。 设置长度单位为in：由主菜单选3D Modeler\Units ，在Set Model Units 对话框中选中in 项。。 2、创建T形波导模型： 创建长方形模型：在Draw 菜单中，点击Box 选项，在Command 页输入尺寸参数以及重命名；在Attribute页我们可以为长方体设置名称、材料、颜色、透明度等参数Transparent（透明度）将其设为0.8。Material（材料）保持为Vacuum。 设置波端口源励：选中长方体平行于yz 面、x=2 的平面；单击右键，选择Assign Excitation\Wave port项，弹出Wave Port界面，输入名称WavePort1；点击积分线(Integration Line) 下的New line ，则提示绘制端口，在绘图区该面的下边缘中部即（2,0,0)处点左键，确定端口起始点，再选上边缘中部即(2,0,0.4)处，作为端口终点。 复制长方体：展开绘图历史树的Model\Vacuum\Tee节点，右键
武大电气工程电磁场仿真实验报告
武汉大学 工程电磁场及高电压综合实验
一、题目 有一极长的方形金属槽，边宽为1cm，除顶盖电位为100sinπxV外，其他三面的电位均为零，试用差分法求槽内电位的分布。 二、解题原理：均匀媒质中的有限差分法 我们在求解场的分布时，当边界形状比较复杂时，解析分析法不再适合了，我们可以采用数值计算的方法，数值计算法的基本思想，是将整体连续的场域划分为若干个细小区域，一般称之为网格或单元，如图1所示，然后用所求的网格交点（一般称为节点或离散点）的数值解，来代替整个场域的真实解。因而数值解，即是所求场域离散点的解。虽然数值解是一种近似解法，但当划分的网格或单元愈密时，离散点数目也愈多，近似解（数值解）也就愈逼近于真实值。 实解。在此处键入公式。 图1场域的剖分，网格节点及步长
(一）、场域的剖分、网格节点及步长 由边界Γ所界定的二维平行平面场（见图1），若采用直角坐标系则可令该场处在xoy 平面内。 所谓场域的剖分就是场域的离散化，即将场域剖分为若干个网格或单元。最常见最简单的剖分为正方形剖分，这种剖分就是在xy 平面上作许多分别与x 轴及y 轴平行的直线，称为网格线。网格线的交点称为节点或离散点，场域内的节点称为内节点，场域边界上的节点称为边界节点。两相邻网格线间距离称为步长，一般以h 表示。若步长相等则整个场域就被剖分为许多正方形网格，这就是正方形剖分。节点（离散点）的布局不一定采用正方形剖分，矩形剖分也常采用，正三角形剖分偶尔也被应用，不过最常见的最简单的仍然是正方形剖分。 （二）、差分与微分 从前面的分析可知，稳恒电、磁场的求解问题，归根到底是求解满足给定边界条件的偏微分方程（泊松方程或拉普拉斯方程）的解的问题所谓差分方法，就是用差商近似代替偏微商，或者说用差分代替微分，从而把偏微分方程转换为差分方程，后者实际上为代数方程。因此这种转化有利于方程的求解。 下面分别对一阶及二阶的差分公式进行推导。首先回顾有关偏导数的定义，有 00(,)(,)(,)(,) lim lim x x f f x x y f x y f x y f x x y x x x →→?+---==? (1) 因此当|x| 充分小时，可近似地用(,)(,)f x x y f x y x +- 或(,)(,) f x y f x x y x -- 代 替 f x ??，所谓差分公式，即是基于上述观点推得的。 设图1所示场域中的位函数为A ，任取一网格节点0，它在xy 平面上的坐标为（x ,i i y ），记节点0的矢量磁位为,i j A ，并把与节点0相邻的其他四个节点1、2、3、4的矢量磁位分别记为1,i j A +、,1i j A +、1,i j A -、,1i j A -，将节点0处函数A 的 一阶偏微商A x ??，用1、0两点函数值的差商1,,i j i j A A h +-近似代替，则有
电磁场与电磁波(第三版)课后答案第9章
第九章习题解答 9.1 设元天线的轴线沿东西方向放置，在远方有一移动接收台停在正南方而收到最大电场强度，当电台沿以元天线为中心的圆周在地面移动时，电场强度渐渐减小，问当电场强 时，电台的位置偏离正南多少度？ 解：元天线（电基本振子）的辐射场为 j k r j θ-=E e 可见其方向性函数为(),sin f θφθ=，当接收台停在正南方向（即090θ=）时，得到最大电场强度。由 s i n θ= 得 045θ= 此时接收台偏离正南方向045±。 9.2 上题中如果接收台不动，将元天线在水平面内绕中心旋转，结果如何？如果接收天线也是元天线，讨论收发两天线的相对方位对测量结果的影响。 解： 如果接收台处于正南方向不动，将天线在水平面内绕中心旋转，当天线的轴线转至沿东西方向时，接收台收到最大电场强度，随着天线地旋转，接收台收到电场强度将逐渐变小，天线的轴线转至沿东南北方向时，接收台收到电场强度为零。如果继续旋转元天线，收台收到电场强度将逐渐由零慢慢增加，直至达到最大，随着元天线地不断旋转，接收台收到电场强度将周而复始地变化。 当接收台也是元天线，只有当两天线轴线平行时接收台收到最大电场强度；当两天线轴线垂直时接收台收到的电场强度为零；当两天线轴线任意位置，接收台收到的电场强介于最大值和零值之间。 9.3 如题9.3图所示一半波天线，其上电流分布为() 11cos 2 2m I I kz z ??=-<< ??? （1）求证：当0r l >>时， 020 cos cos 22sin jkr m z I e A kr πθμπθ -?? ? ??= ? （2）求远区的磁场和电场； （3）求坡印廷矢量； （4）已知22 c o s c o s 20.609sin d π πθθθ ?? ? ?? =? ，求辐射电阻； （5）求方向性系数。 题9.3（1） 图 解：（1）沿z 方向的电流z I 在空间任意一点()0,P r θ产生的矢量磁位为
大学物理演示实验报告：基于电磁学验证流体力学伯努利方程实验
物理演示实验报告物理演示实验自主设计方案
本物理演示实验根据流体流速与压强的关系以及电磁铁的相关性质验证流体力学中伯努利原理 )(2 112111为常数C C gh v p =++ρρ(1)当外界环境被选定后，常数C 可以表示为 gh v p C 2222221ρρ++=(2)将(1)式与(2)式联立，可以得到 gh v p gh v p 22222121112 121ρρρρ++=++(3)这就是我们所说的伯努利方程，下面我们来验证这一原理。 在中学阶段，我们已经知道流体流速越大的地方压强越小这一流体学基本关系。为了验证流速与压强的具体关系，我们不妨选择空气流作为实验流体，大气压强作为外界标准压强，由基本数据可知标准大气的密度ρ=1.29kg/m 3 (温度为0℃，标准大气压p 0=101kpa)，我们只需要测量出流体的某一流速v 以及在该流 速下的压强p 1。进而将p 1,v 代入伯努利方程左右两端，验证等式是否成立。 此时，由于选定的外界是标准大气，故验证的等式为 02121p v p =+ρ(4)下面我们需要清楚流速与该流速下的流体压强的测量原理。 首先我们先测量流速。由于流体是以风的形式存在的，因此我们使用鼓风机作为风的发生装置。我们采取简易风车来测量风速。选择该风车的前提是在无风环境下风车能够静止即处于平衡状态，并且在受到风力时可以较为灵敏地进行转动，即摩擦阻力越小越好。设风车的转动半径为R，风车转动角速度为ω，则根据线速度与角速度的关系有 ωR v =(5) 其中ω可以通过风车的转速n 来测量，即 n πω2=(6) 联立(5)(6)两式，这样我们可以较为准确地得出流速v 的大小为 Rn v π2=(7) 接下来，我们来测量该流速下的压强。该压强的测量需要运用电磁铁以及压一、演示物理原理简介（可以配图说明）
工程电磁场实验报告
实验一 实验目的和要求：学习矢量的定义方法(例A=[1,2,3])，加减运算，以及点积dot(A,B)、叉积cros s (A,B)、求模运算n orm(A)。 实验内容： 1、通过调用函数，完成下面计算【p31，习题1.1】。 给定三个矢量A 、B 和C 如下： 23452x y z y z x z A e e e B e e C e e =+-=-+=- 求（1）A e ；（2）||A B -； （3）A B ?； （4）AB θ （5）A 在B 上的投影 （6）A C ?； （7）()A B C ??和()C A B ??； （8）()A B C ??和()A B C ?? 程序如下： A=[1,2,-3]; B=[0,-4,1]; C=[5,0,-2]; ea=A/norm(A) T2=norm(A-B) T3=dot(A,B) theta=acos(dot(A,B)/(norm(A)*norm(B))) theta*180/pi T5=norm(A)*cos(theta) T6=cross(A,C)
T71=dot(A, cross(B,C)) T72=dot(cross(A,B), C) T81=cross(cross(A,B),C) T82=cross(A,cross(B,C)) 运行如图： 结果如下：
2、三角形的三个顶点位于A(6,-1,2), B(-2,3,-4), C(-3, 1,5)点，求（1）该三角形的面积；（2）与该三角形所在平面垂直的单位矢量。 程序如下： A=[6,-1,2]; B=[-2,3,-4]; C=[-3, 1,5]; n=cross(B-A, C-A); S=1/2*norm(n)
电磁场与电磁波点电荷模拟实验报告
重庆大学 电磁场与电磁波课程实践报告 题目：点电荷电场模拟实验 日期：2013 年12 月7 日 N=28
《电磁场与电磁波》课程实践 点电荷电场模拟实验 1.实验背景 电磁场与电磁波课程内容理论性强，概念抽象，较难理解。在电磁场教学中，各种点电荷的电场线成平面分布，等势面通常用等势线来表示。MATLAB 是一种广泛应用于工程、科研等计算和数值分析领域的高级计算机语言，以矩阵作为数据操作的基本单位，提供十分丰富的数值计算函数、符号计算功能和强大的绘图能力。为了更好地理解电场强度的概念，更直观更形象地理解电力线和等势线的物理意义，本实验将应用MATLAB 对点电荷的电场线和等势线进行模拟实验。 2.实验目的 应用MATLAB 模拟点电荷的电场线和等势线 3.实验原理 根据电磁场理论，若电荷在空间激发的电势分布为V ，则电场强度等于电势梯度的负值，即： E V =-? 真空中若以无穷远为电势零点，则在两个点电荷的电场中，空间的电势分布为： 1 212010244q q V V V R R πεπε=+=+ 本实验中，为便于数值计算，电势可取为
1212 q q V R R =+ 4.实验内容 应用MATLAB 计算并绘出以下电场线和等势线，其中q 1位于(-1,0,0)，q 2位于(1,0,0)，n 为个人在班级里的序号： (1) 电偶极子的电场线和等势线（等量异号点电荷对q 2:q 1 = 1，q 2为负电荷）； (2) 两个不等量异号电荷的电场线和等势线（q 2:q 1 = 1 + n /2，q 2为负电荷）； (3) 两个等量同号电荷的电场线和等势线； (4) 两个不等量同号电荷的电场线和等势线（q 2:q 1 = 1 + n /2）； (5) 三个电荷，q 1、q 2为(1)中的电偶极子，q 3为位于(0,0,0)的单位正电荷。、 n=28 （1） 电偶极子的电场线和等势线（等量异号点电荷对q 2:q 1 = 1，q 2为负电荷）； 程序1： clear all q=1; xm=2.5; ym=2; x=linspace(-xm,xm); y=linspace(-ym,ym); [X,Y]=meshgrid(x,y); R1=sqrt((X+1).^2+Y.^2); R2=sqrt((X-1).^2+Y.^2); U=1./R1-q./R2; u=-4:0.5:4; figure contour(X,Y,U,u,'--'); hold on plot(-1,0,'o','MarkerSize',12); plot(1,0,'o','MarkerSize',12); [Ex,Ey]=gradient(-U,x(2)-x(1),y(2)-y(1));
电磁场实验报告
实验一：静电场的分析与求解 1.求二维标量场u（r）=y^2-x的梯度 [x,y]=meshgrid(-2:.2:2,-2:.2:2); z=y.^2-x; [px,py]=gradient(z,.2,.2); contour(z) hold on quiver(px,py) hold off title('等值线与梯度'); 2.2个等量同号点电荷组成的点电荷系的电势分布图clear v='1./((x-3).^2+y.^2).^0.5+1./((x+3).^2+y.^2).^0.5'; xmax=10; ymax=10; 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'); ylabel('y'); zlabel('电位');
subplot(1,2,2),axis([-xmax xmax -ymax ymax]); cs=contour(x,y,vplot); clabel(cs); hold on quiver(x,y,explot,eyplot) xlabel('x'); ylabel('y'); hold off 3.电偶极子的场（等位线和梯度） clear; clf; q=2e-6; k=9e9; a=1.5; b=-1.5; x=-6:0.6:6; y=x; [X,Y]=meshgrid(x,y); rp=sqrt((X-a).^2+(Y-b).^2); rm=sqrt((X+a).^2+(Y+b).^2); V=q*k*(1./rp-1./rm); [Ex,Ey]=gradient(-V); AE=sqrt(Ex.^2+Ey.^2); Ex=Ex./AE; Ey=Ey./AE; cv=linspace(min(min(V)),max(max(V)),49);
国中学生物理竞赛实验指导书思考题参考答案-电磁学
实验十一 学习使用数字万用表 【思考题参考答案】 1．调节电阻箱的电阻为500Ω和5Ω时，电阻的误差是多大？ 答：以0.1级为例，以每个接点的接触电阻按0.002Ω为例。 对于500Ω电阻从0和9999Ω接线柱输出，误差为 Ω=?+?=?512.0002.06%1.0500R 对于5Ω电阻从0和9.9Ω接线柱输出，误差为 Ω=?+?=?009.0002.02%1.05R 2．电源电压为110V 。是否可以只用一个电阻箱控制，得到0.5A 的电流？ 答：若只用一个电阻箱控制，所需电阻为Ω2205.0110==R 。这需要电阻箱的100?R 档，此档允许电流为0.05A ，实际电流大于额定电流，不能使用。 3．对于一块四位半的数字万用电表的直流电压200mV 量程，可能出现的最大数字是多少？最小分辨率是多少？ 答：最大数字为199.99mV 。最小分辨率为0.01mV 。 4．使用数字万用电表的直流电压2V 量程测量直流电压，测量值为1.5V ，测量误差为多少？如果测量值为0.15V ，测量误差为多少？如果换用200mV 量程测量直流电压0.15V ，误差为多少？ 答：我们以0.5级的三位半表为例，()一个字+±=?x U U %5.0。 2V 量程测量直流电压1.5V 时 ()mV mV V U 5.815.1%5.0±=+?±=? 2V 量程测量直流电压0.15V 时()mV mV V U 8.1115.0%5.0±≈+?±=? 200mV 量程测量直流电压0.15V 时()mV mV mV U 9.01.0150%5.0±≈+?±=? 可见，测量小电压尽量选用低量程档。 5．为什么不宜用数字万用电表的电阻档测量表头内阻？ 答：数字万用电表电阻档内置9V 电池，而微安表头内阻在2000Ω左右。这样测通过表头的电流估计为mA A 5.40045.020009==，这个电流远大于微安表头的满量程电流。 6．为什么不能用数字万用电表的电阻档测量电源内阻？ 答：电阻档的使用条件是被测电阻中无电流通过，或者被测电阻两端无电压。对电源内阻来说，一旦用电阻档测量，电源就为内阻提供了电流，这样容易烧毁电表。 实验十二 制流和分压电路 【思考题参考答案】
电磁场与电磁波第二章课后答案
第二章 静电场 重点和难点 电场强度及电场线等概念容易接受，重点讲解如何由物理学中积分形式的静电场方程导出微分形式的静电场方程，即散度方程和旋度方程，并强调微分形式的场方程描述的是静电场的微分特性或称为点特性。 利用亥姆霍兹定理，直接导出真空中电场强度与电荷之间的关系。通过书中列举的4个例子，总结归纳出根据电荷分布计算电场强度的三种方法。 至于媒质的介电特性，应着重说明均匀和非均匀、线性与非线性、各向同性与各向异性等概念。讲解介质中静电场方程时，应强调电通密度仅与自由电荷有关。介绍边界条件时，应说明仅可依据积分形式的静电场方程，由于边界上场量不连续，因而微分形式的场方程不成立。 关于静电场的能量与力，应总结出计算能量的三种方法，指出电场能量不符合迭加原理。介绍利用虚位移的概念计算电场力，常电荷系统和常电位系统，以及广义力和广义坐标等概念。至于电容和部分电容一节可以从简。 重要公式 真空中静电场方程： 积分形式： ? = ?S S E 0 d εq ?=?l l E 0d 微分形式： 0 ερ= ??E 0=??E 已知电荷分布求解电场强度： 1，)()(r r E ?-?=； ? ' '-'= V V 0 d ) (41)(| r r |r r ρπε ? 2，? ' ''-'-'= V V 3 d |4) )(()(| r r r r r r E πε ρ 3， ? = ?S S E 0 d εq 高斯定律
介质中静电场方程： 积分形式： q S =?? d S D ?=?l l E 0d 微分形式： ρ=??D 0=??E 线性均匀各向同性介质中静电场方程： 积分形式： ε q S = ?? d S E ?=?l l E 0d 微分形式： ε ρ= ??E 0=??E 静电场边界条件： 1，t t E E 21=。对于两种各向同性的线性介质，则 2 21 1εεt t D D = 2，s n n D D ρ=-12。在两种介质形成的边界上，则 n n D D 21= 对于两种各向同性的线性介质，则 n n E E 2211εε= 3，介质与导体的边界条件： 0=?E e n ； S n D e ρ=? 若导体周围是各向同性的线性介质，则 ε ρS n E = ； ε ρ?S n - =?? 静电场的能量：
工程电磁场实验报告
工程电磁场实验报告 姓名： 学号： 联系式： 指导老师：
实验一螺线管电磁阀静磁场分析 一、实验目的 以螺线管电磁阀静磁场分析为例，练习在 MAXWELL 2D 环境下建立磁场模型，并求解分析磁场分布以及磁场力等数据。 二、主要步骤 a) 建立项目：其中包括生成项目录，生成螺线管项目，打开新项目 与运行MAXWELL 2D。 b) 生成螺线管模型：使用MAXWELL 2D 求解电磁场问题首先应该选择求解 器类型，静磁场的求解选择Magnetostatic，然后在打开的新项目中定义画图平面，建立要求尺寸的螺线管几模型，螺线管的组成包括 Core 、Bonnet 、Coil 、Plugnut、Yoke。 c) 指定材料属性：访问材料管理器，指定各个螺线管元件的材料，其中部分 元件的材料需要自己生成，根据给定的BH 曲线进行定义。 图1 元件材料 图2 B-H曲线 d) 建立边界条件和激励源：给背景指定为气球边界条件，给线圈Coil 施加电 流源。 e) 设定求解参数：本实验中除了计算磁场，还需要确定作用在螺线管铁心上 的作用力，在求解参数中要注意进行设定。
f) 设定求解选项：建立几模型并设定其材料后，进一步设定求解项，在对话 框Setup Solution Options 进入求解选项设定对话框，进行设置。 三、实验要求 建立螺线管电磁阀模型后，对其静磁场进行求解分析，观察收敛情况，画各种收敛数据关系曲线，观察统计信息；分析 Core 受的磁场力，画磁通量等势线，分析P lugnut 的材料磁饱和度,画出其B H 曲线。通过工程实例的运行，掌握软件的基本使用法。 四、实验结果 1.螺线管模型 图3 2.自适应求解 图4 收敛数据
电磁场及电磁波实验报告
电磁场与电磁波 实验报告 实验名称：有限差分法解电场边值问题 实验日期：2012年12月8日 姓名：赵文强 学号：100240333 XX工业大学（威海）
问题陈述 如下图无限长的矩形金属导体槽上有一盖板，盖板与金属槽绝缘，盖板电位为U0，金属槽接地，横截面如图所示，试计算此导体槽内的电位分布。 参数说明：a=b=10m, U=100v 实验要求 1)使用分离变量法求解解析解； 2)使用简单迭代发求解，设-10 =100.1,1 x y ε?=?= ，两种情况分别求解数值解； 3)使用超松弛迭代法求解，设-10 =100.1 x y ε?=?= ，确定?（松弛因子）。 求解过程 一、分离变量法求解 因为矩形导体槽在z方向为无限长，所以槽内电位函数满足直 角坐标系中的二维拉普拉斯方程。 22 22 (0,)0,(,)0(0) (,0)0,(,)(0) x y y a y y b x x b U x a ?? ?? ?? ?? += ?? ==≤≤ ==≤≤
根据边界条件可以确定解的形式： 1ππ(,)sin()sinh()n n n x n y x y A a a ?∞ ='=∑ 利用边界条件0(,)x b U ?=求解系数。 01 ππsin( )sinh()n n n x n b A U a a ∞ ='=∑ 01 πsin( )n n n x U f a ∞ ==∑ 0 0041,3,5,2πsin()d π 2,4,6,a n U n n x f U x n a a n ?=? ==??=? ? 011 πππsin()sinh()sin()n n n n n x n b n x A U f a a a ∞ ∞ =='==∑∑ 041,3,5,πsinh(π/) 'πsinh()02,4,6,n n U n f n n b a A n b n a ? =? ==??= ? 01,3,5, 4ππ(,)sin()sinh()πsinh(π/)n U n x n y x y n n b a a a ?∞ == ∑ 简单迭代法求解 二、 有限差分法 有限差分法（Finite Differential Method ）是基于差分原理的一种数值计算法。其基本思想：将场域离散为许多小网格，应用差分原理，将求解连续函数?的泊松方程的问题转换为求解网格节点上?的差分方程组的问题。 泊松方程的五点差分格式 )(4 1 4243210204321Fh Fh -+++=?=-+++?????????? 当场域中,0=ρ得到拉普拉斯方程的五点差分格式
电磁场与电磁波实验指导书要点
电磁场电磁波实验 实验一电磁感应定律的验证 一、实验目的 1、通过电磁感应装置的设计，了解麦克斯韦电磁感应定律的内容 2、了解半波天线感应器的原理及设计方法 3、天线长短与电磁波波长的接收匹配关系 二、预习要求 1、麦克斯韦电磁理论的内容 2、什么是电偶极子？ 3、了解线天线基本结构及其特性 三、实验仪器 HD-CB-IV电磁场电磁波数字智能实训平台：1套 电磁波传输电缆：1套 平板极化天线：1副 半波振子天线：1副 感应灯泡：1个 四、实验原理 麦克斯韦电磁理论经验定律包括：静电学的库仑定律，涉及磁性的定律，关于电流的磁性的安培定律，法拉第电磁感应定律。麦克斯韦把这四个定律予以综合，导出麦克斯韦方程，该方程组系统而完整地概括了电磁场的基本规律，并预言了电磁波的存在。麦克斯韦提出的涡旋电场和位移电流假说的核心思想是：变化的磁场可以激发涡旋电场，变化的电场可以激发涡旋磁场；电场和磁场不是彼此孤立的，它们相互联系、相互激发组成一个统一的电磁场。下面我们通过制作感应天线体，来验证电磁场的存在。 如图示：电偶极子是一种基本的辐射单元，它是一段长度远小于波长的直线电流元，线上的电流均匀同相，一个作时谐振荡的电流元可以辐射电磁波，故又称为元天线，元天线是最基本的天线。电磁感应装置的接收天线可采用多种天线形式，相对而言性能优良，但又容易制作，成本低廉的有半波天线、环形天线、螺旋天线等。
本实验重点介绍其中的一种半波天线。 半波天线又称半波振子，是对称天线的一种最简单的模式。对称天线（或称对称振子）可以看成是由一段末端开路的双线传输线形成的。这种天线是最通用的天线型式之一，又称为偶极子天线。而半波天线是对称天线中应用最为广泛的一种天线，它具有结构简单和馈电方便等优点。 半波振子因其一臂长度为λ /4 ，全长为半波长而得名。其辐射场可由两根单线驻波天线的辐射场相加得到，于是可得半波振子（L= λ /4 ）的远区场强有以下关系式： │ E │ =[60 Im cos( π cos θ /2)]/R 。sin θ=[60 Im/R 。] │ f( θ ) │ 式中，f( θ ) 为方向函数。对称振子归一化方向函数为│ F( θ ) │ = │ f( θ ) │ / fmax=|cos( π cos θ /2)/sin θ | 其中fmax 是f( θ ) 的最大值。由上式可画出半波振子的方向图如下: 半波振子方向函数与ψ无关，故在H 面上的方向图是以振子为中心的一个圆，即为全方性的方向图。在E 面的方向图为8 字形，最大辐射方向为θ = π /2 ，且只要一臂长度不超过0.625 λ，辐射的最大值始终在θ = π /2 方向上；若继续增大L ，辐射的最大方向将偏离θ = π /2 方向。 五、实验步骤 （一）测量电磁波发射频率 1、用N型电缆直接将“输出口1”连接至“功率频率检测口”。 2、在液晶界面上同时显示出发射功率及频率。 3、已知电磁波发射源的频率F,求得波长：λ=F V光，比如，电磁波发射源频率为900MHz,
电磁场与电磁波部分课后答案_郭辉萍版1-6章
第一章 习题解答 1.2解：⑴.A a =A A =149A ++ =（x a +2y a -3z a ）/14 ⑵cos A B θ =A ·B /A B A B θ=135.5o ⑶A ·B =-11, A ?B =-10x a -y a -4z a ⑷A ·（B ?C ）=-42 （A ?B ）·C =-42 ⑸A ?（B ?C ）=55x a -44y a -11z a （A ?B ）?C =2x a -40y a +5z a 1.3有一个二维矢量场F(r) =x a （-y ）+y a (x)，求其矢量线方程，并定性画出该矢量场图 形。 解：由dx/(-y)=dy/x,得2x +2y =c 1.6求数量场ψ=ln （2x +2y +2z ）通过点P （1，2，3）的等值面方程。 解：等值面方程为ln （2x +2 y +2z ）=c 则c=ln(1+4+9)=ln14 那么2x +2 y +2z =14 1.9求标量场ψ（x,y,z ）=62 x 3 y +z e 在点P （2，-1，0）的梯度。 解：由ψ?=x a x ψ??+y a y ψ??+z a z ψ??=12x 3y x a +182x 2y y a +z e z a 得 ψ?=-24x a +72y a +z a 1.10 在圆柱体2 x +2 y =9和平面x=0，y=0,z=0及z=2所包围的区域，设此区域的表面为S: ⑴求矢量场A 沿闭合曲面S 的通量，其中矢量场的表达式为 A =x a 32 x +y a （3y+z ）+z a (3z -x)
错误！未找到引用源。验证散度定理。 解：⑴??s d A =?? 曲+A d S ?? xoz +A d S ?? yoz +A d S ?? 上+A d S ?? 下 A d S ?? 曲 =232 (3cos 3sin sin )z d d ρθρθθρθ++?曲 =156.4 A d S ?? xoz = (3)y z dxdz +? xoz =-6 A d S ?? yoz =- 2 3x dydz ? yoz =0 A d S ?? 上 +A d S ?? 下=(6cos )d d ρθρθρ-?上+cos d d ρθρθ?下 =272π ? ?s d A =193 ⑵dV A V ???=(66)V x dV +?=6(cos 1)V d d dz ρθρθ+?=193 即：??s s d A =dV A V ??? 1.13 求矢量A =x a x+y a x 2y 沿圆周2x +2y =2 a 的线积分，再求A ?? 对此圆周所包围的表 面积分，验证斯托克斯定理。 解：??l l d A =2 L xdx xy dy +? =44a π A ?? =z a 2 y ????S s d A =2S y dS ? =22sin S d d θ ρρρθ? =44a π 即：??l l d A =????S s d A ，得证。 1.15求下列标量场的梯度： ⑴u=xyz+2 x u ?=x a u x ??+y a u y ??+z a u z ??=x a (yz+zx)+y a xz+z a xy ⑵u=42 x y+2 y z -4xz u ?=x a u x ??+y a u y ??+z a u z ??=x a (8xy-4z)+y a (42 x +2yz)+z a (2y -4x) ⑶u ?=x a u x ??+y a u y ??+z a u z ??=x a 3x+y a 5z+z a 5y
电磁场实验报告
电磁场实验报告 姓名：KZY 班级：自动化1405 学号：090114050X 时间：2016年10月23日
实验名称单缝衍射实验、自由空间中电磁波参量的测量 一、实验目的 1、了解电磁波的空间传播特性 2、通过对电磁波波长、波幅和波节的测量进一步了解和认识电磁 波。 3、利用电磁波的干涉原理，研究均匀无耗媒质εr的测量方法。 4、熟悉均匀无耗媒质分界面对电磁波的反射和透射特性。 二、实验仪器设备 1、单缝衍射仪器配置 2、单缝衍射板 3、半透射板 4、全反射板 三、实验原理 1、单缝衍射原理 查阅参考书籍可知，当一平面波入射到一宽度和波长可比拟的狭缝时，就要发生衍射的现象。在缝后面出现的衍射波强度并不是均匀的，中央最强，同时也最宽。在中央的两侧衍射波强度迅速减小，直至出现衍射波强度的最小值，即一级极小，此时衍射角为Фmin=sin-1λ/α。其中λ是波长，α是狭缝宽度。两者取同一长度单位，然后，随着衍射角增大，衍射波强度又逐渐增大，直至出现一级极大值，角
度为：Фmin=sin-1（3/2·λ/α）。 2、迈克尔逊干涉原理 由于两列波存在一定关系的波程差，两列波将发生干涉。而两列波发生干涉，存在合成振幅会出现最大与最小的情况。实验中，为了提高测量波长的精确度，测量多个极小值的位置，设S0为第一个极小值的位置吗，S n为第（n+1）个极小值的位置，L=|S n-S0|，则波长λ=2L/n。 三、实验内容与实验步骤 （1）单缝衍射实验 1、打开DH1121B的电源； 2、将单缝衍射版的缝宽α调整为70mm左右，将其安放在刻度盘上，衍射版的边线与刻度盘上两个90°对齐。
国中学生物理竞赛实验指导书思考题参考答案电磁学
国中学生物理竞赛实验指导书思考题参考答案电磁学
实验十一 学习使用数字万用表 【思考题参考答案】 1．调节电阻箱的电阻为500和5时，电阻的误差是多大？ 答：以0.1级为例，以每个接点的接触电阻按0.002为例。 对于500 电阻从0和9999接线柱输出，误差为 Ω=?+?=?512.0002.06%1.0500R 对于5电阻从0和9.9接线柱输出，误差为 Ω=?+?=?009.0002.02%1.05R 2．电源电压为110V 。是否能够只用一个电阻箱控制，得到0.5A 的电流？ 答：若只用一个电阻箱控制，所需电阻为Ω2205.0110==R 。这需要电阻箱的100?R 档，此档允许电流为0.05A ，实际电流大于额定电流，不能使用。 3．对于一块四位半的数字万用电表的直流电压200mV 量程，可能出现的最大数字是多少？最小分辨率是多少？ 答：最大数字为199.99mV 。最小分辨率为0.01mV 。 4．使用数字万用电表的直流电压2V 量程测量直流电压，测量值为 1.5V ，测量误差为多少？如果测量值为0.15V ，测量误差为多少？如果换用200mV 量程测量直流电压0.15V ，误差为多少？ 答：我们以0.5级的三位半表为例，()一个字+±=?x U U %5.0。 2V 量程测量直流电压1.5V 时 ()mV mV V U 5.815.1%5.0±=+?±=? 2V 量程测量直流电压0.15V 时()mV mV V U 8.1115.0%5.0±≈+?±=? 200mV 量程测量直流电压0.15V 时
()mV mV mV U 9.01.0150%5.0±≈+?±=? 可见，测量小电压尽量选用低量程档。 5．为什么不宜用数字万用电表的电阻档测量表头内阻？ 答：数字万用电表电阻档内置9V 电池，而微安表头内阻在 左右。这样测经过表头的电流估计为mA A 5.40045.020009==，这个电流远大于微安表头的满量程电流。 6．为什么不能用数字万用电表的电阻档测量电源内阻？ 答：电阻档的使用条件是被测电阻中无电流经过，或者被测电阻两端无电压。对电源内阻来说，一旦用电阻档测量，电源就为内阻提供了电流，这样容易烧毁电表。 实验十二 制流和分压电路 【思考题参考答案】 1．在连接分压电路时，有人将电源的 正、负极经过开关分别连到变阻器的一个 固定端和滑动端。这种连接方法对么？会 有什么问题？ 答：电路如图，这种连接方法不对。这种电路负载电阻被短路不会分压。 2．有一分压电路如图（实验的那个电路），负载电阻Ω=k R L 1.5，电压表内阻为Ω=k R V 10，变阻器电阻为0R 。 （1）若希望分压均匀，应选择哪种规格的变阻器？ （a ）A 1,5Ω；（b ）A 5.0,100Ω（c ）A 2.0,1000Ω E K A B C R
北邮电磁场与微波测量实验报告实验五极化实验
北邮电磁场与微波测量实验报告 实验五极化实验 学院：电子工程学院 班号：2011211204 组员： 执笔人： 学号：2011210986
一、实验目的 1．培养综合性设计电磁波实验方案的能力 2．验证电磁波的马吕斯定理 二、实验设备 S426型分光仪 三、实验原理 平面电磁波是横波，它的电场强度矢量E 和波长的传播方向垂直。如果E 在垂直于传播方向的平面沿着一条固定的直线变化，这样的横电磁波叫线极化波。在光学中也叫偏振波。偏振波电磁场沿某一方向的能量有一定关系。这就是光学中的马吕斯定律： 2 0cos I I θ = 式中I 为偏振波的强度，θ为I 与I0间的夹角。 DH926B 型分光仪两喇叭口面互相平行，并与地面垂直，其轴线在一条直线上，由于接收喇叭是和一段旋转短波导连在一起的；在该轴承环的90度围，每隔5度有一刻度，所以接收喇叭的转角可以从此处读到。 四、实验步骤 1．设计利用S426型分光仪验证电磁波马吕斯定律的方案； 根据实验原理，可得设计方案：将S426型分光仪两喇叭口面互相平行，并与地面垂直，其轴线在一条直线上，由于接收喇叭是和一段旋转短波导连在一起的；在该轴承环的90度围，每隔5度有一刻度，接收喇叭课程从此处读取θ（以10度为步长），继而进行验证。 2．根据设计的方案，布置仪器，验证电磁波的马吕斯定律。 实验仪器布置 通过调节，使A1取一较大值，方便实验进行。 然后，再利用前面推导出的θ，将仪器按下图布置。
五、实验数据 I（uA） θ° 0 10 20 30 40 50 60 70 80 90 理论值90 87. 3 79. 5 67. 5 52. 8 37. 2 22. 5 10. 5 2.7 0 实验值90 88 82 69 54 37 20 8 2 0.2 相对误差% 0 0.8 0.6 2.2 2.3 0.5 11. 1 14. 3 25. 9 - 1、数据分析： 由数据可看出，实验值跟理论值是接近的，相对误差基本都很小，在误差允许围，所以可以认为马吕斯定律得到了验证。 2、误差分析： 实验中可能存在仪器仪表误差，人为误差以及各组互相影响造成的误差等。但是角度比较大的时候，相对误差都比较小，也比较精准。角度比较小的时候，由于理论值较小，相对误差会大一点，但是从整体趋势来看，结果也是合理的。所以不影响我们对马吕斯定律进行验证。 六、思考题 1、垂直极化波是否能够发生折射？为什么？给出推导过程。 答：不能。 A1
电磁场与电磁波实验指导书
静电场边值问题实验 对于复杂边界的静电场边值问题，用解析法求解很困难，甚至是不可能的。在实际求解过程中，直接求出静电场的分布或电位又很困难，其精度也难以保证。本实验根据静电场与恒定电流场的相似性用碳素导电纸中形成的恒定电流场来模拟无源区域的二维静电场，从而测出边界比较复杂的无源区域静电场分布。 一、 实验目的： 1、学习用模拟法测量静电场的方法。 2、了解影响实验精度的因素。 二、 实验原理： 在静电场的无源区域中，电场强度E '电位移矢量D '及电位Ф、满足下列方程： ▽×E 、= 0 ▽×D '= 0 D '=ε E 、 E 、 = - ▽φ 、 （1） 式中ε为静电场的介电常数。 在恒定电流场中，电场强度E 、电流密度J 及电位Ф满足下列方程： ▽×E = 0 ▽·J = 0 J = δE E =－▽Φ （2） 式中δ为恒定电流场中导电媒质的电导率。 因为方程组（1）与方程组（2）在形式上完全相似，所以φ、（静电场中的电位分布函数）与Φ（恒定电流场中的电位分布函数）应满足同样形式的微分方程。由方程组（1）和方程组（2）很容易求得： ▽·（ε▽φ、）= 0 （3） ▽·（δ▽Φ）= 0 （4） 式中ε与δ处于相应的位置，它们为对偶量。 若ε与δ在所讨论区域为均匀分布（即其值与坐标无关），则方程（3）、（4）均可简化为拉普拉斯方程： 2?φ'= 0 02=Φ? 电位场解的唯一定理可知：满足相同微分方程的两个电位场，它们具有相同的边界电位值，因此，在保证边界电位值不变的情况下，我们可以用恒定电流场的模型来模拟无源区域的静电场，当静电场中媒质为均匀媒质时，其导电媒质也应为均匀媒质，这样测得的恒定电流场的电位分布就是被模拟的静电场的电位分布，不需要任何改动。
电磁场与电磁波课后习题答案第一章
第一章 给定三个矢量A u r ，B u r ，C u r ： A u r =x a u u r +2y a u u r -3z a u u r B u r = -4y a u u r +z a u u r C u r =5x a u u r -2z a u u r 求：⑴矢量A u r 的单位矢量A a u u r ； ⑵矢量A u r 和B u r 的夹角AB θ； ⑶A u r ·B u r 和A u r ?B u r ⑷A u r ·（B u r ?C u r ）和（A u r ?B u r ）·C u r ； ⑸A u r ?（B u r ?C u r ）和（A u r ?B u r ）?C u r 解：⑴A a u u r =A A u r u r =u r （x a u u r +2y a u u r -3z a u u r ） ⑵cos AB θu r u r =A u r ·B u r /A u r B u r AB θ=135.5o ⑶A u r ·B u r =-11, A u r ?B u r =-10x a u u r -y a u u r -4z a u u r ⑷A u r ·（B u r ?C u r ）=-42 （A u r ?B u r ）·C u r =-42 ⑸A u r ?（B u r ?C u r ）=55x a u u r -44y a u u r -11z a u u r （A u r ?B u r ）?C u r =2x a u u r -40y a u u r +5z a u u r 有一个二维矢量场F(r)r =x a u u r （-y ）+y a u u r (x)，求其矢量线方程，并定性画出该矢量场图形。 解：由dx/(-y)=dy/x,得2x +2y =c 求数量场ψ=ln （2x +2y +2z ）通过点P （1，2，3）的等值面方程。 解：等值面方程为ln （2x +2y +2 z ）=c
电磁场实验指导书及实验报告
CENTRAL SOUTH UNIVERSITY 题目利用Matlab模拟点电荷电场的分布姓名xxxx 学号xxxxxxxxxx 班级电气xxxx班 任课老师xxxx 实验日期2010-10
电磁场理论 实验一 ——利用Matlab 模拟点电荷电场的分布 一．实验目的： 1．熟悉单个点电荷及一对点电荷的电场分布情况； 2．学会使用Matlab 进行数值计算，并绘出相应的图形； 二．实验原理： 根据库伦定律：在真空中，两个静止点电荷之间的作用力与这两个电荷的电量乘积成正比，与它们之间距离的平方成反比，作用力的方向在两个电荷的连线上，两电荷同号为斥力，异号为吸力，它们之间的力F 满足： R R Q Q k F ? 212 = (式1) 由电场强度E 的定义可知： R R kQ E ? 2 = (式2) 对于点电荷，根据场论基础中的定义，有势场E 的势函数为 R kQ U = (式3) 而 U E -?= (式4) 在Matlab 中，由以上公式算出各点的电势U ，电场强度E 后，可以用Matlab 自带的库函数绘出相应电荷的电场分布情况。 三．实验内容： 1. 单个点电荷 点电荷的平面电力线和等势线 真空中点电荷的场强大小是E=kq /r^2 ,其中k 为静电力恒量, q 为电量, r 为点电荷到场点P(x,y)的距离。电场呈球对称分布, 取电量q> 0, 电力线是以电荷为起点的射线簇。以无穷远处为零势点, 点电荷的电势为U=kq /r,当U 取
常数时, 此式就是等势面方程.等势面是以电荷为中心以r 为半径的球面。 平面电力线的画法 在平面上, 电力线是等角分布的射线簇, 用MATLAB 画射线簇很简单。取射线的半径为( 都取国际制单位) r0=, 不同的角度用向量表示( 单位为弧度) th=linspace(0,2*pi,13)。射线簇的终点的直角坐标为: [x,y]=pol2cart(th,r0)。插入x 的起始坐标x=[x; *x].同样插入y 的起始坐标, y=[y; *y], x 和y 都是二维数组, 每一列是一条射线的起始和终止坐标。用二维画线命令plot(x,y)就画出所有电力线。 平面等势线的画法 在过电荷的截面上, 等势线就是以电荷为中心的圆簇, 用MATLAB 画等势 线更加简单。静电力常量为k=9e9, 电量可取为q=1e- 9; 最大的等势线的半径应该比射线的半径小一点 r0=。其电势为u0=k8q /r0。如果从外到里取7 条等势线, 最里面的等势线的电势是最外面的3 倍, 那么各条线的电势用向量表示为: u=linspace(1,3,7)*u0。从- r0 到r0 取偶数个点, 例如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, 其中插入了标题等等。越靠近点电荷的中心, 电势越高, 电场强度越大, 电力线和等势线也越密。