最近在复现一篇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πf0
C = 1 / (ω0² × L)

代入f0 = 85kHz、L = 120μH：
ω0 = 2π × 85000 ≈ 534070 rad/s
C = 1 / (534070² × 120e-6) ≈ 2.92e-8 F

也就是大约29.2nF。这个值听起来不大，实际电路中还要把线圈的分布电容、PCB走线杂散电容考虑进去，仿真初算时可以先用理想值，后面做容差分析时再偏离。

第四，计算耦合系数和互感。两线圈之间的互感M与自感的关系是：

M = k × √(L1 × L2)

如果L1和L2都等于120μH，k取0.2，M就是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欧之间。给电容也加一个很小的ESR，0.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_in

P_in = mean(V_in × I_in)
P_out = mean(V_load × I_load)

注意这里的V_in和I_in必须是逆变器输入端的直流电压电流，而不是电网侧。在读取测量数据时，最好先去除前1ms的启动瞬态，避免效率计算被启动过程污染。

谐波分析中，如果开关频率的整数倍分量很大，可以从SLSPC补偿网络的参数上找问题，比如并联电容太靠近谐振点、串联电感量偏小。高阶网络的优点在这时就会体现出来：相较二阶系统，三次和五次谐波幅值会明显低一截。

最后列一个参数速查表供复现时直接抄作业：

5. 给复现者的一句话经验

真正开始复现这套高阶PT-WPT系统后，我的体会有两点。第一，不要在Simulink里急着把所有模块搭完再验证。先把两侧线圈单独建出来，测出无补偿状态下的耦合特性；再分别加发射补偿和接收补偿，确认各自谐振；最后才组合成PT-WPT闭环。每一层级都用仿真曲线确认一次，最后联调时报错就会少很多。第二，参数不要手写死在模块里。把L、C、M、Rg全部放到初始化脚本里，用变量引用，这样调参、做批量扫描、论文出图都特别顺手。

另外想提醒的是，Simulink的版本差异会影响模块名称和属性结构。我这里用的模块在R2021b之后都兼容，但老版本里Universal Bridge的参数窗口可能叫法不一样。复现时如果找不到某个模块，优先看该模块是否在Simscape Electrical库中，而不是怀疑模型搭错。

最后分享一个小技巧：把耦合系数k和负载电阻R_load设计成MATLAB脚本里的for循环变量，批量跑完工况后，把输出电压、输出功率、效率存成结构体，最后统一绘图。这样复现论文的效率-距离曲线、输出-负载曲线，都能一次性完成，而且图表风格完全可控。这个工作流一旦跑通，后续做参数优化、机器学习代理模型、硬件在环测试，都能直接复用。