MATLAB三相短路分析全攻略：从Simulink建模到短路电流计算

简介：围绕MATLAB/Simulink的电力系统三相短路故障仿真分析资料，面向电力工程专业学生、课程设计者及电气工程师，可系统解决短路计算流程不清晰、保护校验缺乏实例、Simulink建模无从下手等痛点。压缩包共12个文件，核心为MATLAB脚本（.m）与课设报告（PDF），辅以电路原始数据、输出结果及输入输出格式等文本，并配有6张节点阻抗矩阵、节点导纳矩阵过程图；整体仅1.09MB，轻量却覆盖完整分析链路。已有2576人学习下载。内容从网络结构定义、元件参数设定，到Simulink故障建模、仿真运行与结果解读层层推进；课设报告详细说明包括模型建立、故障时间与类型设定、仿真结果分析在内的完整流程，配合可直接运行的脚本，能快速复现三相短路瞬间的电流电压变化，并用于校验保护设备的动作逻辑。节点阻抗/导纳矩阵过程图直观呈现原始网络以及发电机、负载等效接入前后的矩阵演变，txt数据便于对照验算。无论课程设计、毕业设计，还是自学短路电流计算与继电保护整定，这套资料都能提供从理论到仿真的闭环参考。

1. 从一道仿真题说起：三相短路为什么值得花时间做

做电力系统方向的人，几乎都会碰到三相短路分析。不管是课程设计、毕业设计，还是刚进电力行业做保护整定计算的新人，这个课题都是绕不开的入门门槛。原因很简单：三相短路是电力系统中最严重的故障类型，短路电流最大，对设备、电网稳定性的冲击最剧烈，整个系统从规划设计到保护配置都得围着它转。

但说实话，很多人在学校做这个课题时，要么是照着课本把公式推导一遍就收工，要么是拿着Simulink里现成的模块拖几个框图，看个波形就算完成任务。这两种做法都没错，但离“真正理解三相短路分析”还差得远。我自己的体会是，用MATLAB做三相短路分析，意义不在于“跑出一个短路电流数值”，而在于通过这个过程把暂态过程的物理图像、各电气量的变化规律、不同分析方法的适用场景这几件事打通。这篇就把我实际做过的全套流程和踩过的坑整理出来，给正在做类似课题的同学或同行一个参考。

这套分析能帮你解决什么问题？简单说三个：第一，计算短路瞬间各支路的短路电流周期分量和非周期分量，校核设备能否承受；第二，分析短路后系统电压跌落情况，评估对负荷和稳定性的影响；第三，验证保护装置的整定值——没有准确的短路计算，保护就无从谈起。适合的读者包括电气工程专业高年级本科生、电力系统方向研究生，以及刚入门做电力仿真的一线工程师。

2. 三相短路分析到底在算什么——核心物理量拆解

2.1 短路电流的两个分量：周期分量与非周期分量

做仿真之前，建议先把物理概念彻底搞清楚。不然波形出来了一团乱麻，你也分不清哪些是正常现象，哪些是模型搭错了。

三相短路发生时，短路电流由两部分叠加组成。一部分是周期分量（也叫交流分量或稳态短路电流），它由电源电动势和系统总阻抗决定，理论上不衰减，频率就是工频50Hz。另一部分是非周期分量（直流分量），它是电路电感储能导致的自由分量，跟短路时刻的电压相位角有关，按时间常数衰减到零。

非周期分量这东西特别有意思，它只出现在短路后的最初几个周波内，但恰恰是它决定了短路冲击电流的峰值。峰值一般在短路发生后约半个周波（10ms左右）出现，可能达到周期分量幅值的1.8到1.9倍。这个冲击电流就是用来校验断路器开断能力和母线、支柱绝缘子动稳定度的关键指标。

我自己刚开始搭模型时犯过一个典型错误：只盯着周期分量看，忽略非周期分量，结果仿出来的波形和理论手册上的曲线对不上。后来才意识到，非周期分量虽然衰减快，但在前几个周波内对波形形状的影响非常大。如果在MATLAB里观察短路电流波形，你一定会看到开始那一段电流明显“偏”在时间轴一侧，这就是直流分量的作用。

2.2 标幺值系统里的短路计算

三相短路分析里，标幺值系统几乎是必用的工具。原因很实际：电力系统有不同电压等级，各自有功率基准和电压基准，如果都用有名值计算，变压器变比、电抗折算会把人绕晕。用标幺值把所有设备阻抗统一到同一套基准下，计算起来清爽得多。

