电磁场点电荷电场线电势MATLAB仿真中南大学.docxVIP

电磁场理论 实验一 ——利用 Matlab 模拟点电荷电场的分布 一．实验目的： 1．熟悉单个点电荷及一对点电荷的电场分布情况； 2．学会使用 Matlab 进行数值计算，并绘出相应的图形； 二．实验原理： 根据库伦定律：在真空中，两个静止点电荷之间的作用力与这两个电荷的电 量乘积成正比， 与它们之间距离的平方成反比， 作用力的方向在两个电荷的连线 上，两电荷同号为斥力，异号为吸力，它们之间的力 F 满足： Q1Q2 ? F k R2 R (式 1) 由电场强度 E 的定义可知： kQ ? E R2 R (式 2) 对于点电荷，根据场论基础中的定义，有势场 E 的势函数为 U kQ (式 3) R 而 E U (式 4) 在 Matlab 中，由以上公式算出各点的电势 U，电场强度 E 后，可以用 Matlab 自带的库函数绘出相应电荷的电场分布情况 . .实验内容： 单个点电荷 点电荷的平面电力线和等势线 真空中点电荷的场强大小是 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, 其中插入了标题等等 .越靠近点电荷的中心 , 电势越高 , 电场强度越大 , 电力线和等势线也越密 . 单个点电荷的电场线与等势线 0.15 0.1 0.05 0 y -0.05 -0.1 -0.15 -0.2 -0.15 -0.1 -0.05 0 0.05 0.1 0.15 -0.2 x 图 1 点电荷的立体电力线和等势面 立体电力线的画法 先形成三维单位球面坐标 , 绕 z 轴一周有 8 条电力线 [X,Y,Z]=sphere(8), 每维都是 9×9 的网格矩阵 , 将 X 化为行向量 , 就形成各条电力线的终点 x 坐标 x=r0=X(:) ′ , 其他两个坐标也可同样形成终点坐标 y=r0+Y(:) , z=r 0+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 的