基于SLSPC的高阶PT-WPT系统Simulink仿真复现指南
2026/9/9 18:58:37 网站建设 项目流程

你是不是也有这种经历:导师甩给你一篇SCI一区论文,说“把它仿真复现一下”,你满怀信心打开PDF,两页公式看完就懵了。今天说的这个项目——基于SLSPC系列的高阶PT-WPT无线电能传输系统,就是我曾经被“甩”过的一篇。前后折腾了将近三周,把Simulink模型从零搭到能跑,把论文里的补偿网络一步一步翻译成电路模块,最后总算把波形和效率曲线给对上了。这篇文章专门写给正在做无线电能传输仿真复现的人,尤其是刚接触PT对称、高阶补偿这些概念的同学。我会把SLSPC拓扑怎么理解、参数怎么算、Simulink里具体怎么搭、中途踩了哪些坑,原原本本讲清楚。

1. 项目背景与论文思路拆解

1.1 这个课题到底在做什么

很多人第一次看到“基于SLSPC系列的高阶PT-WPT无线电能传输系统”这个标题,第一反应是查缩写。PT-WPT在近几年的无线电能传输文献里,通常指Parity-Time对称无线电能传输,也就是宇称-时间对称WPT。这个思想最早来自量子力学里的PT对称理论,后来被引入到电磁学和电路领域,核心是让系统的增益和损耗刚好达到平衡,使系统在小扰动下表现出非常稳定的输出特性。对WPT来说,这意味着什么?意味着负载变化、线圈偏移、距离改变时,输出功率不用靠复杂的控制算法去补偿,也能基本保持不变。

这个特性对无线电能传输来说太有吸引力了。常规的磁耦合谐振WPT系统,只要负载一变化或者发射线圈和接收线圈稍微错位,输出功率和效率就会剧烈波动,必须额外加阻抗匹配网络、频率跟踪算法或者DC-DC后端来稳压。而PT对称结构相当于把“自适应”做进了电路本身,在特定区间内天然具备恒功率输出能力。这篇论文的贡献就在于,把PT对称机制和高阶补偿网络结合起来,用SLSPC系列拓扑去构造这种有源-无源平衡关系,从而在较宽的耦合系数和负载范围内保持传输特性平稳。

再说SLSPC。这个缩写原论文里给出的是一个拓扑族,逐字母拆开看,S是串联(Series),L是电感支路(Inductor),SPC是串并联补偿(Series-Parallel Compensation),合在一起可以理解为“原边串联加电感支路、副边串并联混合补偿”的一类高阶网络。不同论文里这个缩写可能指代不同结构,复现之前一定要回到原文看拓扑图。我这里讲的是我当时复现时采用的解读,整个Simulink模型也是按这个思路建的。

1.2 为什么值得花时间复现这篇论文

一区论文不一定都复杂,但通常把一个关键问题讲透了。这篇论文的价值不在某个花哨的新元件,而是把两个相对成熟的概念——高阶补偿网络和PT对称原理——有机地结合到了一起。复现这种论文,收获是双重的:既要吃透高阶谐振网络的设计逻辑,又要把PT对称的实现条件在仿真里落地,这比单纯搭一个SS型WPT模型要深入得多。

从实操角度看,这篇论文的可复现性其实不错。它给出的系统结构明确,补偿网络和线圈参数都有据可循,Simulink仿真模型作为论文验证手段,关键波形和效率曲线也比较清晰。我复现的时候,目标是做到两点:一是把论文里提出的SLSPC高阶网络搭成可运行的仿真模型;二是通过参数扫描验证PT对称区间内的恒功率输出特性。如果你正在做类似课题,或者需要复现别的WPT论文,这套方法完全可以平移过去。

整个系统的信号链是这样的:直流电源 → 高频全桥逆变器 → 原边SLSPC补偿网络 → 发射线圈 → 接收线圈 → 副边SLSPC补偿网络 → 整流桥 → 滤波电容 → 负载电阻。和常规WPT系统相比,差别全在补偿网络那一块。二阶系统只有一个谐振电容,高阶系统则有多级LC谐振单元,每个单元对频率、增益、谐波抑制都有各自的贡献。

2. SLSPC拓扑与高阶补偿网络设计要点

2.1 从二阶补偿到高阶补偿,到底多了什么