具体做法是选定基准容量SB（通常取100MVA或系统总容量）和基准电压UB（各电压等级的平均额定电压，比如10.5kV、37kV、115kV、230kV），然后各个设备的电抗标幺值按公式折算：

• 线路电抗标幺值：x* = x（Ω/km）× L（km）× SB / UB²
• 变压器短路电抗标幺值：工程上直接取短路电压百分数UK%换算，x* = UK% / 100 × SB / SN
• 发电机次暂态电抗标幺值：需要按容量折算到统一基准

这里有个常见坑：发电机、变压器铭牌上的电抗都是按自身额定容量给出的标幺值，必须折算到你的统一基准下才能直接相加。很多新手第一步就挂在这里，算出来的短路电流差了好几倍。折算公式很简单：x*(新基准) = x*(额定基准) × SB / SN。

2.3 短路电流计算的两类方法：无穷大电源法与戴维南等值

实际工程里做三相短路计算，根据分析目的不同会选不同的假设条件。两种最常用的方法必须在做MATLAB仿真之前就想清楚用哪个。

第一种是无穷大电源法。假设短路点上游电源容量无穷大、电压恒定不变，短路电流周期分量也就不衰减。这种假设适合估算短路电流的最大可能值，用于设备选型校验。计算简洁，静态分析、保护整定计算里用得最多。

第二种是戴维南等值法。把短路点之前的系统等值为一个电压源串一个阻抗，适用于网络结构复杂、需要精确计算故障后电压分布的场景。在MATLAB里做潮流计算后，用短路点处的节点电压和系统阻抗就可以精确求出短路电流。

如果做的是暂态仿真（比如观察短路全过程的电压电流动态变化），那就要基于微分方程数值积分来做，Simulink正是干这个的。所以做三相短路分析前，先想清楚你要的是“稳态数值”还是“暂态波形”，这将直接决定你的MATLAB实现路线。

3. MATLAB实现路线：Simulink建模与m文件计算怎么选

3.1 两条路线的适用场景对比

用MATLAB做三相短路分析，主流有两条路：一条是Simulink里搭电力系统模型，用Simscape Electrical（早期叫SimPowerSystems）模块库里的三相电源、变压器、线路、断路器等元件搭电网；另一条是纯m文件编程，把系统参数写成矩阵和公式，用数值方法求解。

这两条路我都有过实际使用经历，选哪条主要看你的输出目标是什么。

Simulink建模的好处是直观，模块拖拽就能搭出系统，而且可以直接测量短路电流波形、母线电压动态过程，特别适合做暂态分析、观察断路器动作后的恢复过程。缺点是模型搭好后调试麻烦，有时候仿真步长设置不对，波形会出现数值振荡，很难判断是真实物理现象还是数值问题。

m文件编程的好处是全程可控，每一步都能知道自己算什么、为什么这么算，代码本身就是“可追溯的计算说明书”。缺点是看不到直观的电路结构，纯代码处理复杂网络时比较繁琐。

如果让我给建议：做课程展示、毕业设计仿真演示，优先Simulink，画面感强、答辩时好讲；做深入研究、批量计算不同短路点或多方案对比，强烈建议m文件封装成函数，效率高得多。当然，两者结合起来也很顺手——用m文件算好短路电流的理论值，再放到Simulink里做动态仿真验证，相互印证，结果可信度大幅提高。

3.2 手写代码 vs 工具箱：别小看基础实现

很多人一上手就去找Simulink的现成例子，或者搜“三相短路仿真模型下载”，这当然可以，但我还是建议至少手动写一遍核心计算逻辑，哪怕只是最简模型。

原因是我自己踩过的一次教训：用现成模型下载后跑出来的短路电流用起来没问题，但被老师一问“你这个短路电流用的是什么假设？次暂态还是稳态？冲击系数取了多少？”直接愣住了。这些问题没搞清楚，模型跑得再好看也是虚的。

手动实现一遍以后，你对内部机理的理解会完全不一样。比如手算求短路电流时，你会发现需要先做潮流计算得到短路点正常运行电压，再对短路点以外的网络做戴维南等值求等值阻抗。这个过程会让你真正明白“短路电流的大小取决于故障前电压和故障点等值阻抗”这句话意味着什么。

4. 完整实操：一个10kV配电系统三相短路分析全流程

4.1 系统参数与研究目标设定

下面用一个我实际做过的算例来说明完整流程。算例是一个简化的10kV配电系统，模型如下：

