简介：针对电力系统三相短路故障的建模、仿真与暂态分析需求，这套基于MATLAB/Simulink的配套资料给出了完整可复现的解决方案，适合电气工程专业本科生、研究生及需要快速上手的工程师。资源包共12个文件、压缩后仅1.09MB，包含1个MATLAB脚本（.m）、1份PDF课设报告、4份TXT电路数据与过程输出，以及6张PNG格式的节点阻抗/导纳矩阵与仿真结果图，覆盖从原始网络参数到发电机、负载等效接入的完整建模链条。资料着重演示了Simulink中搭建三相短路模型、设置故障类型与位置、运行仿真并处理暂态数据的方法，同时涉及傅里叶变换分析、牛顿-拉夫森法求解系统平衡方程等关键计算。通过研读报告与代码，读者可快速复现短路电流、电压波形分析流程，理解保护设备动作逻辑，并可直接借鉴到自己的课设或工程验证中。该资源目前已有2576人学习，具有较高的参考价值。
1. 为什么用MATLAB做三相短路分析
1.1 短路分析在电力系统里的地位
做电力系统的人都知道，三相短路是系统中最严重、也最需要优先校核的故障类型。无论是发电厂、变电站还是工业配电网，规划设计阶段都得把短路电流算清楚，否则选断路器、整定保护、校验热稳定和动稳定都会失去依据。我最早接触这个题目是在读研究生的时候，当时导师丢给我一个IEEE 9节点系统，要求手算三相短路电流，那叫一个痛苦——手算对称分量、画序网络、化简等效阻抗，小系统还好，节点一多直接算到怀疑人生。
后来我转到MATLAB/Simulink做仿真，才发现用软件做三相短路分析不是省事一点半点。它不光是能把波形画出来，更重要的是可以把系统模型复用到不同的故障场景里，改个参数就能重新跑一遍。对初学者来说，通过仿真去理解短路的物理过程，比如电压跌落、电流冲击、暂态衰减，比死磕公式要直观得多。对工程师来说，MATLAB脚本加Simulink模型这套组合，既能做快速估算，也能做详细暂态分析，正好覆盖工程设计的两个阶段。
1.2 MATLAB/Simulink 在短路分析中的核心优势
先说计算层面。手算三相短路通常要做一堆假设：忽略负荷、忽略电阻、认为电源电压恒定。这些假设算出来的结果往往偏保守，但有时候会偏离实际情况。Simulink里用三相电压源、变压器、线路、故障模块搭模型，能考虑负荷影响、线路阻抗、变压器励磁特性，仿真结果更贴近真实故障过程。你甚至可以把同步发电机模型放进去，观察次暂态、暂态和稳态电流的变化过程，这是手算完全给不了的。
再说复现性。我自己很烦那种“改一处参数就要重画一张图”的操作，Simulink里改故障电阻、故障时刻、接地方式，点一下运行就出结果，脚本化之后还能批量跑参数扫描。做完一次分析之后，模型文件留着，以后遇到类似的系统改改参数就能复用。这种效率提升，说实话在工程项目里非常实在。
2. 三相短路的数学模型与关键参数
2.1 对称短路的基本方程
先别急着拖模块，得把原理捋清楚。三相短路是对称故障，正序、负序、零序里只有正序存在，所以计算相对简单。假设系统从短路点看过去的等效阻抗为 Z（实际中主要是电抗 X），电源相电压为 E，则短路电流的稳态有效值就是：
I_f = E / (Z + Z_f)
Z_f 是故障点附加阻抗，金属性短路时 Z_f = 0。这个式子看着简单，但工程上有个坑：同步发电机的电抗在短路瞬间不是固定的，会从次暂态电抗 Xd'' 过渡到暂态电抗 Xd'，最后到稳态同步电抗 Xd。所以手算时得看你想算哪个时刻的电流——瞬态的用 Xd''，0.1秒左右的用 Xd'，稳态用 Xd。Simulink里如果用的是简化电源加内部阻抗的模型，这个过渡过程其实是被等效掉的，结果偏向稳态值；如果用详细同步电机模型，就能看到明显的衰减过程。
再往下是冲击电流 i_m，这是短路后半个周波左右的瞬时峰值，用来校核设备动稳定。工程上近似取：
I_m = K_m * sqrt(2) * I_f
K_m 是冲击系数，一般取1.8到1.9，取决于系统电阻与电抗的比值 X/R。X/R越大，K_m越大，衰减越慢。这个参数在短路分析里特别关键，因为实际短路电流是周期分量叠加非周期分量，非周期分量的初始大小跟短路瞬间的电压相位有关。用MATLAB写脚本算冲击电流的时候，我习惯把非周期分量按最恶劣情况考虑，也就是电压过零时短路，这样算出来的峰值最大，用于校验最保险。
2.2 短路容量与保护整定的关系
短路容量 S_f = sqrt(3) * U_N * I_f，这个数据直接决定断路器的开断容量选择。变压器容量、线路长度、系统阻抗都会影响短路容量，在做电力系统规划时，通常需要把最大运行方式下的短路容量算出来，作为开关设备的选型依据。Simulink里可以直接测故障电流有效值，然后用公式换算，也可以直接用Powergui模块里的“Fault Analysis”工具一键计算。我用过几次，那个工具对IEEE节点系统特别方便，能自动生成故障前、故障后的潮流和短路电流，省掉大量手工查表时间。
这里有个经验：Simulink仿真测短路电流，一定要等待波形进入稳态之后再读有效值，尤其是带详细电机模型时，故障后前几个周期电流非常大，那是次暂态分量，跟稳态短路电流不是一个概念。很多新手直接取第一个峰值当稳态值用，结果偏大很多，保护整定值也跟着虚高。
3. Simulink仿真模型搭建实操
3.1 从单机无穷大系统开始搭建
第一次做三相短路仿真，别一上来就搭九节点大网，先从单机无穷大系统练手。所谓无穷大系统，就是用一个理想三相电压源串联等效阻抗去模拟电网，电压和频率在故障前后基本不变。在Simulink里，拖入Three-Phase Source模块，设置线电压有效值比如110kV、频率50Hz，内部串联电阻和电感按系统短路容量折算。举个例子，如果系统短路容量是1000MVA，电压110kV，那等效阻抗大约为：
Z = U_N^2 / S_f = (110k)^2 / 1000M = 12.1 Ω
把这12.1欧姆拆成电阻和电感，比如R=0.1Ω，X=12Ω，L = X/(2pif) ≈ 0.0382H。这样电源模型就具备了实际电网的“支撑强度”。
接着接入一条输电线路，用Three-Phase Series RLC Branch表示，设置正序阻抗和零序阻抗。注意Simulink里默认显示的是R、L、C的值，如果你用的是线路单位长度阻抗，比如每公里0.1Ω + 1mH，那就得把长度乘上去。我最初遇到过单位搞错的情况，把公里数当欧姆填进去，仿真出来电流小得离谱，排查了半天才发现是数量级问题。
3.2 Three-Phase Fault模块的关键设置
故障模块是核心，在Simulink库里的路径是Simscape / Electrical / Specialized Power Systems / Fault Blocks / Three-Phase Fault。双击打开，能看到三相选择框和接地选项。做三相短路，要把A、B、C三相全勾上，同时勾上“Ground”表示三相接地短路。如果你想要三相不接地短路，那就只勾三相不勾Ground，不过工程上最常见的是三相接地故障，波形差别不大，但零序分量会有区别。
故障时序设置里有个Transition times参数，我一般填[0.1 0.3]，意思是在0.1秒投入故障，0.3秒切除。这里要特别注意的是Fault resistance，金属性短路填0，但实际仿真中为了数值稳定，我建议填一个很小的电阻，比如0.001Ω，避免纯电感电路出现数值振荡。真的填0也能跑，但步长要设得很小，仿真时间会翻好几倍，没必要。
还有一个容易忽略的设置是External control，如果勾选了外部控制，故障投切需要额外给控制信号，适合做需要精确控制故障时刻的场景。我一般用内部定时就够了，省事。
3.3 测量模块与示波器配置
为了观察短路前后的电压电流，在故障点前需要加Three-Phase V-I Measurement模块，它可以把三相电压和电流转成Simulink信号，方便接示波器或导入工作区。模块的输出有两组：一组是电压，一组是电流。用Demux或者直接接Scope就能看波形。
示波器配置有讲究。短路电流在故障瞬间可能从正常值跳到几十倍，如果示波器Y轴范围不合适，波形会被削顶。我习惯先把带限幅的magnitude模块放在电压电流后面，显示三相幅值；如果需要看单相波形，再单独把A相电流接出来。另外，仿真时间要设足够长，比如故障0.1秒投入0.3秒切除，那仿真停止时间至少设0.6秒，才能看到恢复过程。采样步长用默认的ode23tb就行，但最大步长要限制在1e-4左右，不然短路瞬间的陡峭波形会失真。
4. 仿真运行与结果解读
4.1 波形里你看不到的东西
运行完以后，别急着截图。先看电流波形，A相电流在0.1秒处突然变大，然后出现一个衰减的直流分量，导致前几个周期波形上下不对称，这是短路电流非周期分量的典型特征。从保护角度讲，这个非周期分量会让电流过零时间推迟，影响断路器开断能力，所以波形上一定要能看到这个现象，才算仿真做对了。如果波形直接就是对称的正弦波，那大概率是故障模块或测量模块接线有问题，或者模型里电机没接对。
电压波形方面，短路瞬间故障点电压会骤降。如果短路是金属性的，故障点残压接近零；如果是经过渡电阻短路，残压会高一些。我在做实验时发现，很多人只盯着电流看，忽略了电压。其实电压跌落的深度直接影响低电压穿越、电压暂降治理等应用，做电能质量分析时电压波形比电流更重要。用FFT工具还能分析谐波含量，不过三相短路是工频故障，谐波不是主要矛盾。
4.2 仿真结果和手算值对比
不管是做课程设计还是工程项目，仿真结果都需要验证。我一般会把手算的稳态短路电流和仿真波形里故障后0.3秒左右的稳态幅值对比。还是刚才的例子，110kV系统，短路容量1000MVA，短路电流 I_f = S_f / (sqrt(3) * U_N) = 1000M / (1.732 * 110k) ≈ 5.25kA。仿真结果在故障后电流幅值应该稳定在5.2kA左右，偏差在5%以内都算正常。如果偏差大，首先检查电源等效阻抗是否正确，其次看线路阻抗有没有填错单位。
还有一个对比项是短路瞬间的冲击电流峰值。如果模型里线路电阻很小，X/R很大，冲击系数取1.8，那峰值大约是 1.8 * sqrt(2) * 5.25 ≈ 13.4kA。仿真波形里第一个峰值如果落在13kA上下，说明模型参数是合理。我第一次做的时候，由于故障模块设置成了0.01秒投入，恰好赶上电压过零点，非周期分量特别大，第一个峰值到了15kA，用这个值校核设备就不对了，因为那是极限情况，不是典型情况。所以对比时要明确你手算用的K_m假设对应的是哪种初相角。
4.3 参数扫描：故障时刻与过渡电阻的影响
仿真模型搭好之后，最有价值的操作就是参数扫描。在MATLAB里写个循环，改变故障模块的Transition times或者Fault resistance，然后把每次仿真的电流峰值导入数组，再画成曲线。我试过把故障电阻从0.001Ω扫到10Ω，短路电流急剧下降，最后趋近于负荷电流。这个曲线在实际工程中很有用，比如评估高阻接地故障对保护灵敏度的影响，或者计算电弧短路时保护能否可靠动作。
故障时刻的扫描则更偏物理过程。把故障投入时刻从电压过零附近逐渐偏移到电压峰值附近，你会发现短路电流初始直流分量的大小完全不同。电压过零时短路，直流分量最大；电压峰值时短路，直流分量接近零。这个现象在书本上叫“合闸相角影响”，看文字没感觉，自己扫一遍就彻底明白了。
5. 常见问题与排查技巧实录
5.1 仿真不收敛，报错“Trouble solving algebraic loop”
这个问题我遇到至少三次，基本都是因为纯电感回路里突然加故障，电流变化率无穷大，导致数值求解器无法处理。解决办法有三个方向：一是给故障模块加一个很小的故障电阻，比如0.001Ω，让回路不再是纯感性；二是调整求解器，把最大步长缩小到1e-5，或者换成ode15s这种刚性求解器；三是在线路旁边并联一个很小的电导。最推荐第一个方向，简单有效，而且对结果影响可以忽略。
5.2 短路前后电流波形不连续，出现尖刺
如果故障瞬间波形出现向上或向下的尖刺，那通常是采样步长太大，没有捕捉到故障切换的准确时刻。Simulink里默认的步长是自适应，但在故障这个突变点上，求解器可能没有精确插值到切换点。把Max step size设为1e-4甚至1e-5，尖刺基本就消失了。另外检查一下故障模块里的Sample time，如果填了-1表示继承，但有时候模块内部默认0.001s，太糙了，改成1e-5也行。
5.3 三相电流波形不完全对称，总有轻微的偏置
这其实不是故障，而是非周期分量的存在，导致前几个周期三相电流的偏置方向不一样。理论上三相互差120度，非周期分量的初值也不同，所以波形看起来不对称。如果你等了0.5秒以后波形还是明显不对称，那就要怀疑故障模块是否真的让三相同时接地了，或者线路参数三相不平衡。检查方法很简单：把故障时间设成0.1到0.2秒的超短故障，看故障后半波是否对称。
5.4 实用避坑清单
- 电源模块里的
Specify impedance using short-circuit level选项，可以直接填短路容量和X/R，省去手工折算阻抗，但要注意基准电压要填对。 - 变压器模块的联结组别会影响零序网络，做三相短路时影响不大，但做单相接地短路时必须复核。
- 测量模块的电压测量是线电压还是相电压，要看模块设置，读波形时别拿线电压峰值去对比相电压有效值。
- 保存波形数据用
To Workspace模块，把时间序列存成timeseries格式，方便后期用MATLAB脚本画图或算峰值。 - 如果你的模型里用了同步电机，记得先把电机的机械功率和励磁设置合理，不然故障后转速飞升，仿真结果没法看。
我做这个题目最大的体会是：仿真软件再方便，也不能替代对短路物理过程的理解。很多时候波形不对，不是软件问题，而是对冲击电流、非周期分量、短路容量的概念没吃透。所以建议大家在跑通仿真的基础上，可以尝试把故障电阻、故障时刻这些参数“极端化”，多观察波形的变化规律，这个过程中获得的手感，远比最后一张完美的图更有价值。后面我还想把这个模型扩展成含分布式电源的配电网短路分析，本质上就是多几个电源和逆变器接口的事，等做完了再跟大家分享。
本文还有配套的精品资源，点击获取