理解SLSPC的前提,是先搞清楚为什么二阶补偿不够用。最基础的SS型WPT系统,原边和副边各一个补偿电容,谐振条件要求工作频率等于线圈自感和补偿电容的谐振频率。它的优点是结构简单、恒流特性好,但缺点也很明显:当耦合系数变化或者负载变化时,系统的传输功率和效率会明显波动。这是因为二阶谐振网络只有一个谐振点,对外部参数变化几乎没有“缓冲”能力。

高阶补偿网络加进去的,本质上是额外的谐振零点和极点。以LCC补偿为例,它在原边线圈前面又插入了一个电感支路和电容支路,这些支路和线圈自身电感、补偿电容共同作用,会在目标频率附近形成更陡峭的阻抗特性,同时让系统的等效输入阻抗在一定范围内对负载变化不敏感。换句话说,高阶网络是用更多的电抗元件,换取系统对参数漂移的容忍度。

把二阶SS、常规LCC-LCC和SLSPC放在一起对比,区别会更直观:

特性二阶SS补偿LCC-LCC补偿SLSPC系列
谐振网络阶数2阶4阶以上高阶,多谐振点
对负载变化的敏感度低,接近PT平衡时几乎恒定
对耦合系数变化的容忍度
恒流/恒压输出能力单一特性可通过参数调节结合串并联结构灵活设计
仿真复杂度较高
谐波抑制能力

2.2 SLSPC网络拓扑是怎么推算出来的

复现的时候,我先把论文给的拓扑图翻译成电路级描述。常见的二阶SS网络,原边就是电感Lp串联电容Cp,副边是Ls串联Cs。SLSPC系列在此基础上做扩展:原边保留串联谐振支路,并在输入端加入一个谐振电感,形成一个L型输入匹配结构;副边不再只是单一串联电容,而是拆成一条串联谐振支路和一条并联补偿支路,两条支路组合后再通过一个电感到整流桥。

用语言描述比较抽象,我拆成三个特征来记:

  • 原边输入端有一个附加电感Lf,和原边补偿电容一起构成输入谐振网络,它的作用是让逆变器的输出电流接近正弦,并抑制高频谐波。
  • 原边线圈和副边线圈各自保留基础串联谐振电容,保证能量能在线圈之间有效耦合。
  • 副边除了串联谐振电容,还并联了一个补偿电容支路,用来调节整流桥入口的等效阻抗,使系统整体更容易落在PT对称的工作区间。

这种结构的好处是,多了两个可调参数维度,原边和副边的阻抗特性不再锁死在单一谐振频率上。我在仿真里验证过,只要Lf和副边并联电容取值合适,系统的输出功率随负载变化会明显变平,这正是PT对称特性的体现。

2.3 关键参数计算与匹配思路

参数设计是复现中最容易翻车的一步,但逻辑其实很固定。先把工作频率定下来,WPT仿真里常用的是50kHz到150kHz这一档,我用的100kHz。然后根据线圈自感算基础谐振电容。以原边线圈Lp=50μH为例:

谐振电容计算公式是C = 1 / (ω²L),其中ω = 2πf。在100kHz下,ω约等于628318.5 rad/s,代入后得到C约等于50.7nF。副边线圈如果自感相同,补偿电容也取50.7nF。

附加谐振电感Lf我取20μH,对应补偿电容Cf也按同一公式算,大约126.7nF。耦合系数k取0.3,互感M = k√(Lp×Ls),也就是0.3×50μH=15μH。负载电阻先按20Ω设置。

我把整套初始参数整理成了一张表,仿真搭建时直接往参数脚本里填就行:

参数符号数值
输入直流电压Uin100 V
开关频率fs100 kHz
原边线圈自感Lp50 μH
副边线圈自感Ls50 μH
耦合系数k0.3
互感M15 μH
原边补偿电容Cp50.7 nF
副边补偿电容Cs50.7 nF
附加谐振电感Lf20 μH
附加谐振电容Cf126.7 nF
负载电阻RL20 Ω

这里要特别说一句:这套参数的核心验证逻辑,是让原边和副边的基础谐振频率都落在100kHz,而附加谐振网络的谐振点稍微偏移,让系统在工作频率附近形成两个相邻的谐振模式。PT对称系统恰好需要这种接近简并的模式结构。如果所有支路都谐振在同一个频率,PT特性反而出不来。