• 110kV无穷大电源接入，110/10kV降压变压器，容量SN = 31.5MVA，短路电压UK% = 10.5%，变比110±8×1.25%/10.5kV
• 10kV出线为架空线路，长度L = 8km，正序电抗x1 = 0.35Ω/km，正序电阻r1 = 0.15Ω/km
• 短路点设在10kV线路末端
• 分析目标：计算短路点的三相短路电流周期分量起始值I''、冲击电流峰值iimp、短路容量Sk，并观察短路电流的暂态波形

这个算例覆盖了从变压器到线路的完整阻抗链，既可以手动校验，也能看MATLAB仿真结果是否吻合。

4.2 手算理论值：先算准，再仿真

按标幺值法手算。取基准容量SB = 100MVA，基准电压UB = 10.5kV（对应10kV电压等级的平均额定电压）。

变压器电抗标幺值：

xT* = UK% / 100 × SB / SN = 0.105 × 100 / 31.5 = 0.3333

线路长度8km，归算到10.5kV基准下的电抗标幺值：

xL* = x1 × L × SB / UB² = 0.35 × 8 × 100 / 10.5² = 2.54

线路电阻标幺值：

rL* = r1 × L × SB / UB² = 0.15 × 8 × 100 / 10.5² = 1.09

系统总等值阻抗标幺值（忽略电源内阻抗时）：

z* = sqrt((xT* + xL*)² + rL*²) = sqrt((0.3333+2.54)² + 1.09²) ≈ 3.08

短路点电压近似取平均额定电压，U* = 1.0，则短路电流周期分量起始值标幺值：

I* = U* / z* = 1.0 / 3.08 ≈ 0.325

折算为有名值：

I'' = I* × SB / (√3 × UB) = 0.325 × 100 / (1.732 × 10.5) kA ≈ 1.79kA

冲击电流按冲击系数Kimp = 1.8计算：

iimp = Kimp × √2 × I'' = 1.8 × 1.414 × 1.79 ≈ 4.56kA

短路容量：

Sk = √3 × UB × I'' = 1.732 × 10.5 × 1.79 ≈ 32.6MVA

注意：实际工程中计算线路阻抗折算到低压侧时，需要留意变压器变比到底取哪个电压值，10.5kV是变压器分接头的额定电压。如果取10kV，线路电抗标幺值会略有差异，但工程上误差在可接受范围内。

手算值出来后，后面仿真才有“对答案”的依据。

4.3 Simulink建模的核心步骤与参数设置

Simulink建模流程并不复杂，但有几个关键参数容易设置错，这里把我的操作过程拆开写一下。

第一步，先建一个空白模型，从Simscape Electrical库中拖出以下元件：Three-Phase Source（三相电压源）、Three-Phase Transformer（三相变压器）、Three-Phase Series RLC Branch（线路阻抗，也可以拆成RL分离的两个元件）、Three-Phase Fault（三相短路故障模块）、Three-Phase V-I Measurement（电压电流测量模块），以及一个Power GUI模块。

第二步，设置电源参数。110kV电压源采用无穷大电源等效时，内阻抗设一个极小值（如0.001Ω），幅值按110kV有效值设置，频率50Hz。这里有一个需要注意的点：如果电压源被设置为“Ideal 3-Phase Source”（理想三相电源）而没有任何内阻抗，短路电流中的非周期分量衰减时间常数可能异常，数值上可能出现不收敛的情况。建议给电源串个极小阻抗，既不影响短路计算结果，又能避免数值问题。

第三步，变压器参数设置。额定功率31.5MVA，频率50Hz，一次绕组额定电压110kV，二次绕组额定电压10.5kV，短路阻抗百分比UK%填10.5%，接线组别选Delta-Wye（Δ-Yg），这是国内10kV系统最典型的接线方式。

第四步，线路用Series RLC Branch，电阻0.15×8=1.2Ω，电感按电抗算：XL = 0.35×8=2.8Ω，L = XL/(2π×50) ≈ 0.0089H。注意这里用的是有名值（因为Simulink中可以选择用有名值或标幺值），和手算的标幺值计算过程正好能互相对照。

第五步，Three-Phase Fault模块设置。把故障类型选为Three-Phase（三相短路），切换时间为[0.1 0.3]秒，即0.1秒时发生短路，0.3秒时故障切除。这个时间设置是为了能同时看到短路前稳态、短路中暂态、短路切除后恢复三个阶段的波形。

第六步，给测量模块加上示波器，重点观察三相短路电流、短路点母线三相电压和线路有功功率。

仿真参数建议这样设置：仿真时间0.5秒（覆盖短路前后完整的动态过程），求解器选变步长ode23tb，相对容差1e-4。这里特别注意——如果选固定步长或ode45默认参数，非周期分量的衰减过程可能失真，波形尾部会出现毛刺。

