地铁牵引供电系统Simulink仿真：24脉波整流与故障复现全解析

简介：面向轨道交通电气系统设计、维护与研究人员的Matlab/Simulink地铁牵引供电系统建模与仿真复现资料，聚焦24脉波整流、接触网PI模型、SVPWM控制及故障仿真等关键技术。资源基于成都地铁实际参数，构建了含外部电源、整流变电站、直流母线、接触网与钢轨大地模型的完整供电系统，并搭建牵引传动系统及牵引负荷模型，可复现空载、负载与故障工况仿真，并与现场故障录波对比验证模型有效性。

内容预览部分提供了详细可运行代码与解释，涵盖24脉波整流器建模（两组相位差15°变压器实现）、接触网分布参数特征阻抗计算、空间矢量PWM逆变控制及系统级仿真框架，方便读者从零搭建仿真模型并理解各环节原理。整个资源为1个PDF文档，压缩包约671KB，文件小巧但思路完整，适合用于系统设计优化、运行维护和故障分析。

目前已有244人学习，尤其适合需要快速掌握地铁牵引供电系统数字化建模方法的工程师与研究人员；其分层验证体系和多时间尺度仿真策略可推广至其他城市轨道交通线路，为后续数字孪生与AI故障预测研究奠定基础。

1. 地铁牵引供电系统仿真不只是画拓扑：24脉波整流与故障复现的关键拆解

做牵引供电仿真的工程师都知道，Simulink 里把系统框图画出来不难，测点波形和现场录波对得上才是真门槛。这份论文复现资源用 Matlab/Simulink 搭了完整的地铁直流牵引供电系统，涵盖 24 脉波整流变电站、直流母线、接触网 PI 模型、钢轨大地回流网络和牵引传动系统，并给出了短路故障注入、录波对比和误差分析的完整代码路径。它解决的核心问题是：理论模型怎么落到可运行的 Simulink 框图，整流脉波数怎么定、接触网分布参数怎么设、SVPWM 控制和故障仿真参数怎么调。适合做牵引供电设计、整流机组选型和故障分析的工程师，也适合在校研究生直接拿它当课题起点。

2. 系统架构与模型选型：为什么24脉波整流和PI接触网是这套仿真的核心

2.1 从外部电源到钢轨大地的子系统拆解

论文模型的拓扑关系是：外部电源经过电缆接入整流变电站，变电站在 Simulink 里输出直流母线电压，直流母线经馈电线连接接触网，列车通过受电弓从接触网取流，电流经钢轨和大地回流至变电站负极。整个回路必须全部建模，不能只搭供电侧而忽略回流网络，否则短路电流和电位分布都算不准。

原文把系统拆成了五个明确子系统，每个子系统都有各自的简化逻辑，这也是复现时最容易直接照搬的部分：

外部电源用戴维南等效电路表示，是刻意简化。地铁供电系统的短路容量主要由电网侧决定，线路长度和变压器阻抗才是关键，戴维南等效的电压源串阻抗完全能覆盖这个特性。

整流变电站是整个模型里最影响波形精度的模块。论文用两组移相变压器配合两组三相全桥整流，输出合成直流电压，理论脉波数为 24。这个设计的选型理由是：在地铁直流牵引中，24 脉波整流是谐波性能和设备成本之间的主流平衡点，比 6 脉波和 12 脉波的直流纹波更小，比 48 脉波的变压器结构简单。直流母线用 RLC 集中参数建模，主要反映母线电缆和滤波电抗的储能特性。接触网采用 PI 型链式电路，用单位长度电阻、电感、电容和电导四个参数描述分布特性。钢轨大地采用四层土壤模型，用 Carson 公式计算地回路阻抗。

这套分层的价值在于每个子系统的参数可以独立标定，调试时可以分块验证。我复现时通常先把整流变电站单独跑一遍，确认空载直流电压在规定范围内，再接接触网和负载，避免问题交叉。实际模型里子系统之间用电气节点连接，Simulink 的 Specialized Power Systems 库提供了现成的连接方式，不需要自己写接口。

2.2 24脉波整流：移相15度背后的谐波与电压逻辑

