CN107546994B - 一种多电平逆变器并网系统及方法 - Google Patents

Abstract

本发明公开了一种多电平逆变器并网系统及方法。该系统包括主功率电路、控制电路和检测电路，其中主功率电路包括NPC三相三电平逆变器和L型低通滤波器，检测电路包括入网电流检测电路、电网电压检测电路和直流侧中点电压检测电路。方法为：首先引入校正参数λ，依据定频原则对多电平逆变器的输出电压矢量表进行校正；然后根据校正后的输出电压矢量表，对多电平逆变器的直流侧中点电压和交流侧输出电流进行两步预测控制。本发明能够在开关频率恒定的条件下，对多电平逆变器进行有限集模型预测并网控制，实现直流侧中点电压的强烈抑制和交流侧输出电流的精确跟踪。

Description

一种多电平逆变器并网系统及方法

技术领域

本发明涉及电能变换装置的直流-交流变换器技术领域，特别是一种多电平逆变器并网系统及方法。

背景技术

NPC(Neutral Point Clamped)三相三电平并网逆变器是一种广泛应用于中大容量的分布式并网发电系中的多电平逆变器，其控制策略主要有：PI控制、PR控制、滞环控制和有限集模型预测控制。其中有限集模型预测控制具有强鲁棒性、快速性、易于数字实现和考虑逆变器非线性等优点，广泛应用于多电平逆变器的并网控制之中。

然而，有限集模型预测控制本身属于非线性控制策略，应用在多电平逆变器并网控制时，具有不能形成恒定的开关频率的缺点，这不仅使得后期的逆变器并网滤波器的设计变的十分困难，而且会带来电磁兼容等问题。

发明内容

本发明的目的在于提供一种基于矢量校正式定频有限集模型预测的多电平逆变器并网系统及方法，以实现具有恒定开关频率的多电平逆变器有限集模型预测并网控制。

实现本发明目的的技术解决方案为：一种多电平逆变器并网系统，包括主功率电路、控制电路和检测电路，所述主功率电路包括输入电压源V_dc、NPC三相三电平逆变器、L型低通滤波器和三相电网e_a、e_b、e_c，其中输入电压源V_dc与NPC三相三电平逆变器输入端连接，NPC三相三电平逆变器输出端与L型低通滤波器输入端连接，L型低通滤波器输出端与三相电网e_a、e_b、e_c连接；

所述的检测电路包括入网电流检测电路、电网电压检测电路、直流侧中点电压检测电路，其中入网电流检测电路输入端与三相电网连接，入网电流检测电路输出端与控制电路第一输入端连接，电网电压检测电路输入端与三相电网连接，电网电压检测电路输出端与控制电路第二输入端连接，直流侧中点电压检测电路输入端与逆变器直流侧电容连接，直流侧中点电压检测电路输出端与控制电路第三输入端连接，控制电路的输出端接入主功率电路；

所述控制电路包括Clark转换模块和参考电流计算模块。

进一步地，所述控制电路采用DSP芯片TMS320F28335。

一种多电平逆变器并网方法，该方法基于矢量校正式定频有限集模型，采用矢量校正后的输出电压矢量表进行两步预测控制，包括以下步骤：

步骤1，校正：引入校正参数λ，依据定频原则对多电平逆变器的输出电压矢量表进行校正，得到校正电压矢量表；

步骤2，最优输出：输出当前，即第k个采样周期最优矢量S_opt(k)控制逆变器，S_opt(k)由上一个，即第k-1个采样周期计算得到；

步骤3，采样及变换：对入网电流进行检测得到a、b、c三相入网电流i_a(k)、i_b(k)、i_c(k)，并对i_a(k)、i_b(k)、i_c(k)进行Clark变换得到i_α(k)、i_β(k)；对电网电压进行检测得到a、b、c三相电网电压e_a(k)、e_b(k)、e_c(k)，并对e_a(k)、e_b(k)、e_c(k)进行Clark变换得到e_α(k)、e_β(k)；根据变换所得e_α(k)、e_β(k)，计算入网参考电流，得到αβ坐标系下的入网参考电流i^* _α(k)、i^* _β(k)；对电容电压进行检测得到直流侧中点电压Δv_c(k)；

步骤4，相角补偿：对e_α(k)、e_β(k)做一个采样周期T_s的相角补偿得到e_α(k+1)、e_β(k+1)，对i^* _α(k)、i^* _β(k)做两个采样周期2T_s的相角补偿得到

