企业官网建设流程全解析

电偶极子这事儿，我在 MATLAB 里折腾了很久，踩了不少坑，也攒了不少经验。这个内容其实特别适合电磁场课程的作业、考研复习的辅助，或者像我一样非要可视化才觉得“学明白了”的强迫症患者。这篇博客我直接把思路和能跑的源码都拿出来，照着敲一遍就能出图。

1. 电偶极子模型梳理与整体设计思路

1.1 先想清楚电偶极子到底在算什么

电偶极子最经典的模型就是两个带电量为 +q 和 -q 的点电荷，距离为 d，放在空间里。我们关心的是这个电荷系统产生的电场分布和电势分布。很多同学一开始喜欢直接堆公式，但用 MATLAB 做场模拟，核心不是套公式，而是把“空间离散化”这件事想明白。

空间里的每一个点 (x, y)，都对应一个电势值 V(x, y)，而电场是电势的负梯度，即E = -grad(V)。如果所有点都手算，这不现实，MATLAB 强就强在能用矩阵并行处理几千几万个点。我的思路很简单：先建网格，再算每个点的电势，最后求梯度得到电场分量。

这里有个容易懵的地方：电势是标量场，电场是矢量场。用 MATLAB 做场模拟，这两样东西都要画出来才能体现出“场”的感觉。电势用等值线或伪彩图，电场用箭头图，两者叠在一张图上才是完整的电偶极子场分布图。

1.2 为什么选电势作为计算的中间量

我在第一次写代码时，直接算电场分量 Ex 和 Ey，发现代码长且容易出错。后来改用“先算电势 V，再梯度求电场”的思路，代码短了一半，逻辑也顺很多。这背后的道理是：电势是标量，对一个点只用做一次加法运算；而电场是矢量，对一个点要做两次除法运算，还得处理方向符号。

从数值角度看，通过电势求梯度还有一个额外的好处：梯度操作在 MATLAB 里本身就是针对矩阵设计的，gradient()函数一次调用就能完成整场计算，性能高，代码可读性强。

2. 核心源码拆解与参数选取逻辑

2.1 网格生成与电荷参数设置

我先贴一段核心代码，然后逐段解释为什么这么写：

% 参数设置 q = 1e-9; % 电荷量（库仑） d = 0.02; % 电荷间距（米） x = linspace(-0.05, 0.05, 201); y = linspace(-0.05, 0.05, 201); [X, Y] = meshgrid(x, y);

这里最关键的选择是电荷量单位。如果你直接把 q 设为 1，画出来的电势数值会非常大（因为库仑常数 k = 9e9），等值线图看起来会很拥挤，反而不利于观察电场线走势。我习惯用纳库（1e-9）量级，算出的电势在几百伏量级，等值线分布均匀，观感最好。

网格密度 201 x 201 是性能和精度的平衡点。太密（比如 501 x 501）导致 grad 运算变慢，操作起来有迟滞感；太疏（比如 51 x 51）则等值线会出现明显的折线，不光滑。对于演示和课程作业来说，201 x 201 够用。

2.2 电势计算与奇点处理

深入一点说，电偶极子在某点的电势公式是：

V = kq / r1 - kq / r2

其中 r1 和 r2 分别是该点到正电荷和负电荷的距离。直接用这个公式，会掉进一个经典的坑：电荷所在位置处的电势是无穷大。如果不对这些点做处理，算出来的矩阵里会有 Inf，画图时会出现难看的空白或异常色块。

我的做法是给距离加一个很小的偏移量，相当于在物理上把点电荷换成一个小球体，等效于给距离做正则化处理：

r1 = sqrt((X - d/2).^2 + Y.^2 + 0.002^2); r2 = sqrt((X + d/2).^2 + Y.^2 + 0.002^2); k = 9e9; V = k * q ./ r1 - k * q ./ r2;

加的 0.002 米偏移量，相当于把点电荷理解为半径 2 毫米的带电小球。这个“球体近似”在远离电荷的区域误差很小，但规避了无穷大的逻辑错误。

算电场：

[Ex, Ey] = gradient(-V);

对，就这一行。gradient 函数返回的是 V 在 x 和 y 方向的导数的负数，正好就是电场分量。不过要注意，gradient 默认按数据点间距为 1 计算。如果我们的 x 和 y 不是以 1 为间隔，需要告诉 gradient 真实的间距：

hx = x(2) - x(1); hy = y(2) - y(1); [Ex, Ey] = gradient(-V, hx, hy);

