电磁场与电磁波 双电荷静电系统 电力线 等位面（MATLAB模拟）

电磁场与电磁波 第一阶段大作业（XD）

题目：采用计算机编程绘制下列双电荷静电系统的电力线分布和等位面分布（二维xy平面分布）

该程序模拟了四种情况
（默认q2=e）
（1）q1=q2
（2）q1=-q2
（3）q1=2q2
（4）q1=-2q2
如图所示四种情况模拟的结果

以下列出q1=q2情况下模拟的源代码

clear;
radius = 1; %点电荷的半径
x = 20; y = 20;
k = 9e+9; %k = 1/(4πε)
limitdist = 1; %避免出现无限大情况
a = 6; %点电荷的位置
angles = linspace(-pi, pi, 20);
starty = radius*sin(angles);

[X, Y] = meshgrid(-x : 0.5 : x, -y : 0.5 : y);
r1 = sqrt((X + a).^2 + Y.^2 + limitdist);
r2 = sqrt((X - a).^2 + Y.^2 + limitdist);

q = 1.602e-19;
q1 = q;
q2 = q;
V1 = q1 * (k./r1) + q2 * (k./r2);
hold on;

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

%绘制电力线
[Ex, Ey] = gradient(-V1);
startlx = -a + radius*cos(angles);
startrx = a - radius*cos(angles);
streamline(X, Y, Ex, Ey, startlx, starty);
streamline(X, Y, Ex, Ey, startrx, starty);

%标记点电荷
if(q1>0)
    plot(-a, 0, 'r+', 'MarkerSize', 5);
elseif(q1<0)
    plot(-a, 0, 'ro', 'MarkerSize', 5);
end
if(q2>0)
    plot(a, 0, 'r+', 'MarkerSize', 5);
elseif(q2<0)
    plot(a, 0, 'ro', 'MarkerSize', 5);
end
title('q1 = q2');