24 脉波整流的核心是两组三相整流桥的输出在直流侧并联，而两组桥的交流侧电压相位错开 15°。代码里对应的关键行是：

% 24脉波整流器建模 function [output_voltage] = rectifier_24pulse(input_voltage) % 参数设置 V_primary = 35e3; % 35kV交流输入 V_secondary1 = 1180; % 两组次级绕组电压 V_secondary2 = 1180; phase_shift = 15; % 移相角度(度) % 创建三相变压器组(带移相) transformer1 = threePhaseTransformer(... 'PrimaryVoltage', V_primary, ... 'SecondaryVoltage', V_secondary1, ... 'PhaseShift', 0); transformer2 = threePhaseTransformer(... 'PrimaryVoltage', V_primary, ... 'SecondaryVoltage', V_secondary2, ... 'PhaseShift', phase_shift); % 整流桥配置(4个三相全桥) rectifier1 = universalBridge('NumDiodes', 6, 'Configuration', 'Three-phase bridge'); rectifier2 = universalBridge('NumDiodes', 6, 'Configuration', 'Three-phase bridge'); % 并联输出直流 output_voltage = (rectifier1.output + rectifier2.output)/2; end

逻辑说明：变压器 1 的相移为 0°，变压器 2 的相移为 15°，两组次级绕组电压同为 1180V。每台变压器接一组三相全桥二极管整流器，输出通过并联取平均值合成直流电压。两组桥的换相时刻彼此错开，使得合成直流电压的纹波频率提高到工频的 24 倍，幅值显著降低。

参数说明：V_primary=35kV 是地铁中压环网常用电压等级；V_secondary=1180V 决定了空载直流电压的理论值，按三相全桥整流关系，空载直流电压约为线电压有效值的 1.35 倍，即约 1593V，这和论文案例中 1650V 的实测值接近，差值留给线路压降和变压器漏抗解释；phase_shift=15 是本代码实现 24 脉波的关键参数，如果改成 0，两组桥的换相时刻完全同步，等效为一个 12 脉波整流器，直流纹波和网侧谐波都显著增大。

在 Simulink 里复现时，移相变压器不是用普通三相双绕组变压器直接改相位，而是用专门的移相变压器模块，设置二次侧绕组连接组别，或者通过绕组延边方式实现移相。论文代码里的 threePhaseTransformer 是封装函数，实际落地时对应 Simulink 的 Three-Phase Transformer (Two Windings) 或 Specialized Power Systems 库里的对应元件。我一般会在二次侧加一组理想电压源做验证，把输出空载直流电压和 1.35 倍关系核对一遍再往下接，这个习惯帮我避开了不少变压器绕组接反的低级错误。

2.3 接触网PI模型与四层大地回流模型的参数边界

接触网的分布参数在仿真里不是直接拖一个 RLC 元件，而是按单位长度参数乘以线路长度，然后用 PI 型链式电路近似。论文代码给出了四个基础参数：

% 接触网分布参数模型 function [Z_contact] = contact_line_model(freq) % 某城市地铁线路参数 R_per_km = 0.02; % 单位电阻(Ω/km) L_per_km = 1.2e-3; % 单位电感(H/km) C_per_km = 12e-9; % 单位电容(F/km) G_per_km = 1e-6; % 单位电导(S/km) % 计算特征阻抗 Z_series = R_per_km + 1j*2*pi*freq*L_per_km; Y_shunt = G_per_km + 1j*2*pi*freq*C_per_km; Z_contact = sqrt(Z_series/Y_shunt); end

逻辑说明：这段代码计算的是接触网的串联阻抗和并联导纳，再求特征阻抗。Z_series 由单位电阻和感抗组成，Y_shunt 由单位电导和容抗组成。特征阻抗用于后续判断线路的阻抗匹配和反射特性。

参数说明：R_per_km=0.02Ω/km 对应接触网导线的直流电阻水平，这是一个偏小的值，适合大截面铜导线；L_per_km=1.2mH/km 包含接触线和回流回路的等效电感，这个数量级在直流牵引网中常见；C_per_km=12nF/km 反映接触线对地的分布电容；G_per_km=1μS/km 代表绝缘泄漏。