3. Simulink仿真建模全流程实操

3.1 环境准备与模块选型

我用的是MATLAB R2021b以上版本,Simscape Electrical库里的Specialized Power Systems子库,老版本叫SimPowerSystems,界面略有不同,但元件用法基本一致。搭建WPT功率回路,核心用到的模块就四类:电源和测量、开关器件、无源元件、耦合线圈。

  • 直流电压源用DC Voltage Source。
  • 逆变器我建议直接用Universal Bridge模块,设成两桥臂(2 arms),开关器件选MOSFET。这样省去自己搭四个开关管的麻烦,而且内部自带缓冲电路,仿真不容易发散。
  • 无源元件用Series RLC Branch,把电阻值设成0或者一个很小的寄生电阻,设成电感或电容模式。
  • 耦合线圈用Mutual Inductance,双绕组模式,直接填Lp、Ls、M三个参数。

所有参数不要手填在模块对话框里,改成用MATLAB脚本统一初始化,模型变量名引用脚本变量。这样后面做参数扫描,不用反复改模块参数,非常省事。初始化脚本内容如下:

% 系统参数初始化 fs = 100e3; % 开关频率 100 kHz w = 2*pi*fs; % 角频率 Lp = 50e-6; % 原边线圈自感 Ls = 50e-6; % 副边线圈自感 k = 0.3; % 耦合系数 M = k*sqrt(Lp*Ls); % 互感 Cp = 1/(w^2*Lp); % 原边补偿电容 Cs = 1/(w^2*Ls); % 副边补偿电容 Lf = 20e-6; % 附加谐振电感 Cf = 1/(w^2*Lf); % 附加谐振电容 RL = 20; % 负载电阻 Uin = 100; % 输入直流电压

3.2 功率主回路搭建步骤

搭建顺序我建议从后往前:先搭负载和整流,再搭副边补偿,然后是线圈,接着是原边补偿,最后才是逆变器和电源。这种从后往前的顺序有个好处,每一级模块插上后,都能有一个明确的测试点,方便定位问题。

整流部分比较简单,用四个二极管搭全桥整流,后面加一个滤波电容和负载电阻。滤波电容一般取10μF到100μF之间,太小了输出电压纹波大,太大了启动电流冲击明显。我测试下来20μF在这个功率等级下比较合适,纹波和启动速度都能接受。

副边补偿网络按之前说的结构连:接收线圈的一端接串联谐振电容Cs,然后接整流桥;接收线圈的另一端接一个并联补偿电容支路,再汇合到整流桥的另一端。之后在整流桥入口处串入附加谐振电感Lf。这个电感的位置很关键,它同时承担谐振和滤波两个功能,不能省略。

线圈部分用Mutual Inductance模块,原边绕组和副边绕组没有电气连接,各自引出两个端子。设置的时候,自感分别填Lp和Ls,互感填M,注意同名端符号。同名端接反会导致耦合振荡相位完全不对,输出功率会异常低落,这是非常容易踩的坑。

原边补偿网络稍微简单一些,一端接发射线圈,另一端接逆变器输出的一个桥臂中点。原边补偿电容Cp与电感Lp串联,附加输入端再加一个Lf就行。需要说明的是,这里Lf和Cp的节点位置如果接错,输入阻抗特性会完全改变,判断标准是逆变器输出电流应该接近正弦,如果看到明显的尖峰电流,多半就是网络形态不对。

逆变器用Universal Bridge,直流侧接100V电压源,交流侧两个输出端分别接原边补偿网络的输入端和电源地。

3.3 PWM驱动与死区配置

WPT系统里的逆变器工作在固定频率方波模式,不需要SPWM或SVPWM,用50%占空比方波驱动全桥即可。关键是两个桥臂的驱动信号要相差180度,对角桥臂同时导通。

最简单的方法是用两个Pulse Generator,一个初始相位0度,另一个初始相位180度,频率都设成fs。脉冲宽度设为50%。如果用的是Universal Bridge,四个开关管的门极信号在模型里要通过一个Demux拆开,或者使用内部的PWM Generator模块。

死区时间不能省。虽然纯仿真里开关管是理想器件,但如果你后面要接硬件,没有死区习惯的话很容易出问题。在Universal Bridge里可以直接设死区时间,我通常设0.5μs到1μs,对应100kHz开关频率,死区占开关周期的5%到10%,是合理的。

