CN104701845A - 基于虚拟磁链的不平衡电网下pwm整流器直接功率控制方法 - Google Patents

Abstract

本发明公开了一种基于虚拟磁链的不平衡电网下PWM整流器直接功率控制方法，属于电力电子技术领域。所述控制方法包括以下步骤：采集三相电网电流和直流母线电压，并结合开关管信号计算出虚拟磁链；对虚拟磁链和电网电流进行坐标变换，并分离出正负序分量；根据电网电流和虚拟磁链计算出瞬时功率和瞬时功率误差；计算整流器参考电压，并进行坐标反变换；对整流器参考电压进行限幅后利用空间矢量调制方法，产生PWM整流器的开关信号。本发明在有效抑制系统有功功率和直流侧电压波动、减小电网电流畸变程度的同时，省去了电网电压传感器，降低了控制系统硬件成本，提高了控制系统可靠性，且在功率内环省去了多个PI调节器，降低了系统的复杂程度。

Description

基于虚拟磁链的不平衡电网下PWM整流器直接功率控制方法

技术领域

本发明涉及一种基于虚拟磁链的不平衡电网下PWM整流器直接功率控制方法，属于电力电子技术领域。

背景技术

PWM整流器具有越来越好的发展前景，其性能要求也越来越严，这就对PWM整流器的控制技术提出了更高的要求。

近年来，直接功率控制以其控制算法简单，动态响应快等优点得到了许多学者的关注。直接功率控制分成基于滞环控制和基于空间矢量调制两种方式。传统的直接功率控制采用滞环控制，它的优点是控制算法简单，但具有开关频率不固定，采样频率要求高等缺点。基于空间矢量调制的直接功率控制具有采样频率低、动态响应好、损耗小等优点，目前已经成为研究的热点。

在实际情况中，电网电压不平衡会使得PWM整流器网侧电流发生畸变，输出直流母线电压和瞬时功率产生波动，造成运行性能恶化。因此，需要发明一种简单有效的直接功率控制，提高系统运行效率，满足实际运行需求。

发明内容

本发明的目的是针对电网电压不平衡的情况下，提出了一种基于虚拟磁链的不平衡电网下PWM整流器直接功率控制方法，该方法可以有效抑制电网电压不平衡时直流母线电压和瞬时功率的波动，减小电网电流谐波，以满足电网运行要求。

一种基于虚拟磁链的不平衡电网下PWM整流器直接功率控制方法，包括下列步骤：

(1)采集三相电网电流和直流母线电压，建立两相静止坐标系，按下列公式计算虚拟磁链的α轴和β轴分量，