4.4 m文件实现的代码逻辑与结果对比

m文件实现的核心思路是直接计算短路电流的理论值，并画出三相短路的电流波形。这里写一个简化版的核心逻辑，方便理解。

% 三相短路分析 - 基于戴维南等值法 % 系统参数（标幺值，SB=100MVA, UB=10.5kV） UB = 10.5e3; % 基准电压 SB = 100e6; % 基准容量 xT = 0.3333; % 变压器电抗标幺值 xL = 2.54; % 线路电抗标幺值 rL = 1.09; % 线路电阻标幺值 % 短路点等值阻抗 Zeq = sqrt((xT+xL)^2 + rL^2); I_sc_pu = 1.0 / Zeq; % 短路电流标幺值 % 折算有名值 I_base = SB / (sqrt(3)*UB); I_sc = I_sc_pu * I_base; % 周期分量有效值，单位A % 冲击电流 Kimp = 1.8; i_imp = Kimp * sqrt(2) * I_sc; % 输出结果 fprintf('短路电流周期分量起始值: %.2f kA\n', I_sc/1000); fprintf('冲击电流峰值: %.2f kA\n', i_imp/1000); fprintf('短路容量: %.2f MVA\n', sqrt(3)*UB*I_sc/1e6);

这段代码的逻辑很简单，但它是整个短路计算的核心骨架。我在实际项目中会把它封装成函数，输入是系统阻抗参数，输出是短路电流、冲击电流、短路容量，然后批量遍历不同短路点、不同运行方式下的计算结果，效率非常高。

用Simulink仿真跑完后的结果和手算值高度一致。周期分量有效值约1.79kA，冲击电流约4.5kA，短路容量约32.6MVA，和4.2节的手算理论值完全对得上。这说明两种方法交叉验证成功，模型搭得没有遗漏。

5. 波形怎么看——三步读懂短路暂态规律

5.1 第一步：看短路电流的整体形状

打开Simulink的Scope窗口，先别急着截图。第一件事是看短路电流波形在短路瞬间的形状变化。

正常运行时电流幅值很小（负载电流），0.1秒短路发生后，电流瞬间跳变到十几千安的幅值，而且注意观察波形相对时间轴是否对称。三相短路电流在暂态初期明显偏在轴的一侧，这就是非周期分量叠加的效果。偏得越明显，说明短路初相角越靠近电压过零时刻，电感储能初始条件越极端。

大约0.1秒到0.11秒之间，电流峰值会达到一个明显的最高点，这个最高点就是冲击电流。用手动光标工具量一下，数值应该接近4.5kA。如果量出来比理论值高，比如超过5kA，通常有两种可能：一是冲击系数取得过大，二是变压器短路阻抗百分比设置得和实际不一致。

5.2 第二步：看电压跌落与恢复

同一个Scope里同时观察母线电压。短路瞬间，母线电压会从10.5kV相对值1.0迅速掉到一个很低的水平，三相对称短路时三相电压同时跌落，不存在负序或零序分量，所以三个相电压波形仍然是对称的。

电压跌落深度取决于短路点到母线之间的阻抗大小。如果母线就在变压器出口，电压几乎跌到0；如果距离短路点经过一段线路，电压跌落程度会“打折”。这就是为什么保护装置需要用“距离”或“方向”来区分故障位置——电压不同，对应故障点的远近也不同。

故障在0.3秒切除后，母线电压重新恢复，但恢复过程不是瞬间的，会有一次小幅度的过冲然后进入稳态。如果过冲幅度很大，说明发电机的励磁调节系统在起作用，这可能引出一个新的研究课题——暂态稳定分析。

5.3 第三步：对比不同短路时刻的影响

这里分享一个很多人忽略的细节：故障发生的时刻不同，暂态波形形状差异非常大。

可以在Simulink里把Three-Phase Fault的切换时间改成0.103秒（改变短路初相角），重新仿真一次。你会发现短路瞬间的电流波形“偏移”方向和幅度都不一样了。原因是短路电流非周期分量的初始值等于短路瞬间电流周期分量瞬时值与故障前负载电流瞬时值之差，初相角不同，这个差值就不同。

这个细节在实际工程中很重要。保护装置的整定计算一般按最严重情况考虑（冲击电流取最大可能的相位角），但如果做短路电流实测记录比对，就必须记录故障的真实发生时刻，否则波形对不上会误判设备状态。

6. 常见问题与避坑指南——我踩过的那些坑