步骤5，计算入网电流反馈值：结合S_opt(k)在校正电压矢量表中的矢量坐标值

和步骤3的采样及变换结果，根据逆变器的校正离散数学模型，计算出第k+1个采样周期入网电流反馈值i_α(k+1)、i_β(k+1)；

步骤6，计算直流侧中点电压反馈值：结合S_opt(k)在校正电压矢量表中的相开关函数状态

和步骤3的采样及变换结果，根据逆变器的校正离散数学模型，计算出第k+1个采样周期直流侧中点电压反馈值Δv_c(k+1)；

步骤7，反Clark变换：对第k+1个采样周期入网电流反馈值i_α(k+1)、i_β(k+1)进行反Clark变换，得到i_a(k+1)、i_b(k+1)、i_c(k+1)；

步骤8，遍历计算：结合i_α(k+1)、i_β(k+1)、e_α(k+1)、e_β(k+1)、Δv_c(k+1)、i_a(k+1)、i_b(k+1)、i_c(k+1)、逆变器交流输出侧校正电压矢量和校正电压矢量表中各个校正电压矢量对应的相开关函数状态，根据逆变器的校正离散数学模型，计算出第k+2个采样周期的入网电流i_α(k+2)(i)、i_β(k+2)(i)和直流侧中点电压Δv_c(k+2)(i)；

步骤9，建立目标函数g，作为选出校正电压矢量表中下一个采样周期，即第k+1个采样周期最优矢量S_opt(k+1)的依据；

步骤10，存储目标函数g的最小值对应的矢量，该矢量即最优矢量S_opt(k+1)，将S_opt(k+1)对应的相开关函数状态做为第k+1个采样周期的最优状态输出；

步骤11，等待本次采样周期时间结束，返回步骤2，进入下一个寻优周期。进一步地，步骤1所述校正：引入校正参数λ，依据定频原则对多电平逆变器的输出电压矢量表进行校正，得到校正电压矢量表，具体如下：

(1.1)将逆变器输出电压矢量表中的矢量分为零矢量组、小矢量组、中矢量组、大矢量组四组；其中：

零矢量组中的矢量对应的相开关函数状态S_aS_bS_c为：-1-1-1、000、111，

小矢量组中的矢量对应的相开关函数状态S_aS_bS_c为：100、0-1-1、110、00-1、010、-10-1、011、-100、001、-1-10、101、0-10，

中矢量组中的矢量对应的相开关函数状态S_aS_bS_c为：10-1、01-1、-110、-101、0-11、1-10，

大矢量组中的矢量对应的相开关函数状态S_aS_bS_c为：1-1-1、11-1、-11-1、-111、-1-11、1-11；

(1.2)引入校正参数λ，在一个采样周期，即开关周期内对矢量进行分组校正，其中，执行时间

的矢量为辅助矢量，属于校正矢量，执行时间λT_s的矢量为主矢量，属于被校正矢量；依据定频原则，具体校正如下：

零矢量校正：000执行

-1-1-1执行λT_s、000再执行

构成校正矢量表中的1号矢量；-1-1-1执行

000执行λT_s、-1-1-1再执行

构成校正矢量表中的2号矢量；000执行

111执行λT_s、000再执行

构成校正矢量表中的3号矢量；

小矢量校正：-1-1-1执行100执行λT_s、-1-1-1再执行

构成校正矢量表中的4号矢量；111执行

0-1-1执行λT_s、111再执行构成校正矢量表中的5号矢量；-1-1-1执行110执行λT_s、-1-1-1再执行

构成校正矢量表中的6号矢量；111执行

00-1执行λT_s、111再执行

构成校正矢量表中的7号矢量；-1-1-1执行

010执行λT_s、-1-1-1再执行构成校正矢量表中的8号矢量；111执行

-10-1执行λT_s、111再执行构成校正矢量表中的9号矢量；-1-1-1执行

011执行λT_s、-1-1-1再执行

构成校正矢量表中的10号矢量；111执行

-100执行λT_s、111再执行

构成校正矢量表中的11号矢量；-1-1-1执行001执行λT_s、-1-1-1再执行

构成校正矢量表中的12号矢量；111执行

-1-10执行λT_s、111再执行

构成校正矢量表中的13号矢量；-1-1-1执行

101执行λT_s、-1-1-1再执行

