单相统一功率因数控制：不平衡d-q坐标系与Simulink仿真实践

1. 从“统一功率因数”说起一个被低估的控制目标在电力电子和电机驱动的圈子里聊到变流器控制大家张口闭口都是矢量控制、直接转矩控制或者各种五花八门的现代控制理论。但今天我想聊一个听起来有点“复古”但实际上在单相系统中极具工程价值且实现起来颇有门道的控制策略——统一功率因数Unity Power Factor, UPF控制。你可能会问功率因数校正PFC不是老生常谈了吗Boost PFC电路满大街都是。没错但那通常是针对整流器前端目标是让输入电流正弦且与电压同相减少谐波。而我们今天讨论的“统一功率因数”模式其内涵更广。它针对的是单相电压源型逆变器VSI或交直变流器其核心思想是无论变流器是工作在整流模式从电网取电还是逆变模式向电网馈电都强制使其交流侧呈现为一个纯电阻特性。这意味着从电网侧看进去变流器的等效阻抗是纯阻性的其输入电流与电压不仅同相位而且波形为正弦。这就实现了在任何功率流向下的“单位功率因数”运行。为什么这个目标有价值想象一下一个单相并网逆变器比如小型的光伏逆变器或储能变流器。如果它只能以固定的功率因数比如0.8滞后运行当它需要从电网吸收少量功率时比如夜间自用电其无功电流可能会恶化本地电网的电能质量。而UPF模式让它无论在充电还是放电时都对电网“最友好”只交换有功功率不产生无功负担。再比如在一些精密实验电源或电力电子负载模拟器中UPF模式可以确保被测设备或电网不受变流器本身引入的无功和谐波干扰。实现这一目标的技术核心就是标题中提到的在不平衡d-q坐标系下的同步参考帧控制。对于三相系统d-q变换是标准操作但对于单相系统我们只有一个物理量如何构造出正交的“两相”信号来进行d-q变换这就是第一个技术门槛。而“不平衡d-q控制”则是为了应对单相系统本身不对称的特性对传统三相d-q理论进行的适配和修正。最后正弦PWM是实现控制指令的最终执行者。整个技术链条从信号构造、坐标变换、控制器设计到PWM生成环环相扣在Simulink中搭建模型进行仿真是理解和验证这套控制逻辑最高效、最直观的途径。这篇文章我就结合自己的仿真实践把这套技术的里里外外、坑坑洼洼给你一次讲透。2. 单相系统的“两相”魔术如何构造正交信号在深入控制器之前我们必须解决单相系统的先天不足它只有一相电压和电流信号比如Va Ia而d-q变换需要α-β两相静止坐标系下的信号作为输入。因此我们的第一个任务就是从一个实时的单相信号中实时地生成与其正交的虚拟信号。2.1 正交信号生成器的几种玩法这不是什么新问题工程师们想出了不少办法各有优劣。方法一基于传输延迟的T/4延时法这是最直观的方法。对于一个基波频率为f的正弦信号其正交信号理论上就是延迟四分之一周期T/4后的信号。在数字控制中我们可以用一个深度为N/4的环形缓冲区FIFO来实现其中N是一个基波周期内的采样点数。% 伪代码示例离散实现T/4延时 persistent buffer; % 声明一个持久化数组作为缓冲区 if isempty(buffer) buffer zeros(N/4, 1); % 初始化缓冲区 end beta buffer(1); % 缓冲区最“老”的数据就是延迟了T/4的信号 % 更新缓冲区将新数据放入末尾最老的数据移出 buffer [buffer(2:end); alpha];注意这种方法简单但在频率变化时固定的延时长度会导致相位差不再是精确的90度从而在d-q变换中引入误差。它适用于电网频率非常稳定的场合。方法二基于二阶广义积分器的SOGI-PLL这是目前更主流、性能更优的方案。SOGISecond-Order Generalized Integrator本质上是一个具有频率自适应能力的正交信号发生器。它的结构是一个带中心频率ω的谐振器当输入信号频率等于ω时它的两个输出恰好是正交的。 其传递函数可以表示为Vα(s) (kω‘ s) / (s² kω‘ s ω‘²) * V(s) Vβ(s) (kω‘ ω‘) / (s² kω‘ s ω‘²) * V(s)其中V(s)是输入单相信号Vα和Vβ是生成的正交输出k是阻尼系数影响带宽和响应速度。 在Simulink中你可以用基本的积分器、增益和加法器搭建它也可以直接搜索“SOGI”相关的用户提交模块或使用Simscape Electrical库中的PLL模块其内部可能包含类似结构。方法三基于Hilbert变换的数字滤波器Hilbert变换可以在理论上完美地产生一个信号的90度相移版本。在离散时间域可以通过设计一个FIR或IIR滤波器来近似实现Hilbert变换。MATLAB的hilbert函数返回解析信号其虚部就是正交分量很方便但它是非因果的适用于离线分析。实时实现需要使用因果的FIR滤波器这会引入固定的群延时需要在控制环路中予以补偿。我的选择与考量在追求高精度和频率自适应的仿真中我强烈推荐使用SOGI-PLL方案。它不仅生成了高质量的正交信号Vα, Vβ其内部的PLL还能同时提取出电网电压的相位角θ和频率ω这个θ正是我们后续进行同步坐标变换所必需的关键信息。一个模块解决了两个核心问题性价比极高。在Simulink里我通常会把它封装成一个子系统清晰明了。2.2 从α-β到d-q同步旋转坐标变换一旦我们有了静止坐标系下的两相正交信号Xα, Xβ就可以通过派克Park变换将它们转换到同步旋转坐标系d-q下。变换公式如下Xd Xα * cosθ Xβ * sinθ Xq -Xα * sinθ Xβ * cosθ这里的θ就是上一步SOGI-PLL锁定的电网电压相位角。这个变换的物理意义非常美妙对于一个与电网电压同步旋转的正弦量在d-q坐标系下它被转换成了直流量。例如与电网电压同相位的分量会体现在d轴上而正交的分量体现在q轴上。对于我们的单相统一功率因数控制目标而言d轴电流Id代表有功电流分量。Id 0 表示变流器从电网吸收有功功率整流模式Id 0 表示变流器向电网输送有功功率逆变模式。q轴电流Iq代表无功电流分量。Iq 0 就是我们追求的统一功率因数状态意味着没有无功交换。因此控制目标简化为通过调节变流器的开关状态使得交流侧电流的q轴分量Iq恒为零同时d轴分量Id跟踪其有功功率参考值Id_ref。Id_ref由外环如直流母线电压控制环给出。3. “不平衡”d-q控制的内涵与PI调节器设计标题中提到的“不平衡d-q控制”可能会让人困惑。在三相平衡系统中d轴和q轴的控制环路通常是对称的使用相同的PI参数。但在单相系统中经过上述正交构造和变换后d轴和q轴所控制的物理对象在动态特性上并非完全对称尤其是在考虑数字控制延时、PWM延时以及系统参数偏差时。3.1 为什么需要区别对待d轴和q轴信号来源的非理想性我们通过SOGI等算法构造的β分量Vβ, Iβ并非直接测量得到其动态响应和噪声特性可能与原始的α分量略有不同。这种非理想性在变换到d-q轴后可能会在两个轴上表现出微小的差异。控制目标的优先级差异在UPF控制中q轴电流控制维持为零是保证功率因数的核心其动态性能要求极高需要快速抑制任何无功扰动。而d轴电流控制主要跟踪功率指令其动态可以稍慢但稳态精度要求高。为q轴设计带宽更高、响应更快的PI调节器有时是必要的。系统耦合与解耦即使在旋转坐标系下d轴和q轴之间仍然存在由旋转电动势引起的交叉耦合项ωL * Iq 和 ωL * Id。标准的解耦控制会将这些项作为前馈补偿加入。但在参数不准确电感L估算误差或存在扰动时两个环路的实际耦合程度可能不同因此对PI参数的鲁棒性要求也可能不同。因此“不平衡控制”更准确的理解是根据d轴和q轴的实际控制特性与性能要求分别独立地设计和整定其PI调节器参数而不是机械地使用同一套参数。3.2 电流内环PI参数设计实战电流内环是控制性能的关键。我们通常将其建模为一阶惯性环节加纯延时环节。被控对象是变流器的交流侧电感L和电阻R。Gp(s) 1 / (Ls R)考虑数字计算延时一个采样周期Ts和PWM等效延时通常近似为0.5个开关周期总的小惯性时间常数TΣ ≈ 1.5 * TswTsw为开关周期。采用典型的“零极点对消”法设计PI调节器PI传递函数Gc(s) Kp Ki/s Kp * (1 1/(Ti s)) 其中 Ti Kp / Ki。设计思路让PI调节器的零点-1/Ti去对消被控对象的大惯性极点-R/L。即令 Ti L / R。确定Kp对消后开环传递函数简化为 Kp / (s * L)。为了获得期望的闭环带宽ωc根据幅频特性在穿越频率ωc处开环增益应为1。因此有 |Kp / (jωc * L)| 1 推导出Kp ωc * L。参数调整这是一个理论起点。在Simulink中你需要进行扫参仿真。通常我会先为d轴和q轴设定相同的初始参数基于上述计算。然后通过观察阶跃响应重点关注q轴给q轴电流参考一个阶跃扰动从0到一个小值再回0观察其跟踪速度和超调。为了确保无功的快速抑制可以适当提高q轴的Kp增加带宽或者减小Ti使零点频率更高响应更快。验证d轴给d轴电流参考一个阶跃观察其跟踪和抗负载扰动的能力。d轴更注重稳态无差可以保持或微调参数。耦合测试改变d轴参考观察是否引起q轴的波动反之亦然。如果耦合严重可能需要检查解耦前馈项的准确性而非一味调整PI。在Simulink中搭建这个双闭环控制结构时清晰的信号命名和子系统封装至关重要。我会把d轴和q轴的PI调节器做成两个独立的模块方便分别调整参数。4. Simulink仿真建模从理论到可视化的关键步骤理论分析完毕是时候在Simulink里“搭积木”了。一个可靠的仿真模型是验证想法、调试参数、预测性能的必备工具。4.1 主电路建模要点单相全桥拓扑这是最常用的单相VSI/变流器拓扑。使用Simscape Electrical库中的“Mosfet”或“IGBT”模块搭建注意要并联续流二极管。直流侧接一个电容和直流电压源或负载交流侧通过一个L或LCL滤波器连接到单相交流电压源。参数设置务求真实开关频率fsw设为10kHz或20kHz这是一个常见的范围。滤波电感L根据开关纹波电流要求计算。公式ΔI (Vdc * Duty) / (L * fsw)。通常选取使纹波在额定电流的20%-30%以内。直流母线电容Cdc其大小影响直流电压的纹波和动态响应。可根据允许的电压纹波ΔVdc来计算C ≈ (P * Δt) / (Vdc * ΔVdc)其中P为功率Δt为工频半周期内电容放电时间。电网电压源设置为你所在地区的标准电压和频率如220Vrms, 50Hz。测量与传感准确测量交流侧电压Va、电流Ia以及直流母线电压Vdc。使用“Voltage Sensor”和“Current Sensor”模块。非常重要的一点为了模拟真实传感器的延迟和噪声可以在测量信号后添加一个一阶低通滤波器截止频率远高于基频但低于开关频率和小幅度的高斯白噪声这能让你的控制算法更具鲁棒性。4.2 控制算法子系统搭建这是模型的核心建议分层封装Level 1: 信号处理与坐标变换层封装SOGI-PLL子系统输入Va输出sinθ, cosθ, ω。封装Clark Park变换子系统。输入Ia, 从SOGI-PLL来的sinθ/cosθ。输出Id, Iq。同样对电网电压Va进行变换得到Vd, Vq主要用于前馈和解耦计算。Level 2: 电流环控制层分别搭建d轴和q轴的PI调节器。输入Id_ref, Iq_ref0与反馈值Id, Iq的误差。实现解耦和前馈这是提升性能的关键。计算电压前馈项Ud_ff Vd - ω * L * Iq Uq_ff Vq ω * L * Id将PI调节器的输出Ud_pi, Uq_pi分别与这些前馈项相加得到最终的d-q坐标系下的电压指令Ud和Uq。注意这里的电感L值必须与主电路实际值一致否则解耦不彻底会导致控制环路耦合振荡。Level 3: 反变换与PWM生成层进行反Park变换将Ud*, Uq转换回α-β静止坐标系下的电压指令Uα, Uβ*。Uα* Ud* * cosθ - Uq* * sinθ Uβ* Ud* * sinθ Uq* * cosθ正弦PWMSPWM实现将Uα*因为单相通常只用Uα*作为调制波与一个三角载波进行比较生成驱动全桥上下桥臂的PWM信号。注意设置正确的死区时间Dead Time这是硬件保护的必须项在Simulink的PWM Generator模块中可以设置。对于单相全桥需要两路互补的PWM信号。Level 4: 外环直流电压环这是一个更外层的控制环用于稳定直流母线电压。测量Vdc与参考值Vdc_ref比较误差经过一个PI调节器带宽远低于电流内环通常为10-20Hz其输出作为内环d轴电流的参考值Id_ref。这个环决定了变流器是吸收还是发出有功功率。4.3 仿真配置与调试技巧求解器选择对于这种包含快速开关动作的电力电子仿真使用变步长Variable-step求解器是更高效的选择。推荐使用ode23tb或ode15s刚性方程求解器它们能更好地处理系统刚性问题。将最大步长Max step size限制在开关周期的1/50到1/100例如对于10kHz开关频率设为2e-6秒以保证PWM细节的准确性。调试流程第一步开环验证。先让控制器输出固定的占空比检查主电路功率流向是否正确传感器测量是否正常。第二步仅投入电流内环。将Id_ref和Iq_ref设为固定值断开外环。观察Id, Iq能否快速、准确地跟踪参考值。这是整定d轴和q轴PI参数的关键阶段。使用“PID Tuner”工具或手动扫参。第三步投入外环。整定外环PI参数观察直流电压的建立过程和动态响应如突加负载时的电压跌落与恢复。第四步统一功率因数验证。在稳态和动态过程中观察电网电压Va和电流Ia的波形。使用“Powergui”模块的FFT分析工具计算总谐波畸变率THD和功率因数。真正的UPF模式下Va和Ia应该同相位且THD很低5%。关键波形观测除了Va和Ia一定要在示波器上同时观测Id, Iq, Vd, Vq, 调制波以及PWM信号。通过观察Id和Iq的稳态值Iq应为0和动态响应你能最直接地判断控制器是否工作正常。5. 仿真结果分析与常见问题排坑经过一番调试模型应该能跑起来了。我们来看看成功的仿真应该是什么样子以及那些年我踩过的坑。5.1 成功的仿真波形解读运行一个完整的仿真比如让变流器先整流稳定直流电压然后直流侧负载突增观察动态过程。稳态波形电网侧Va和Ia应该是同频率、同相位的完美正弦波Ia的幅值由负载决定。功率因数表应显示为1.0或无限接近1。d-q坐标系Id稳定在一个正值整流或负值逆变Iq应该在零值附近有非常微小的纹波主要由开关谐波引起其平均值应为0。直流侧Vdc应稳定在参考值有微小的开关频率纹波和两倍工频纹波对于单相系统这是固有的无法完全消除只能通过增大电容减小。动态过程波形负载突增直流电压Vdc会瞬间有一个跌落外环PI迅速响应增大Id_ref。内环电流控制器快速动作使Id跟踪新的参考值电网电流Ia的幅值随之增大。在整个动态过程中应能观察到Iq始终被牢牢地控制在零附近没有明显的偏移这证明了q轴控制器的快速性和鲁棒性。Va和Ia的相位始终保持同步。5.2 高频振荡与不收敛问题排查清单如果你的仿真出现了振荡、发散或者功率因数很差请按以下清单排查PI参数问题最常见现象Id或Iq持续振荡频率较高接近或高于电流环带宽。排查首先将d轴和q轴的PI参数都设小比如Kp0.1, Ki1观察系统是否稳定。如果稳定再缓慢增大Kp。特别注意q轴的积分系数Ki过大的Ki会引起超调和振荡。可以先整定比例部分Kp待响应快速无超调后再加入较小的Ki消除静差。解耦前馈不准确现象改变Id_ref时Iq出现明显的扰动或者改变运行点功率大小时控制性能变差。排查检查解耦公式中的电感L值是否与主电路模型中的电感值完全一致。在仿真中你可以尝试将解耦项暂时去掉只用PI控制观察耦合是否严重。如果严重则说明解耦是必要的并需核对参数。SOGI-PLL失锁或动态慢现象电网电压轻微畸变或频率变化时d-q变换后的信号出现低频振荡导致控制紊乱。排查检查SOGI-PLL的输出相位θ是否平滑增长360度/周期。调整SOGI的增益k和PLL的PI参数。k值影响带宽k越大滤波效果越差但响应越快PLL的带宽应低于电网频率但远高于其可能的变化率。传感器与PWM延时未补偿现象系统在理论计算下应该稳定但仿真中在较高带宽时出现不稳定。排查在控制算法中你是否考虑了数字控制的一个采样周期延时一个简单的补偿方法是在电压前馈指令Ud*, Uq*上乘以一个相位超前补偿因子 e^(s * Td)其中Td是总延时通常为1.5倍开关周期。在离散域这可以通过预测控制或状态观测器来实现。在初步仿真中可以尝试略微降低电流环的带宽来规避这个问题。死区时间的影响现象在轻载或电流过零点附近电流波形畸变THD增大功率因数偏离1。排查死区时间会导致实际输出电压与理想调制波之间产生误差这个误差与电流方向有关。在追求高性能的UPF控制中可能需要加入死区补偿算法。在仿真中你可以先尝试减小或取消死区时间观察波形是否改善以确认是否是死区引起的问题。搭建和调试这样一个完整的单相UPF变流器Simulink模型是对电力电子和控制理论的一次深度实践。它强迫你去思考每一个环节的物理意义和数学本质从信号构造的近似处理到坐标变换的物理映射再到多环控制器的参数耦合与解耦。当你在示波器上看到电网电压和电流完美同相的正弦波时那种成就感是看多少篇论文都无法替代的。这个模型也是一个强大的试验平台你可以在此基础上尝试将PI控制器替换为PR比例谐振控制器、重复控制器以改善对特定谐波的抑制或者尝试模型预测控制等更先进的方法探索性能的边界。