有一点要特别注意:Pulse Generator是连续时间模块,而Simscape Electrical里很多元件用离散求解器,两者混用会导致数据采样不匹配。解决方法是把所有PWM信号通过Rate Transition模块统一到同一个采样时间,推荐设置成1μs,也就是开关周期的十分之一,既能分辨占空比和死区,又不会让仿真速度慢到难以接受。

3.4 求解器配置与波形提取

Simscape Electrical的仿真步长配置直接影响模型能不能稳定跑起来。我试过默认的ode45,效果很差,开关瞬间的电流尖峰会让求解器频繁报错。后来改成ode23tb,问题少了很多。这个求解器对刚性系统有更好的适应性,电力电子电路正属于典型的刚性系统。

具体设置上,相对误差1e-4,最大步长限制为开关周期的百分之一,也就是1e-7秒。如果不限制最大步长,仿真会在开关切换时刻跳过关键细节,波形看起来平滑,实际上全是错的。仿真时间先跑2ms,大概200个开关周期,足够系统进入稳态。

波形测量用Current Measurement和Voltage Measurement模块,分别接在需要观测的支路上。输出功率不要用瞬时功率直接看,那样波形毛刺很大,建议把电压电流信号送到To Workspace,然后在MATLAB里用rms函数计算有效值,再用有效值乘积计算输出功率。效率计算则是在输入端加一个功率测量模块,把直流侧输入功率和交流侧输出功率相除。

最后处理的时候,我习惯把节点电压、原边电流、副边电流、输出电压这四路信号存下来,然后画在一张图里,和论文里的波形图对比。论文里的波形通常是在稳态下截取的,所以仿真数据也要取后半段的稳态区间,不要包含启动暂态。

4. 踩坑记录与问题排查实录

4.1 波形毛刺和系统发散

仿真最容易遇到的现象是波形前面一段剧烈震荡,然后求解器报错停止,或者波形布满密密麻麻的毛刺。这种情况九成是步长问题。高频开关电路的最大步长必须足够小,我之前用1μs最大步长跑100kHz系统,波形还能看,但一提高到200kHz就崩了。后来设置最大步长不超过0.1μs,稳定多了。

另一种常见原因是电感和电容的初始状态。Simscape里的电感默认初始电流是0,这没问题,但电容要特别注意,如果整流后的滤波电容初始电压是0,启动瞬间会有一个很大的充电电流冲击,可能导致求解器发散。可以在电容模块里设置初始电压,也可以接受这个启动暂态,只取稳态数据。我的做法是模型里加一个很小的预充电路电阻,启动后1ms断开,仿真更稳定。

4.2 输出功率远低于预期值的排查顺序

如果系统跑通了,但输出功率比论文低很多,不要急着调参数,按这个顺序排查:

先看逆变器输出端波形,确认输出电压是±100V方波,没有异常压降。再看原边电流,正常应该是接近正弦的交流波形,如果电流很小,说明原边补偿网络没有谐振起来。接着看副边电流和整流后电压,如果所有这些都正常但输出功率低,问题多半在耦合线圈或负载参数。

互感同名端接反是很多人的盲区。Mutual Inductance模块如果同名端接反,两个线圈的感应电压方向相反,耦合效果几乎为0,输出功率会掉到正常值的十分之一以下。判断方法是看原边电流和副边电流的相位关系,正常应该接近同相或反相,如果相差90度以上,多半同名端不对。

4.3 谐振频率和驱动频率不匹配的问题

很多人觉得“谐振频率100kHz,开关频率100kHz,两个一样不就行了”。实际不是这么简单。SLSPC高阶网络里有多个谐振支路,它们的谐振频率并不完全一样,而是分布在100kHz左右形成一个工作带。如果所有支路都严格谐振在100kHz,系统可能会表现出非常陡峭的频率选择性,稍微偏离就掉功率。反之,如果各支路谐振点稍微错开,系统在100kHz附近的响应会更平缓,PT对称特性更容易出现。

参数设置里,Lf和Cf是调节这个“频带宽度”的关键。Lf取20μH、Cf取126.7nF时,谐振频率其实还是100kHz,但由于和线圈支路参数不同,整个网络形成两个相邻模式,实际工作点处于两个模式之间的“最佳平衡点”。这个平衡点需要通过参数扫描来找到,不是一眼能看出来的。