这里有一个容易忽略的坑：接触网在直流牵引系统里传输的是直流电流，但故障暂态过程和高频振荡分析必须考虑分布参数。论文把接触网模型独立封装并计算特征阻抗，实际 Simulink 实现时是把单位参数乘以线路长度，再用多个 PI 小节串联。线路长度超过 2 公里时，建议至少拆成 4 到 5 个 PI 小节，否则电压分布和短路电流波形会失真。

钢轨大地回流网络的建模原文用了四层土壤模型：

% 钢轨-大地阻抗计算 function [Z_rail] = rail_impedance_model(freq) % 某城市地铁线路参数 soil_layers = [0.5 1.0 2.0 Inf]; % 各层深度(m) resistivity = [100 80 50 30]; % 电阻率(Ω·m) rail_params = struct('R', 0.02, 'L', 1.5e-3); % Carson公式计算大地回路阻抗 Z_earth = carson_equation(freq, soil_layers, resistivity); % 总阻抗计算 Z_rail = rail_params.R + 1j*2*pi*freq*rail_params.L + Z_earth; end

逻辑说明：钢轨自身用集中电阻和电感描述，大地回流路径用 Carson 公式计算。Carson 公式适用于以大地为回路的线路阻抗计算，在地铁场景中用于评估钢轨对地泄漏和回流路径的地阻抗。

参数说明：soil_layers 的四层深度分别对应地表土壤、浅层回填层、深层土壤和无限深基层；resistivity 从 100Ω·m 递减到 30Ω·m，反映了土壤电阻率随深度的变化。这个参数在短路和人身安全分析中直接影响地电位升和杂散电流分布。

这个模型复现时的难点是 Carson 公式在 Simulink 里没有现成模块，我一般会把它写成一个 MATLAB Function 模块，输入频率、层深和电阻率，输出地阻抗的实部和虚部，然后在 Simulink 里用受控电压源或阻抗元件接入回流回路。

3. 在Simulink里搭出可运行模型：关键模块实现、参数表与联合仿真设置

3.1 整流变电站模块：从封装函数到元件落位

论文代码里的 rectifier_24pulse 是功能函数，到了 Simulink 里要落成具体元件。推荐的做法是新建一个子系统，内部包含两个移相变压器、两个三相二极管整流桥、直流电抗器和输出滤波器。

子系统输入输出定义：输入是三相 35kV 交流电压源，输出是直流母线电压和整流侧三相电流。我用 Simscape Electrical Specialized Power Systems 库里的 Three-Phase Source 和 Three-Phase Transformer (Two Windings) 来实现。变压器参数设置时注意把绕组连接组设为 Dyn11，二次侧额定电压设为 1180V，移相角按 0° 和 15° 分别配置。

直流母线侧接一个平波电抗器，电感取 5mH 左右，然后是直流母线电容。电容值按负载功率的 1/10 到 1/5 设计。原文给了 LCL 滤波器的计算方法：

% LCL滤波器设计 function [filter] = design_lcl_filter() % 设计参数 P_nom = 1500e3; % 额定功率 V_dc = 1500; % 直流电压 f_sw = 2e3; % 开关频率 % 计算参数 L1 = 0.1 * V_dc^2 / (P_nom * 2*pi*f_sw); C_f = 0.05 * P_nom / (V_dc^2 * 2*pi*f_sw); L2 = 0.5 * L1; end

逻辑说明：L1、C_f、L2 三个参数构成了一个三阶低通滤波器，截止频率设计在开关频率的 1/10 左右，用于衰减 IGBT 开关产生的谐波。额定功率 1500kW、直流电压 1500V、开关频率 2kHz 与论文的牵引系统规模一致。

参数说明：L1 的表达式是按电感上的纹波电流不超过额定电流的 10% 来反推的；C_f 的表达式是按无功功率不超过额定功率的 5% 来约束的；L2 取 L1 的一半是为了在满足滤波效果的前提下减小磁性元件体积。

这个滤波器在 Simulink 里很好验证：用 LTI Viewer 或直接跑一次满载阶跃，看直流母线电压纹波和电流 THD。我调试时习惯把 L1 初始值放大两倍再逐步缩小，因为 L1 偏大会牺牲动态响应，偏小了纹波电流超标，需要找一个平衡点。