这一步不写，画出来的箭头方向会乱，大小也失真。我最初就栽在这里，箭头指向看着没问题，但长度完全不对。

2.3 完整的电偶极子模拟脚本

把上面的片段拼起来，再配上可视化部分，一个能直接跑的完整脚本如下：

% 电偶极子场模拟完整脚本 clear; clc; % 常数设置 k = 8.99e9; % 库仑常数 q = 1e-9; % 电荷量 1nC d = 0.02; % 距离 2cm % 空间网格 x = linspace(-0.06, 0.06, 301); y = linspace(-0.06, 0.06, 301); [X, Y] = meshgrid(x, y); % 计算电势（含奇点修正） R1 = sqrt((X - d/2).^2 + Y.^2 + 0.002^2); R2 = sqrt((X + d/2).^2 + Y.^2 + 0.002^2); V = k * q ./ R1 - k * q ./ R2; % 计算电场 hx = x(2) - x(1); hy = y(2) - y(1); [Ex, Ey] = gradient(-V, hx, hy); % 可视化 figure('Color', 'w'); % 画等势线 contour(X, Y, V, 30, 'LineWidth', 0.8); hold on; % 画电场矢量 quiver(X, Y, Ex, Ey, 1.5, 'Color', [0.8 0.1 0.1]); % 标出电荷位置 plot(d/2, 0, 'ro', 'MarkerSize', 10, 'MarkerFaceColor', 'r'); plot(-d/2, 0, 'bo', 'MarkerSize', 10, 'MarkerFaceColor', 'b'); axis equal; xlabel('x (m)'); ylabel('y (m)'); title('电偶极子电场与等势线分布'); hold off;

这段代码直接粘贴就能跑通。上面的quiver 的缩放因子 1.5是视觉调参后的结果。太小箭头看不清方向，太大会相互遮挡。如果改用 1 或 2，视觉效果都会差一些，这个参数值得多试几次找到最适合自己图的。

3. 可视化策略与绘图细节优化

3.1 不同图形之间的取舍

代码写出来是第一步，图好不好看直接影响报告和论文的质量。我试过几种不同的可视化方法，各有特点：

我最推荐的是等值线 + 矢量箭头叠加图，这是几乎所有电磁学教材里电偶极子图的呈现方式，信息量大且对比清晰。

3.2 解决箭头过密或过稀的问题

很多时候网格是 301 x 301，如果直接把所有点的电场画出来，图上一片红，什么都看不清。解决思路是“抽样画箭头，全量算场”：

% 每 15 个点取一个箭头 step = 15; xs = X(1:step:end, 1:step:end); ys = Y(1:step:end, 1:step:end); Exs = Ex(1:step:end, 1:step:end); Eys = Ey(1:step:end, 1:step:end); quiver(xs, ys, Exs, Eys, 1.2, 'Color', [0.7 0.2 0.2], 'LineWidth', 0.6);

这一步看起来简单，但实际效果差别非常大。没有抽样的图是“红色马赛克”，抽样后才是干净清晰的箭头流场。建议 step 值在 10~20 之间试，找到和你图形大小匹配的密度。

3.3 隐藏电偶极子附近的极端电场

电偶极子附近的电场强度远高于远处，如果用颜色标度统一显示，远处的场几乎看不见。我处理办法是截断色标范围，具体做法是把 V 的数值在两倍标准差附近截断，或者用 caxis 控制显示范围（新版本用 clim ）：

clim([-500 500]);

这样设置后，远离电荷区域也能看到清晰的电势梯度变化，不会只有两个亮点。这个细节在做伪彩图时尤其关键。

4. 常见报错、异常排查与避坑经验

4.1 为什么画出来的全是 NaN 或 Inf

最常见的原因是 R1 或 R2 出现了 0。网格点恰好落在电荷位置时，距离为 0，除以 0 就是 Inf。前面说的加偏移量 0.002 就是最直接的解决办法。还有一种情况是网格中心出现在坐标原点，而电荷放在正负 d/2，这时候原点处电势的确是 0，但个别线画到电荷附近依然会有异常，加偏移量同样能解决。

4.2 箭头方向看起来是对的但大小不对

这是绝大多数人忽略的问题。gradient(V) 如果不传 hx 和 hy，默认假设数据点间距为 1。但我们的实际坐标单位是米，间距是 0.0004 米左右，导数算出来会被放大 2500 倍。箭头方向不受影响，但长度全错。加了 hx 和 hy 参数后，量级才正确。我建议无论何时都写上这两个参数，不要嫌麻烦。

