MATLAB小白也能看懂的电场仿真：手把手教你用代码画三电荷电场线与等势面-开发者社区

MATLAB零基础实战：三电荷系统电场可视化全解析

刚接触电磁场仿真的同学往往会被复杂的公式和编程吓退，但今天我要分享的这套方法，能让没有任何MATLAB基础的小白也能轻松绘制出专业级的电场分布图。我们以经典的正三角形三电荷系统为例，从原理到代码逐行解析，带你体验物理与编程碰撞出的火花。

1. 环境准备与基础概念

在开始之前，确保你的电脑已经安装了MATLAB R2018a或更高版本。打开MATLAB后，我们会使用最基本的绘图功能，不需要任何额外工具箱。先来理解几个核心概念：

电场线：描述电场方向的假想曲线，切线方向表示该点电场方向
等势面：电势相同的点构成的面，在二维情况下表现为等势线
叠加原理：多个电荷产生的总电势/电场等于各电荷单独产生的电势/电场的矢量和

我们的目标是绘制出三个点电荷（+1C, -1C, +1C）在正三角形顶点分布时的电场线和等势线。这三个电荷的位置坐标分别为(5,0)、(-5,0)和(0,5)毫米。

提示：MATLAB中所有长度单位保持一致即可，本文使用厘米为单位，1cm=10mm

2. 参数设置与网格生成

首先我们需要定义计算区域和精度。在MATLAB命令窗口输入以下代码：

% 基本参数设置 number = 3; % 电荷数量 charge = [1, 0.5, 0; % 第一个电荷：+1C，位置(0.5,0)cm 1, -0.5, 0; % 第二个电荷：+1C，位置(-0.5,0)cm -1, 0, 0.866]; % 第三个电荷：-1C，位置(0,0.866)cm k = 9e9; % 静电力常数，单位N·m²/C² % 创建计算网格 x = -2.5:0.01:2.5; % x轴范围-2.5到2.5cm，步长0.01cm y = -2.5:0.01:2.5; % y轴相同设置 [X,Y] = meshgrid(x,y); % 生成501×501的网格点

这段代码建立了我们的计算坐标系。三个电荷位于边长为1cm的正三角形顶点（转换了单位）。网格精度设为0.01cm，确保图像足够精细。

3. 电势计算与等势线绘制

电势计算是绘制等势线的基础。根据叠加原理，空间某点总电势等于各电荷在该点产生电势的代数和。在MATLAB中继续输入：

% 初始化电势矩阵 U = zeros(size(X)); % 计算各点电势 for j = 1:number U = U + k * charge(j,1) ./ sqrt((X-charge(j,2)).^2 + (Y-charge(j,3)).^2); end % 处理电荷位置处的无穷大问题 for a = 1:number for b = -10:10 for c = -10:10 if abs(U(round(charge(a,3)*100)+251+b, round(charge(a,2)*100)+251+c)) > 1e10 U(round(charge(a,3)*100)+251+b, round(charge(a,2)*100)+251+c) = sign(U(round(charge(a,3)*100)+251+b, round(charge(a,2)*100)+251+c)) * 1e10; end end end end % 绘制等势线 contour(X,Y,U,100); % 绘制100条等势线 hold on; % 保持图形，准备叠加电场线

关键点解析：

U = U + k*q/r实现了电势叠加公式
电荷位置处电势理论上是无穷大，我们将其限制在1e10以内
contour函数自动选择100个电势值绘制等势线

4. 电场线计算与绘制

电场线从正电荷发出，终止于负电荷。我们采用数值方法逐步追踪电场线：

% 电场线参数设置 total_charge = sum(abs(charge(:,1))); % 总电荷量 px = cell(number,50*number); % 存储电场线x坐标 py = cell(number,50*number); % 存储电场线y坐标 % 计算每条电场线轨迹 for a = 1:number lines_per_charge = round(abs(charge(a,1))/total_charge * 50 * number); for b = 1:lines_per_charge % 设置起始点（电荷表面均匀分布） theta = b * 2 * pi / lines_per_charge; px{a,b}(1) = charge(a,2) + 0.05*sqrt(2) * cos(theta); py{a,b}(1) = charge(a,3) + 0.05*sqrt(2) * sin(theta); % 逐步追踪电场线 i = 1; while true Ex = 0; Ey = 0; % 计算当前点的电场强度 for c = 1:number Rx = px{a,b}(i) - charge(c,2); Ry = py{a,b}(i) - charge(c,3); E = k * charge(c,1) / (Rx^2 + Ry^2); Ex = Ex + E * Rx / sqrt(Rx^2 + Ry^2); Ey = Ey + E * Ry / sqrt(Rx^2 + Ry^2); end % 确定步长方向 step_size = 0.001; if charge(a,1) < 0 % 负电荷电场线反向 step_size = -step_size; end % 计算下一个点 norm_E = sqrt(Ex^2 + Ey^2); px{a,b}(i+1) = px{a,b}(i) + step_size * Ex/norm_E; py{a,b}(i+1) = py{a,b}(i) + step_size * Ey/norm_E; i = i + 1; % 检查终止条件 stop = false; for c = 1:number if ((px{a,b}(i) - charge(c,2))^2 + (py{a,b}(i) - charge(c,3))^2) < 0.005 stop = true; break; end end if stop || px{a,b}(i) < -2.5 || px{a,b}(i) > 2.5 || py{a,b}(i) < -2.5 || py{a,b}(i) > 2.5 break; end end end end % 绘制电场线 for a = 1:number for b = 1:size(px,2) if ~isempty(px{a,b}) plot(px{a,b}, py{a,b}, 'b'); end end end

5. 图像美化与结果分析

最后添加标签和装饰，使图像更专业：

% 添加标签和装饰 xlabel('X position (cm)'); ylabel('Y position (cm)'); title('Electric Field Lines and Equipotential Surfaces'); text(1.8, 2.3, 'Charge unit: C'); text(2.1, -2.35, 'Coordinate unit: cm'); % 标记电荷位置和电量 for a = 1:number text(charge(a,2)-0.04, charge(a,3), num2str(charge(a,1))); end % 设置坐标轴比例一致 axis equal; grid on;

运行完整代码后，你将看到一幅清晰的电场分布图，其中：