简介：面向轨道交通电气系统设计、维护及相关研究人员的MATLAB/Simulink建模仿真资源，围绕地铁牵引供电系统展开，涵盖外部电源、24脉波整流变电站、直流母线、接触网与钢轨大地模型，以及含逆变电路和牵引负荷的牵引传动系统。基于成都地铁实际参数进行空载、负载与故障工况仿真，并与现场故障录波对比验证模型有效性，可帮助工程师理解建模方法、优化系统设计并开展故障分析。资源为单个PDF文件，共671KB，内容包含论文复现所需的详细可运行代码及逐段中文解释，涉及24脉波整流移相变压器组、接触网PI型分布参数模型、SVPWM逆变控制及系统级仿真框架等关键技术，并体现了分层验证体系与多时间尺度仿真策略。已有244人学习下载，适合需要对照论文落地仿真、掌握地铁牵引供电建模与故障仿真方法的工程技术及研究人员。

1. 地铁牵引供电仿真：为什么值得啃下这篇论文复现

地铁牵引供电系统是城市轨道交通的能量大动脉，而基于MATLAB/Simulink的建模与仿真，一直是硕博论文和工程验证里出现频率最高的方向之一。这篇标题对应的复现工作，核心不是把别人的图重新画一遍，而是用Simulink把牵引变电所、整流机组、接触网、钢轨回流和列车负载串成一个能跑的动态系统，再人为注入短路、整流器故障等工况，从波形里读出保护配合逻辑和故障特征。对做供电专业的学生来说，它是打通电路理论、电力电子和继电保护三块知识的最佳实践载体；对现场工程师而言，一套可调参的仿真模型能提前验证保护定值和运行方式变更的影响，远比直接动真设备安全。

我最早做这个方向时，被Simulink里整流桥的数值震荡折磨过，也因忽略了钢轨散流模型导致故障电流和实测对不上。这篇文章会把建模流程、关键参数、故障注入方法和常见坑一次讲透，代码段都能直接复制到你的Simulink环境里改参数跑通，希望帮你在复现这条路上少走几天弯路。

2. 从牵引变电所到接触网：Simulink模型拓扑与最小可跑通方案

2.1 牵引供电系统的基本组成与仿真模型选型

地铁牵引供电系统的主干结构并不复杂：高压交流电源（通常110kV或35kV）进入主变电所，降压后供给牵引变电所；牵引变电所内通过整流机组把交流变成直流1500V（部分线路用750V），输出到接触网（或第三轨）；列车通过受电弓取流，电流经钢轨和回流线流回牵引变电所。模型里必须把这几个环节都放进去，否则故障分析时回路不完整，短路电流路径会少一条。

在MATLAB/Simulink里做这个系统，选型上有三条路：一是用Simscape Electrical（旧称SimPowerSystems）里的基础电气元件自己搭整流桥、变压器和直流母线，灵活度最高；二是用Simscape Electrical的Specialized Power Systems库里的“Universal Bridge”和“Three-Phase Transformer”模块，搭起来快，内部开关逻辑已经封装好；三是用Simscape基础元件库做物理域建模，但和传统Simulink控制信号交互稍麻烦。我一般选第二种，原因是论文复现时整流机组的脉波数、变压器联结组别都需要显式设置，Universal Bridge可以直接指定桥臂数和开关器件类型，最贴近技术论文里的描述。

模型拓扑我建议按“交流电源—整流变压器—整流桥—直流母线—接触网—列车负载—钢轨回流”的单线图来搭。为了后续能注入故障，直流母线正极和接触网之间加一个分段开关，接触网末端预留故障点。仿真步长不用全局固定，Simscape Electrical会自动变步长，只需要把最大步长限制在10微秒以内，避免整流桥换相时数值突变引起步长收缩过慢。

2.2 用Simscape Electrical搭建整流机组与直流母线