4.3 为什么等值线在电荷之间“打架”

这其实不是错误，而是物理本质。在电偶极子中轴线上，正负电荷之间的电势梯度很大，等值线会很密集地挤在一起。如果嫌图不美观，可以在 contour 里限制等级数：

contour(X, Y, V, [-800:50:800])

用明确指定的等值线值代替“自动分成 30 条”，可以避免过密堆积，也能针对自己关心的电势范围做定制。

4.4 3D 视角看电势表面的经验

如果要把 V 画成三维曲面，我推荐用 surf 加一个简单的光照调整：

surf(X, Y, V); shading interp; colormap(jet); view([-30, 30]);

这样能从侧面看到“两个尖峰一个低谷”的形状，正电荷处是高峰，负电荷处是深谷，直观地呈现了势能景观。如果觉得尖峰太高影响观察，也可以用 zlim 截断。

5. 从静电场到动态演示的进阶扩展

5.1 让电荷动起来，看场如何变化

静态图看够了，可以模拟旋转偶极子的场变化。这个进阶玩法我强烈建议试试，用for 循环改变电荷角度，不断更新电势并重绘等值线，就能得到“场随电荷移动而变化”的动画效果：

% 定义旋转角度序列 theta = 0:0.05:2*pi; R = 0.02; % 旋转半径 figure('Color', 'w'); for t = 1:length(theta) % 计算当前电荷位置 xp = R * cos(theta(t)); yp = R * sin(theta(t)); xn = -xp; yn = -yp; % 重新计算距离和电势 RP = sqrt((X - xp).^2 + (Y - yp).^2 + 0.002^2); RN = sqrt((X - xn).^2 + (Y - yn).^2 + 0.002^2); Vt = k * q ./ RP - k * q ./ RN; % 重绘图 clf; contourf(X, Y, Vt, 20, 'LineStyle', 'none'); hold on; plot(xp, yp, 'ro', 'MarkerSize', 8, 'MarkerFaceColor', 'r'); plot(xn, yn, 'bo', 'MarkerSize', 8, 'MarkerFaceColor', 'b'); axis equal; title(sprintf('旋转电偶极子角度: %.1f°', rad2deg(theta(t)))); drawnow; end

这段代码跑起来后你会看到等势线随着电荷旋转而变化，比静态图有意思得多。这也方便在场模拟演示课上当场展示“电场是瞬时的，电荷一动场就跟随着变”的物理图像。

5.2 从二维扩展到三维视角

再进一步，可以把二维网格扩展成三维，计算空间的电势并画切片：

[x3, y3, z3] = meshgrid(linspace(-0.05, 0.05, 51)); R1_3d = sqrt((x3 - d/2).^2 + y3.^2 + z3.^2 + 0.002^2); R2_3d = sqrt((x3 + d/2).^2 + y3.^2 + z3.^2 + 0.002^2); V3 = k * q ./ R1_3d - k * q ./ R2_3d; % 用 slice 画三个切面 slice(x3, y3, z3, V3, 0, 0, 0);

三维版本的计算量明显增加，网格点数建议控制在 51 以下，否则内存会吃紧。切片图适合展现场的对称性，能看出电偶极子场关于 z 轴旋转对称的结构。

5.3 隐藏的高阶玩法：矢势和辐射场

如果想再进一步，可以计算偶极子的矢势，甚至加上时间因子变成辐射场。这里可以做一个简化的偶极子辐射场模拟，用sin(theta) / r来描述远场方向图，画出偶极子的辐射花瓣图。不过这个属于天线工程的内容了，和纯静电场模拟的数值方法差别较大，感兴趣的话可以自己往前研究。

6. 总结一点实操后的心得

上面这些经验总结下来，其实最核心的思想就是“先算标量，再推矢量”。电偶极子的场模拟，不管你是本科课程设计，还是研究生科研项目的前期验证，掌握这个思路后可以迁移到任意电荷系统，比如四极子、线电荷，甚至平行板电容器。MATLAB 的优势在于快速验证，10 分钟写出代码，5 分钟调出图，这是手算画图没法比的效率。

最后分享一个小技巧：调试时先用 51 x 51 的粗网格把逻辑跑通，再加密到 301 x 301 渲染最终图。这样每次迭代只需几秒钟，不用等待漫长的矩阵运算。我早期就是一把梭直接上 501 x 501，改一个参数等半分钟，一天下来效率低得让人想放弃。粗网格调试、细网格出图，这个习惯帮我至少省了一个星期的时间。