简介面向轨道交通电气系统设计、维护及相关研究人员的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电路参数也需要注意默认Rs1e5Ω、Csinf这相当于把开关理想化。若要更真实的故障电流需要把Snubber电路的电容设为有限值例如10nF电阻降到10Ω但这会显著增加仿真时长建议只在最终验证阶段开启。3.2 列车牵引/制动特性曲线怎么进模型列车负载的动态特性直接影响直流母线电压的波动范围。地铁列车的运行工况分为牵引、巡航、惰行、制动四象限其中牵引时列车从接触网吸取功率制动时一部分能量通过再生制动回馈电网。Simulink里最直观的建模方式是用受控电压源和受控电流源组合牵引时电流方向为正制动时为负电压始终由直流网侧决定。我把列车模型拆成两部分功率计算和电流注入。功率计算根据列车速度、加速度和质量按简化力学公式算出瞬时牵引力再乘以车速得到机械功率除以效率后得到电功率电流注入则把这个电功率除以直流母线电压得到参考电流送给受控电流源。这个模型能体现恒功率特性——母线电压下降时电流会上升和实际列车特性一致。下面是一个最小实现function I_ref train_current(v, a, m, eta, Vdc) % 列车牵引电流计算函数 % v 车速 m/sa 加速度 m/s^2m 质量 kgeta 效率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时理论上牵引力公式里v0会导致功率为0但实际列车启动时会有静止牵引力所以需要在模型里加一个启动逻辑当v0.1m/s且F0时按恒力矩处理即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模块求导然后输出去找最大值。这个值受直流回路电感影响仿真中如果发现上升率比论文给的大很多去看看平波电抗器参数是否设置正确、接触网线路电感是否被忽略了。电感量每减少1mHdi/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空载直流电压应该约1590V1180×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牵引供电建模与故障分析这条路上少一点玄学多一点可复现的确定性。如果你跑通了自己的第一个模型再去对比论文的波形、保护动作时序和故障特征表你会发现那些看似复杂的图其实都是同一套原理在不同参数下的投影祝搭建顺利。本文还有配套的精品资源点击获取 
长河君 · 2026/10/21 16:11:43 · 阅读 1.2k 地铁牵引供电系统Simulink建模与故障仿真分析
简介面向轨道交通电气系统设计、维护及相关研究人员的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电路参数也需要注意默认Rs1e5Ω、Csinf这相当于把开关理想化。若要更真实的故障电流需要把Snubber电路的电容设为有限值例如10nF电阻降到10Ω但这会显著增加仿真时长建议只在最终验证阶段开启。3.2 列车牵引/制动特性曲线怎么进模型列车负载的动态特性直接影响直流母线电压的波动范围。地铁列车的运行工况分为牵引、巡航、惰行、制动四象限其中牵引时列车从接触网吸取功率制动时一部分能量通过再生制动回馈电网。Simulink里最直观的建模方式是用受控电压源和受控电流源组合牵引时电流方向为正制动时为负电压始终由直流网侧决定。我把列车模型拆成两部分功率计算和电流注入。功率计算根据列车速度、加速度和质量按简化力学公式算出瞬时牵引力再乘以车速得到机械功率除以效率后得到电功率电流注入则把这个电功率除以直流母线电压得到参考电流送给受控电流源。这个模型能体现恒功率特性——母线电压下降时电流会上升和实际列车特性一致。下面是一个最小实现function I_ref train_current(v, a, m, eta, Vdc) % 列车牵引电流计算函数 % v 车速 m/sa 加速度 m/s^2m 质量 kgeta 效率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时理论上牵引力公式里v0会导致功率为0但实际列车启动时会有静止牵引力所以需要在模型里加一个启动逻辑当v0.1m/s且F0时按恒力矩处理即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模块求导然后输出去找最大值。这个值受直流回路电感影响仿真中如果发现上升率比论文给的大很多去看看平波电抗器参数是否设置正确、接触网线路电感是否被忽略了。电感量每减少1mHdi/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空载直流电压应该约1590V1180×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牵引供电建模与故障分析这条路上少一点玄学多一点可复现的确定性。如果你跑通了自己的第一个模型再去对比论文的波形、保护动作时序和故障特征表你会发现那些看似复杂的图其实都是同一套原理在不同参数下的投影祝搭建顺利。本文还有配套的精品资源点击获取