第一步是在Simulink里新建模型，拖入以下模块：Three-Phase Source（三相交流源）、Three-Phase Transformer（整流变压器）、Universal Bridge（整流桥）、Series RLC Branch（直流平波电抗器）、DC Voltage Measurement和Scope。交流源设为35kV/50Hz，变压器二次侧线电压设为1180V左右——这个值是为了让空载直流电压达到1500V的1.35倍关系，即约2025V空载电压，带负载后自然跌落到1500V左右。整流桥选“Universal Bridge”，Power Electronic Device选“Diode”，桥臂数选“3”，这样就是典型的三相6脉波整流。

变压器联结组别是关键，工程上常用Dd0或Dy11。6脉波整流用Dd0和Dy11均可，但12脉波必须用两组变压器分别接Dd0和Dy11，利用相位差30度来抵消5、7次谐波。如果论文里讨论24脉波，就要在变压器二次侧加延边三角形移相绕组。复现时先跑通6脉波，再扩展到12脉波，验证步骤更清晰。

下面是最小模型的参数脚本，可以写在模型初始化回调里：

% 牵引供电系统最小模型参数初始化 % 交流电源部分 V_line_rms = 35e3; % 电源线电压有效值，单位V f_grid = 50; % 电网频率，单位Hz R_source = 0.5; % 电源等效内阻，单位ohm L_source = 5e-3; % 电源等效电感，单位H % 整流变压器参数 V_primary = 35e3; % 一次侧线电压，单位V V_secondary = 1180; % 二次侧线电压，单位V，空载直流电压约2025V R_leak = 0.02; % 变压器短路电阻标幺值或折算到二次侧的欧姆值 L_leak = 0.15e-3; % 变压器短路电感，单位H % 整流桥与直流侧 cnv_pulses = 6; % 整流脉波数，6脉波基础模型 L_smooth = 5e-3; % 直流平波电抗器，单位H R_smooth = 0.01; % 平波电抗器内阻，单位ohm C_dc = 10e-3; % 直流支撑电容，单位F，论文中有时不装，按需开启

这段脚本的关键是变压器二次侧电压的取值逻辑：三相桥式整流电路空载直流电压是二次侧线电压有效值的1.35倍，所以要让空载电压约2025V，二次侧线电压就是2025/1.35≈1500V？这里其实有个常见误解。工程上牵引变压器二次侧额定电压取1180V~1200V，是考虑负载压降和整流管压降后的设计值，空载直流电压约1600V上下，满载时直流电压约1500V。不同论文取值有差异，复现时最好先按论文给定的空载电压折算回去。如果论文没给，就用二次侧1200V起步，跑负载看直流电压是否落在1500V附近，再微调。

2.3 列车负载与线路参数的模拟方法

列车负载在牵引供电仿真里是最容易“凑合”的部分，但它直接决定直流电流波形。最简单的做法是恒定电阻负载，用直流电压除以期望电流算电阻值；更接近实际的用受控电流源或Simscape的DC-DC变换器模型，模拟列车牵引时的恒功率特性和制动时的回馈特性。论文复现建议分步走：第一步用恒定电阻跑通系统，第二步换成可变电阻模拟列车启动—加速—惰行—制动四象限工况。

线路参数不能只给一个纯电阻，接触网和钢轨都有单位长度电阻、电感和对地电容。1500V直流牵引网通常用走行轨回流，钢轨对地有泄漏电导，这个参数在故障分析中影响站间距离与短路电流的对应关系。在Simulink里可以用“PI Section Line”模块或者多段串联的RLC分支来模拟。每段线路长度按站间距取值，例如1.5km一段，每段电阻0.01Ω/km、电感0.8mH/km、对地电导0.5S/km（潮湿隧道环境取值会更高）。注意钢轨回流路径必须和接触网正极一起构成回路，否则故障时仿真会提示“Ground connection”问题。

