简介：这份资源是一份基于Matlab模拟线电荷电场分布的仿真实验报告，面向学习电磁场与电磁波课程、需要完成仿真实验或课程设计的高校学生。报告围绕库仑定律、电场强度与势函数展开，通过设置线电荷分布、构建二维坐标网格、循环累加各点电势，并调用surf、contour与quiver函数绘制电势三维图、等位线图和电场矢量图，完整呈现从原理推导到代码实现再到结果分析的实验流程。压缩包内为1个docx文档，约227KB，包含实验目的、原理说明、程序流程、完整Matlab代码及运行结果图，并附有实验结论与感悟，便于读者对照复现和修改参数。目前已有3041人学习，适合作为电磁场仿真入门参考与实验报告撰写模板。
1. 从一张“糊成一片”的电位图说起：线电荷电场仿真到底要交什么
第一次把这段代码跑通的人，十有八九会盯着surf出来的那张图发愣：中间一道深沟，两边翘起，等位线挤成一团，矢量箭头密密麻麻像刺猬。图是出来了，可报告里要写的“电场分布规律”却无从下笔。问题不在 Matlab 用得熟不熟，而在于线电荷电场仿真这件事，本质是把连续分布的电荷做离散化近似，再用数值梯度反推场强，中间每一步都有可调参数和隐含假设。这份实验报告要交付的，是一套能自圆其说的仿真流程：电荷怎么离散、网格怎么划、电势怎么叠加、场强怎么由电势求得、图怎么画才读得出物理意义。适合正在做电磁场课程设计、Matlab 仿真实验报告，或者想拿一个最小可复现案例入门数值场计算的人。下面按“需求—设计—实现—避坑—验收”的顺序，把这份报告拆成可以直接照着做的方案。
2. 需求与物理模型：把连续线电荷拆成有限个点电荷
2.1 实验要回答的三个问题
这份仿真实验报告表面上是“画两张图”，实际要回答三个层次的问题。第一层是数值层面：给定电荷量、线电荷长度和网格范围，Matlab 能不能稳定算出每个网格点的电势，不出现NaN、不出现维度报错。第二层是物理层面：等位线是否围绕线电荷呈对称分布，场强矢量是否垂直于等位线、由高电位指向低电位。第三层是方法层面：离散电荷个数nr取多少合适，网格步长0.3是粗是细，这些参数变化对结果有什么影响。
很多报告翻车就翻在只做了第一层，图一贴、结论一写“电场呈对称分布”就交差。合格的报告要把三层都覆盖到，尤其是第三层，它才是“仿真”区别于“画图”的地方。通常的做法是：固定其他参数，单独改变nr或网格步长，各出一组图做对比，用图说话。
2.2 库仑定律到离散叠加的推导链
物理模型这条链要写清楚，否则代码里的每一行都成了黑匣子。真空中点电荷q在距离r处产生的电势为：
U = q / (4 * pi * e0 * r)其中e0是真空电容率。线电荷不是单个点电荷，而是电荷沿一条线段连续分布。严格求解要对线电荷做线积分，但数值仿真里更常用的做法是“以点代线”：把长度为L的线电荷等分成nr个小段，每段近似为一个点电荷，总电势就是所有点电荷电势的标量叠加。这就是“利用微积分原理模拟线电荷”的落地方式——积分变求和，连续变离散。
这里有个容易被忽略的细节：离散后每个小段代表的电荷量应该是总电荷除以段数，而不是每段都用同一个q。原程序里每段都取q=1.6e-19，相当于总电荷随nr线性增长，这在物理上对应“电荷线密度固定、线电荷随段数变长”的情形。报告里必须把这一点讲明白，否则读者会误以为nr只是精度参数。我一般会在报告里明确写一句：本仿真采用等电荷量离散，nr增大时总线电荷同步增大，因此对比图时只看分布形态、不看绝对幅值。
2.3 为什么选 Matlab 而不是别的工具
选型理由要写，但不能写成软件广告。对这个实验而言，Matlab 的优势集中在三点：矩阵化运算让meshgrid生成的网格可以直接参与逐元素计算，省掉双重循环；surf、contour、quiver三个绘图函数刚好对应电势三维图、等位线图、场强矢量图三种表达；gradient函数能直接从电势矩阵数值求梯度，免去手写差分。换成 Python 的 NumPy 加 Matplotlib 同样能做，但对课程实验来说，Matlab 的“一个脚本出三张图”更贴合实验报告的交付节奏。选型不是比谁强，是比谁在这个场景下改动成本低。
3. 总体设计：网格、离散、叠加、求梯度四步走
3.1 程序流程的四个阶段
整套仿真可以拆成四个阶段，每个阶段对应一段可独立调试的代码。第一阶段建网格：用meshgrid在[-40,40]×[-40,40]范围内按0.3步长生成坐标点，得到X、Y两个矩阵。第二阶段布电荷：用linspace在y=0的直线上均匀放置nr+1个点电荷，横坐标从-25附近展开。第三阶段算电势：循环每个点电荷，算出它到所有网格点的距离矩阵，套电势公式，累加到总电势U。第四阶段求场强并绘图：对-U求梯度得到Ex、Ey，归一化后画矢量图，同时画等位线和三维曲面。
这个流程的顺序不能乱。电势必须先算完再求梯度，因为gradient作用在整个矩阵上；矢量图必须在等位线之后画，否则hold on的叠加顺序会让图例混乱。报告里的“实验步骤”一节，就按这四个阶段写，每阶段配一句“输入是什么、输出是什么”，比笼统写“设置参数、计算、绘图”要清楚得多。
3.2 关键参数表与取值依据
参数不是随手填的，每个都要能说出理由。下面这张表是我做这类实验时习惯整理的参数清单，报告里可以直接引用：
| 参数 | 取值 | 含义 | 取值依据 |
|---|---|---|---|
nr | 50 | 线电荷离散段数 | 段数过少分布不均，过多循环耗时且图线过密 |
e0 | 1e-9/(36*pi) | 真空电容率 | 与1/(4*pi*e0)=9e9的量级对应 |
q | 1.6e-19 | 单段电荷量 | 取基本电荷量级，便于观察 |
| 网格范围 | -40:0.3:40 | 计算区域 | 覆盖电荷分布区并留出足够边界 |
| 等位线条数 | 30 | contour层数 | 太少看不出层次，太多线条糊在一起 |
| 矢量缩放 | 0.5 | quiver缩放因子 | 控制箭头长度，避免重叠 |
网格步长0.3是个折中：步长1时等位线出现明显折角，步长0.1时矩阵规模到801×801，循环 51 次后内存和耗时都上去了。报告里可以补一句“步长减半后图形平滑度提升有限，但计算量翻四倍”，这就是参数讨论的价值。
3.3 电势叠加的矩阵化写法
原程序用for循环逐段累加，写法直观但有个隐患：U的初始化。代码里写U=[0]，第一次循环时U是标量0，加上Uk（矩阵）后自动扩展成矩阵，Matlab 能跑通，但这是依赖隐式扩展的写法，换个语言或换个版本就可能报维度错误。更稳妥的写法是先U=zeros(size(X))，让累加从一开始就在同维度矩阵上进行。这一点在报告的“实验结论”里被点到了——“矩阵求和维度必须一致，累加时初始矩阵赋值为 0”——但代码里没落实，属于典型的“结论写得对、代码没改”。我在方案里会直接把初始化改成zeros，并在注释里说明原因。
4. 核心实现：从电势矩阵到三张图的完整代码
4.1 电势计算主循环
下面这段是整理后的主计算部分，相比原始代码做了三处改动：初始化改为同维零矩阵、距离计算加极小量防止除零、电荷位置用变量统一管理。
clear; clc; close all; nr = 50; % 线电荷离散段数 e0 = 1e-9/(36*pi); % 真空电容率 q = 1.6e-19; % 每段电荷量 k_e = 1/(4*pi*e0); % 库仑常数，约 9e9 [X, Y] = meshgrid(-40:0.3:40, -40:0.3:40); % 坐标网点 U = zeros(size(X)); % 同维初始化，避免隐式扩展 x0 = linspace(-25, 25, nr+1); % 线电荷横坐标，关于 y 轴对称分布 y0 = 0; % 线电荷位于 y=0 直线上 for k = 1:numel(x0) Rk = sqrt((X - x0(k)).^2 + (Y - y0).^2); Rk = Rk + 1e-12; % 防止网格点与电荷重合导致除零 Uk = k_e * q ./ Rk; % 单段点电荷电势 U = U + Uk; % 标量叠加 end这段代码的关键在x0的取法。原程序用k-25让电荷从-24排到26，中心偏了 1 个单位，等位线看起来“歪了一点”却说不清原因。改成linspace(-25,25,nr+1)后电荷严格关于原点对称，等位线左右对称，报告里“对称分布”的结论才站得住。1e-12这个小量是血泪经验：网格步长0.3时，电荷点恰好落在网格点上的概率不低，一旦Rk=0，Uk变Inf，整张图报废，加个极小量就能绕过去。
4.2 由电势求场强：gradient 的方向问题
电场强度是电势的负梯度，E = -∇U。Matlab 的gradient返回的是数值差分，默认对矩阵第二维（列方向，即x）和第一维（行方向，即y）分别求偏导，所以调用顺序是[Ex, Ey] = gradient(-U)，注意先取负再求梯度，或者求完梯度再取负，两者等价但别漏。
[Ex, Ey] = gradient(-U); % 负梯度得到场强分量 AE = sqrt(Ex.^2 + Ey.^2); % 场强幅值 Ex = Ex ./ (AE + 1e-12); % 归一化，只保留方向 Ey = Ey ./ (AE + 1e-12);归一化这一步是为了画矢量图。不归一化的话，靠近电荷处场强幅值极大，箭头长度差异悬殊，quiver自动缩放后远处的箭头几乎看不见。归一化后所有箭头等长，只表达方向，配合等位线看“箭头是否垂直穿过等位线”，物理意义更直观。报告里要说明：归一化牺牲了幅值信息，幅值信息由三维曲面图surf承担，两张图分工不同。
4.3 三张图的绘制与读图要点
绘图部分分三张图，每张图对应一个物理结论。第一张surf(X,Y,U)看电势的空间起伏，读图要点是“电荷所在直线附近电势最高、向两侧单调衰减”。第二张contour加电荷位置标记，看等位线形状，读图要点是“等位线在电荷附近密集、远处稀疏，且关于y=0对称”。第三张quiver叠加在等位线上，读图要点是“箭头方向由高电位指向低电位，且与等位线正交”。
figure(1); surf(X, Y, U); shading interp; % 平滑着色，避免网格线干扰 xlabel('X'); ylabel('Y'); zlabel('U'); title('线电荷电势三维分布'); figure(2); Vmin = min(U(:)); Vmax = max(U(:)); CV = linspace(Vmin, Vmax, 30); % 30 条等位线 contour(X, Y, U, CV, 'LineWidth', 1); hold on; plot(x0, y0*ones(size(x0)), 'r.', 'MarkerSize', 8); % 标出电荷位置 quiver(X, Y, Ex, Ey, 0.5, 'g'); % 矢量场叠加 hold off; xlabel('X'); ylabel('Y'); title('线电荷等位线与场强方向');shading interp是容易被省掉的一行，不加的话surf默认带网格线，三维图看起来像铁丝网，报告里放这种图显得粗糙。Vmin、Vmax用U(:)取全局极值，比原程序的min(min(U))更通用，矩阵维度变了也不用改。等位线取 30 条是经验值，少于 15 条层次不够，多于 50 条线条挤在一起反而看不清疏密变化。
5. 避坑与排查：这类仿真最容易翻车的五个点
5.1 现象：图能出但等位线明显偏心
原因通常是电荷位置没有关于原点对称。原程序k-25的写法让电荷中心落在x=1附近，等位线整体右移，肉眼能看出但说不出所以然。解决办法是用linspace显式指定对称区间，并在报告里写明电荷分布范围。排查时可以直接plot(x0, y0, 'r*')单独画一张电荷位置图，一眼就能看出偏没偏。
5.2 现象：运行报错“矩阵维度必须一致”
原因有两类。一类是U初始化成标量[0]，某些 Matlab 版本或后续改动后不再自动扩展；另一类是循环里某个变量被意外转置。解决办法是统一用zeros(size(X))初始化，并在循环内所有涉及网格的运算都用.*、.^、./逐元素运算符。排查技巧：在循环里临时加一句size(U)打印，看第一次累加前后维度是否一致。
5.3 现象：图上出现NaN或整片空白
原因是网格点与电荷位置重合，Rk=0导致除零。步长0.3和电荷坐标-25:1:25都是 0.3 的倍数关系时，重合概率不低。解决办法是在Rk上加极小量1e-12，或者把电荷坐标整体偏移半个步长。这个坑不显眼，但一旦踩中，surf会直接画出空洞，contour可能报错，属于典型的“玄学”故障，加个极小量就根治。
5.4 现象：矢量箭头方向看起来“反了”
原因是gradient求的是梯度，电场是负梯度，漏掉负号就会让箭头由低电位指向高电位。排查方法很简单：在等位线图上找一条线，看箭头是穿向更高电位一侧还是更低电位一侧。另一个容易混淆的点是gradient的返回顺序，[Ex,Ey]对应的是列方向导数和行方向导数，如果网格的X、Y定义顺序反了，箭头会整体旋转 90 度。解决办法是固定meshgrid(x,y)的参数顺序，并在注释里写清楚。
5.5 现象：nr改大后图反而更糊
原因是电荷段数增加后，总线电荷同步增大，近场电势幅值飙升，contour的等位线间隔被近场主导，远处细节被压缩。这不是 bug，是“等电荷量离散”模型的固有特性。解决办法有两个方向：要么改成“总线电荷固定、每段电荷量随nr反比缩小”，要么在报告里明确说明本模型是等电荷量离散、对比时只看形态。我一般选后者，因为改动小、结论也诚实。报告里把这一点写进“讨论”，反而是加分项。
6. 验收与答辩：怎么证明这套仿真“算对了”
验收这类仿真实验，不能只看“图好不好看”，要有可核对的判据。我通常准备三组对照，答辩时按顺序摆出来。第一组是解析解对照：在y轴上取几个点，用无限长线电荷的解析公式U = (λ/(2πε₀))·ln(r₀/r)估算电势，和仿真值比量级和衰减趋势，允许离散化带来的偏差，但趋势必须一致。第二组是参数敏感性对照：固定其他参数，把nr从 20 改到 80，把网格步长从 0.5 改到 0.2，各出一张等位线图，说明“段数增加后近场更平滑、步长减小后等位线折角消失”，证明自己知道每个参数在干什么。第三组是物理自洽性检查：场强矢量是否与等位线正交、是否由高电位指向低电位、电荷两侧是否对称，这三条用肉眼就能验，但必须主动指出来。
答辩时最容易被问的问题是“你这个线电荷为什么不是理想线电荷”。标准回答是：离散化本身就是近似，nr有限导致电荷分布有间隙，近场会出现周期性起伏；改进方向是增大nr或改用连续积分解析式做基准。另一个高频问题是“电势参考点选在哪”，回答：点电荷电势公式默认无穷远处为零，离散叠加后仍以无穷远为零参考，因此绝对幅值有物理意义，但近场因离散化偏大。把这两个问题答清楚，基本就稳了。
还有个小技巧：报告里放图时，每张图下面配一句“读图结论”，而不是只写图号。比如等位线图下面写“等位线在电荷附近密集，说明场强在近场更强；等位线关于 y=0 对称，验证了电荷分布的对称性”。评委翻报告时先看图注，图注写得好，印象分直接上去。
从那以后我每次做数值仿真实验，都强制自己先跑一遍“参数敏感性对照”再写结论——图能出来只说明代码没报错，参数一动结论还成立，才说明模型站得住。这套流程放在这份线电荷电场仿真报告里，从需求、模型、代码到验收是一条完整的线，照着改参数、换电荷分布就能复用到其他静电场仿真上。希望帮到你。
本文还有配套的精品资源，点击获取
