点电荷分布matlab仿真,利用Matlab-模拟点电荷电场的分布的仿真 实验报告.pdf

利用Matlab-模拟点电荷电场的分布的仿真 实验报告

利用 Matlab模拟点电荷电场的分布

一、实验目的：

1．熟悉单个点电荷及一对点电荷的电场分布情况；

2．学会使用 Matlab 进行数值计算，并绘出相应的图形；

二、实验原理：

根据库伦定律：在真空中，两个静止点电荷之间的作用力与这两个电荷的电

量乘积成正比，与它们之间距离的平方成反比，作用力的方向在两个电荷的连线

上，两电荷同号为斥力，异号为吸力，它们之间的力 F 满足：

 Q Q

1 1

ˆ

F k R                        (式 1)

R 2

由电场强度 E 的定义可知：

 kQ ˆ

E R                        (式 2)

R 2

对于点电荷，根据场论基础中的定义，有势场 E 的势函数为

kQ

U                        (式 3)

R



而                             E U                         (式 4)

在 Matlab 中，由以上公式算出各点的电势 U，电场强度 E 后，可以用 Matlab

自带的库函数绘出相应电荷的电场分布情况。

三、实验内容

1、画单个点电荷的平面电场线与等势线

等势线就是以电荷为中心的圆簇，用Matlab画等势线更加简单。静电力常量

为k=9*e9，电量可取为q=1e‐19；最大的等势线的半径应该比射线的半径小一点，

r 0.1 。其电势为u k *q / r 。如果从外到里取7条等势线，最里面的等势线的

0 0 0

电势是最外面的3倍，那么各条线的电势用向量表示为：u linspace(1,3,7)*u0 。

从‐r0到r0取偶数个点，例如100个点，使最中心点的坐标绕过0，各点的坐标可用

向量表示：x linspace (r ,r ,100) ，在直角坐标系中可形成网格坐标：

0 0

[X ,Y] meshgrid(x) 。各点到原点的距离为：r X .^ 2 Y .^ 2 ，在Matlab中进行乘

方运算时，乘方号前面要加点，表示对变量中的元素进行乘方计算。各点的电势

为U k.* q./ r ；同样地，在进行除法运算时，除号前面也要加点，同样表示对

变量中的元素进行除法运算。用等高线命令contour X ,Y ,U ,u 即可画出等势线。

如图1所示：

图 1 单个点电荷的平面电场线与等势线

其代码如下：

%单个点电荷的平面电场线与等势线 u1=linspace(1,3,7)*u;

a=0.2; x=linspace(-0.1,0.1,100);

k=9*10^9;