3.2 接触网与回流网络：把单位参数换算成PI小节

接触网在 Simulink 里的落法是：先计算线路总长度，比如单线 5 公里，然后分成 5 段，每段 1 公里。每段用串联阻抗和并联导纳组成 PI 结构。具体数值计算：

% 接触网PI小节参数计算 total_length = 5; % 接触网总长度(km) num_sections = 5; % PI小节数量 section_length = total_length / num_sections; % 每段长度(km) R_section = 0.02 * section_length; % 每段串联电阻(Ω) L_section = 1.2e-3 * section_length; % 每段串联电感(H) C_section = 12e-9 * section_length; % 每段并联电容(F) G_section = 1e-6 * section_length; % 每段并联电导(S)

逻辑说明：这个换算的本质是把单位长度参数乘以 PI 小节的长度，得到该段的集中参数。串联阻抗放在该段的支路里，并联导纳放在该段两端对地。

参数说明：section_length 是 PI 小节长度。对于 5 公里的接触网，分成 5 段意味着每一段的谐振频率远高于关注的暂态频率区间，仿真结果就能收敛到稳定值。实际项目里我会先按 5 段跑通，再用分段数翻倍做一次对比，如果电压分布误差小于 1%，就不需要继续加密。

回流网络这边，钢轨自身阻抗和大地回路阻抗串在回流支路里。我在 Simulink 里用 Series RLC Branch 封装，把钢轨电阻和电感固定，地阻抗通过 MATLAB Function 模块实时计算。连接关系是钢轨支路一端接列车受电弓回流点，另一端接变电站负极母线，大地支路从钢轨对地节点接入无穷远地电位参考点。

3.3 牵引传动系统：SVPWM控制与FOC的落法

牵引传动系统由前端 LCL 滤波器、IGBT 三相逆变桥、永磁同步电机或异步电机模型组成。论文给了 SVPWM 和 FOC 两套控制算法。SVPWM 的核心是扇区判断和占空比计算：

% 空间矢量PWM逆变控制 function [PWM] = svpwm_control(V_ref, V_dc, f_sw) % 输入参数： % V_ref - 参考电压矢量 % V_dc - 直流母线电压 % f_sw - 开关频率 % 坐标系变换(Clark变换) V_alpha = V_ref(1); V_beta = V_ref(2); % 扇区判断 theta = atan2(V_beta, V_alpha); sector = floor(theta/(pi/3)) + 1; % 计算占空比 T1 = sqrt(3)*abs(V_ref)*sin(pi/3 - mod(theta, pi/3))/V_dc; T2 = sqrt(3)*abs(V_ref)*sin(mod(theta, pi/3))/V_dc; T0 = 1 - T1 - T2; % 生成PWM波形 PWM = zeros(1,6); switch sector case 1 PWM = [T1+T2+T0/2, T2+T0/2, T0/2, T0/2, T0/2, T1+T0/2]; % 其他扇区类似处理... end end

逻辑说明：SVPWM 的本质是把参考电压矢量在 α-β 平面分解，判断它落在哪个扇区，再用该扇区相邻两个基本电压矢量来合成。T1 和 T2 是两个基本矢量的作用时间，T0 是零矢量作用时间。

参数说明：V_ref 通常由电流环或速度环输出；V_dc 是直流母线电压，直接决定了最大可输出电压幅值，代码里 T1、T2 计算公式的分母是 V_dc，反映的是调制比的概念；f_sw 影响 PWM 的载波周期，最终影响电流纹波和开关损耗。

这里有个容易出错的地方：扇区判断用 atan2 得到的角度范围是 -pi 到 pi，而 sector 计算只考虑了 0 到 2pi 的情况。当 theta 为负值时，floor(theta/(pi/3))+1 会给出 0 或负值扇区。我复现时的处理办法是先把 theta 归一化到 0 到 2pi 范围再做扇区判断，后面避坑章节会展开讲。

FOC 控制代码：

