电磁学实验党福利：用Python模拟静电场与磁场的完整教程（附代码）

最新推荐文章于 2026-05-31 07:10:15 发布

原创

最新推荐文章于 2026-05-31 07:10:15 发布 · 474 阅读

标签

#Python模拟 #电磁学 #静电场 #磁场可视化

电磁学实验党福利：用Python模拟静电场与磁场的完整教程（附代码）

还记得大学物理课上，面对那些抽象的电场线、磁感线图，你是否也曾感到一丝困惑？公式推导固然严谨，但总感觉少了点什么——一种直观的、可以亲手“把玩”的体验。电磁学的美妙之处，恰恰在于其规律背后那幅可以“看见”的物理图景。今天，我们不再满足于纸笔计算，而是拿起Python这个强大的工具，将库仑定律、高斯定理、安培环路定理等核心原理，转化为屏幕上动态可视的电场与磁场分布图。这不仅是学习方式的革新，更是理解深度的跃迁。无论你是正在啃电磁学教材的学生，还是对物理模拟充满好奇的技术爱好者，或是希望将理论应用于仿真项目的工程师，这篇教程都将带你从零开始，构建一套属于自己的电磁场可视化实验平台。我们将使用NumPy进行高效的数值计算，并用Matplotlib绘制出媲美教科书的精美图像，更重要的是，整个过程完全由你掌控，参数随心调，现象即时现。

1. 环境搭建与基础原理回顾

工欲善其事，必先利其器。我们的“数字实验室”建立在Python的科学计算生态之上。首先，确保你的环境中已安装以下核心库：

pip install numpy matplotlib scipy

NumPy：负责所有数组运算和数学计算的核心，没有它，我们寸步难行。
Matplotlib：我们的“画笔”和“画布”，用于将数值结果转化为直观的图形。
SciPy：虽然不是必须，但其scipy.integrate等模块在处理复杂积分时会非常方便，可以作为备选工具。

安装完成后，在代码开头导入它们，这是我们的标准起手式：

import numpy as np
import matplotlib.pyplot as plt
# 为了显示中文和负号，通常需要额外设置
plt.rcParams['font.sans-serif'] = ['SimHei', 'DejaVu Sans'] # 用来正常显示中文标签
plt.rcParams['axes.unicode_minus'] = False # 用来正常显示负号

在开始编码前，我们快速回顾两个基石性的概念，这决定了我们模拟的“世界观”。

静电场的本质是电荷对周围空间产生的影响。一个点电荷Q在距离r处产生的电场强度E，由库仑定律给出： E = k * Q / r^2，方向沿径向。多个电荷产生的总电场，是每个电荷产生电场的矢量叠加。而电势V则是电场的标量积分，反映了单位正电荷在某点的势能。

稳恒磁场则由运动的电荷（电流）产生。一段电流元Idl在空间某点产生的磁感应强度dB，由毕奥-萨伐尔定律描述，其大小与电流元、距离和夹角有关。对于具有高度对称性的电流分布（如无限长直导线、螺线管），使用安培环路定理计算磁场会简便得多：磁场沿任意闭合环路的线积分，等于穿过该环路所围面积的电流代数和乘以一个常数。

我们的模拟，就是将上述物理定律离散化、数值化。我们将空间网格化，在每个网格点上计算电场或磁场的矢量，然后绘制出来。

2. 静电场可视化：从点电荷到复杂体系

让我们从最简单的系统开始：一个孤立的点电荷。这是理解一切复杂电场的基础。

2.1 单个与多个点电荷的电场

我们首先在二维平面上定义一个计算区域，比如从-2米到2米。在这个区域内生成网格点，然后计算每个点上的电场矢量。

def electric_field_point_charge(q, pos_q, x, y):
    """
    计算点电荷在网格点(x, y)上产生的电场。
    参数:
        q: 电荷量（库仑），可正可负。
        pos_q: 电荷位置 (x0, y0)。
        x, y: 网格坐标数组。
    返回:
        Ex, Ey: 电场在x和y方向的分量。
    """
    k = 8.9875517873681764e9  # 真空介电常数，单位 N·m²/C²
    # 计算网格上每点到电荷的距离向量分量
    dx = x - pos_q[0]
    dy = y - pos_q[1]
    r = np.sqrt(dx**2 + dy**2)
    r3 = r**3
    r3[r == 0] = np.inf  # 避免除以零，电荷所在点电场无穷大

    Ex = k * q * dx / r3
    Ey = k * q * dy / r3
    return Ex, Ey

现在，让我们绘制一个正电荷的电场。电场线从正电荷发出，我们使用plt.streamplot来绘制流线，它能很好地表现矢量的方向和连续性。

# 创建计算网格
x = np.linspace(-2, 2, 50)
y = np.linspace(-2, 2, 50)
X, Y = np.meshgrid(x, y)

# 定义正电荷 (+1nC) 位于原点
charge_q = 1e-9
charge_pos = (0, 0)

# 计算电场
Ex, Ey = electric_field_point_charge(charge_q, charge_pos, X, Y)

# 绘制
fig, ax = plt.subplots(figsize=(7, 7))
# 绘制电场线
ax.streamplot(X, Y, Ex, Ey, color='blue', linewidth=1, density=1.5, arrowsize=1.5)
# 标记电荷位置
ax.scatter(*charge_pos, c='red', s=200, marker='o', label=f'+{charge_q/1e-9:.0f} nC')
ax.set_xlabel('x (m)')
ax.set_ylabel('y (m)')
ax.set_title('孤立正点电荷的电场线分布')
ax.legend()
ax.set_aspect('equal')
plt.grid(True, alpha=0.3)
plt.show()

运行这段代码，你将看到一个标准的辐射状电场线图。将charge_q改为负值（如-1e-9），重新运行，观察电场线方向如何变为指向电荷中心。