聊到单相 PWM 整流器仿真很多人第一反应是拿三相那套 dq 双闭环直接改。真搭过一次就会发现单相系统少了一相信息锁相和坐标变换都得重新想办法。最近整理了一个 MATLAB 2022a 下的单相整流器双闭环仿真项目把单极性调制和单极性倍频调制两种 PWM 方案都做了实现和对比还顺手把模型降版到低版本 MATLAB 也能打开的程度。这篇文章就把整个实现思路、控制参数设计、Simulink 搭建细节和调试踩坑都写透适合正在做电力电子课程设计、毕业设计或者想快速验证整流器控制算法的朋友参考。1. 项目到底做什么单相整流器的控制目标与方案选型1.1 单相 PWM 整流器解决了什么问题传统不控整流桥后面挂个大电容交流侧电流是尖峰脉冲谐波大、功率因数低晶闸管相控整流虽然能调压但会从电网吸入大量无功和谐波。单相 PWM 整流器用全控开关管IGBT/MOSFET替代二极管或晶闸管通过高频斩波让交流侧电流主动跟踪电压波形实现三件事直流输出电压稳定可调、输入电流正弦化、功率因数近似为 1。这个电路的拓扑本质上是 boost 型变换器直流侧电压必须高于交流电压峰值才能正常工作。以 220Vrms/50Hz 电网为例峰值约 311V工程上直流母线一般取 360~400V。控制上如果只调占空比开环输出电压会随负载剧烈波动电流波形也很难看所以才需要闭环。1.2 为什么选“电压外环电流内环”双闭环单相整流器的控制策略大致分两类间接电流控制和直接电流控制。间接电流控制通过控制调制波的幅值和相位来间接控制电流结构简单但动态响应差电网电压稍微波动电流就失控。直接电流控制把交流侧电流作为反馈量用电流闭环把实际电流快速拉到指令值附近抗扰动能力和动态性能都好一个量级。双闭环是直接电流控制里最经典的架构电压外环采样直流母线电压和参考值做误差通过 PI 调节器输出电流幅值指令 id_ref。本质上是功率平衡环——负载功率变化时母线电压先跌电压环检测到误差后调节有功电流指令让输入功率跟上输出功率。电流内环采集电感电流经过锁相和坐标变换得到 dq 轴分量用 PI 调节器快速跟踪 id_ref 和 iq_ref通常 iq_ref0输出调制波指令。电压环响应慢、电流环响应快这种级联结构天然带限幅保护把电压环输出限制在允许电流范围内负载突变或启动瞬间电流不会失控。相比单环直接控占空比双闭环把“慢”的电压控制和“快”的电流控制解耦参数设计也更直观。1.3 单极性调制和单极性倍频调制为什么要做对比标题里同时出现单极性调制和单极性倍频调制说明这个项目不只是做一个能跑的整流器而是要做两种 PWM 方案的对比验证。单极性调制实现简单、开关损耗低单极性倍频调制等效开关频率翻倍电流纹波和滤波要求都更低但四个开关管都高频工作损耗和死区影响需要权衡。这两种方案在 MATLAB 里用同一个双闭环框架就能切换主电路和控制部分完全共用只换 PWM 生成模块非常适合做横向对比。这也是为什么整篇文章我会把两个 PWM 子系统的搭建差异作为重点。2. PWM 调制策略拆解双极性、单极性、单极性倍频到底差在哪2.1 三种调制方式的开关状态与输出波形双极性 PWM 是最基础的方案四个开关管在载波周期内同时切换H 桥对角线上的两个管子一起导通输出电压 Uab 直接在高电平 Ud 和低电平 -Ud 之间跳变。好处是逻辑简单坏处是电压跳变幅度大电流纹波大四个管子都是硬开关。单极性 PWM 的思路是让 H 桥一个桥臂高频斩波另一个桥臂按工频切换。以正半周为例左桥臂上管高频开关右桥臂下管常通此时 Uab 只在 0 和 Ud 之间切换负半周左桥臂下管高频开关右桥臂上管常通Uab 只在 0 和 -Ud 之间切换。电压跳变幅度从 ±Ud 减半到 Ud但输出脉冲频率仍然是开关频率 fs因为每个工频半周内只能产生一种极性的脉冲。单极性倍频 PWM 把两个桥臂都变成高频开关但两个桥臂使用相位相差 180° 的三角载波。这样在一个开关周期内Uab 的脉冲被拆分成了两个等效输出纹波频率变成 2fs。正半周 Uab 在 0 和 Ud 之间以 2fs 的频率交替负半周在 0 和 -Ud 之间交替。2.2 单极性倍频的“倍频”原理倍频效果从波形上看得最清楚。普通单极性调制里载波频率是 fs每个载波周期 Uab 产生一个脉冲脉冲宽度由调制波与三角载波的交点决定。倍频调制里两路三角载波错开半个周期左桥臂和右桥臂分别在各自的载波交点处动作两个桥臂中点电压一叠加Uab 在每个载波周期里就出现了两个脉冲等效脉冲频率自然翻倍。这就好比两个人接力搬砖一个人搬一整天的效率和两个人错开时间轮流搬的效率完全不同——倍频不是让某个人干得更快而是让两个人交替干活整体节奏翻倍。2.3 三种方案的关键性能对比做一个直观的对比表对比项双极性 PWM单极性 PWM单极性倍频 PWMUab 电平数±Ud 两电平0、±Ud 三电平0、±Ud 三电平等效纹波频率fsfs2fs高频开关管数量4个2个4个相同条件下电流纹波最大中间最小约为双极性的1/2交流侧 THD 表现差中优滤波器体积大中小开关损耗高低高死区影响中中较大单极性倍频的电流纹波之所以小是因为电感电流的纹波峰峰值与开关周期成正比等效开关频率翻倍意味着每个脉冲周期缩短一半纹波幅值也几乎减半。同样的电流纹波指标下倍频调制可以把交流侧电感选得更小或者用同样的电感获得更好的电流质量。代价是四个开关管都高频开关损耗增加而且两路载波错相之后死区时间对波形的影响更敏感这个在后续调试时要特别注意。3. 双闭环控制的数学原理与 PI 参数设计3.1 单相系统里的虚拟正交信号与 dq 坐标变换三相系统可以直接用自然坐标系变换到 dq 同步旋转坐标系单相系统没这么方便——只有一相电压电流没法直接做 Park 变换。工程上有两种常用办法把实际信号当作 α 轴分量用延时四分之一基波周期20ms/45ms构造虚拟的 β 轴分量或者用 SOGI二阶广义积分器生成正交信号。SOGI 的频率适应性更好但结构稍复杂延时法实现最简单代码就一行对固定 50Hz 电网完全够用。仿真项目里两种都可以我习惯先用延时法快速搭通整个闭环再换成 SOGI 看频率波动下的表现。这个项目里为了简洁和低版本兼容用延时法足够了。构造出 vα、vβ 和 iα、iβ 后就可以用标准 Park 变换把交流量转成直流量d 轴分量 cosθ × α 轴分量 sinθ × β 轴分量 q 轴分量 -sinθ × α 轴分量 cosθ × β 轴分量θ 是电网电压的相位由单相锁相环PLL提供。锁相环的作用就是把 θ 锁到电网电压的相位上让 d 轴和电网电压同相。这样 d 轴电流就对应有功电流q 轴电流对应无功电流控制 q 轴电流为 0 就能实现单位功率因数。3.2 电流内环 PI 参数的计算过程电流内环的被控对象是交流侧电感 L 和它的寄生电阻 R写成传递函数是G_i(s) 1 / (L × s R)这个传递函数是个典型的一阶惯性环节时间常数 T L/R。用 PI 控制器进行极点对消把 PI 的零点放在被控对象的极点位置也就是令积分时间常数 Ti L/R这样闭环就变成一个干净的纯惯性环节PI 参数只剩下一个比例系数需要设计。工程上电流环带宽 w_ci 通常取开关频率的 1/20 到 1/10。以开关频率 10kHz 为例电流环带宽取 500Hz换算成角频率约 3140 rad/s。比例系数 Kp 按以下公式估算Kp L × w_ci Ki Kp / Ti w_ci × R代入实际参数L 4mHR 0.1Ω得到 Kp ≈ 12.6Ki ≈ 314。这里要注意 Ki 的量纲是“每秒积分增益”在 Simulink 的离散 PI 控制器里对应的是积分增益别和连续域混着设。3.3 电压外环 PI 参数与限幅设计电压外环的被控对象是直流母线电容加负载。电容本身是个积分环节负载决定系统的阻尼但负载会变化所以电压环不适合像电流环那样精确建模。工程上一般取电压环带宽为电流环带宽的 1/10 到 1/5让电压环明显慢于电流环这样才能保证级联稳定。我习惯先给 Kp_v 一个较小的初值 0.5~2Ki_v 取 20~50然后在仿真里看母线电压的启动波形微调。如果电压超调大、振荡就降低 Kp_v、抬高 Ki_v 同时配合限幅如果响应太慢就适当加 Kp_v。电压环 PI 的输出就是电流幅值指令 id_ref必须加限幅。参考这个例子交流 220V、直流 400V、额定功率 2kW输入电流峰值约 18.2A留 20% 裕量限幅设为 ±20A。这个限幅直接保护了开关管也是软启动能力的一部分——启动瞬间电压误差大电压环饱和在限幅值上电流被限制在 20A 以内不会出现恐怖的冲击电流。3.4 完整控制链路把整个信号流串一遍母线电压采样和参考值做差进入电压环 PI输出限幅后作为 id_ref电网电压进入锁相环得到相位 θ电感电流经虚拟正交和 Park 变换得到 id、iqd 轴电流环 PI 跟踪 id_refq 轴电流环 PI 跟踪 iq_ref0两个 PI 输出 vd、vq反 Park 变换得到调制波 m调制波进入 PWM 发生器产生四路开关信号。整个链路里唯一需要手动标定的就是两个 PI 和限幅值其他都是标准模块。4. Simulink 模型搭建实操从主电路到 PWM 子系统4.1 参数初始化脚本不要直接在 Simulink 模块里敲参数一定要用 .m 脚本统一初始化。这样换 MATLAB 版本、调参、给别人复现都方便。下面这个脚本是整套模型的基准参数%% 单相PWM整流器仿真参数 % 电网参数 Vg_rms 220; % 交流电压有效值 220V f_grid 50; % 电网频率 50Hz Vg_peak Vg_rms*sqrt(2); % 峰值 311V % 主电路参数 L 4e-3; % 交流侧电感 4mH R_L 0.1; % 电感寄生电阻 0.1Ω C 4700e-6; % 直流母线电容 4700uF Udc_ref 400; % 直流母线电压参考 400V P_out 2000; % 额定功率 2kW R_load Udc_ref^2/P_out; % 负载电阻 80Ω % 调制与采样 fs_sw 10000; % 开关频率 10kHz Ts_ctrl 1e-6; % 控制周期/仿真步长 1us Vtri 1; % 三角载波幅值(归一化) % PI参数(基于极点对消法整定) w_ci 2*pi*500; % 电流环带宽 500Hz Kp_i L*w_ci; % ≈12.57 Ki_i w_ci*R_L; % ≈314.16 w_cv w_ci/10; % 电压环带宽约为电流环的1/10 Kp_v 1.0; % 先给初值再微调 Ki_v 30; id_max 20; % 电流幅值限幅 % 死区时间 T_dead 2e-6; % 2微秒把这段脚本放在模型的 InitFcn 回调里或者直接在命令行跑一遍模型里的变量就都有值了。4.2 主电路搭建主电路用 Simulink 的 Simscape Electrical Specialized Power Systems 库这套库在各版本 MATLAB 里都很稳定兼容性最好。需要拉的模块AC Voltage Source设置峰值 311V、频率 50Hz或直接填有效值 220V 并勾选“Peak amplitude”选项区分。Series RLC Branch当作电感 L 用R 填 R_LL 填 LC 设 0。Universal Bridge选单相桥Single-phase bridge开关器件选 IGBT/Diodes这个模块自带反并联二极管模拟整流器再合适不过。Series RLC Branch当负载电阻用C 填 0R 填 R_load。直流母线电容用另一个 Series RLC BranchR 填 inf 或者很大L 填 0C 填 C。注意电容初始电压可以设成 400这样仿真启动直接进入稳态附近调试控制参数时非常省时间。测量模块Current Measurement 串在电感和桥臂之间Voltage Measurement 测量母线电压和交流电源电压。别忘了放一个 Powergui 块连续仿真选 Continuous离散仿真选 Discrete。我建议用离散离散模式加固定步长 1e-6速度和稳定性都好控制。4.3 控制子系统分层建议把控制部分拆成三个子系统PLL和坐标变换、双闭环控制器、PWM发生器。这样能单独调试每一级也方便替换调制策略做对比。PLL 部分用延时法构造正交信号把电网电压测量值同时给到两条支路一条直接输出作为 vα另一条进一个 Transport Delay 块设延迟 5ms输出作为 vβ。然后用锁相环常用的闭环结构或者直接用 Simscape 库里的 Three-Phase PLL 改成单相配置。更省事的方法是自己搭一个简单的单相 PLL把 vα、vβ 做 Park 变换q 轴电压通过 PI 调节器得到角频率修正再积分得到相位 θ。锁相环的 PI 参数我一般给 Kp_pll50、Ki_pll1000然后微调。电流环部分把电感电流延迟得到正交 β 分量和原信号一起做 Park 变换得到 id、iq。反 Park 变换得到调制波时公式是调制波 cosθ × vd - sinθ × vq这里的调制波 m 范围在 [-1,1]和三角载波幅值 Vtri1 直接比较。4.4 单极性倍频 PWM 子系统的搭建细节这是整个项目里最值得细看的部分。单极性倍频调制用同一个调制波但两个桥臂各自配一个相位相差 180° 的三角载波。在 Simulink 里用手动比较器实现最直观也不依赖特定库模块用两个 Repeating Sequence 块生成三角波频率设 fs10kHz幅值 0 到 1第二个块的波形偏移设成半周期也就是偏移向量和时间向量都平移 0.5/fs 秒这样两个载波天然错相 180°。调制波 m 经过归一化0 到 1 范围也可以但建议用 [-1,1] 配合幅值 [-1,1] 的载波后一路和载波1比较输出左桥臂上管的开关信号另一路也和载波2比较输出右桥臂上管的开关信号。上管信号取反得到下管信号然后各自加死区。死区用 On-Off Delay 块实现对上升沿延迟 T_dead下降沿不延迟。这样保证互补的两只管不会同时导通。四路信号进 Universal Bridge 的门极输入。用 MATLAB Function 块写逻辑也可以而且更容易看懂function [S1, S4, S3, S2] pwm_bipolar_multiple(m, c1, c2, Tdead) % m: 调制波 [-1,1] % c1, c2: 两路反相三角载波 S1_raw double(m c1); S3_raw double(m c2); % 下管取反 S4_raw 1 - S1_raw; S2_raw 1 - S3_raw; % 死区逻辑简化为上升沿延时(理想开关仿真下可以直接输出) S1 S1_raw; S4 S4_raw; S3 S3_raw; S2 S2_raw; end如果只做理想开关仿真验证死区那一行可以直接跳过如果要更接近真实把 Tdead 用 On-Off Delay 处理上升沿。需要说明的是两路载波错相后每个开关周期内两个桥臂的动作点错开了半个周期Uab 的等效脉冲频率就是 2fs仿真波形里观察电感电流的纹波频率能明显看到这个效果。4.5 单极性 PWM 子系统的实现单极性调制和倍频差别只在右桥臂的驱动逻辑。我的做法是保留左桥臂的高频斩波右桥臂用工频切换正半周调制波 m 大于 0让右桥臂下管常通把 B 点接到直流母线负极负半周让右桥臂上管常通把 B 点接到直流母线正极。判断极性用一个简单的 Sign 或者比较器m 0就能实现。这个方案里右桥臂两个管子只有一只在工频动作开关损耗几乎为零。对比做出来以后倍频方案在电流波形质量上的优势会一目了然但如果你算损耗单极性方案的低频管优势也很明显。5. 仿真验证与波形分析怎么判断模型算对了5.1 正常启动与稳态波形判读仿真时长设 0.5s时钟用固定步长 1e-6 秒。母线电容初始电压设 400V 的话启动第一个周期电压环就直接在稳态附近工作电流跟随指令建立起来。把示波器接四路信号母线电压 Udc、电感电流 i_s、电网电压 v_s、调制波 m。判断整流器是否正常工作先看母线电压是否稳定在 400V纹波峰峰值应该在几伏以内再看电感电流波形理想情况下是光滑正弦波相位和电网电压基本对齐功率因数接近 1。如果电流相位明显滞后或超前电压锁相环相位没调对优先检查 θ 和电网电压的对应关系。把电感电流和电网电压放大看电流过零点应该和电压过零点几乎重合。如果看到电流有反向尖峰或者明显削顶多半是调制波 m 超出载波幅值进入了过调制区这时要回头检查电流环 PI 输出限幅和母线电压参考值是否合理。5.2 单极性调制和倍频调制的波形对比两种方案在同一个模型里切换后最直观的差异是电感电流的纹波。同样工况下单极性调制的电流纹波呈明显的锯齿状峰峰值比如是 1.2A切到倍频调制后电流纹波的锯齿频率变成两倍但峰峰值会掉到 0.6A 左右。原因是等效开关频率翻倍后电感在每个脉冲周期里的电流变化量减半。这个差异也可以用 FFT 工具Powergui 里的 FFT Analysis量化电流的频谱里单极性调制在 10kHz 及其整数倍处有明显的谐波尖峰倍频调制后 10kHz 处的尖峰大幅减小主要谐波被推到 20kHz 附近。看到这个现象就可以确认倍频 PWM 逻辑搭对了。5.3 动态响应验证负载突变实验在 0.3s 时给负载电阻 R_load 做一个阶跃变化比如从 160Ω 跳到 80Ω等效功率从 1kW 跳到 2kW。观察母线电压的响应理想情况是母线电压短暂下跌几伏到十几伏然后电压环迅速增大 id_ref大约几十毫秒内母线电压回到 400V。这个过程中电流内环应该稳稳地跟踪住升高的指令电流幅值平滑过渡不该出现振荡或过冲。动态实验里最容易暴露两个问题电压环 PI 太激进会导致母线电压振荡电压环 PI 太慢会导致负载突变后母线电压长时间偏低。我调试时习惯先记录母线电压的最大跌落量和恢复时间再对照这两个指标微调 Kp_v 和 Ki_v。电流环一般不要动太多因为按极点对消法算出来的参数已经比较准了。5.4 一个容易忽略的判断点调制波的波形调制波 m 在稳态下应该是一个平滑的正弦波和电网电压同频幅值远小于载波幅值。如果看到 m 里叠加了高频毛刺或明显的畸变说明电流反馈信号里有严重的开关噪声或者坐标变换里的相位有抖动。这种时候可以在电流采样环节加一个低通滤波器截止频率设 2kHz 左右能滤掉大部分开关纹波而不影响基波跟踪。但要注意滤波器会引入相位延迟加完之后电流环的相位裕度会变差PI 参数可能需要微调。6. 常见问题与排查技巧实录从死区到版本兼容6.1 仿真发散和振荡的五张“症状对照表”现象常见原因排查方向仿真直接发散、数值爆炸步长过大或 PI 参数过大缩小固定步长到 1e-6 以内先调小 Kp 再试母线电压稳态振荡电压环带宽过高或相位裕度不足降低 Kp_v必要时把电压环带宽取到电流环的 1/15电流波形上下不对称锁相环相位偏差检查虚拟正交信号的延时是否为 5msθ 波形是否连续光滑电流纹波异常大电感太小或 PWM 逻辑错误核对比较器逻辑确认倍频下 Uab 是否出现 2fs 脉冲序列调制波 m 频繁饱和顶到 1电流环饱和或限幅过小增大电流限幅或降低负载功率排查时有个原则我记得很清楚先开环后闭环。把电流环和电压环全部断开给调制波 m 直接赋一个固定幅值的正弦波看主电路能不能正常工作。开环波形正常再逐个闭合电流环、电压环这样能快速缩小问题范围不会一上来就面对一堆互相耦合的失控信号。6.2 PI 参数调试的先后顺序新手最常见的错误是一上来就同时调两个 PI。正确顺序是先把电压环断开给定 id_ref 为常数比如 5A调电流环电流环跟踪良好之后再接上电压环调电压环参数。电流环调好的标准是 id 快速无超调地跟踪指令iq 稳态为 0电压环调好的标准是负载突变时母线电压恢复平滑、无振荡。我实际用的粗调方法是Kp_i 从计算值的一半开始如果电流跟踪慢就往大调如果出现高频振荡就往小调Ki_i 决定稳态误差如果 id 和 id_ref 之间有固定偏差就加大 Ki_i。电压环同理但 Kp_v 的调整幅度要更保守每次改动不超过 50%。6.3 MATLAB 版本兼容2022a 项目如何“降至更低版本”这个项目的标题特别提到“可将至更低版本”说明模型需要在旧版 MATLAB 里也能打开运行这是很多同学在实际答辩、迁移时会遇到的硬需求。核心原则有三个第一建模时尽量避开新版专属模块。Simscape Electrical Specialized Power Systems 里的 Universal Bridge、Series RLC Branch、Voltage Measurement、Current Measurement 这些模块在 R2015a 时代就有跨版本兼容性最好。新出的 Simscape 生态里很多模块在旧版根本没有对应物能用老库解决的尽量用老库。第二用保存为旧版本的选项。Simulink 的 Save As 对话框里可以选择保存格式版本比如 R2021b、R2020a、R2019b 等。操作路径文件菜单 → Save As → 在“保存类型”或“版本”下拉菜单里选目标版本。但要注意这个操作只对全部由兼容模块组成的模型有效如果模型里用了新版本才加入的模块旧版打开时会报错或提示替换。第三把参数全部脚本化。模型里的所有参数都引用工作区变量参数定义统一放在 .m 脚本里。这样换版本时模型文件本身不含任何版本特有的参数存储格式用户只需要在新版本里跑一次脚本模型就能直接运行。6.4 死区与离散化的工程经验倍频调制下死区的影响比普通单极性大。两路载波错相后上下桥臂的切换动作非常密集2 微秒的死区在这个频率下已经能造成输出电压的微小畸变。仿真时可以先不加死区验证控制逻辑最后再叠加死区观察波形。如果加了死区后母线电压出现低频振荡不要急着调 PI先检查死区方向和死区时间死区时间过大超过开关周期的 5%会让交流电流在过零点附近发生明显畸变这种情况优先把死区降到 1~2 微秒再考虑调整补偿策略。离散化方面控制周期和开关周期最好同步。开关频率 10kHz 时控制周期可以取 100 微秒但仿真步长要远小于控制周期取 1 微秒。如果仿真速度太慢可以把开关频率降到 5kHz再用倍频获得 10kHz 的等效纹波频率仿真时间能缩短一半左右。结尾一点个人的实操体会这个仿真项目做完我最深刻的体会是单相整流器不是三相整流器的“简化版”而是有自己的坑。单相的锁相和坐标变换需要额外构造正交信号这是第一个容易卡壳的地方等双闭环跑通了PWM 从单极性切到倍频时又会在波形细节上发现一堆新问题。建议复现这个项目的朋友不要急着追求一步到位先把开环跑通再逐级闭合电流环、电压环最后再折腾两种调制的对比。所有 PI 参数都用脚本管理所有波形都用 FFT 和示波器双确认。按照这个流程走下来不仅这个项目能顺利出结果后续无论是扩展到三相整流器、无桥 PFC还是把控制代码移植到 DSP 上你都会有一个非常扎实的仿真基础。最后再分享一个小技巧模型里把 Powergui 的模式设为离散步长 1e-6仿真 0.5s 的时长足够看稳定和动态两种工况还不会因为连续求解器把电脑卡到风扇狂转。