% 牵引电机矢量控制系统实现 function [PWM] = foc_control(id_ref, iq_ref, theta_elec) % 参数设置 Ld = 8e-3; % d轴电感(H) Lq = 12e-3; % q轴电感(H) Rs = 0.2; % 定子电阻(Ω) % 电流环PI控制器 Kp_d = 0.5; Ki_d = 10; Kp_q = 0.8; Ki_q = 15; % 坐标变换 id = ia*cos(theta_elec) + ib*sin(theta_elec); iq = -ia*sin(theta_elec) + ib*cos(theta_elec); % PI控制输出 Vd = Kp_d*(id_ref - id) + Ki_d*integral(id_ref - id); Vq = Kp_q*(iq_ref - iq) + Ki_q*integral(iq_ref - iq); % 反Park变换生成PWM PWM = svpwm_generator(Vd, Vq, theta_elec); end

逻辑说明：FOC 把三相定子电流变换到 d-q 旋转坐标系，d 轴电流控制磁通，q 轴电流控制转矩，两个电流环 PI 控制器输出 d-q 轴电压，再通过反 Park 变换生成 SVPWM 的输入参考矢量。

参数说明：Ld=8mH、Lq=12mH 是典型的凸极电机参数，d 轴电感小于 q 轴电感，这是内置式永磁电机的特征。Kp_d=0.5 和 Ki_d=10 是 d 轴电流环的 PI 参数，Kp_q=0.8 和 Ki_q=15 是 q 轴电流环的 PI 参数。这组参数在仿真里可以实现约 100ms 以内的电流阶跃响应。

调试这种方法时要注意：FOC 依赖转子位置角度 theta_elec，Simulink 里要用电机模块输出的机械角度乘以极对数得到电角度，然后把电角度送给 Park 变换和反 Park 变换。这个角度同步一旦出错，转矩方向会反掉，电机表现为失速。如果仿真里电机转速上不去，先检查电角度极性，再查 PI 输出限幅。

3.4 主仿真脚本：负荷曲线、仿真时长与数据集管理

主仿真脚本把各个子系统串起来。原文给的框架：

% 主仿真脚本 function main_simulation() % 1. 初始化参数 sim_time = 10; % 仿真时长(s) load_profile = [0 0; 2 1500; 5 3000; 8 1500; 10 0]; % 负荷曲线(kW) % 2. 构建子系统 power_supply = build_power_supply(); % 供电系统 contact_net = build_contact_net(); % 接触网 train_model = build_train_model(); % 车辆模型 % 3. 运行仿真 sim_out = sim('metro_system.slx', 'StopTime', num2str(sim_time)); % 4. 结果分析 plot_results(sim_out, load_profile); end

逻辑说明：主脚本的作用是统一管理仿真时长、负荷曲线和子系统构建。power_supply 对应整流变电站，contact_net 对应接触网与回流网络，train_model 对应牵引传动系统。sim 命令调用 Simulink 模型文件 metro_system.slx。

参数说明：sim_time=10 秒覆盖了列车从启动、巡航到制动的完整过程。load_profile 是一个两列矩阵，第一列是时间点，第二列是负荷功率。这个曲线在 0-2 秒空载，2-5 秒负荷 1500kW，5-8 秒负荷 3000kW，8-10 秒回落到 1500kW。这个动态负荷切换是仿真最有价值的部分，可以观察直流母线电压随负载的波动。

我在实际复现时，会在 load_profile 基础上再加一个列车运行阻力模型：

% 列车运行阻力与负荷动态特性 load_profile = struct(... 'Acceleration', [0 1.0; 10 0.8; 20 0], % 加速度曲线(m/s²) 'MaxSpeed', 80, % 最高速度(km/h) 'TrainMass', 280000, % 列车质量(kg) 'Resistance', @(v) 5.2 + 0.045*v + 0.00086*v^2); % 运行阻力公式(N)

逻辑说明：这个结构体把列车运行特性封装起来。加速度曲线分为三个阶段，最高速度 80km/h，列车质量 280 吨。阻力公式是 Davis 方程的简化形式，包含滚动阻力、线性阻力和空气阻力三项。

