最近在复现一篇SCI一区论文里基于SLSPC系列的高阶PT-WPT无线电能传输系统Simulink仿真这部分前前后后花了接近三周。老实说刚拿到标题的时候一度看不明白SLSPC到底指什么补偿网络PT-WPT跟普通的感应耦合无线电能传输系统差别在哪为什么一定要用高阶拓扑。把这些概念一个个啃透、再把模型一段段搭起来之后才意识到这套系统真正有意思的地方。这篇就把从项目拆解到Simulink搭建、参数计算、调参和报错排查的完整过程记录下来给准备复现SCI电路论文或者做无线电能传输研究的同行做个参考。整个过程最适合三类读者一是正在做WPT方向的研究生想快速把PT对称原理变成可运行的仿真模型二是做高频电力电子变换器的工程师想了解高阶补偿网络怎么建模三是刚接触Simulink仿真的同学想完整地看一个高级电路仿真项目是如何从零组织起来的。注意下文中的参数设计和实现细节是基于“常见PT-WPT系统和SLSPC高阶补偿网络”的通用工程实践补全的不是对某篇具体论文的逐字复刻。复现任何论文前一定要先对照原文的电路拓扑和参数表把方向和框架确定后再动手改参数。1. 项目到底在做什么SLSPC与PT-WPT的核心逻辑1.1 PT-WPT传输不再是“调准”而是“锁定”普通磁耦合WPT系统大家都知道发射线圈和接收线圈通过感应耦合传递能量系统工作点是靠发射端和接收端同时谐振在一个固定频率上实现的。问题也出在这里耦合系数一变、负载一变、间距一变谐振状态就会跑偏严重时出现“频率分裂”输出功率明显掉下来。传统解决办法是加阻抗匹配网络、加频率跟踪PLL控制起来挺麻烦。PT-WPT这个名字里的PT指的是“Parity-Time对称”也就是“宇称-时间对称”。这个概念最早来自量子力学中的非厄米系统后来被引入到电路和WPT领域。放到无线电能传输系统里看PT对称的核心做法是在发射端构造一个等效负电阻让增益和接收端的损耗平衡。这样一来系统的工作频率就不再依赖两端精确调谐而是会在满足PT条件时自动锁定到某个实数特征频率上。我用大白话解释一下普通系统像两个人手拉手跑步必须步调一致才能一起跑PT-WPT则像是给领跑者加了节拍器只要一方能量给得够整个系统会自动合拍。这个特性在给电动汽车、植入医疗设备、水下设备无线供电时非常有用因为线圈距离和负载经常变传统系统很难保证稳定输出。PT-WPT带来的直接好处有三个一是抑制频率分裂过耦合区域不再出现双峰二是对耦合系数变化不敏感线圈对准差一点也能维持效率三是系统具有“恒功率”输出的倾向负载在一定范围内变化时输出先保持稳定这正好匹配很多无线充电场景。1.2 SLSPC系列高阶拓扑的定位SLSPC这个缩写不同论文里可能也有不同定义。按我这套复现用的常见结构理解S代表发射端串联电感L是谐振电感SPC是串联和并联组合电容。把这些元件按多级结构连接起来就形成一个比SS、LCC更复杂的“高阶补偿网络”。为什么要用高阶网络两层原因。第一层是滤波。二阶系统的LC只在一个频率点附近呈现低阻抗受到逆变器谐波影响后线圈电流中会混入大量高频分量。高阶网络相当于在路上加了多道“滤网”谐波抑制能力明显更强。第二层是让系统在宽负载范围内维持接近零相位角ZPA。补偿网络设计得当时逆变器输出电流和电压几乎同相开关管的无功负担小软开关也更容易实现。我把SLSPC系列理解为一系列可扩展的拓扑族二阶是最基础的串联补偿三阶加入一个并联谐振支路四阶再增加一级串联电感电容。复现时不需要从零推导网络函数关键是把每一级的固有谐振频率对齐到系统工作频率然后通过交流阻抗分析校验总输入阻抗的相频特性。1.3 Simulink在复现工作中的不可替代性复现论文里的电能传输波形用数学软件做矩阵方程能算出效率但看不到开关过程、谐波和瞬态。Simulink加上Simscape Electrical的优势在于多物理域建模电源、开关管、变压器、补偿电容、整流桥、负载可以全部放同一张图里并且能把开关损耗、导通内阻、磁饱和这些非理想因素也加进去。另一个实用点是参数化管理。把系统频率、线圈自感、补偿电容、耦合系数全部定义成MATLAB工作区变量Simulink模型里直接引用。改一次参数所有模块同步更新。配合信号查看器和脚本绘图可以快速画出频率扫描曲线、负载扫描曲线和FFT频谱这对验证PT-WPT的恒功率特性特别重要。2. 参数设计与架构拆解2.1 系统整体架构一个完整的PT-WPT无线电能传输系统仿真模型从电源侧到负载侧可以分成六段直流电源提供系统能量电压大小决定功率等级高频逆变器把直流变成高频交流通常用全桥或半桥结构发射端补偿网络SLSPC高阶结构负责阻抗匹配和谐波抑制线圈耦合环节发射线圈和接收线圈通过互感传递能量接收端补偿网络与发射端对称或互补实现谐振和阻抗变换整流滤波与负载把接收到的交流电变成直流供给负载。PT-WPT的特殊之处在逆变器和耦合环节之间逆变器不只是输出固定PWM而是在控制器的配合下“模拟”出发射端等效负电阻。负载侧线圈和补偿网络呈现正电阻损耗发射端只要提供大小相等、方向相反的等效负电阻增益整个系统的谐振模式就从共轭频率变成实数频率频率分裂随之消失。2.2 频率、电感、补偿电容的计算参数计算是整个仿真项目里最容易被忽略的一步但也是决定成败的一步。我复现时先从系统工作频率出发向后推元件值第一确定工作频率。工业无线充电常用85kHz高频实验室常用6.78MHz或13.56MHz。为了兼顾Simulink仿真速度和开关器件模型精确度我选了85kHz。第二确定线圈自感。自感值由线圈尺寸、匝数、磁芯形状决定。仿真里可以先给定值比如发射和接收线圈自感都取120微亨。第三计算补偿电容。如果使用最基本的串联谐振谐振角频率满足ω0 2πf0C 1 / (ω0² × L)代入f0 85kHz、L 120μHω0 2π × 85000 ≈ 534070 rad/sC 1 / (534070² × 120e-6) ≈ 2.92e-8 F也就是大约29.2nF。这个值听起来不大实际电路中还要把线圈的分布电容、PCB走线杂散电容考虑进去仿真初算时可以先用理想值后面做容差分析时再偏离。第四计算耦合系数和互感。两线圈之间的互感M与自感的关系是M k × √(L1 × L2)如果L1和L2都等于120μHk取0.2M就是0.2 × 120μH 24μH。这个互感会随线圈距离剧烈变化我在做PT-WPT验证时会把k设为可扫描变量从0.05一直扫到0.4。高阶网络的参数不能一个个独立算必须用整体阻抗匹配的视角来解。我的做法是先用基础公式估算串联谐振电容再在各谐振支路并联电容处做一次“频率对齐”最后看仿真得到的输入电压电流相位角。如果相位角不为零微调并联支路的电容值直到逆变器输出接近单位功率因数。2.3 PT对称条件与等效负电阻设计PT-WPT的灵魂在“增益等于损耗”这个平衡条件上。把电路简化成两侧谐振器来看发射端有等效负电阻-Rg接收端包括线圈内阻、整流等效输入电阻在内的总正电阻为Rr。系统存在实数特征频率的条件是Rg ≥ Rr等号是临界PT对称点大于号则进入PT对称区。实际仿真时如果Rg远大于Rr系统失去PT特性会重新出现频率分裂如果Rg小于Rr系统又可能进入欠增益状态输出功率不足。所以负电阻不是越大越好而是应该严格追着总损耗走。在Simulink里等效负电阻有两种实现方式。第一种是理想受控源实现用电流传感器把发射端电流采出来乘以负增益再输入到受控电压源形成v -Rg × i。这个办法简单适合验证PT理论但容易出现数值振荡必须配合一阶低通滤波器使用。第二种是自激振荡逆变器实现通过正反馈回路让逆变器的开关频率跟着谐振网络自动走逆变器前端等效为负阻。这个办法更接近论文实际电路但建模复杂。我在这套复现模型里先用理想负阻模型把PT-WPT的整体特性验证通过再把理想负阻替换成全桥逆变器加反馈控制这样每一步的变量都控制在最小范围内出问题时很容易定位。3. Simulink建模一步步搭起来3.1 模型顶层结构打开Simulink后不要急着拖模块。先把顶层结构画出来直流源、全桥逆变器、发射补偿、发射线圈、接收线圈、接收补偿、整流桥、负载、测量模块。我的习惯是每个功能块用Subsystem包装起来Subsystem之间用物理信号线或Electrical信号连接。这样后面出问题时可以快速定位到某个子系统不用在一张图上找半天。Simscape Electrical中的电气网络需要参考地所以在直流电源负端和整流桥输出负端都要接“Electrical Reference”模块。如果漏接地仿真会报“unconnected port”或者电压显示成浮地状态这是新手最容易踩的一个坑。3.2 高频逆变桥与PWM控制高频逆变器我优先用“Universal Bridge”模块桥臂数设为2也就是单相全桥。器件类型根据工作频率选MOSFET把默认的缓冲电阻snubber resistance关掉或者设置一个很大的值比如1e6欧避免默认吸收电阻影响谐振回路等效阻抗。控制信号方面如果只是验证传递功率可以用PWM Generator给固定占空比的互补方波加入死区时间0.5微秒。对于PT-WPT固定PWM不是最终方案因为系统需要负阻抗效应。更接近原理的做法是加一个频率跟随控制器采样发射端线圈电流相位用锁相环锁定到谐振频率再生成开关信号。Simulink里的PLL模块直接可用但参数要微调否则在仿真初期容易失锁。如果只为了快速看到PT特性负阻模型更稳把全桥逆变器替换为Controlled Voltage Source输入为-Rg × i周围留一个测量环节。两种模型我都会写在测试脚本里用标志位切换这种可配置结构对复现论文中的多工况对比非常有帮助。3.3 耦合线圈与互感建模线圈建模是无线电能传输仿真中的关键。Simulink里面最常规的做法是用“Mutual Inductance”模块分别设置Primary自感、Secondary自感和互感M。如果版本里没有这个模块也可以用两个电感线圈配合互感参数实现但要注意同名端。线路连接时必须严格区分同名端。发射端电流流入同名端的正方向要和接收端感应电压极性配合一旦接反互感会变成负值传输功率方向颠倒仿真结果直接错掉。我建议把M定义成变量而不是直接写数字。比如在初始化脚本里写L 120e-6; k 0.2; M k * L;然后在Mutual Inductance模块的互感栏填写“M”。之后做耦合系数扫描时只需要改k模型不用动。这个习惯救了很多次后续的批量仿真。如果想把漏感也体现出来可以拆成耦合系数加理想变压器方式建模初级侧总电感L1、次级侧总电感L2再加一个理想变压器变比√(L1/L2)。这种做法在分析较高阶补偿网络时更方便也更容易看漏感对PT区宽度的影响。3.4 补偿网络与杂散参数处理SLSPC高阶补偿网络在Simulink里实现并不难无非是R、L、C模块按拓扑结构连接。难在元件的“非理想性”。如果全用理想电容和理想电感谐振Q值会接近无穷仿真中会出现极高电压尖峰数值上很容易发散这和实际电路也不符合。建模时我给每个电感串一个小电阻阻值用线圈的ESR等效串联电阻一般在0.05欧到0.5欧之间。给电容也加一个很小的ESR0.01欧级别。这样系统在谐振时的阻抗特性才接近真实硬件PT平衡条件才有实际的求解意义。高阶拓扑的具体连接方式我建议先用交流扫描法验证在Simulink里用“AC Sweep”或者手动进行线性分析看输入阻抗模值和相位。理想情况下在工作频率处输入阻抗相位为零模值接近预期的等效负载。若相位不正返回去调整补偿电容值而不是盲目改仿真步长。3.5 整流、负载、后处理模块接收端用二极管整流桥模块负载可以接纯电阻、电池模型或者受控电流源。复现论文的恒功率平台时我会在负载端接一个电阻然后通过脚本改变阻值记录输出功率。整流桥后的大电容会对PT特性产生明显影响电容太小输出电压纹波大电容太大整流桥非线性负载会让接收端电阻折算变得复杂。一般先取输出滤波电容为10μF到100μF之间以电压纹波小于5%为准。测量用“Current Sensor”和“Voltage Sensor”分别放到逆变器输出侧和负载侧信号通过“To Workspace”输出到MATLAB工作区。建议把保存格式设为“Timeseries”后面可以直接用mean、rms、fft这些函数做分析。4. 仿真中遇到的坑与排查技巧4.1 频率分裂怎么观察和修复我在早期用理想负阻模型时输出功率和频率曲线出现两个峰按说PT-WPT不应该有频率分裂。排查后发现负电阻增益Rg比发射端损耗大了三倍系统已经跳出PT区从“实数特征频率解”退化到“共轭特征频率解”。判断方法很简单把输入阻抗的相频特性画出来如果只在一个频率处过零说明PT条件成立如果看到两个过零点就是频率分裂。修复方法用公式计算总损耗把Rg调到略大于损耗而不是大很多。实际操作中Rg一般取发射端损耗的1.05到1.2倍能够在宽耦合范围内稳定。然后扫描耦合系数k看输出功率是否会像WPT文献里那样平顶化。我的结果里k从0.15到0.35之间时负载电压波动不到5%说明系统确实进入了PT区。4.2 求解器配置与仿真时间过长85kHz的开关频率周期约11.76微秒。如果采用全桥开关模型仿真步长必须能解析开关边缘一般最大步长设为开关周期的1/200也就是大约0.06微秒。想仿真100毫秒的物理时间步数会非常庞大直接卡死。我的解决方案分三层第一层全模型跑5ms到10ms只看稳态结果把Stop time控制住没必要跑太长时间。第二层把全桥开关模型换成平均模型或负阻受控源模型这样步长可以放到1微秒以上仿真速度提升十倍。第三层用Powergui求解器工具栏把模式设为Continuous或Phasor。对稳态功率验证用Phasor求解器最快但看不到瞬态验证波形时再用Continuous。实际处理里我用连续模式加ode23tb最大步长设2e-6秒多数工况都能在几分钟内跑完。4.3 电压电流尖峰处理高阶谐振网络在启动瞬间容易产生很大的浪涌电流尤其是电容初始电压为零而直流源直接接入时瞬态电流可能达到稳态的几十倍。这不一定是物理问题但会让仿真变慢甚至中断。处理办法有两个一是给直流源加缓慢启动斜坡比如通过一个Step模块控制受控电压源使母线电压在0.2ms内从0升到设定值二是在逆变器桥臂上并联RC吸收支路吸收电容0.1μF、吸收电阻10欧这同样是实际电路中常见的做法。仿真里出现数值振荡时不要急着改小步长而是先看是不是负阻模型产生了正反馈。理想受控源v -Rg × i的响应带宽要限制我用一阶低通传递函数把增益在高频处滚降模型稳定性提升非常明显。4.4 结果验证与FFT分析复现论文光看波形还不够还得量化传给效率和谐波含量。Simulink里可以用power_fft_scope打开FFT分析窗口或者在脚本里用内置的fft函数处理To Workspace保存的数据。效率计算定义η P_out / P_inP_in mean(V_in × I_in)P_out mean(V_load × I_load)注意这里的V_in和I_in必须是逆变器输入端的直流电压电流而不是电网侧。在读取测量数据时最好先去除前1ms的启动瞬态避免效率计算被启动过程污染。谐波分析中如果开关频率的整数倍分量很大可以从SLSPC补偿网络的参数上找问题比如并联电容太靠近谐振点、串联电感量偏小。高阶网络的优点在这时就会体现出来相较二阶系统三次和五次谐波幅值会明显低一截。最后列一个参数速查表供复现时直接抄作业参数推荐取值备注开关频率 f085 kHz也可根据论文取6.78MHz线圈自感 L120 μH两线圈对称设计耦合系数 k0.15~0.35PT区典型工作区间初级补偿电容 C约29.2 nF按公式计算后微调负电阻增益 Rg总损耗的1.05~1.2倍过大容易失去PT特性仿真最大步长2e-6 s全开关模型建议小于开关周期/200输出滤波电容10~100 μF以纹波小于5%为准5. 给复现者的一句话经验真正开始复现这套高阶PT-WPT系统后我的体会有两点。第一不要在Simulink里急着把所有模块搭完再验证。先把两侧线圈单独建出来测出无补偿状态下的耦合特性再分别加发射补偿和接收补偿确认各自谐振最后才组合成PT-WPT闭环。每一层级都用仿真曲线确认一次最后联调时报错就会少很多。第二参数不要手写死在模块里。把L、C、M、Rg全部放到初始化脚本里用变量引用这样调参、做批量扫描、论文出图都特别顺手。另外想提醒的是Simulink的版本差异会影响模块名称和属性结构。我这里用的模块在R2021b之后都兼容但老版本里Universal Bridge的参数窗口可能叫法不一样。复现时如果找不到某个模块优先看该模块是否在Simscape Electrical库中而不是怀疑模型搭错。最后分享一个小技巧把耦合系数k和负载电阻R_load设计成MATLAB脚本里的for循环变量批量跑完工况后把输出电压、输出功率、效率存成结构体最后统一绘图。这样复现论文的效率-距离曲线、输出-负载曲线都能一次性完成而且图表风格完全可控。这个工作流一旦跑通后续做参数优化、机器学习代理模型、硬件在环测试都能直接复用。