列车负载模拟我常用一个受控电流源，输入信号是列车牵引电流曲线。这个电流曲线可以用Lookup Table导入实测或论文给定的数据点。比如某列车启动电流从0平缓升到1500A，维持到速度拐点后线性下降。生成这个曲线的脚本如下：

% 生成列车牵引电流曲线，存入Lookup Table t = [0 5 10 15 20 25 30]; % 时间轴，单位s I_train = [0 600 1500 1500 900 300 0]; % 牵引电流，单位A % 在Simulink中创建Lookup Table时直接使用这两组数据 plot(t, I_train, 'o-'); xlabel('时间 (s)'); ylabel('牵引电流 (A)'); title('列车牵引电流参考曲线'); grid on;

电流源的输出接在接触网正极和钢轨回流之间，方向是从接触网流向钢轨，正号代表牵引取流，负号代表制动回馈。Simulink里的受控电流源需要并联一个大电阻（比如1MΩ）避免数值孤立，这个并联电阻不影响正常工况，但在短路故障瞬间会产生一个微小分流，分析故障电流时需要减掉这部分误差。列车位置的移动性可以用多个电流源加开关在不同时间接入来实现，但初期模型里固定一个位置就够了，论文复现的重点是波形的“形态”，不是精确的时空分布。

3. 整流机组与列车负载：两个决定仿真精度的关键技术

3.1 整流机组等效模型与脉波数选择

整流机组是地铁牵引供电系统里最核心的电力电子环节，它的脉波数直接决定了直流侧电压脉动和注入交流系统的谐波含量。6脉波整流是最基础的结构，直流电压纹波较大，5次、7次谐波突出；12脉波由两组6脉波整流器通过移相30度并联或串联组成，5、7次谐波被大幅抵消，直流纹波也明显变小；24脉波则是在12脉波基础上再加移相绕组，工程上用于大型主变电所。

在Simulink里实现12脉波，常见做法是使用两台三绕组变压器，一次侧并联，二次侧一组Y接、一组Δ接，分别带各自的整流桥，然后直流侧串联输出。注意两套整流桥的直流输出电压必须极性一致且幅度接近，否则会形成环流。参数上，两组变压器的变比必须严格一致，短路阻抗误差尽量控制在±5%以内。我用Universal Bridge分别搭两套整流桥，桥臂数和器件类型保持一致，只改变变压器联结组别。

为了在论文里定量说明脉波数的影响，可以在模型里加一个FFT Analysis工具，稳态运行后测量直流电压纹波系数和谐波频谱。纹波系数的定义是直流电压交流分量有效值与直流平均值之比。6脉波理论上纹波系数约4.2%，12脉波降到1.0%左右，24脉波进一步降到0.26%。如果你复现的论文给出了纹波系数对比表，仿真值应该落在理论值附近，差太多就检查整流桥触发角是否为零或变压器漏抗是否过大。

整流机组的等效模型还有一层：IGBT或晶闸管开关的动态特性。在故障暂态分析里，开关器件的关断延迟和反向恢复过程会影响短路电流的初始上升率。Universal Bridge默认的Snubber电路参数也需要注意，默认Rs=1e5Ω、Cs=inf，这相当于把开关理想化。若要更真实的故障电流，需要把Snubber电路的电容设为有限值，例如10nF，电阻降到10Ω，但这会显著增加仿真时长，建议只在最终验证阶段开启。

3.2 列车牵引/制动特性曲线怎么进模型

列车负载的动态特性直接影响直流母线电压的波动范围。地铁列车的运行工况分为牵引、巡航、惰行、制动四象限，其中牵引时列车从接触网吸取功率，制动时一部分能量通过再生制动回馈电网。Simulink里最直观的建模方式是用受控电压源和受控电流源组合：牵引时电流方向为正，制动时为负，电压始终由直流网侧决定。