4.4 与论文数据对比时的注意事项

复现论文最忌讳的是拿着仿真结果追求和论文图一模一样,这不现实,也没有必要。论文里的参数可能经过归一化处理,波形可能有测量误差。正确做法是看趋势一致性:论文说负载从5Ω变到50Ω时输出功率恒定,你的仿真曲线也应该是基本水平的;论文说某个频率点效率最高,你的效率曲线峰值位置应该能对上,具体数值差两个百分点问题不大。

如果趋势都不对,那问题就大了。这时候检查论文的补偿网络结构是否和你搭建的一致,特别是原边和副边的串并联关系。有的论文把副边分成两个支路,其中一个支路经整流桥后接负载,另一个支路是空闲谐振器,看漏一个支路整个模型就废了。

常见现象可能原因处理方法
波形毛刺多,求解器报错最大步长过大设置最大步长为开关周期的1/100
启动冲击电流过大滤波电容初始电压为0允许启动暂态,只取稳态数据
输出功率异常低互感同名端接反检查Mutual Inductance同名端
输出功率随负载变化剧烈高阶网络参数失配查看Lf、Cf参数,重新扫描匹配
输出电压纹波很大滤波电容偏小增大滤波电容
PWM驱动信号不完整采样时间不匹配用Rate Transition统一采样时间

5. 扩展思路与实操心得

5.1 用参数扫描验证PT对称特性

复现论文的高阶PT-WPT系统,仿真模型跑通只是第一步,真正验证PT特性需要做参数扫描。最直观的实验是扫描负载电阻,从5Ω逐步增加到50Ω,记录每个负载下的输出功率。普通二阶SS系统在这个范围内,输出功率曲线会明显下坠;而SLSPC系统如果PT对称条件满足,输出功率应该呈现出一个“平台”,波动控制在10%以内。

扫描耦合系数也有意义。固定负载20Ω,把k从0.1扫到0.5,观察输出功率。PT对称区间内,耦合系数在一定范围内变化时输出功率基本不变,超出区间后功率迅速下降,形成一个明显的工作窗口。这个窗口的宽度直接反映了系统抗偏移能力,也是论文里最常见的对比曲线。

5.2 从Simulink仿真走向实验验证

仿真模型的一个重要用途是给硬件实验指明方向。比如Lf和Cf的取值可以直接对应到实际电感电容选型,逆变器的开关频率和死区时间可以直接照搬到驱动板。我自己做硬件时,仿真里的稳态电流有效值用来选功率器件和导线线径,电容电压峰值用来确定器件耐压等级。

但仿真和硬件之间有明显的鸿沟。仿真里可以忽略线圈的交流电阻和电容的ESR,硬件里这些寄生参数对效率有显著影响。我做实验时,线圈用利兹线绕制,工作频率100kHz下的Q值明显高于普通漆包线。电容也尽量选C0G或薄膜电容,避免陶瓷电容的直流偏压特性导致容值偏移。

5.3 复现论文的经验与通用方法

这套复现过程走下来,我最大的体会是:高阶WPT系统的仿真搭建,核心不在于会连电路图,而在于知道每一级网络在干什么。先搭二阶SS系统,跑通,看波形;然后加入附加电感Lf,看波形有什么变化;再加入副边并联补偿电容,继续观察。每加一级,系统行为都会变,只有每一级的变化都在预期之内,最后整个SLSPC网络才能稳定工作。

我在实际复现中还养成了一个习惯:每次修改参数前,先把当前模型的波形截图保存,写上参数和现象,做成一份简单的实验记录。SLSPC这种多参数系统,调参空间很大,靠记忆很容易乱。有了记录,回头排查问题时一眼就能看出是哪次改动引入了问题。

后面再复现类似的LCC-LCC、LCL-LCL或者其他高阶WPT论文时,这套方法可以直接平移。万变不离其宗,先把基础谐振单元理解透,再谈高阶拓扑的组合,你也能在一个星期内把一篇一区论文的模型搭得明明白白。

需要专业的网站建设服务?

联系我们获取免费的网站建设咨询和方案报价,让我们帮助您实现业务目标

立即咨询