参数说明：TrainMass=280000kg 对应一列 6 节编组的 B 型地铁列车满载质量。阻力公式的三个系数 5.2、0.045、0.00086 分别对应滚动阻力、线性阻力和空气阻力系数，速度 v 单位是 km/h，阻力结果单位是 N。

把这个阻力和牵引力一起折算成功率，再映射到 load_profile，比直接用固定负荷曲线更符合实际。这也是论文验证时能把电压误差控制在 2% 以内的前提。我一般会在主脚本里预留一个 debug 开关，控制是否输出中间变量到工作区，方便事后排查。

4. 故障工况仿真与验证：短路注入、录波对比和误差边界

4.1 短路故障注入方法：从set_param到故障电阻的选取

故障仿真的核心是在指定时刻给系统注入一个低阻抗路径。论文代码提供了两种方式：一种是直接修改模型中故障模块的参数，另一种是动态插入短路支路。

第一种方式：

% 短路故障仿真 function fault_simulation() % 设置故障参数 fault_time = 5; % 故障发生时间(s) fault_duration = 0.1;% 持续时间(s) fault_R = 0.01; % 短路电阻(Ω) % 修改Simulink模型参数 set_param('metro_system/Fault', 'Time', num2str(fault_time)); set_param('metro_system/Fault', 'Duration', num2str(fault_duration)); % 运行仿真 sim_out = sim('metro_system.slx'); % 绘制故障电流波形 plot(sim_out.tout, sim_out.fault_current); title('短路故障电流特性'); end

逻辑说明：set_param 命令在仿真开始前修改故障模块的触发时间和持续时间。故障模块是一个受控开关并联一个 0.01Ω 电阻，开关在 5 秒闭合，0.1 秒后断开。

参数说明：fault_time=5s 意味着系统在稳定运行后进入故障状态；fault_duration=0.1s 模拟保护装置在故障后 100ms 内动作切除故障；fault_R=0.01Ω 是金属性短路过渡电阻的典型值，也可以改成 0.1Ω 或 0.5Ω 来模拟经过渡电阻的短路。

第二种方式是在 Simulink 模型里动态插入短路支路：

% 极间短路故障注入 function inject_pole_fault(modelname, t_fault) % 获取系统所有支路 branches = find_system(modelname, 'Type', 'Branch'); % 在故障点插入短路电阻(10mΩ) add_block('powerlib/Elements/Series RLC Branch', ... [modelname '/Fault_Path'], ... 'Resistance', '0.01', ... 'Inductance', '0', ... 'Capacitance', 'inf'); % 配置故障时间控制 set_param([modelname '/Fault_Control'], ... 'EventTime', num2str(t_fault)); end

逻辑说明：add_block 动态在模型里添加一个短路支路，电感设 0、电容设 inf，相当于一个纯阻性故障路径，然后通过故障控制模块的 EventTime 参数指定动作时刻。

参数说明：Resistance=0.01 对应 10mΩ 的金属性短路电阻；Capacitance=inf 表示该支路没有电容效应，只在故障瞬间提供一个直流短路通道。这个做法的好处是每跑一次仿真不用重新打开模型手动改参数，适合批量扫故障类型和过渡电阻。

故障位置的选择也要注意。论文验证场景里短路点设在接触网中部，这是最不利的位置之一，因为距离整流变电站和列车都有一定距离，故障电流要经过线路阻抗，波形特征和变电站出口短路完全不同。我复现时会同时设置三个故障备选位置：变电站出口、接触网中点、列车受电弓处，分别观察电压跌落和电流峰值的变化。

4.2 典型工况对比：空载、负载与故障的误差分析

论文验证的核心是三类工况：空载、满载和短路故障。对比表：

三组数据的误差特性说明了模型的适用边界。空载电压误差 0.9% 来自变压器漏抗和二极管导通压降的模型简化，这个精度已经足够用于设备选型。满载压降误差 1.1% 说明接触网电阻和钢轨回流阻抗的参数取值得到了很好的校准。短路电流误差 3.2% 在可接受范围内，但对于断路器开断能力校核，建议把仿真值加上 1.2 倍裕量再用于工程判断。