我把列车模型拆成两部分：功率计算和电流注入。功率计算根据列车速度、加速度和质量，按简化力学公式算出瞬时牵引力，再乘以车速得到机械功率，除以效率后得到电功率；电流注入则把这个电功率除以直流母线电压得到参考电流，送给受控电流源。这个模型能体现恒功率特性——母线电压下降时电流会上升，和实际列车特性一致。

下面是一个最小实现：

function I_ref = train_current(v, a, m, eta, Vdc) % 列车牵引电流计算函数 % v 车速 m/s，a 加速度 m/s^2，m 质量 kg，eta 效率，Vdc 直流母线电压 V % 基本阻力简化为 A + B*v + C*v^2 A = 4500; % 滚动阻力系数 N B = 90; % 线性阻力系数 N·s/m C = 15; % 空气阻力系数 N·s^2/m^2 F = m * a + A + B*v + C*v^2; % 所需牵引力 N if F > 0 P_mech = F * v; % 机械功率 W P_elec = P_mech / eta; % 牵引时电功率大于机械功率 else P_mech = F * v; % 制动力，F为负则P_mech为负 P_elec = P_mech * eta; % 回馈时电功率小于机械功率，且为负 end I_ref = P_elec / Vdc; % 直流侧电流参考值 A end

在Simulink里用MATLAB Function模块调用这个函数，输入车速、加速度和直流电压测量值，输出I_ref给受控电流源。注意初始时刻车速为0时，理论上牵引力公式里v=0会导致功率为0，但实际列车启动时会有静止牵引力，所以需要在模型里加一个启动逻辑：当v<0.1m/s且F>0时，按恒力矩处理，即I_ref用F乘以一个等效半径换算的常数代替。忽略这个逻辑会导致仿真启动瞬间电流参考值跳变到0，列车根本没法“动起来”，这个坑我踩过两次。

3.3 直流侧短路故障的Simulink仿真实现

直流侧短路是地铁供电系统最严重的故障类型，接触网正极对钢轨直接短接，短路电流可达数千安到数十千安。仿真中实现短路最简单的方式是用一个断路器（Three-Phase Breaker或Ideal Switch）加一个很小电阻，在指定时间闭合。为了模拟不同短路位置的过渡电阻，可以用Breaker串联一个可变电阻。

故障时序设定：先让系统空载或带负载稳定运行到0.1s，再在0.1s时闭合短路支路，持续0.15s后由保护动作断开。保护动作在仿真里用Breaker的打开时间模拟，也可以用Simulink的Step模块触发开关。关键是要记录短路发生前后直流电压和电流的波形，观察冲击电流峰值、电压跌落幅度和恢复到正常的时间。

注入故障时，我建议检查三个细节：第一，短路支路的闭合电阻不能设置为0，否则Simulink会出现代数环或数值奇异，一般取1e-3Ω；第二，故障点与整流机组之间必须有线路阻抗，否则短路电流峰值会远超真实值，因为变压器漏抗和直流电抗器是限制短路电流的主要元件，如果线路阻抗被忽略了，峰值电流可能高估30%以上；第三，故障后的恢复过程要保留平波电抗器，它决定短路电流上升率。

下面是一个简单的故障控制脚本，可以用作Model Callback，也可以放在Constant模块的Signal Builder里：

% 直流侧短路故障控制 % 使用Simulink Signal Builder或Step模块，本脚本生成故障开关时间序列 t_step = [0 0.1 0.1 0.25 0.25]; % 时间节点 Brk_state = [0 0 1 1 0]; % 0表示打开，1表示闭合 % 对应的断路器控制信号，在Simulink中通过Repeating Sequence模块导入 % 短路电阻取0.005ohm，闭合后经历0.15s的保护动作时间

实际运行时会发现，短路瞬间直流电压从1500V骤降到接近0，直流电流迅速攀升到峰值，随后由于交流侧阻抗和直流电感的作用逐渐衰减到稳态短路电流。这个波形的上升沿斜率主要受直流回路电感限制，而稳态值受回路电阻限制。论文里如果给出了短路电流的解析计算式，仿真峰值和解析值对不上时，优先排查故障点位置和直流侧总电感是否正确。