构成校正矢量表中的14号矢量；111执行

0-10执行λT_s、111再执行

构成校正矢量表中的15号矢量；

中矢量校正：0-10执行

10-1执行λT_s、0-10再执行构成校正矢量表中的16号矢量；010执行

10-1执行λT_s、010再执行

构成校正矢量表中的17号矢量；100执行

01-1执行λT_s、100再执行

构成校正矢量表中的18号矢量；-100执行

01-1执行λT_s、-100再执行

构成校正矢量表中的19号矢量；00-1执行

-110执行λT_s、00-1再执行

构成校正矢量表中的20号矢量；001执行

-110执行λT_s、001再执行

构成校正矢量表中的21号矢量；010执行

-101执行λT_s、010再执行

构成校正矢量表中的22号矢量；0-10执行

-101执行λT_s、0-10再执行

构成校正矢量表中的23号矢量；-100执行

0-11执行λT_s、-100再执行

构成校正矢量表中的24号矢量；100执行

0-11执行λT_s、100再执行

构成校正矢量表中的25号矢量；001执行

1-10执行λT_s、001再执行构成校正矢量表中的26号矢量；00-1执行1-10执行λT_s、00-1再执行

构成校正矢量表中的27号矢量；

大矢量校正：000执行

1-1-1执行λT_s、000再执行

构成校正矢量表中的28号矢量；000执行

11-1执行λT_s、000再执行

构成校正矢量表中的29号矢量；000执行

-11-1执行λT_s、000再执行

构成校正矢量表中的30号矢量；000执行-111执行λT_s、000再执行

构成校正矢量表中的31号矢量；000执行

-1-11执行λT_s、000再执行

构成校正矢量表中的32号矢量；000执行

1-11执行λT_s、000再执行

构成校正矢量表中的33号矢量。

进一步地，步骤4所述相角补偿：对e_α(k)、e_β(k)做一个采样周期T_s相角补偿得到e_α(k+1)、e_β(k+1)，对i^* _α(k)、i^* _β(k)做两个采样周期2T_s相角补偿得到

具体公式如下：

其中，ω为电网电压基波角频率。

进一步地，步骤5所述计算入网电流反馈值：结合S_opt(k)在校正电压矢量表中的矢量坐标值

和步骤3的采样及变换结果，根据逆变器的校正离散数学模型，计算出第k+1个采样周期入网电流反馈值i_α(k+1)、i_β(k+1)，具体公式如下：

其中，L为滤波电感的感值，R为逆变器桥臂电阻和滤波电感电阻折合后等效电阻的阻值；i_α(k)、i_β(k)为第k个采样周期，三相入网电流的采样值经过Clark变换后的值；e_α(k)、e_β(k)为第k个采样周期，电网电压的采样值经过Clark变换后的值。

进一步地，步骤6所述计算直流侧中点电压反馈值：结合S_opt(k)在校正电压矢量表中的相开关函数状态

和步骤3的采样及变换结果，根据逆变器的校正离散数学模型，计算出第k+1个采样周期直流侧中点电压反馈值Δv_c(k+1)，具体公式如下：

进一步地，步骤8所述遍历计算：结合i_α(k+1)、i_β(k+1)、e_α(k+1)、e_β(k+1)、Δv_c(k+1)、i_a(k+1)、i_b(k+1)、i_c(k+1)、逆变器交流输出侧校正电压矢量和校正电压矢量表中各个校正电压矢量对应的相开关函数状态，根据逆变器的校正离散数学模型，计算出第k+2个采样周期的入网电流i_α(k+2)(i)、i_β(k+2)(i)和直流侧中点电压Δv_c(k+2)(i)，具体公式如下：

其中，i＝1、2……33，u_α(i)、u_β(i)为第i组逆变器交流输出侧校正电压矢量值，

为第i组逆变器交流输出侧校正电压矢量值对应的相开关函数状态。

进一步地，步骤9所述建立目标函数g，作为选出校正电压矢量表中下一个采样周期，即第k+1个采样周期最优矢量S_opt(k+1)的依据，具体公式如下：

其中，i＝1、2……33，λ_dc为直流侧中点电压权重系数。

本发明与现有技术相比，其显著优点在于：(1)仅需对逆变器交流输出侧电压矢量做校正，即可以实现恒定开关频率的有限集模型预测控制，算法简单，易于数字实现；(2)方法可靠性高。

附图说明