注意：短路电流误差对过渡电阻参数非常敏感。用 0.01Ω 仿真得到 12.1kA，把电阻改成 0.05Ω，短路电流会掉到 10kA 以下。所以现场对比时分清录波是金属性短路还是经过渡电阻短路，如果不是金属性短路，不要强行调模型参数去匹配录波峰值。

我复现时的经验是先把空载和满载两组数据校准干净，再跑故障工况。空载电压差了就调整流变压器漏抗比例；满载压降差了就调接触网电阻和钢轨回流阻抗；最后再让短路电流对齐。按这个顺序来，三组数据能同时满足误差要求的时间最短。

4.3 三级验证体系与多时间尺度仿真策略

论文提出了三级验证方法论：

每个层级都有独立的合格标准。单元测试针对整流器、滤波器和控制器，比如整流器空载电压要和 1.35 倍的线电压理论值对标。集成测试看子系统联手后的动态性能，比如满载阶跃时直流母线电压的恢复时间。系统验证是把整个模型和现场录波直接对比。

我复现时的习惯是先跑单元测试再做系统联调。整流变电站单独仿真时，把输出空载直流电压记录到工作区，和 1593V 的理论值对比，误差超过 5% 就回头检查变压器参数和整流桥配置。接触网单独测试时，注入一个阶跃电流，观察线路末端的电压响应和理论上的一阶惯性特性是否一致。

多时间尺度仿真策略也是一个有实用价值的点。牵引供电系统的暂态过程跨越几个数量级：IGBT 开关周期是微秒级，保护动作是毫秒级，列车运行工况是秒级，全天运行分析是小时级。一个固定步长仿真无法同时覆盖这些尺度。原文的思路是分层验证，稳态用大步长跑，故障暂态用变步长求解器重点捕捉。我一般会为故障仿真单独设置一个求解器配置，把最大步长限制在 10μs，绝对容差调到 1e-6，这样可以精确捕捉短路电流的峰值和上升率。

5. 复现避坑指南：五个高频问题与排查路径

复现这套模型最折磨人的不是原理理解，而是 Simulink 的工程细节。我前后跑通了三次，每一次踩的坑都不一样，但归纳起来集中在五类问题上。下面按现象、原因、解决的顺序展开，基本覆盖了从模型搭建到故障仿真的完整链路。

5.1 现象：仿真步长与开关频率不匹配导致波形发散

现象：Simulink 模型跑 SVPWM 整流器或逆变器时，波形出现严重的高频振荡，或者仿真直接报 Solver encountered an error。明明参数是按论文设置的，换个电脑就发散。

原因：开关频率 f_sw=2kHz 对应开关周期 500μs，如果固定步长求解器把步长设到 100μs 以上，一个开关周期内只有不到 5 个采样点，IGBT 开关沿附近的电压电流跳变无法被准确捕捉，数值算法就会出现振荡。

解决：把仿真步长降到开关周期的 1/10 到 1/20。对于 2kHz 开关频率，固定步长建议不大于 20μs。如果是变步长求解器，设最大步长也是 20μs，相对容差设 1e-4。每次修改 f_sw 之后，同步检查求解器步长设置。

5.2 现象：移相变压器设置不当导致脉波数不足

现象：直流侧波形还呈现 6 脉波或 12 脉波的纹波特征，说好的 24 脉波完全看不出来，谐波分析结果也不对。

原因：两个变压器的移相角没有正确配置，或者配置了移相角但在 Simulink 里没有生效。最常见的错误是用普通三相变压器模块改了 PhaseShift 参数，但绕组连接组没设对，移相没有实际作用到二次侧线电压上。

解决：复现时用带移相功能的专用变压器模块，明确设置二次侧连接组别为延边三角形或 Dyn11，并检查二次侧线电压的相角。我每次都会在整流变电站输出端接一个 PLL 锁相环，直接测量两组整流桥交流侧线电压的相位差，确认是 15° 而不是 0° 或 30°。

5.3 现象：接触网分布参数没有按链式电路落

现象：仿真出来的直流电压沿线路分布是均匀的，但实际上距离变电站越远电压降越大。对比实测电压分布时偏差很大。