4. 故障分析：把短路电流波形和继电保护动作讲清楚

4.1 整流器内部短路故障分析

整流器内部故障常见有整流二极管击穿短路、同一桥臂上下管直通、交流侧相间短路等。在Simulink里模拟二极管击穿，最简单的做法是在被击穿的二极管两端并联一个受控电压源或小电阻开关，但这种建模太理想化，不能反映故障对交流侧的影响。更接近实际的建模方法是在整流桥外面并联故障支路，故障点设在整流桥直流输出端，模拟逆变失败导致的直流母线短路。

另一种常见故障是整流变压器二次侧匝间短路。由于变压器模型在Simscape Electrical里是基于理想变压器的，匝间短路需要把变压器拆成多绕组模型，在二次绕组上并联一个小电阻来模拟局部短路。复现论文里讨论整流器故障时，重点看的是故障后的特征：交流侧电流出现严重不平衡，直流电压跌落且伴有较大谐波，直流电流出现周期性尖峰。

我在仿真整流器内部故障时发现一个需要注意的现象：整流桥内部故障如果不改变触发条件，仿真器会自动报“Singularity”错误，因为故障支路把电源直接短路了。解决办法是在故障支路串入一个0.01Ω小电阻，并把仿真器最大步长放宽到1e-5s。如果仍然不收敛，把整流桥的Snubber电路电阻调大，比如从1e5Ω调到1e6Ω。

4.2 接触网短路与保护动作配合

接触网短路是地铁系统中最需要保护配合的故障，因为列车不断移动，故障点位置不同导致短路电流大小不同。仿真里通常设置几个典型位置：牵引变电所出口、线路中点、末端。出口短路时短路电流最大，末端短路时电流最小，保护需要通过电流大小和距离信息来区分故障区段，避免越级跳闸。

Simulink里实现多位置短路很简单，把接触网分成若干段，每段末端接一个断路器到钢轨，不同时刻闭合不同的断路器就行。但要注意，各短路点不应该同时闭合，否则会重合短路。合理的做法是用一个故障序列控制逻辑，例如0.1s在第一个点短路，0.3s断开后又切换第二个点。

保护动作配合仿真可以用简单的继电器逻辑实现：测量直流牵引网电流，超过设定值后延时动作，断开牵引变电所出口断路器。为了模拟整定值的配合关系，可以分别设置两段保护——瞬时速断和定时过流。在Simulink里用Relay模块加延时环节搭，输出触发断路器的信号。关键是把保护动作时间设定得比故障持续时间短，否则断路器不会断开，仿真会一直维持短路状态。

比较有指导意义的仿真实验是调整保护定值和线路参数，观察保护动作的覆盖范围。比如末端短路电流为1500A，如果速断保护定值设在3000A，那么末端故障只能靠过流保护（定值1200A，延时0.5s）动作，导致故障切除时间更晚。论文里讨论保护配合不合理时，通常就是拿这种仿真结果作对比，说明不同定值下保护范围的变化。

4.3 故障波形怎么看：电压、电流、功率的异常特征

故障分析最终要落到波形判读上。直流电压、电流波形是最直观的，但功率波形、谐波频谱、交流侧电流也要一起看。整流器出口短路时，直流电压瞬间塌陷，电流冲高后回落，功率波形会出现一个正向尖峰后持续低值；交流侧电流则出现明显的不对称，某相电流幅值异常增大。

我还有几个读波形的习惯：第一，不要只看瞬时值，用RMS转换模块计算故障期间的电流有效值，注意Simulink里到的RMS是滑动窗口值，窗口长度要设置为工频周期20ms，否则波形平滑度不够，无法和工程上的保护动作值对应；第二，用Power Spectral Density或FFT工具看故障前后谐波变化，整流器内部故障通常会产生明显的偶次谐波，而外部短路主要产生直流分量和低频分量；第三，把仿真数据和短路电流解析公式的计算结果画在同一张图里，验证模型是否在可接受误差范围内。