图1是NPC三相三电平并网逆变器主功率电路示意图。

图2是本发明多电平逆变器并网系统的结构示意图。

图3是矢量校正前NPC三相三电平并网逆变器交流侧输出电压矢量分布示意图。

图4是矢量校正后NPC三相三电平并网逆变器交流侧输出电压矢量分布示意图。

图5是矢量校正后NPC三相三电平并网逆变器两步预测法流程图。

图6是实施例中直流侧中点电压仿真结果图。

图7是实施例中入网有功功率仿真结果图。

图8是实施例中入网无功功率仿真结果图。

图9是实施例中实际入网电流和给定入网电流仿真结果图。

图10是实施例中逆变器逆变器交流输出侧(a)相对直流侧中点(O)的电压u_ao的频谱分布图。

具体实施方式

以下结合附图及具体实施例对本发明做进一步详细说明。

1、NPC三相三电平并网逆变器基本离散数学模型

图1是NPC三相三电平并网逆变器主功率电路，系统采用如图2所示的三相三线制接法，本发明设定：直流侧电容C₁＝C₂＝C，且C足够大，直流侧电容电压

三相滤波电感的感值L_a＝L_b＝L_c＝L，并且逆变器交流侧的等效阻值R_a＝R_b＝R_c＝R。

定义：相开关函数

其中：i＝a、b、c。

则逆变器交流输出侧(a、b、c)相对直流侧中点(O)的电压：

利用Clark变换：

得到αβ坐标系下逆变器交流输出侧(a、b、c)相对直流侧中点(O)的电压矢量分布，即三电平逆变器的输出电压矢量表，如图3所示。三相三电平逆变器有3³＝27个开关状态，对应输出27个电压矢量，19种不同的电压矢量，图3中“01-1”表示S_a＝0，S_b＝1，S_c＝-1，其余依此类推。

根据基尔霍夫电压定律可得到逆变器交流输出侧电压平衡方程：

其中，v_no为电网电压中性点(n)相对直流侧中点(O)的电压。对(4)式两端进行Clark变换，可以得到αβ坐标系下逆变器交流输出侧电压平衡方程：

取采样周期(即开关周期)为T_s，利用一阶前向差分方程(6)式：

得到αβ坐标系下，αβ坐标系下入网电流的离散数学模型：

式中，i_α(k)、i_β(k)为第k个采样周期，三相入网电流的采样值经过Clark变换后的值，u_α(k)、u_β(k)为第k个采样周期所采用的电压矢量在αβ坐标系下的坐标值，e_α(k)、e_β(k)为第k个采样周期，电网电压的采样值经过Clark变换后的值，i_α(k+1)、i_β(k+1)为αβ坐标系下第k+1个采样周期入网电流的预测计算值。

根据基尔霍夫电流定律可得直流侧中点(O)电流平衡方程：

i_o＝i_c1-i_c2 (8)

其中：

令直流侧中点电压Δv_c＝v_c1-v_c2，同时将(9)、(10)、(11)式带入(8)式，可以得到：

取采样周期为T_s，利用一阶前向差分方程(13)式：

得到自然坐标系下，直流侧中点电压的离散数学模型：

式中，Δv_c(k)为自然坐标系下第k个采样周期直流侧中点电压采样值，Δv_c(k+1)为自然坐标系下第k+1个采样周期直流侧中点电压预测计算值，S_a(k)、S_b(k)、S_c(k)为第k个采样周期所采用的相开关函数状态值。

2、校正矢量与NPC三相三电平并网逆变器校正离散数学模型

引入校正参数λ，依据定频原则对三电平逆变器的输出电压矢量表进行校正，得到校正电压矢量表如图4所示。

第一步：将逆变器输出电压矢量表中的矢量进行分组，零矢量组、小矢量组、中矢量组、大矢量组；其中：

零矢量组中的矢量对应的相开关函数状态S_aS_bS_c为：-1-1-1、000、111，

小矢量组中的矢量对应的相开关函数状态S_aS_bS_c为：100、0-1-1、110、00-1、010、-10-1、011、-100、001、-1-10、101、0-10，

中矢量组中的矢量对应的相开关函数状态S_aS_bS_c为：10-1、01-1、-110、-101、0-11、1-10，

大矢量组中的矢量对应的相开关函数状态S_aS_bS_c为：1-1-1、11-1、-11-1、-111、-1-11、1-11；

