SVPWM仿真波形详解：马鞍形调制波与六阶梯波的Matlab实现

最近在Matlab2019B里搭SVPWM仿真模型折腾到后半夜示波器里同时出现了6阶梯波和马鞍形调制波那一刻还挺有成就感的。很多人第一次看到这两个波形第一反应是“模型搭错了”——马鞍形怎么看都不像正弦6阶梯又怎么看都不像常见的高频PWM。其实这两个波形恰好是SVPWM区别于SPWM最直观的两个特征前者说明调制波里注入了零序分量后者说明电压利用率已经被拉满。这篇就把这套仿真模型怎么搭、波形怎么认、Matlab2019B里哪些地方容易翻车一次说清楚。适合正在做逆变器仿真、电机控制或者被课程设计、毕业设计里的SVPWM折磨的读者。1. 先解释清楚马鞍形和6阶梯波到底从哪里来1.1 SPWM的电压利用率瓶颈要理解SVPWM为什么要搞出个马鞍形得先看看SPWM的短板。SPWM的套路很简单三相正弦波分别和同一个高频三角载波比较比较结果驱动六个开关管。正弦波峰值一旦超过载波峰值就会出现削顶所以正弦调制波幅值最多只能到载波幅值。这样算下来SPWM在线性区内相电压基波幅值最大只能是 Vdc/2线电压基波幅值最大是 (√3/2)Vdc约等于 0.866Vdc。这意味着什么假设直流母线是600VSPWM能输出的线电压基波峰值最多519.6V而理想情况下你明明有600V的母线有将近14%的电压被“浪费”了。在电机控制里这就直接限制了电机能跑到的最高转速。想跑更高要么抬高母线电压要么换更高压的器件成本嗖嗖往上走。所以业内一直在琢磨怎么把这个电压榨出来。1.2 SVPWM的空间矢量视角SVPWM换了条思路它不直接看三相正弦波而是把三个相电压看成一个空间矢量。三相桥臂只有八种开关状态六个有效矢量和两个零矢量。六个有效矢量在复平面上把一个360°圆周切成六个扇区每个扇区60°。只要参考电压矢量落在某个扇区里就用这个扇区相邻的两个有效矢量去合成再穿插零矢量调整作用时间。采样的那一刻算出两个有效矢量各作用多久下个采样周期再重新算一遍。这样磁链轨迹就不是SPWM那种“正弦推出来”的形状而是用六个方向的矢量一点一点逼近圆形。这就是SVPWM的核心不是盲目提高正弦调制波的幅值而是改变调制策略本身让电压矢量能跑满整个六边形内切圆。1.3 马鞍形的本质零序分量注入那马鞍形调制波是哪来的现在很多教材会给一个结论SVPWM等效于在SPWM的三相正弦波里注入一个零序分量。这个零序分量一般取三相参考电压最大值和最小值的平均值再取负写成公式就是U0 -0.5 × (Umax Umin)把U0加到三相正弦波上每一相的调制波就从纯正弦变成了马鞍形。因为它把正弦波峰值附近往里“压”了一些峰值降下来了原来那些超过载波幅值的部分被拉回线性区内所以基波幅值可以继续往上加直到和SPWM过调制前的极限相比线电压基波幅值能提高到 Vdc直接比SPWM多了15.47%。这个零序分量有个特点三相里加的完全一样所以线电压里它会被抵消。换句话说马鞍形里多出来的那些3次、9次、15次谐波在相电压里看得到在线电压和线电流里基本不存在。很多新手看到相电压FFT里一堆3次谐波就慌了其实这是SVPWM正常工作的铁证。1.4 6阶梯波的真正含义很多人以为SVPWM的线电压波形直接就是一个六阶梯拿来示波器一接发现全是高频脉冲还以为模型坏了。这里要花点功夫说清楚如果把示波器直接接在逆变器输出端你看到的永远是高频PWM脉冲串不是光滑的阶梯。所谓6阶梯波是指把这些脉冲按开关周期平均之后或者直接在控制器里重构等效线电压调制波得到的电压轮廓。为什么平均之后是阶梯状因为在一个扇区里三相桥臂的开关组合是固定的几个状态线电压的电平在60°范围内基本维持同一个水平到下一个扇区切换成另一个固定电平。六个扇区走完线电压平均轮廓就成了六个台阶每个台阶恰好对应60°。所以正确的理解方式是SVPWM实际输出的瞬时电压是PWM脉冲但它的基波等效电压是六阶梯波马鞍形是相电压调制波六阶梯是等效线电压轮廓。两者一个是“因”一个是“果”。2. Matlab2019B仿真模型的整体结构与核心模块2.1 顶层架构从直流母线到三相桥臂我用的版本是Matlab2019BSimulink里搭的模型整体分两大部分主电路和控制系统。主电路很简单直流电压源600V接三相全桥逆变器输出接星形负载每相电阻10Ω、电感5mH。逆变器我直接用了Simscape Electrical里的Universal Bridge选IGBT/Diode桥臂三相桥臂结构这样省去自己一个一个搭IGBT的麻烦。2019B的Universal Bridge里需要手动指定Ron、Snubber电阻电容这些参数Snubber我习惯设成默认值只是仿真速度会稍慢一点。控制系统负责生成6路PWM先产生三相正弦参考波频率50Hz幅值280V送进SVPWM算法模块算法输出三相占空比占空比再和三角载波比较最终生成IGBT的6路门极信号。这里面最容易忽略的是Universal Bridge的门极输入需要的是Simulink信号但比较器出来的往往是boolean类型2019B的电力电子模型对输入信号类型很敏感类型不对直接报错或者仿真卡死。我在门极前统一加了Data Type Conversion全部转成double再送进去。2.2 扇区判断与矢量作用时间核心算法实现SVPWM算法这块我在模型里用了一个MATLAB Function模块直接写函数实现扇区判断和作用时间计算。第一步先把三相参考电压做Clarke变换得到αβ轴分量然后计算三个中间变量X、Y、Z再根据αβ分量的符号判断扇区。中间变量的好处是六个扇区的T1、T2查表可以直接复用不用每个扇区都重新算三角函数。function [sector, T1, T2] svpwm_sector_time(ua, ub, uc, Vdc, Ts) % Clarke变换等幅值变换 ualpha (2/3) * (ua - 0.5*ub - 0.5*uc); ubeta (2/3) * (sqrt(3)/2*ub - sqrt(3)/2*uc); % 中间变量X、Y、Z X sqrt(3) * ubeta * Ts / Vdc; Y (1.5*ualpha sqrt(3)/2*ubeta) * Ts / Vdc; Z (-1.5*ualpha sqrt(3)/2*ubeta) * Ts / Vdc; % 扇区判断N值映射到1~6扇区 B1 double(ubeta 0); B2 double((sqrt(3)/2*ualpha - 0.5*ubeta) 0); B3 double((-sqrt(3)/2*ualpha - 0.5*ubeta) 0); N B1 2*B2 4*B3; switch N case 3, sector 1; case 1, sector 2; case 5, sector 3; case 4, sector 4; case 6, sector 5; case 2, sector 6; end end扇区N的映射关系不同教材写法可能不一样关键是你的扇区判断逻辑要和后面的T1、T2查表严格对应否则波形一定会乱。这块我踩过一次换了个教材的扇区判断公式忘了改查表出来的相电流波形整个是歪的。2.3 七段式发波与三相占空比分配算出T1、T2之后要分配给三相这里用最常见的七段式SVPWM。七段式的思路是把零矢量拆成两段一段放在周期开头一段放在周期中间这样三相开关动作对称谐波集中且幅值较小。先算出零矢量时间 T0 Ts - T1 - T2然后定义三个切换点Tcm1 T0 / 4 Tcm2 T0 / 4 T1 / 2 Tcm3 T0 / 4 T1 / 2 T2 / 2不同的扇区要往A、B、C三相分配不同的切换点。我常用的切换点分配表如下扇区A相切换点B相切换点C相切换点ITcm1Tcm2Tcm3IITcm2Tcm1Tcm3IIITcm3Tcm1Tcm2IVTcm3Tcm2Tcm1VTcm2Tcm3Tcm1VITcm1Tcm3Tcm2得到三相切换点之后把它归一化到0到1的占空比再和三角载波比较当载波值小于切换点时输出高电平否则输出低电平。注意这里的比较方向决定了输出波形是对齐还是反相我习惯用载波小于切换点时输出1。如果你不想这么麻烦其实SVPWM在载波比较视角下有个等效做法直接在MATLAB Function里算出三相正弦波后注入零序分量生成马鞍形参考波再和三角载波比较。两种方法结果等价但查表法更容易理解扇区切换零序注入法更容易直接观察到马鞍形。我在模型里两种都做了用Mux把两条路的波形拉到同一个示波器上对比互相佐证。3. 波形怎么看正确波形和错误波形的差别3.1 马鞍形调制波的观测点这是新手最容易搞混的地方。马鞍形调制波不是让你在逆变器输出端测出来的——输出端全是硬邦邦的PWM脉冲。马鞍形是控制器内部“等效调制波”的概念必须从算法模块里拉出来看。我用两个方法观测马鞍形方法一查表法里重构。在SVPWM函数外面把六个扇区对应的三相切换点还原成等效调制波即每相切换点相对载波峰值的比值。这个方法稍微绕一点但能明确看到每个扇区切换时刻的波形变化。方法二零序注入法。直接搭一个简单的数学模块组用MinMax模块取三相参考的最大值和最小值相加取负再乘0.5得到U0然后把三路正弦分别加上U0送到Scope。这种方法最直观看到的波形就是经典的马鞍形。我自己实际观测时是把A相的马鞍形调制波和三角载波放在同一个Scope里看。载波频率10kHz基波50Hz一个基波周期内能看到200个载波周期马鞍形的包络非常清晰。观察时间轴要放到20ms以上如果你只放1ms看到的就是满屏三角波和窄窄的脉冲根本看不出马鞍形状。马鞍形正常的特征波峰比纯正弦平缓波谷也往上收整体像一个“压扁的馒头”两边对称。如果只看到上下不对称先检查注入的零序分量公式符号是不是取反了。3.2 6阶梯波怎么观察6阶梯波是等效线电压轮廓所以最直接的做法是把三相马鞍形参考波两两相减得到线电压参考。因为马鞍波里的零序分量在相减时抵消了剩下的就是只含基波和6的倍数次谐波的线电压参考在示波器上呈现六个台阶。具体操作在模型里复制一份马鞍形调制波的信号加一个Subtract模块用A相减B相接到Scope。波形应该是一个周期内出现六个平台每个平台宽度大约3.33ms对应60°。平台之间的跳变刚好对应扇区切换。有朋友问能不能直接看逆变器输出的线电压VAB。可以但你看到的是高频脉冲不是阶梯。如果你在Scope里打开均值滤波或者把仿真数据导到MATLAB工作区里做滑动平均窗口设一个开关周期100us也能恢复出六阶梯轮廓。这个方法虽然多一步但好处是顺便验证了主电路和驱动信号是否正确——如果主电路接线错了重构波形再漂亮也没用。3.3 正确波形特征的量化对照光说“看起来像”不行仿真调试要有点量化的标准。我给自己列了一个波形特征检查表检查项正常表现异常表现可能原因马鞍形对称性上下半周期对称顶部平、底部尖零序注入符号错误或载波比较方向不一致马鞍波峰值280V参考时峰值约300V接近正弦峰值346V占空比归一化错误或Ts与载波周期不匹配六阶梯平台宽度六个台阶均匀各3.33ms平台宽度不均扇区判断逻辑顺序错或N值映射表不对线电压脉冲脉冲宽度按载波周期均匀调制脉冲成片缺失死区过大或驱动信号极性接反我这里马鞍波峰值“约300V”是因为注入了零序后波形峰值比纯正弦小具体数值取决于调制比不必死记关键是和理论计算对得上。4. Matlab2019B调参陷阱与避坑记录4.1 代数环警告SVPWM模块最容易报的错Simulink模型跑之前出现Algebraic Loop警告这个问题几乎每个搭SVPWM的人都会撞上。原因是SVPWM计算模块的输出依赖当前时刻的输入而在连续系统中Simulink会尝试在每个仿真步内同时求解输入输出关系形成环。2019B遇到这个情况并不会直接拒绝仿真但会严重影响速度有时波形还抖得厉害。我的处理办法是在SVPWM算法模块前加一个Unit Delay把三相参考输入延迟一个采样周期。开关频率10kHz延迟一个周期100us对50Hz基波的相位影响微乎其微但一下子把代数环破掉了。更干净的方案是把SVPWM算法模块设成离散模块采样时间设成Ts1e-4这样MATLAB Function在每个固定步长运行一次天然不存在代数环问题。但要注意离散模块的输出类型要保持一致2019B里如果输出信号有维数变化记得用Squeeze或者Reshape调整。4.2 求解器步长和载波周期的匹配这一条对波形观感影响极大。很多人搭好模型跑出来的马鞍形毛刺多得没法看第一反应是算法错了其实只是仿真步长太粗。三角载波频率10kHz载波周期100us。要比较器正确解析出载波形状仿真步长至少要小于载波周期的1/100也就是1us。我习惯把求解器设成离散固定步长1e-6这样每个载波周期内能采100个点三角波看起来平滑马鞍形包络也非常干净。如果把步长设成1e-510us一个载波周期只能采10个点三角载波变成折线比较器翻转的瞬间经常被跳过马鞍波边缘全是锯齿六阶梯平台宽窄也不均匀。这种问题不是算法错误纯粹是步长不够。我在做FFT分析时会把仿真数据用To Workspace存到MATLAB工作区设置成Structure with time格式再在MATLAB脚本里处理。存数据时采样间隔和仿真步长一致1e-6否则FFT的频率分辨率不够3次谐波根本分不出来。4.3 数据类型不匹配导致的驱动信号异常2019B对boolean和double的区分很严格。三角载波比较器输出的通常是boolean而Universal Bridge的门极输入期望的是double如果你直接连轻则warning重则仿真直接终止。另外IGBT门极信号如果是“互补”的上下桥臂不能同时导通。我一开始图省事用两个Relation Operator分别比较“大于”和“小于”结果在零序分量过零附近上下桥臂切换瞬间出现了一小段同时导通区母线电流直接冲高。后来统一改成NOT逻辑取反把互补PWM的关断延迟通过Transport Delay模块加死区才把尖峰压下去。提醒一点Simulink里Transport Delay模拟死区时间很直观但延迟模块的初始输出要设置成0否则仿真一开始会出现莫名其妙的误导通。4.4 死区时间对马鞍形对称性的影响加死区之后马鞍形调制波会出现一个不太起眼但很烦人的现象波峰和波谷附近的小台阶被抹掉波形边缘发圆有时候看起来左右不对称。原理很简单死区时间内上下桥臂全关输出端被负载钳位到某个不确定电平相当于在理想PWM波形上叠加了一个与电流方向相关的误差脉冲。这个误差脉冲在高频下呈现固定宽度反映到平均电压上就是一个和开关频率相关的电压偏差。我实测下来10kHz开关频率下死区时间从2us加到8us线电压基波幅值会往下掉1%到3%马鞍形波峰处明显出现“削顶”感。如果只是验证马鞍形和六阶梯建议先不要加死区把算法波形验证通过之后再补死区环节否则两个问题混在一起很难定位是算法错还是死区错。5. 用基波幅值和FFT反推算法正确性5.1 基波幅值对比SVPWM与SPWM的差距仿真波形“看着对”还不够我一向会用数学结果给自己一个确凿的证据。最硬的指标是基波幅值。前面说过SVPWM在线性区内相电压基波峰值最大Vdc/√3线电压基波峰值最大Vdc。而SPWM相电压基波峰值最大Vdc/2线电压基波峰值最大0.866Vdc。我在同一个模型里搭了一个SPWM对照组用纯正弦波和三角载波比较两套系统共用同一个直流母线600V和同一个负载。FFT分析结果显示SPWM线电压基波幅值约519VSVPWM线电压基波幅值约600V两者比值1.1547和理论值分毫不差。这个数据一出比看十个波形都管用。马鞍形波峰降低了但基波幅值反而上去了——SVPWM省出来的那14.5%电压就在这里。5.2 FFT分析零序分量去了哪里这是我最想让大家亲手做一次的验证。把相电压UA和线电压UAB分别存下来用MATLAB的fft函数或者2019B里的powergui FFT Analysis工具做频谱分析。结果应该是相电压UA基波50Hz占主导同时有明显3次、9次、15次谐波3次谐波含量通常能到基波的10%~20%这正是马鞍形里零序分量的体现。线电压UAB3次、9次谐波几乎消失基波占绝对主导。为什么因为0序分量在三相里大小相等、方向相同线电压是两相相减零序直接抵消。马鞍形里那些“不平滑”的部分在线电压里根本不造成污染。这个验证做完你对SVPWM的理解会上一个台阶——马鞍形不是畸变而是注入的“好谐波”。5.3 扇区切换时刻核对最后一个小技巧把扇区判断模块输出的N信号也拉进Scope和马鞍形调制波放在一起看。正常情况下N信号应该严格对应参考电压矢量所在的扇区切换时刻正好落在六阶梯波每个平台的跳变沿上。基波50Hz一个周期20ms六个扇区每个停留3.33ms节拍均匀不多不少。如果某个扇区停留时间明显偏长或偏短那说明Clarke变换的系数、或者扇区判断里的阈值符号出了问题优先检查ualpha和ubeta的符号是否和电机旋转方向约定一致。6. 一点个人体会与后续可做的事仿真跑通不等于真的理解了SVPWM我自己的感觉是波形能复现只是完成了“验算”真正理解是在把马鞍形和六阶梯的成因讲清楚之后。建议拿到本文模型后把直流母线电压、参考幅值、载波频率各改一次每次都用FFT和基波幅值核验一遍比单纯看波形快得多。如果还想往下深挖可以试试把七段式改成五段式发波。五段式两个扇区之间会有一相钳位不动作马鞍形会变成底部贴地的“钳位形”开关损耗直接降三分之一但波形里会多出较强的6次谐波成分。系统是否真的“更优”取决于负载工况和EMI需求值得在仿真里对比一把。另外把开环的SVPWM接上转速闭环或者磁链观测器又是另一个大工程。建议每一步都保留波形快照和FFT结果不然等闭环调PI参数的时候算法环节的风吹草动全都会被PI背锅那才是真正的噩梦。