6.1 仿真不收敛怎么办

Simulink仿真三相短路最常见的坑是数值不收敛。症状是仿真跑到一半报错，或者波形发散成锯齿状。我遇到这种情况的排查顺序是：先看Power GUI的仿真类型，必须在偏好设置里选择“Continuous”（连续系统仿真），选成“Discrete”（离散）且步长设置不合理时，短路电流的非周期分量可能被严重扭曲。

再看求解器，不要用默认的ode45，换成ode23tb或ode15s这类刚性求解器。三相短路系统的阻尼角频率较高，短路瞬间又是一个剧烈的阶跃变化，属于典型的刚性系统，默认求解器容易逼不出来。

最后看元件参数，特别是线路的电阻值。如果线路电阻设得太小（理想情况模拟），非周期分量衰减时间常数过长，仿真收敛性也会变差。如果只是为了快速验证电路拓扑，可以在线路里串一个小电阻（比如0.001Ω）兜底。

6.2 仿真结果与手算对不上

这种情况我至少遇到三次，每次都花了不少时间排查。最常见的原因是变压器接线方式设置不对，特别是Δ-Yg接线的中性点接地方式没有选对，导致零序网络和故障边界条件对不上。

第二个常见原因是参数单位混乱。Simulink里既可以填有名值，也可以选标幺值，一旦某个元件用了标幺值而另一个用了有名值，计算结果的误差会非常离谱。建议全模型统一使用有名值，或者统一使用标幺值，不要混用。

第三个原因更隐蔽——变压器分接头变比设置。Simulink变压器默认的额定变比是基于一次和二次额定电压的，如果你在二次侧设置的是10.5kV还是10kV，会在潮流计算时影响短路点的初始电压，进而影响短路电流。我习惯在变压器分接头设置里把变比调成110/10.5，和手算基准电压保持一致。

6.3 波形看起来“不对劲”的三种典型情况

波形异常的情况各有各的特征，我给你列一个速查表，照着对就能快速定位问题：

另外补充一个学术论文里特别常见的坑：很多人把短路故障时间设置在0秒，结果仿真从零时刻开始就是暂态过程，初始潮流的建立和故障同时发生，波形一片混乱。建议所有暂态仿真的故障时间都设在0.05秒甚至0.1秒之后，确保系统先进入稳态再发生短路故障，这样分析波形时才能分清哪个是短路暂态，哪个是启动暂态。

7. 从三相短路再往前走一步：这个方法的扩展空间

做完一个基础的三相短路分析，这套方法完全可以扩展出去。我在实际项目中接着做过几个方向，简单提一下供你参考。

一个是扩展到不对称短路分析。把三相短路模块换成单相接地短路、两相短路或两相接地短路，用对称分量法拆解正序、负序、零序网络，在MATLAB里分别计算各序网络阻抗，再按故障类型组合求解。你会发现三相短路虽然是最严重的故障状态，但不一定是最常出现的故障，实际电网中单相接地短路占了故障总数的大头。

另一个扩展方向是结合保护装置的整定计算。短路电流计算结果出来后，可以直接用MATLAB写脚本自动计算限时电流速断保护、过电流保护的整定值，每次修改系统参数后一键更新整定结果。我做过一个小工具，输入系统阻抗参数，输出保护整定值表和灵敏度校验结果，后期配合批量仿真，工作效率提升非常明显。

再一个方向是短路电流的动态过程用于设备选型校核。比如断路器开断能力校核，需要知道的是断路器开断瞬间的短路电流值，而不是短路开始时的最大冲击电流。如果把仿真波形导出后做数据后处理，结合断路器分闸时间，就能算出实际开断电流，这个数据才是设备选型的直接依据。

写在最后的实操体会

做三相短路分析，工具只是手段，理解物理过程才是核心。MATLAB帮我们省去了大量手算的繁琐，但前提是你得先知道自己的计算目标和理论边界在哪里。建议第一次做这个课题的朋友，哪怕你的任务只是Simulink仿真出一张漂亮的波形图，也务必先手动算一遍短路电流的理论值，用理论值去校准仿真结果，而不是仿真结果说什么就信什么。这个“先理论后仿真”的习惯，比任何工具技巧都值钱。

另外一个小经验：仿真做完后，记得把关键波形用MATLAB的exportgraphics函数导出成高分辨率图片，把仿真参数和理论计算结果做成一张汇总表。做课程设计或毕业设计时，这套“理论推导+仿真验证+参数汇总”三件套能让你的工作量直接翻倍呈现，答辩时心里也踏实。

本文还有配套的精品资源，点击获取