第二步：引入校正参数λ，在一个采样周期T_s(即是开关周期)中对三电平逆变器的输出电压矢量表中的矢量进行分组校正，其中，执行时间

的矢量为辅助矢量，属于校正矢量，执行时间λT_s的矢量为主矢量，属于被校正矢量；依据定频原则，具体校正如下：

零矢量校正：000执行

-1-1-1执行λT_s、000再执行构成校正矢量表中的1号矢量；-1-1-1执行000执行λT_s、-1-1-1再执行构成校正矢量表中的2号矢量；000执行111执行λT_s、000再执行

构成校正矢量表中的3号矢量；

小矢量校正：-1-1-1执行

100执行λT_s、-1-1-1再执行

构成校正矢量表中的4号矢量；111执行

0-1-1执行λT_s、111再执行

构成校正矢量表中的5号矢量；-1-1-1执行

110执行λT_s、-1-1-1再执行

构成校正矢量表中的6号矢量；111执行

00-1执行λT_s、111再执行

构成校正矢量表中的7号矢量；-1-1-1执行

010执行λT_s、-1-1-1再执行构成校正矢量表中的8号矢量；111执行

-10-1执行λT_s、111再执行构成校正矢量表中的9号矢量；-1-1-1执行

011执行λT_s、-1-1-1再执行构成校正矢量表中的10号矢量；111执行-100执行λT_s、111再执行构成校正矢量表中的11号矢量；-1-1-1执行

001执行λT_s、-1-1-1再执行构成校正矢量表中的12号矢量；111执行

-1-10执行λT_s、111再执行

构成校正矢量表中的13号矢量；-1-1-1执行

101执行λT_s、-1-1-1再执行

构成校正矢量表中的14号矢量；111执行

0-10执行λT_s、111再执行

构成校正矢量表中的15号矢量；

中矢量校正：0-10执行

10-1执行λT_s、0-10再执行

构成校正矢量表中的16号矢量；010执行

10-1执行λT_s、010再执行

构成校正矢量表中的17号矢量；100执行

01-1执行λT_s、100再执行

构成校正矢量表中的18号矢量；-100执行

01-1执行λT_s、-100再执行

构成校正矢量表中的19号矢量；00-1执行-110执行λT_s、00-1再执行

构成校正矢量表中的20号矢量；001执行

-110执行λT_s、001再执行

构成校正矢量表中的21号矢量；010执行

-101执行λT_s、010再执行

构成校正矢量表中的22号矢量；0-10执行

-101执行λT_s、0-10再执行

构成校正矢量表中的23号矢量；-100执行

0-11执行λT_s、-100再执行

构成校正矢量表中的24号矢量；100执行0-11执行λT_s、100再执行

构成校正矢量表中的25号矢量；001执行

1-10执行λT_s、001再执行构成校正矢量表中的26号矢量；00-1执行

1-10执行λT_s、00-1再执行构成校正矢量表中的27号矢量；

大矢量校正：000执行1-1-1执行λT_s、000再执行

构成校正矢量表中的28号矢量；000执行

11-1执行λT_s、000再执行

构成校正矢量表中的29号矢量；000执行

-11-1执行λT_s、000再执行

构成校正矢量表中的30号矢量；000执行

-111执行λT_s、000再执行

构成校正矢量表中的31号矢量；000执行

-1-11执行λT_s、000再执行

构成校正矢量表中的32号矢量；000执行

1-11执行λT_s、000再执行

构成校正矢量表中的33号矢量。

经过上述校正步骤，可以将校正电压矢量表中的矢量分为两类，一类校正矢量的辅助矢量是零矢量(1号校正矢量到15号校正矢量以及28号校正矢量到33号校正矢量)，另一类校正矢量的辅助矢量是冗余小矢量(16号校正矢量到27号校正矢量)。

由于逆变器交流输出侧的数学模型属于线性数学模型，根据基于线性合成的伏秒平衡原理可知，校正电压矢量对应的逆变器交流输出侧数学模型形式不变，即：

其中，

为校正电压矢量表中的校正矢量的坐标值。

对于直流侧非线性数学模型，校正后的数学模型：

其中，为第i号校正矢量对应的校正相开关函数状态。

3、目标函数与预测计算

模型预测控制属于最优控制范畴，定义与被控变量相关的目标函数g作为最优选择的依据。本发明需要控制三相入网电流i_a、i_b、i_c的Clark变换值i_α、i_β跟踪入网参考电流i^* _α、i^* _β以及抑制直流侧中点电压Δv_c，定义目标函数g如下：