如果论文中提到了“直流牵引网短路电流上升率”这一指标，你需要在仿真里测量的就是di/dt。Simulink里可以对电流信号用Derivative模块求导，然后输出去找最大值。这个值受直流回路电感影响，仿真中如果发现上升率比论文给的大很多，去看看平波电抗器参数是否设置正确、接触网线路电感是否被忽略了。电感量每减少1mH，di/dt可能上升近一倍，参数敏感性很高。

5. 避坑指南：MATLAB/Simulink牵引供电仿真的5个常见坑与排查

5.1 整流桥仿真不收敛，报“Singularity”错误

现象：仿真跑到整流桥换相瞬间突然报错，提示“Derivative of state is infinite”或“Singularity”，步长反复收缩到极小值后崩溃。

原因：整流桥开关器件切换时，理想开关两端电感电流突变，Simulink的变步长求解器无法处理不连续状态；或者故障支路电阻设为了0，直接把电容或电压源短路。

解决：给所有开关器件和故障支路串入一个0.001~0.01Ω的小电阻，并检查直流侧是否串有平波电抗器；同时在求解器设置里把“Max Step Size”设为1e-5s，把“Relay”模块的切换容差调大一些。我一般还会把Universal Bridge的Snubber电阻从1e5Ω改成100Ω，这能大幅缓解突变。

5.2 直流电压偏低，空载值只有1000V而不是1500V

现象：模型搭建完成后，交流侧参数都对，但直流电压测量值远低于设计值，而且波形的脉波数量不对。

原因：变压器二次侧电压设置错误，或者整流桥桥臂数选错。比如应该选3桥臂（三相全桥）却选了1桥臂；更常见的是把变压器二次侧线电压和相电压搞混，计算空载直流电压时应该用线电压乘1.35而不是相电压。

解决：先断开所有负载和故障支路，单独测试整流机组空载输出。确认变压器二次侧额定电压是否正确，用万用表模块分别测量二次侧相电压和线电压。如果线电压是1180V，空载直流电压应该约1590V（1180×1.35），在正常范围。若实测值只有约1400V，检查整流桥是否选成了“半控”模式或触发电平有问题。

5.3 故障电流峰值比论文小很多，甚至只有一半

现象：用断路器模拟直流侧短路，理论估算短路电流峰值约9kA，仿真测量只有4.5kA。

原因：常见原因有三个：一是直流侧平波电抗器电感值设置过大，限制了电流上升率；二是交流电源等效内阻和变压器漏抗被隐式放大了，比如变压器励磁电感默认值过小；三是故障点位置离整流机组太远，线路阻抗分掉了大量压降。

解决：把平波电抗器电感从5mH降到2mH看峰值是否上升；再检查变压器“Leakage inductance”参数是否与实际短路阻抗匹配。还有一个隐性参数是“Three-Phase Source”内部的“Source internal resistance”，默认0.01Ω偏大时也会限制短路电流。把这个参数改成0.001Ω，同时确认故障断路器的导通电阻不超过1mΩ。

5.4 仿真速度极慢，跑一秒要十分钟

现象：模型加上整流桥、列车负载、线路分段后，仿真速度急剧下降，步长始终小于1e-7s。

原因：电气系统中存在极小的时间常数，比如大电容和极小电阻串联，或开关器件的缓冲电路阻容时间常数太小，迫使求解器采用极小步长。另一个常见原因是过度使用了理想开关，开关频率太高。

解决：检查所有串联电阻和并联电容组合，算出时间常数，尽量让它们不小于1e-6s；把Universal Bridge的Snubber电路电容设为Inf（即无缓冲）可以显著提速，前提是建模重点是稳态和保护配合，而不是开关过程精确波形。另外，把线路模型从分布式参数改成集中参数，用RLC串联分支代替PI Section Line，速度能提升数倍。