原因：接触网直接拖了一个 RLC 串联支路，没有拆成多个 PI 小节。单位长度参数乘以总长度得到的集中参数模型无法反映线路的分布特性，负载变化时电压分布曲线失真。

解决：把接触网按线路长度分段建模，每段一个 PI 小节。线路 5 公里就至少分成 4 到 5 段，每段长度不超过 1.5 公里。分段数越多分布特性越精确，但仿真速度会下降。我复现时先按 5 段跑通，再根据电压分布精度决定是否加密。

5.4 现象：故障注入后直流侧电流不突变

现象：设置了故障模块，也定义了故障时间和持续时间，但仿真结果里短路电流和正常电流几乎一样。

原因：故障模块的连接可能没有真正接入直流母线回路。在 Simulink 里，故障模块如果只是从元件库拖到画布而没有用线连接，或者连接到了错误的节点，set_param 改参数不会产生任何作用。

解决：故障注入前先建立模型快照，检查故障模块的输入输出端口与直流母线正负极是否正确连接。我用 trace 检查所有连接到故障两端节点的支路，确保故障闭合时电流有通路。另一个快速验证方法是临时把故障电阻改成 1Ω，看直流电流是否立即下降。

5.5 现象：SVPWM扇区判断在边界处出错

现象：PWM 波形在某个角度出现不连续，输出电流过零点附近有畸变。看起来像是调制算法在扇区切换时跳变。

原因：扇区判断用 atan2 函数返回的角度范围是 -pi 到 pi，而占空比计算里的 theta 直接用 mod(theta, pi/3) 处理，当 theta 是负值时，归一化不一致，T1 和 T2 的计算结果错误，PWM 波形在扇区边界处跳变。

解决：先把 theta 规范到 0 到 2pi 范围再计算扇区。具体做法是对 theta 加 2pi 直到落在 0 到 2pi 区间，再做扇区判断和占空比计算。这个改动解决了波形畸变，也让调制比在边界处连续变化。

6. 从复现到扩展：把模型改造成多车运行与数字孪生实验台

复现完整的单列车模型之后，最有价值的扩展是把它改造为多车运行场景。做法是把 train_model 复制为多份，分别用不同的启动时间和负荷曲线驱动，在直流母线上并行接入。子系统的复制不用改内部结构，只需要把每份的 load_profile 时间参数错开。

多车运行仿真的关键在于负荷曲线时间错开。真实线路里不同位置的列车加速、巡航和制动状态不同，如果多列车同时处于加速阶段，直流母线电压会跌得很快。我在一个模拟项目里仿真了 3 列车错峰运行，母线电压最低点出现在两列车同时加速的重叠区间，这个结果直接用于牵引所容量校核。具体做法是把第 2 列车和第 3 列车的启动时间分别延迟 20 秒和 40 秒，然后在同一张图里叠加三条负荷曲线：

% 多车错峰运行仿真参数 sim_time = 60; train1_load = [0 0; 5 1500; 15 3000; 25 1500; 30 0]; train2_load = [20 0; 25 1500; 35 3000; 45 1500; 50 0]; train3_load = [40 0; 45 1500; 55 3000; 60 1500; 65 0]; sim_out = sim('metro_system_multi.slx', 'StopTime', num2str(sim_time)); plot(sim_out.tout, sim_out.dc_bus_voltage);

第二个扩展方向是数字孪生应用。把 Simulink 模型封装成 FMU 或者生成 C 代码，接入实时数据平台，用实测录波数据做在线参数校准。论文里提到的 AI 故障预测也能在这个模型上做数据增强：用故障仿真批量生成短路、过载和整流器失效的样本，训练故障分类模型。我做过一个简单的尝试，用 200 组不同短路电阻和位置的仿真数据训练一个分类器，准确率能到八成以上，但样本覆盖不足时边界情况容易误判。

我个人的习惯做法是：跑一套完整的系统仿真前，先跑通单模块的单元测试，再联调，然后是故障注入、录波对比。从那以后，我每次拿到新的线路参数，都要强制走一遍整流器空载验证、接触网分段数校核、短路电流灵敏度扫描这三步，没有这三步，我不会往报告里放仿真结论。希望帮到你。

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