其中，λ_dc为直流侧中点电压权重系数，λ_dc越大，对Δv_c抑制效果越好，电流跟踪效果相对变差；λ_dc越小，对、i_b、抑制效果越差，电流跟踪效果相对变好。

因为实际采样和计算存在延时，本发明采用两步预测法对延时进行补偿，基本原理如图5所示，在k时刻(即当前时刻)，第一步：应用最优矢量S_opt(k)控制并网逆变器，S_opt(k)是上一个采样周期的存储值；第二步：采样i_a、i_b、i_c、e_a、e_b、e_c、Δv_c以及做Clark变换，并对电流做两个采样周期补偿，对电网电压做一个采样周期补偿；第三步：基于电流i_α、i_β与直流侧中点电压Δv_c的校正离散预测模型和S_opt(k)，预测计算被控变量对应的在t_k+1时刻的预测值；第四步：将t_k+1时刻的计算值作为反馈值，结合三相三电平并网逆变器的33种校正矢量与校正离散预测模型，计算出被控变量在t_k+2时刻的33个预测值，并带入目标函数

g的值越小，则被控变量越接近其给定值，因此使目标函数值最小的校正矢量即为最优，将该最优矢量作为t_k+1时刻的控制输出，便能够使t_k+2时刻的被控变量最接近于t_k+2时刻的被控变量给定。

实施例1

实施例的仿真结果如图所示：图6直流侧中点电压、图7入网有功功率、图8入网无功功率、图9实际入网电流和参考入网电流、图10逆变器交流输出侧(a)相对直流侧中点(O)的电压u_ao的频谱分布，仿真参数如表1所示。

从实施例的仿真结果可以看出，本发明提出的多电平逆变器并网系统及方法，能够在保持开关频率恒定的条件下，很好地抑制直流侧中点电压以及控制实际入网电流跟踪给定参考电流，从而实现对入网有功功率和无功功率的控制。

表1仿真参数

电网电压e<sub>a</sub>(＝e<sub>b</sub>＝e<sub>c</sub>)幅值	311(V)
		V<sub>dc</sub>	800(V)
C<sub>1</sub>(＝C<sub>2</sub>)	500e-6(F)
		λ<sub>dc</sub>	1
T<sub>s</sub>	25e-6(s)
		R	1(Ω)
L	4e-3(H)
		P<sup>*</sup><sub>0</sub>	10(kW)
Q<sup>*</sup><sub>0</sub>	500(Var)

综上所述，本发明所多电平逆变器并网系统及方法，能够在几乎不增加算法复杂度的情况下，很好的解决传统模型预测控制开关频率不恒定的问题的同时，实现对直流侧中点电压的强烈抑制和交流侧电流侧精确跟踪。

Claims (7)

1.一种多电平逆变器并网方法，其特征在于，该方法基于矢量校正式定频有限集模型，采用矢量校正后的输出电压矢量表进行两步预测控制，包括以下步骤：

步骤1，校正：引入校正参数λ，依据定频原则对多电平逆变器的输出电压矢量表进行校正，得到校正电压矢量表；

步骤2，最优输出：输出当前，即第k个采样周期最优矢量S_opt(k)控制逆变器，S_opt(k)由上一个，即第k-1个采样周期计算得到；

步骤3，采样及变换：对入网电流进行检测得到a、b、c三相入网电流i_a(k)、i_b(k)、i_c(k)，并对i_a(k)、i_b(k)、i_c(k)进行Clark变换得到i_α(k)、i_β(k)；对电网电压进行检测得到a、b、c三相电网电压e_a(k)、e_b(k)、e_c(k)，并对e_a(k)、e_b(k)、e_c(k)进行Clark变换得到e_α(k)、e_β(k)；根据变换所得e_α(k)、e_β(k)，计算入网参考电流，得到αβ坐标系下的入网参考电流i^* _α(k)、i^* _β(k)；对电容电压进行检测得到直流侧中点电压Δv_c(k)；

步骤4，相角补偿：对e_α(k)、e_β(k)做一个采样周期T_s的相角补偿得到e_α(k+1)、e_β(k+1)，对i^* _α(k)、i^* _β(k)做两个采样周期2T_s的相角补偿得到