$\left\{\begin{array}{c}{\mathrm{\ψ}}_{\mathrm{\α}}={\∫u}_{\mathrm{\α}}\mathrm{dt}+{\mathrm{Li}}_{\mathrm{\α}}\\ {\mathrm{\ψ}}_{\mathrm{\β}}={\∫u}_{\mathrm{\β}}\mathrm{dt}+{\mathrm{Li}}_{\mathrm{\β}}\end{array}\right.$

式中

$\left\{\begin{array}{c}{u}_{\mathrm{\α}}=\frac{2}{3}{U}_{\mathrm{DC}}(S_a-\frac{1}{2}(S_b+S_c))\\ u_{\mathrm{\β}}=\frac{1}{\sqrt{3}}{U}_{\mathrm{DC}}(S_b-S_c)\end{array}\right.$

其中，u_α、u_β为整流器参考电压α轴和β轴分量，i_α和i_β为采集的电网电流α轴和β轴分量，U_DC为采集的直流母线电压，S_a、S_b、S_c为PWM整流器开关信号；

(2)虚拟磁链和电网电流正、负序分量计算

建立正、负序双旋转坐标系，将虚拟磁链和电网电流正、负序分量分别定向在正、负序旋转坐标系的d轴上，d^p和q^p为正序旋转坐标系dq轴，以角速度ω逆时针旋转；dⁿ和qⁿ为负序旋转坐标系dq轴，以角速度ω顺时针旋转；为正序虚拟磁链d轴分量；为负序虚拟磁链d轴分量；和为电网正序电流d轴和q轴分量；和为电网负序电流d轴和q轴分量；θ^p和θⁿ为正、负序旋转坐标系d轴与A轴夹角；

将步骤(1)采集的电网电流和计算的虚拟磁链进行正、负序旋转坐标系变换，即可得到虚拟磁链和电网电流的正、负序分量；

(3)瞬时功率计算

按下列公式分别计算瞬时有功功率P和瞬时无功功率Q，

$\left\{\begin{array}{c}P=P_0+P_{c2}\cos ({\mathrm{\θ}}^p+{\mathrm{\θ}}^n)+P_{s2}\sin ({\mathrm{\θ}}^p+{\mathrm{\θ}}^n)\\ Q=Q_0+Q_{c2}\cos ({\mathrm{\θ}}^p+{\mathrm{\θ}}^n)+Q_{s2}\sin ({\mathrm{\θ}}^p+{\mathrm{\θ}}^n)\end{array}\right.$

式中

$\left\{\begin{array}{c}{P}_0=\frac{3}{2}(-{\mathrm{\ψ}}_q^p{i}_d^p+{\mathrm{\ψ}}_d^p{i}_q^p+{\mathrm{\ψ}}_q^n{i}_d^n-{\mathrm{\ψ}}_d^n{i}_q^n)\mathrm{\ω}\\ P_{c2}=\frac{3}{2}({\mathrm{\ψ}}_q^n{i}_d^p-{\mathrm{\ψ}}_d^n{i}_q^p-{\mathrm{\ψ}}_q^p{i}_d^n+{\mathrm{\ψ}}_d^p{i}_q^n)\mathrm{\ω}\\ P_{s2}=\frac{3}{2}(-{\mathrm{\ψ}}_d^n{i}_d^p-{\mathrm{\ψ}}_q^n{i}_q^p-{\mathrm{\ψ}}_d^p{i}_d^n-{\mathrm{\ψ}}_q^p{i}_q^n)\mathrm{\ω}\\ Q_0=\frac{3}{2}({\mathrm{\ψ}}_d^p{i}_d^p+{\mathrm{\ψ}}_q^p{i}_q^p-{\mathrm{\ψ}}_d^n{i}_d^n-{\mathrm{\ψ}}_q^n{i}_q^n)\mathrm{\ω}\\ Q_{c2}=\frac{3}{2}(-{\mathrm{\ψ}}_d^n{i}_d^p-{\mathrm{\ψ}}_q^n{i}_q^p+{\mathrm{\ψ}}_d^p{i}_d^n+{\mathrm{\ψ}}_q^p{i}_q^n)\mathrm{\ω}\\ Q_{s2}=\frac{3}{2}(-{\mathrm{\ψ}}_q^n{i}_d^p+{\mathrm{\ψ}}_d^n{i}_q^p-{\mathrm{\ψ}}_q^p{i}_d^n+{\mathrm{\ψ}}_d^p{i}_q^n)\mathrm{\ω}\end{array}\right.$

其中，P₀、Q₀为网侧瞬时有功、无功功率直流分量；P_c2、Q_c2为以余弦规律变化的网侧瞬时有功、无功功率二次脉动分量幅值；P_s2、Q_s2为以正弦规律变化的网侧瞬时有功、无功功率二次脉动分量幅值；ω为虚拟磁链角速度；

(4)瞬时功率误差计算

将期望得到的直流母线电压平方值与步骤(1)采集到的直流母线电压平方值进行做差，将该差值经PI调节器得到瞬时有功功率P₀的参考值；P_c2、P_s2和Q₀的参考值为零，则按下列公式计算瞬时功率误差为：

$\left\{\begin{array}{c}{\mathrm{\ΔP}}_0=P_0^{*}-P_0\\ {\mathrm{\ΔQ}}_0=Q_0^{*}-Q_0\\ {\mathrm{\ΔP}}_{s2}=P_{s2}^{*}-P_{s2}\\ {\mathrm{\ΔP}}_{c2}=P_{c2}^{*}-P_{c2}\end{array}\right.$

其中，ΔP₀、ΔQ₀、ΔP_s2和ΔP_c2分别为P₀、Q₀、P_s2和P_c2的参考值与实际值的误差；

(5)整流器参考电压计算

根据步骤(2)计算的虚拟磁链和电网电流正、负序分量和步骤(4)计算的瞬时功率误差，按下式计算整流器参考电压：

$\left\{\begin{array}{c}{u}_d^p=\frac{2L}{3\mathrm{\ω}{T}_s}\frac{{\mathrm{\ΔP}}_{s2}-{\mathrm{\ΔQ}}_0}{[{\left({\mathrm{\ψ}}_d^p\right)}^2+{\left({\mathrm{\ψ}}_d^n\right)}^2]}-{\mathrm{\ω\ψ}}_q^p+{\mathrm{\ωLi}}_q^p\\ u_d^n=\frac{2L}{3\mathrm{\ω}{T}_s}\frac{{\mathrm{\ΔP}}_{s2}+{\mathrm{\ΔQ}}_0}{[{\left({\mathrm{\ψ}}_d^p\right)}^2+{\left({\mathrm{\ψ}}_d^n\right)}^2]}+{\mathrm{\ω\ψ}}_q^n-{\mathrm{\ωLi}}_q^n\\ u_q^p=\frac{2L}{3\mathrm{\ω}{T}_s}\frac{{\mathrm{\ΔP}}_{c2}+{\mathrm{\ΔP}}_0}{[{\left({\mathrm{\ψ}}_d^p\right)}^2-{\left({\mathrm{\ψ}}_d^n\right)}^2]}+{\mathrm{\ω\ψ}}_d^p-{\mathrm{\ωLi}}_d^p\\ u_q^n=\frac{2L}{3\mathrm{\ω}{T}_s}\frac{{\mathrm{\ΔP}}_{c2}+{\mathrm{\ΔP}}_0}{[{\left({\mathrm{\ψ}}_d^n\right)}^2-{\left({\mathrm{\ψ}}_d^p\right)}^2]}-{\mathrm{\ω\ψ}}_d^n+{\mathrm{\ωLi}}_d^n\end{array}\right.$

$\left\{\begin{array}{c}{u}_{\mathrm{dq}}^p=u_d^p+{\mathrm{ju}}_q^p\\ u_{\mathrm{dq}}^n=u_d^n+{\mathrm{ju}}_q^n\end{array}\right.$

式中，和为整流器正序参考电压d、q轴分量；和为整流器负序参考电压d、q轴分量；和为两相旋转坐标系整流器正、负参考电压合成矢量；T_s为系统采样周期；L为电网侧输入电感值；

(6)对步骤(5)计算的整流器正、负序参考电压进行坐标反变换，并对坐标反变换后的整流器参考电压进行限幅，计算公式如下：

$\left\{\begin{array}{c}{u}_{\mathrm{\α\β}}^P=e^{j{\mathrm{\θ}}^p}{u}_{\mathrm{dq}}^p\\ u_{\mathrm{\α\β}}^n=e^{-j{\mathrm{\θ}}^n}{u}_{\mathrm{dq}}^n\end{array}\right.$

$\left\{\begin{array}{c}{u}_{\mathrm{\α}}=u_{\mathrm{\α}}^p+u_{\mathrm{\α}}^n\\ u_{\mathrm{\β}}=u_{\mathrm{\β}}^p+u_{\mathrm{\β}}^n\end{array}\right.$

式中，和为整流器正序参考电压α轴和β轴分量；和为整流器负序参考电压α轴和β轴分量；和为两相静止坐标系整流器正、负参考电压合成矢量；u_α和u_β为两相静止坐标系整流器参考电压α轴和β轴分量；

(7)采用空间矢量调制算法对限幅后的整流器参考电压进行调制后产生PWM整流器开关信号，并将开关信号通过驱动电路驱动功率开关器件。

综上所述，在电网电压不平衡的情况下，本发明在有效抑制直流母线电压和瞬时有功功率波动、减小电网电流谐波的同时，省去了电网电压传感器，增强了系统可靠性，降低了系统硬件成本，并且在功率内环省去了多个PI调节器，简化了系统结构，降低了控制系统的复杂程度。

附图说明

图1为三相PWM整流器主电路拓扑结构图；

图2为系统整体控制框图；

图3为虚拟磁链定向矢量图；

图4为采用传统控制方法时在电网不平衡情况下三相电网电流、直流母线电压仿真结果；

图5为采用本发明所提的控制方法时在电网不平衡情况下三相电网电流、直流母线电压仿真结果；

具体实施方式

下面结合附图和实施例对本发明作进一步说明。

图1为三相PWM整流器主电路拓扑结构。如图所示，e_a、e_b、e_c为电网电压；i_a、i_b、i_c为电网侧相电流；u_a、u_b、u_c为整流器参考相电压；U_DC为直流母线电压；L为电网侧滤波电感；R为电网侧输入电阻；R_L为负载电阻。

图2为系统整体控制框图，该控制方法在图1的三相PWM整流器主电路拓扑结构上按照如下步骤依次实现：

1数据采集

控制系统通过电流霍尔传感器采集三相电网电流值i_a、i_b和i_c，通过电压霍尔传感器采集到直流母线电压U_DC。其中，通过三相电网电流和直流母线电压估算出电网侧虚拟磁链，并且计算出PWM整流器瞬时功率，通过直流母线电压平方值与直流母线参考电压平方值做差后经PI调节器输出瞬时有功功率参考值。模拟量采集需选用专用的AD转换芯片。在本实施例中，选用飞思卡尔公司的MC56F8346控制芯片，使用该控制芯片内部的AD转换接口采集模拟量。

采集三相电网电流和直流母线电压，建立两相静止坐标系，按下列公式计算虚拟磁链的α轴和β轴分量：

$\left\{\begin{array}{c}{\mathrm{\ψ}}_{\mathrm{\α}}={\∫u}_{\mathrm{\α}}\mathrm{dt}+{\mathrm{Li}}_{\mathrm{\α}}\\ {\mathrm{\ψ}}_{\mathrm{\β}}={\∫u}_{\mathrm{\β}}\mathrm{dt}+{\mathrm{Li}}_{\mathrm{\β}}\end{array}\right.$

式中

$\left\{\begin{array}{c}{u}_{\mathrm{\α}}=\frac{2}{3}{U}_{\mathrm{DC}}(S_a-\frac{1}{2}(S_b+S_c))\\ u_{\mathrm{\β}}=\frac{1}{\sqrt{3}}{U}_{\mathrm{DC}}(S_b-S_c)\end{array}\right.$

其中，u_α、u_β为整流器参考电压α轴和β轴分量，i_α和i_β为采集的电网电流α轴和β轴分量，U_DC为采集的直流母线电压，S_a、S_b、S_c为PWM整流器开关信号；

2虚拟磁链和电网电流正、负序分量计算

建立正、负序双旋转坐标系，将虚拟磁链和电网电流正、负序分量分别定向在正、负序旋转坐标系的d轴上，则系统虚拟磁链和电网电流矢量图如图3所示。在图3中，A、B、C为三相静止坐标轴；d^p和q^p为正序旋转坐标系dq轴，以角速度ω逆时针旋转；dⁿ和qⁿ为负序旋转坐标系dq轴，以角速度ω顺时针旋转；为正序虚拟磁链d轴分量；为负序虚拟磁链d轴分量；和为电网正序电流d轴和q轴分量；和为电网负序电流d轴和q轴分量；θ^p和θⁿ为正、负序旋转坐标系d轴与A轴夹角；

将虚拟磁链正负分量分别定向于正负旋转坐标系的d轴，则 和为常量。

3瞬时功率计算

根据虚拟磁链和电网电流正负序分量，计算瞬时有功功率P和瞬时无功功率Q：

$\left\{\begin{array}{c}P=P_0+P_{c2}\cos ({\mathrm{\θ}}^p+{\mathrm{\θ}}^n)+P_{s2}\sin ({\mathrm{\θ}}^p+{\mathrm{\θ}}^n)\\ Q=Q_0+Q_{c2}\cos ({\mathrm{\θ}}^p+{\mathrm{\θ}}^n)+Q_{s2}\sin ({\mathrm{\θ}}^p+{\mathrm{\θ}}^n)\end{array}\right.$

式中

$\left\{\begin{array}{c}{P}_0=\frac{3}{2}(-{\mathrm{\ψ}}_q^p{i}_d^p+{\mathrm{\ψ}}_d^p{i}_q^p+{\mathrm{\ψ}}_q^n{i}_d^n-{\mathrm{\ψ}}_d^n{i}_q^n)\mathrm{\ω}\\ P_{c2}=\frac{3}{2}({\mathrm{\ψ}}_q^n{i}_d^p-{\mathrm{\ψ}}_d^n{i}_q^p-{\mathrm{\ψ}}_q^p{i}_d^n+{\mathrm{\ψ}}_d^p{i}_q^n)\mathrm{\ω}\\ P_{s2}=\frac{3}{2}(-{\mathrm{\ψ}}_d^n{i}_d^p-{\mathrm{\ψ}}_q^n{i}_q^p-{\mathrm{\ψ}}_d^p{i}_d^n-{\mathrm{\ψ}}_q^p{i}_q^n)\mathrm{\ω}\\ Q_0=\frac{3}{2}({\mathrm{\ψ}}_d^p{i}_d^p+{\mathrm{\ψ}}_q^p{i}_q^p-{\mathrm{\ψ}}_d^n{i}_d^n-{\mathrm{\ψ}}_q^n{i}_q^n)\mathrm{\ω}\\ Q_{c2}=\frac{3}{2}(-{\mathrm{\ψ}}_d^n{i}_d^p-{\mathrm{\ψ}}_q^n{i}_q^p+{\mathrm{\ψ}}_d^p{i}_d^n+{\mathrm{\ψ}}_q^p{i}_q^n)\mathrm{\ω}\\ Q_{s2}=\frac{3}{2}(-{\mathrm{\ψ}}_q^n{i}_d^p+{\mathrm{\ψ}}_d^n{i}_q^p-{\mathrm{\ψ}}_q^p{i}_d^n+{\mathrm{\ψ}}_d^p{i}_q^n)\mathrm{\ω}\end{array}\right.$

其中，P₀、Q₀为网侧瞬时有功、无功功率直流分量；P_c2、Q_c2为以余弦规律变化的网侧瞬时有功、无功功率二次脉动分量幅值；P_s2、Q_s2为以正弦规律变化的网侧瞬时有功、无功功率二次脉动分量幅值；ω为虚拟磁链角速度；

4瞬时功率误差计算

将期望得到的直流母线电压平方值与采集到的直流母线电压平方值进行做差，将该差值经PI调节器得到瞬时有功功率P₀的参考值；P_c2、P_s2和Q₀的参考值为零，则按下列公式计算瞬时功率误差为：

$\left\{\begin{array}{c}{\mathrm{\ΔP}}_0=P_0^{*}-P_0\\ {\mathrm{\ΔQ}}_0=Q_0^{*}-Q_0\\ {\mathrm{\ΔP}}_{s2}=P_{s2}^{*}-P_{s2}\\ {\mathrm{\ΔP}}_{c2}=P_{c2}^{*}-P_{c2}\end{array}\right.$

其中，ΔP₀、ΔQ₀、ΔP_s2和ΔP_c2分别为P₀、Q₀、P_s2和P_c2的参考值与实际值的误差；

5整流器参考电压计算

根据计算的虚拟磁链和电网电流正、负序分量和瞬时功率误差，按下式计算整流器参考电压：

$\left\{\begin{array}{c}{u}_d^p=\frac{2L}{3\mathrm{\ω}{T}_s}\frac{{\mathrm{\ΔP}}_{s2}-{\mathrm{\ΔQ}}_0}{[{\left({\mathrm{\ψ}}_d^p\right)}^2+{\left({\mathrm{\ψ}}_d^n\right)}^2]}-{\mathrm{\ω\ψ}}_q^p+{\mathrm{\ωLi}}_q^p\\ u_d^n=\frac{2L}{3\mathrm{\ω}{T}_s}\frac{{\mathrm{\ΔP}}_{s2}+{\mathrm{\ΔQ}}_0}{[{\left({\mathrm{\ψ}}_d^p\right)}^2+{\left({\mathrm{\ψ}}_d^n\right)}^2]}+{\mathrm{\ω\ψ}}_q^n-{\mathrm{\ωLi}}_q^n\\ u_q^p=\frac{2L}{3\mathrm{\ω}{T}_s}\frac{{\mathrm{\ΔP}}_{c2}+{\mathrm{\ΔP}}_0}{[{\left({\mathrm{\ψ}}_d^p\right)}^2-{\left({\mathrm{\ψ}}_d^n\right)}^2]}+{\mathrm{\ω\ψ}}_d^p-{\mathrm{\ωLi}}_d^p\\ u_q^n=\frac{2L}{3\mathrm{\ω}{T}_s}\frac{{\mathrm{\ΔP}}_{c2}+{\mathrm{\ΔP}}_0}{[{\left({\mathrm{\ψ}}_d^n\right)}^2-{\left({\mathrm{\ψ}}_d^p\right)}^2]}-{\mathrm{\ω\ψ}}_d^n+{\mathrm{\ωLi}}_d^n\end{array}\right.$

$\left\{\begin{array}{c}{u}_{\mathrm{dq}}^p=u_d^p+{\mathrm{ju}}_q^p\\ u_{\mathrm{dq}}^n=u_d^n+{\mathrm{ju}}_q^n\end{array}\right.$

式中，和为整流器正序参考电压d、q轴分量；和为整流器负序参考电压d、q轴分量；和为两相旋转坐标系整流器正、负参考电压合成矢量；T_s为系统采样周期；L为电网侧输入电感值；

6对整流器正、负序参考电压进行坐标反变换，计算公式如下：

$\left\{\begin{array}{c}{u}_{\mathrm{\α\β}}^P=e^{j{\mathrm{\θ}}^p}{u}_{\mathrm{dq}}^p\\ u_{\mathrm{\α\β}}^n=e^{-j{\mathrm{\θ}}^n}{u}_{\mathrm{dq}}^n\end{array}\right.$

$\left\{\begin{array}{c}{u}_{\mathrm{\α}}=u_{\mathrm{\α}}^p+u_{\mathrm{\α}}^n\\ u_{\mathrm{\β}}=u_{\mathrm{\β}}^p+u_{\mathrm{\β}}^n\end{array}\right.$

式中，和为整流器正序参考电压d、q轴分量；和为整流器负序参考电压d、q轴分量；和为两相旋转坐标系整流器正、负参考电压合成矢量；T_s为系统采样周期；L为电网侧输入电感值；

为防止计算出的u_α、u_β超出整流器输出最高限定电压，因此需要对坐标反变换后的整流器参考电压进行限幅。

7采用空间矢量调制算法对限幅后的整流器参考电压进行调制后产生PWM整流器开关信号，并将开关信号通过驱动电路驱动功率开关器件。

本发明提出控制方法的有效性可以通过对比图4和图5所示的仿真结果得出，二者仿真环境完全相同，区别仅在于图4采用的是传统控制策略，而图5采用的是本发明提出的控制方法。从图4看出，当电网电压不平衡时，网侧电流发生严重畸变，直流母线电压产生二次谐波分量；而从图5看出，采用本文提出的控制方法后，网侧电流正弦化程度好，直流母线电压输出恒定。

Claims (1)

1.一种基于虚拟磁链的不平衡电网下PWM整流器直接功率控制方法，其特征在于，包括以下步骤：

(1)采集三相电网电流和直流母线电压，建立两相静止坐标系，按下列公式计算虚拟磁链的α轴和β轴分量，

$\left\{\begin{array}{c}{\mathrm{\ψ}}_{\mathrm{\α}}=\∫u_{\mathrm{\α}}\mathrm{dt}+{\mathrm{Li}}_{\mathrm{\α}}\\ {\mathrm{\ψ}}_{\mathrm{\β}}=\∫u_{\mathrm{\β}}\mathrm{dt}+{\mathrm{Li}}_{\mathrm{\β}}\end{array}\right.$

式中

$\left\{\begin{array}{c}{u}_{\mathrm{\α}}=\frac{2}{3}{U}_{\mathrm{DC}}(S_a-\frac{1}{2}(S_b+S_c))\\ u_{\mathrm{\β}}=\frac{1}{\sqrt{3}}{U}_{\mathrm{DC}}(S_b-S_c)\end{array}\right.$

其中，u_α、u_β为整流器参考电压α轴和β轴分量，i_α和i_β为采集的电网电流α轴和β轴分量，U_DC为采集的直流母线电压，S_a、S_b、S_c为PWM整流器开关信号；

(2)虚拟磁链和电网电流正、负序分量计算

建立正、负序双旋转坐标系，将虚拟磁链正、负序分量分别定向在正、负序旋转坐标系的d轴上，d^p和q^p为正序旋转坐标系dq轴，以角速度ω逆时针旋转；dⁿ和qⁿ为负序旋转坐标系dq轴，以角速度ω顺时针旋转；为正序虚拟磁链d轴分量；为负序虚拟磁链d轴分量；和为电网正序电流d轴和q轴分量；和为电网负序电流d轴和q轴分量；θ^p和θⁿ为正、负序旋转坐标系d轴与A轴夹角；

将步骤(1)采集的电网电流和计算的虚拟磁链进行正、负序旋转坐标系变换，即可得到虚拟磁链和电网电流的正、负序分量；

(3)瞬时功率计算

按下列公式分别计算瞬时有功功率P和瞬时无功功率Q，

$\left\{\begin{array}{c}P=P_0+P_{c2}\cos ({\mathrm{\θ}}^p+{\mathrm{\θ}}^n)+P_{s2}\sin ({\mathrm{\θ}}^p+{\mathrm{\θ}}^n)\\ Q=Q_0+Q_{c2}\cos ({\mathrm{\θ}}^p+{\mathrm{\θ}}^n)+Q_{s2}\sin ({\mathrm{\θ}}^p+{\mathrm{\θ}}^n)\end{array}\right.$

式中

$\left\{\begin{array}{c}{P}_0=\frac{3}{2}(-{\mathrm{\ψ}}_q^p{i}_d^p+{\mathrm{\ψ}}_d^p{i}_q^p+{\mathrm{\ψ}}_q^n{i}_d^n-{\mathrm{\ψ}}_d^n{i}_q^n)\mathrm{\ω}\\ P_{c2}=\frac{3}{2}({\mathrm{\ψ}}_q^n{i}_d^p-{\mathrm{\ψ}}_d^n{i}_q^p-{\mathrm{\ψ}}_q^p{i}_d^n+{\mathrm{\ψ}}_d^p{i}_q^n)\mathrm{\ω}\\ P_{s2}=\frac{3}{2}(-{\mathrm{\ψ}}_d^n{i}_d^p-{-\mathrm{\ψ}}_q^n{i}_q^p-{\mathrm{\ψ}}_d^p{i}_d^n-{\mathrm{\ψ}}_q^p{i}_q^n)\mathrm{\ω}\\ Q_0=\frac{3}{2}({\mathrm{\ψ}}_d^p{-}_d^p+{\mathrm{\ψ}}_q^p{i}_q^p-{\mathrm{\ψ}}_d^n{i}_d^n-{\mathrm{\ψ}}_q^n{i}_q^n)\mathrm{\ω}\\ Q_{c2}=\frac{3}{2}(-{\mathrm{\ψ}}_d^n{i}_d^p-{\mathrm{\ψ}}_q^n{i}_q^p+{\mathrm{\ψ}}_d^p{i}_d^n+{\mathrm{\ψ}}_q^p{i}_q^n)\mathrm{\ω}\\ Q_{s2}=\frac{3}{2}(-{\mathrm{\ψ}}_q^n{i}_d^p+{\mathrm{\ψ}}_d^n{i}_q^p-{\mathrm{\ψ}}_q^p{i}_d^n+{\mathrm{\ψ}}_q^p{i}_d^n+{\mathrm{\ψ}}_d^p{i}_q^n)\mathrm{\ω}\end{array}\right.$

其中，P₀、Q₀为网侧瞬时有功、无功功率直流分量；P_c2、Q_c2为以余弦规律变化的网侧瞬时有功、无功功率二次脉动分量幅值；P_s2、Q_s2为以正弦规律变化的网侧瞬时有功、无功功率二次脉动分量幅值；ω为虚拟磁链角速度；

(4)瞬时功率误差计算

将期望得到的直流母线电压的平方值与步骤(1)采集到的直流母线电压的平方值进行做差，将该差值经PI调节器得到瞬时有功功率P₀的参考值；P_c2、P_s2和Q₀的参考值为零，则按下列公式计算瞬时功率误差为：

$\left\{\begin{array}{c}{\mathrm{\ΔP}}_0=P_0^{*}-P_0\\ {\mathrm{\ΔQ}}_0=Q_0^{*}-Q_0\\ {\mathrm{\ΔP}}_{s2}=P_{s2}^{*}-P_{s2}\\ {\mathrm{\ΔP}}_{c2}=p_{c2}^{*}-P_{c2}\end{array}\right.$

其中，ΔP₀、ΔQ₀、ΔP_s2和ΔP_c2分别为P₀、Q₀、P_s2和P_c2的参考值与实际值的误差；

(5)整流器参考电压计算

根据步骤(2)计算的虚拟磁链和电网电流正、负序分量和步骤(4)计算的瞬时功率误差，按下式计算整流器参考电压：

$\left\{\begin{array}{c}{u}_d^p=\frac{2L}{3\mathrm{\ω}{T}_s}\frac{{\mathrm{\ΔP}}_{s2}-{\mathrm{\ΔQ}}_0}{[{\left({\mathrm{\ψ}}_d^p\right)}^2+{\left({\mathrm{\ψ}}_d^n\right)}^2]}-{\mathrm{\ω\ψ}}_q^p+{\mathrm{\ωLi}}_q^p\\ u_d^n=\frac{2L}{3\mathrm{\ω}{T}_s}\frac{{\mathrm{\Δp}}_{s2}+{\mathrm{\ΔQ}}_0}{[{\left({\mathrm{\ψ}}_d^p\right)}^2+{\left({\mathrm{\ψ}}_d^n\right)}^2]}+{\mathrm{\ω\ψ}}_q^n-{\mathrm{\ωLi}}_q^n\\ u_q^p=\frac{2L}{{3\mathrm{\ωT}}_s}\frac{{\mathrm{\ΔP}}_{c2}+{\mathrm{\ΔP}}_0}{[{\left({\mathrm{\ψ}}_d^p\right)}^2-{\left({\mathrm{\ψ}}_d^n\right)}^2]}{+\mathrm{\ω\ψ}}_d^p-{\mathrm{\ωLi}}_d^p\\ u_q^n=\frac{2L}{{3\mathrm{\ωT}}_s}\frac{{\mathrm{\ΔP}}_{c2}+{\mathrm{\ΔP}}_0}{[{\left({\mathrm{\ψ}}_d^n\right)}^2-{\left({\mathrm{\ψ}}_d^p\right)}^2]}-{\mathrm{\ω\ψ}}_d^n+{\mathrm{\ωLi}}_d^n\end{array}\right.$

$\left\{\begin{array}{c}{u}_{\mathrm{dq}}^p=u_d^p+{\mathrm{ju}}_q^p\\ u_{\mathrm{dq}}^n=u_d^n+{\mathrm{ju}}_q^n\end{array}\right.$

式中，和为整流器正序参考电压d、q轴分量；和为整流器负序参考电压d、q轴分量；和为两相旋转坐标系整流器正、负参考电压合成矢量；T_s为系统采样周期；L为电网侧输入电感值；

(6)对步骤(5)计算的整流器正、负序参考电压进行坐标反变换，并对坐标反变换后的整流器参考电压进行限幅，计算公式如下：

$\left\{\begin{array}{c}{u}_{\mathrm{\α\β}}^P=e^{{\mathrm{j\θ}}^P}{u}_{\mathrm{dq}}^p\\ u_{\mathrm{\α\β}}^n=e^{-{\mathrm{j\θ}}^n}{u}_{\mathrm{dq}}^n\end{array}\right.$

$\left\{\begin{array}{c}{u}_{\mathrm{\α}}=u_{\mathrm{\α}}^p+u_{\mathrm{\α}}^n\\ u_{\mathrm{\β}}=u_{\mathrm{\β}}^p+u_{\mathrm{\β}}^n\end{array}\right.$

式中，和为整流器正序参考电压α轴和β轴分量；和为整流器负序参考电压α轴和β轴分量；和为两相静止坐标系整流器正、负参考电压合成矢量；u_α和u_β为两相静止坐标系整流器参考电压α轴和β轴分量；

(7)采用空间矢量调制算法对限幅后的整流器参考电压进行调制后产生PWM整流器开关信号，并将开关信号通过驱动电路驱动功率开关器件。