多点电荷系统电力线和等位面计算机仿真

本博文依据“电磁场与波”课程设计，写篇博文，记录一下自己的心路历程---

---

前言

电磁场与电磁波是信号传输处理领域的基础学科之一，了解和掌握电磁波领域的相关知识对于信号处理而言至关重要。 电磁场与电磁波研究基础之一在于点电荷问题。电磁场领域中的诸多问题都可以看作点电荷问题的扩展----即多点电荷系统，因此了解和掌握点电荷以及多点电荷系统性质对于电磁场问题的学习至关重要。

---

本篇博文以电磁场与波课程设计为依托，对于点电荷以及多点电荷系统进行研究，并借助于matlab编程，实现问题可视化，辅助问题理解与学习。

一、点电荷

点电荷，物理学上把本身的线度比相互之间的距离小得多的带电体叫做点电荷。相当于运动学的“质点”模型。 电荷都是有体积，有大小的。电荷之间存在相互作用，同种电荷相互推斥，异种电荷相互吸引。在定量地研究电荷之间相互作用的时候，发现有些电荷的大小对所研究问题的结果带来的影响微不足道，这个时候就完全可以把电荷的体积和大小忽略掉，把*电荷看做只有电量，没有大小的电荷*，这就是点电荷模型。 点电荷模型是电磁场问题研究中的基本模型，对于点电荷的理解和掌握对于电磁场领域的学习是至关重要的。

二、多点电荷系统

1.多点电荷系统简介

双点电荷系统与点电荷模型类似，同样是电磁场问题中的基本模型，可抽象为电磁场问题中的“电偶极子模型”。 根据库仑定律： 在真空中，两个静止点电荷之间的作用力与这两个电荷的电量乘积成正比，与它们之间距离的平方成反比，作用力的方向在两个电荷的连线上，两电荷同号为排斥力，异号为吸引力，它们之间的作用力F满足：

其中$ϵ$ 被称为介电常数，电场强度E的计算公式：


由电场强度计算公式可以推导出电势的计算公式：
在Matlab编程中，先求出各点的电势φ 后，可以用 Matlab 自带的库函数求出E并绘出相应电荷的电场分布情况。

2.多点电荷系统仿真编码

代码如下（示例）：

map_size=[20,20];    % 图窗大小
a=[-4,0;4,0];        % 电荷位置 eg.[-4,0;4,0]
Q=[-1,1];            % 每个电荷带电量 eg.[-1,1]
% 程序编写设计----可自由更改电荷的数量和位置


radius = 1;          % 点电荷的半径
k = 9e+9;            % k = 1/(4πε)
limitdist = 1;       % 避免出现无限大情况
q = 1.602e-19;       % 元电荷电量初始化
angles = linspace(-pi, pi, 20);
x=map_size(1);
y=map_size(2);
[X, Y] = meshgrid(-x : 2*x/80 : x, -y : 2*y/80 : y); % 张成二维网格平面
ls=size(a);
l=ls(1);
r={};
for i=(1:l)
    r{1,i}=sqrt((X - a(i,1)).^2 + ((Y-a(i,2)).^2 + limitdist)); % 网格平面中的每个点到每个电荷之间的距离
end

hold on;
V1=zeros(size(X));
for i=(1:l)
    V1=V1+Q(i)*q*(k./r{1,i});  % 每个电荷对（x，y）点的电势均有影响，累加求得电势
end

%绘制等势线
min1 = min(min(V1));
max1 = max(max(V1));
range = linspace(min1, max1, 100);
contour(X, Y, V1, range);

%绘制电力线
[Ex, Ey] = gradient(-V1); % 电势的负梯度求解电场强度
for i=(1:l)
    if Q(i)>0
        factor=1;
    else
        factor=-1;
    end
    startlx = a(i,1) + radius*cos(angles);
    startrx = a(i,1) + radius*cos(angles);
    startly = a(i,2) + radius*sin(angles);
    startry = a(i,2) + radius*sin(angles);
    streamline(X, Y, factor*Ex, factor*Ey, startlx, startly);
    streamline(X, Y, factor*Ex, factor*Ey, startrx, startry);
end

for i=(1:l)
    %标记点电荷
    q1=Q(i);
    locx=a(i,1);
    locy=a(i,2);
    if(q1>0)
        plot(locx, locy, 'r+', 'MarkerSize', 5);
    elseif(q1<0)
        plot(locx, locy, 'ro', 'MarkerSize', 5);
    end
end
legend('等位线','电场线')
title('多点电荷电力线与等位线示意图');

3.多点电荷系统仿真结果展示

“电偶极子”电力线和等位面计算机仿真：

二维平面平行板电力线及等位面计算机仿真：

总结

研究多点电荷系统电力线和等位面分布情况对于理解电磁场基本问题具有非常重要的意义。本文通过matlab编程将多点电荷周围电力线和等位面情况做可视化分析，以清晰直观的方式展现多点电荷系统自身性质以及周围场的分布性质，对于电磁场问题的理解和掌握有非常大的帮助。