步骤5，计算入网电流反馈值：结合S_opt(k)在校正电压矢量表中的矢量坐标值

和步骤3的采样及变换结果，根据逆变器的校正离散数学模型，计算出第k+1个采样周期入网电流反馈值i_α(k+1)、i_β(k+1)；

步骤6，计算直流侧中点电压反馈值：结合S_opt(k)在校正电压矢量表中的相开关函数状态

和步骤3的采样及变换结果，根据逆变器的校正离散数学模型，计算出第k+1个采样周期直流侧中点电压反馈值Δv_c(k+1)；

步骤7，反Clark变换：对第k+1个采样周期入网电流反馈值i_α(k+1)、i_β(k+1)进行反Clark变换，得到i_a(k+1)、i_b(k+1)、i_c(k+1)；

步骤8，遍历计算：结合i_α(k+1)、i_β(k+1)、e_α(k+1)、e_β(k+1)、Δv_c(k+1)、i_a(k+1)、i_b(k+1)、i_c(k+1)、逆变器交流输出侧校正电压矢量和校正电压矢量表中各个校正电压矢量对应的相开关函数状态，根据逆变器的校正离散数学模型，计算出第k+2个采样周期的入网电流i_α(k+2)(i)、i_β(k+2)(i)和直流侧中点电压Δv_c(k+2)(i)；

步骤9，建立目标函数g，作为选出校正电压矢量表中下一个采样周期，即第k+1个采样周期最优矢量S_opt(k+1)的依据；

步骤10，存储目标函数g的最小值对应的矢量，该矢量即最优矢量S_opt(k+1)，将S_opt(k+1)对应的相开关函数状态做为第k+1个采样周期的最优状态输出；

步骤11，等待本次采样周期时间结束，返回步骤2，进入下一个寻优周期。

2.根据权利要求1所述的多电平逆变器并网方法，其特征在于，步骤1所述校正：引入校正参数λ，依据定频原则对多电平逆变器的输出电压矢量表进行校正，得到校正电压矢量表，具体如下：

(1.1)将逆变器输出电压矢量表中的矢量分为零矢量组、小矢量组、中矢量组、大矢量组四组；其中：

零矢量组中的矢量对应的相开关函数状态S_aS_bS_c为：-1-1-1、000、111，

小矢量组中的矢量对应的相开关函数状态S_aS_bS_c为：100、0-1-1、110、00-1、010、-10-1、011、-100、001、-1-10、101、0-10，

中矢量组中的矢量对应的相开关函数状态S_aS_bS_c为：10-1、01-1、-110、-101、0-11、1-10，

大矢量组中的矢量对应的相开关函数状态S_aS_bS_c为：1-1-1、11-1、-11-1、-111、-1-11、1-11；

(1.2)引入校正参数λ，在一个采样周期，即开关周期内对矢量进行分组校正，其中，执行时间

的矢量为辅助矢量，属于校正矢量，执行时间λT_s的矢量为主矢量，属于被校正矢量；依据定频原则，具体校正如下：

零矢量校正：000执行

-1-1-1执行λT_s、000再执行

构成校正矢量表中的1号矢量；-1-1-1执行

000执行λT_s、-1-1-1再执行

构成校正矢量表中的2号矢量；000执行

111执行λT_s、000再执行

构成校正矢量表中的3号矢量；

小矢量校正：-1-1-1执行100执行λT_s、-1-1-1再执行

构成校正矢量表中的4号矢量；111执行

0-1-1执行λT_s、111再执行

构成校正矢量表中的5号矢量；-1-1-1执行

110执行λT_s、-1-1-1再执行构成校正矢量表中的6号矢量；111执行

00-1执行λT_s、111再执行

构成校正矢量表中的7号矢量；-1-1-1执行010执行λT_s、-1-1-1再执行

构成校正矢量表中的8号矢量；111执行

-10-1执行λT_s、111再执行

构成校正矢量表中的9号矢量；-1-1-1执行

011执行λT_s、-1-1-1再执行

构成校正矢量表中的10号矢量；111执行

-100执行λT_s、111再执行

构成校正矢量表中的11号矢量；-1-1-1执行

001执行λT_s、-1-1-1再执行

构成校正矢量表中的12号矢量；111执行

-1-10执行λT_s、111再执行

构成校正矢量表中的13号矢量；-1-1-1执行101执行λT_s、-1-1-1再执行

构成校正矢量表中的14号矢量；111执行