5.5 三相交流侧电流波形严重畸变，谐波含量异常

现象：交流侧电流波形不是预期的近似正弦波，而是尖峰密集的畸波，FFT结果里谐波高达20%以上，与论文的仿真图完全不符。

原因：整流变压器二次侧绕组联结错误，或者整流桥触发逻辑用错了PWM信号。牵引供电整流是不控整流，二极管应该自然换相，如果你错误接入了触发脉冲，等于给二极管强行加上了开关信号，导致电流断续。

解决：检查Universal Bridge的“Power Electronic Device”是否确实选为“Diode”，而不是“IGBT/Diode”。IGBT混合模型默认允许外部触发信号，未接线时可能处在错误的门控状态。另一个原因是交流侧电源内阻太小且无滤波电抗器，整流器换相时产生尖峰。在交流侧每个相线上串入一个小电感（约0.1mH）来模拟实际变压器漏抗，波形会立即平滑下来。

6. 进阶：把仿真结果变成论文复现证据的三种做法

到这里模型已经能跑，但复现的关键是让仿真曲线和论文结果“对上话”。我的习惯是先做“参数敏感性验证”，再做“故障特征对比表”，最后生成“三图合一”的复合波形。

先做一个简单的脚本，批量修改平波电抗器电感值，记录直流电压纹波系数和短路电流峰值，生成对比表。这个表直接放进论文里，能说明模型对参数的响应符合物理规律，比单张波形更有说服力。脚本如下：

% 批量扫描平波电抗器电感，记录关键指标 L_values = [1e-3, 2e-3, 5e-3, 8e-3]; results = zeros(length(L_values), 3); for i = 1:length(L_values) set_param('traction_sim/L_smooth', 'L', num2str(L_values(i))); sim('traction_sim'); results(i, 1) = L_values(i); results(i, 2) = Vdc_ripple_ratio; % 从工作区读取纹波系数 results(i, 3) = Isc_peak; % 短路电流峰值 end disp(table(results(:,1), results(:,2), results(:,3), ... 'VariableNames', {'L_H', 'ripple_ratio', 'Isc_peak'}));

第二种做法是验证相对误差。找到论文里给出的稳态直流电压、短路电流峰值或谐波含量数据，把仿真值算成相对误差，控制在±10%以内即可认为模型可信。如果超过这个范围，优先怀疑模型参数与论文不一致，特别是变压器短路阻抗、线路单位长度参数和列车负载特性曲线。

第三种做法是把仿真波形、解析计算波形和论文截图（如果你有版权许可）放在同一张图里对比。没有论文原图时，至少要把仿真波形和解析公式曲线对比，比如空载直流电压、短路稳态电流都可以手算。这样审查老师会认为你的模型经过了双向验证，不是“黑匣子”跑出来的。

我做这类复现的一个深刻教训是：仿真波形“像不像”论文是次要的，“为什么像”才是重点。某次我在整流器故障分析里得到一条和论文几乎一致的电压波形，后来仔细检查才发现是因为我在整流桥输出端并了一个大电容，恰好把这个故障特征“滤”得和论文一样了。这个电容在真实系统里不存在，纯粹是我为了波形好看加上去的。后来我删掉电容，波形露出了原本的尖峰和振荡，反而更像真实故障。现在我的习惯是：每加一个元件前先问自己，真实系统里有没有这个东西，如果没有，那就说明模型建错了，而不是论文画错了。

希望这篇笔记能让你在MATLAB/Simulink牵引供电建模与故障分析这条路上少一点玄学，多一点可复现的确定性。如果你跑通了自己的第一个模型，再去对比论文的波形、保护动作时序和故障特征表，你会发现那些看似复杂的图，其实都是同一套原理在不同参数下的投影，祝搭建顺利。

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