0-10执行λT_s、111再执行

构成校正矢量表中的15号矢量；

中矢量校正：0-10执行

10-1执行λT_s、0-10再执行

构成校正矢量表中的16号矢量；010执行

10-1执行λT_s、010再执行

构成校正矢量表中的17号矢量；100执行

01-1执行λT_s、100再执行

构成校正矢量表中的18号矢量；-100执行

01-1执行λT_s、-100再执行

构成校正矢量表中的19号矢量；00-1执行-110执行λT_s、00-1再执行

构成校正矢量表中的20号矢量；001执行-110执行λT_s、001再执行

构成校正矢量表中的21号矢量；010执行-101执行λT_s、010再执行

构成校正矢量表中的22号矢量；0-10执行

-101执行λT_s、0-10再执行

构成校正矢量表中的23号矢量；-100执行

0-11执行λT_s、-100再执行

构成校正矢量表中的24号矢量；100执行

0-11执行λT_s、100再执行构成校正矢量表中的25号矢量；001执行

1-10执行λT_s、001再执行

构成校正矢量表中的26号矢量；00-1执行1-10执行λT_s、00-1再执行

构成校正矢量表中的27号矢量；

大矢量校正：000执行

1-1-1执行λT_s、000再执行

构成校正矢量表中的28号矢量；000执行

11-1执行λT_s、000再执行

构成校正矢量表中的29号矢量；000执行

-11-1执行λT_s、000再执行

构成校正矢量表中的30号矢量；000执行

-111执行λT_s、000再执行

构成校正矢量表中的31号矢量；000执行

-1-11执行λT_s、000再执行

构成校正矢量表中的32号矢量；000执行1-11执行λT_s、000再执行构成校正矢量表中的33号矢量。

3.根据权利要求1所述的多电平逆变器并网方法，其特征在于，步骤4所述相角补偿：对e_α(k)、e_β(k)做一个采样周期T_s相角补偿得到e_α(k+1)、e_β(k+1)，对i^* _α(k)、i^* _β(k)做两个采样周期2T_s相角补偿得到具体公式如下：

其中，ω为电网电压基波角频率。

4.根据权利要求1所述的多电平逆变器并网方法，其特征在于，步骤5所述计算入网电流反馈值：结合S_opt(k)在校正电压矢量表中的矢量坐标值

和步骤3的采样及变换结果，根据逆变器的校正离散数学模型，计算出第k+1个采样周期入网电流反馈值i_α(k+1)、i_β(k+1)，具体公式如下：

其中，L为滤波电感的感值，R为逆变器桥臂电阻和滤波电感电阻折合后等效电阻的阻值；i_α(k)、i_β(k)为第k个采样周期，三相入网电流的采样值经过Clark变换后的值；e_α(k)、e_β(k)为第k个采样周期，电网电压的采样值经过Clark变换后的值。

5.根据权利要求1所述的多电平逆变器并网方法，其特征在于，步骤6所述计算直流侧中点电压反馈值：结合S_opt(k)在校正电压矢量表中的相开关函数状态

和步骤3的采样及变换结果，根据逆变器的校正离散数学模型，计算出第k+1个采样周期直流侧中点电压反馈值Δv_c(k+1)，具体公式如下：

6.根据权利要求1所述的多电平逆变器并网方法，其特征在于，步骤8所述遍历计算：结合i_α(k+1)、i_β(k+1)、e_α(k+1)、e_β(k+1)、Δv_c(k+1)、i_a(k+1)、i_b(k+1)、i_c(k+1)、逆变器交流输出侧校正电压矢量和校正电压矢量表中各个校正电压矢量对应的相开关函数状态，根据逆变器的校正离散数学模型，计算出第k+2个采样周期的入网电流i_α(k+2)(i)、i_β(k+2)(i)和直流侧中点电压Δv_c(k+2)(i)，具体公式如下：

其中，i＝1、2……33，u_α(i)、u_β(i)为第i组逆变器交流输出侧校正电压矢量值，

为第i组逆变器交流输出侧校正电压矢量值对应的相开关函数状态。

7.根据权利要求1所述的多电平逆变器并网方法，其特征在于，步骤9所述建立目标函数g，作为选出校正电压矢量表中下一个采样周期，即第k+1个采样周期最优矢量S_opt(k+1)的依据，具体公式如下：

其中，i＝1、2……33，λ_dc为直流侧中点电压权重系数。

Images