CN110297426B - 采用周期干扰差分抑制策略的并联型半周期重复控制方法 - Google Patents

Abstract

本发明公开了一种采用周期干扰差分抑制策略的并联型半周期重复控制方法。根据齐次/偶次谐波干扰信号的半周期负/正对称特性，分别构造出相应的等效干扰信号；将等效干扰嵌入到吸引律，构造e/v信号转换模块，其输出信号用于并联型半重复控制器的修正量；继而计算出其输出信号作为逆变器对象的输入，使逆变器系统跟随参考信号变化。本发明提供的时域设计的离散并联型半周期重复控制器是一种快速收敛、高控制精度、消除齐次谐波和偶次谐波干扰信号、有效抑制分数次谐波和其它非周期干扰信号的控制器。

Description

采用周期干扰差分抑制策略的并联型半周期重复控制方法

技术领域

本发明涉及一种采用周期干扰差分抑制策略的并联型半周期重复控制方法，该方法适于逆变电源，也适用于工业控制中的周期运行过程。

背景技术

多年来，周期信号的跟踪和干扰抑制补偿问题一直是众多学者关注的课题。重复控制是一种适用于周期系统的控制技术，现有的重复控制技术主要集中于基于内模原理的频域设计方法。内模原理的本质是采用延迟时间为T的延迟环节的正反馈形式1/(1-e^-Ts)来构造周期为T的周期信号内模，并将其嵌入稳定的闭环系统中，内模输出就会对输入信号逐周期累加形成控制作用，从而解决周期性参考信号的跟踪或周期性干扰信号的抑制问题。这种控制技术已经广泛应用于电机伺服系统、电力电子逆变器、硬盘/光盘伺服系统及其它重复运行过程。

实际工程中采用计算机控制技术，控制系统多是以离散时间方式实现。离散重复控制器设计主要有两种途径：一种是通过对连续重复控制器离散化得到；另一种是直接针对离散时间系统进行控制器设计。取采样周期T_s，使得参考信号周期为采样周期的整数倍，记每个周期中的采样点个数为N，即T＝NT_s。这样，离散周期信号内模为1/(1-Z^-N)。减少控制器的内存占用量和内模的响应速度是重复控制实现时需要考虑的问题。专利号US7265932B2的美国发明专利针对满足的半周期对称信号x(t)＝-x(t-N/2)，提出了一种半周期重复控制器，比整周期重复控制器的内存占用量减少一半，且系统在半周期后就开始内模响应。Costa-Castello等(Costa-Castello R,Grino R,FossasE.Odd-Harmonicdigital repetitive control of a single-phase current active filter.IEEETransactions on Power Electronics,19(4):1060-1068,2004)提出一种齐次谐波重复控制方法也减少一半内存占用量。但对于偶次谐波是关于半周期正对称，仅用齐次谐波重复控制器，甚至会放大偶次谐波干扰。Keliang Zhou等(Keliang Zhou,Danwei Wang,etal.Plug-in dual-mode-structure repetitive controller for CVCF PWM inverters,IEEE Transactions on Industrial Electronics,56(3):784-791,2009)设计了一种双模结构重复控制器，可同时对齐次谐波和偶次谐波进行抑制，但并未考虑分数次谐波干扰的抑制问题。上述重复控制器的设计是在频域进行的，然而信号对称性表现在时域中，对于更为复杂的对称性信号并不能进行有效处理。因此，仍然有必要继续深入地研究重复控制技术。

发明内容

为了克服已有重复控制方法的未考虑干扰信号的对称性质和抑制干扰信号缓慢、动态品质较差以及控制时滞过长等问题，本发明提供了一种采用周期干扰差分抑制策略的并联型半周期重复控制方法，可同时抑制齐次谐波和偶次谐波，且设计过程是在时域进行，更有利于结合现有的干扰观测器技术。采用这种并联型半周期重复控制技术能够实现周期参考信号跟踪任务，同时可快速消除齐/偶次谐波干扰。

本发明解决上述技术问题采用的技术方案为：

采用周期干扰差分抑制策略的并联型半周期重复控制方法，包括以下步骤：

步骤1：给定参考信号r_k是周期为N的正弦信号，满足

r_k＝Asin(2πk/N)，r_k＝r_k-N (1)

其中，r_k,r_k-N分别为k，k-N时刻的给定参考信号，A为给定参考信号r_k的幅值；

步骤2：齐次谐波m_k满足半周期的负对称特性：

m_k＝-m_k-N/2 (2)

其中，m_k,m_k-N/2分别为第k和k-N/2时刻的齐次谐波；偶次谐波n_k满足半周期的正对称特性：

n_k＝n_k-N/2 (3)

其中，n_k,n_k-N/2分别为第k和k-N/2时刻的偶次谐波；根据齐次/偶次谐波的半周期对称特性，分别构造两种等效干扰。对于式(2)，齐次谐波的等效干扰d_k,m为

d_k,m＝w_k+w_k-N/2 (4)

其中，w_k,w_k-N/2分别为k和k-N/2时刻的干扰信号；d_k,m为第k时刻的齐次谐波等效干扰；

对于式(3)，偶次谐波的等效干扰d_k,n为

d_k,n＝w_k-w_k-N/2 (5)

其中，d_k,n为第k时刻的偶次谐波等效干扰；

步骤3：构造下述离散吸引律：

e_k+1＝(1-ρ)e_k-εΦ(e_k,δ) (6)

其中，Φ(e_k,δ)为改进型非线性幂次函数，具有“大误差幂次、小误差放大”的特性，其具体形式如下：

δ为改进型非线性幂次函数的边界层参数，e_k+1,e_k分为k+1，k时刻的跟踪误差，且e_k＝r_k-y_k；0＜ρ＜1,δ＞0,0＜λ＜1为用于调整吸引速度的参数；

为了提高系统对齐次和偶次谐波干扰的抑制能力，将周期干扰抑制补偿项嵌入到吸引律(6)中，形成如下理想误差动态

e_k+1＝(1-ρ)e_k-εΦ(e_k,δ)+k₁(d_k,m-d_k+1,m)+k₂(d_k,n-d_k+1,n) (8)

其中，0≤k₁,k₂≤1，k₁+k₂＝1；d_k+1,m和d_k+1,n可分别抑制齐次谐波和偶次谐波干扰信号；d_k,m和d_k,n分别用于补偿d_k+1,m和d_k+1,n，实现对分数次谐波和非周期干扰信号的抑制。

步骤4：依据吸引律(8)，并联型半周期重复控制器的表达式为

u_k＝k₁u_k,m+k₂u_k,n (9)

其中，u_k,m，u_k,n分别为第k时刻的齐次谐波和偶次谐波半周期重复控制器，具体形式如下：

其中，u_k,m，u_k-1,m，u_k-N/2,m，u_k-1-N/2,m分别为第k，k-1，k-N/2，k-1-N/2时刻的齐次谐波半周期重复控制器；u_k,n，u_k-1,n，u_k-N/2,n，u_k-1-N/2,n分别为第k，k-1，k-N/2，k-1-N/2时刻的偶次谐波半周期重复控制器；r_k+1分别为k+1时刻的给定参考信号；y_k,y_k-1,y_k-1-N/2,y_k-N/2,y_k+1-N/2分别为第k，k-1，k-1-N/2，k-N/2，k+1-N/2时刻的输出信号；式(10)和式(11)中的a₁,a₂,b₁,b₂为逆变器的系统参数，逆变器的数学模型如下：

y_k+1+a₁y_k+a₂y_k-1＝b₁u_k+b₂u_k-1+w_k+1 (12)

其中，y_k+1,y_k,y_k-1为逆变器k+1，k，k-1时刻的输出信号，u_k,u_k-1表示逆变器k，k-1时刻的控制输入信号，a₁,a₂,b₁,b₂为系统参数；w_k+1＝m_k+1+n_k+1+Δw_k+1为k+1时刻的系统干扰信号，包括齐次谐波m_k+1、偶次谐波n_k+1以及其他间谐波和参数摄动Δw_k+1；

步骤5：将u_k作为逆变器的控制输入信号，可量测获得逆变器系统输出信号y_k，跟随参考信号r_k变化，且系统跟踪误差的动态行为由式(8)表征；

步骤6：进一步，为表征吸引律的吸引过程，本发明给出单调减区域、绝对吸引层以及稳态误差带这3个指标的表达式；这3个指标可用于指导控制器参数整定，其中单调减区域、绝对吸引层、稳态误差带的定义如下：

1)单调减区域Δ_MDR

2)绝对吸引层Δ_AAL

|e_k+1|＜|e_k|，当|e_k|＞Δ_AAL (14)

3)稳态误差带Δ_SSE

|e_k+1|≤Δ_SSE，当|e_k|≤Δ_SSE (15)

这里，Δ_MDR为单调减区域边界，Δ_AAL为绝对吸引层边界，Δ_SSE为稳态误差带边界。

在并联型半周期重复控制器(13)作用下，且干扰补偿误差满足

|k₁(d_k,m-d_k+1,m)+k₂(d_k,n-d_k+1,n)|≤Δ (16)

时，各指标的表达式如下：

1)单调减区域Δ_MDR表示为：

Δ_MDR＝max{Δ_MDR1,Δ_MDR2} (17)

式中，Δ_MDR1，Δ_MDR2为正实数，由下式确定

其中，Δ为干扰补偿误差的上确界；

2)绝对吸引层Δ_AAL表示为：

Δ_AAL＝max{Δ_AAL1,Δ_AAL2} (19)

式中，Δ_AAL1，Δ_AAL2为正实数，且满足

3)稳态误差带Δ_SSE表示为：

Δ_SSE＝max{Δ_SSE1,Δ_SSE2,Δ}

(21)

式中，Δ_SSE1，Δ_SSE2为正实数，由下式确定

其中，当ξ＞δ时，

当ξ≤δ时，

本发明的技术构思为：提出一种采用周期干扰差分抑制策略的并联型半周期重复控制方法。根据齐次谐波的负半周期对称和偶次谐波的正半周期对称，分别设计出齐次谐波和偶次谐波半周期重复控制器，将其以并联方式组合成并联型半周期重复控制器。依据该控制器可同时抑制齐次谐波和偶次谐波，实现对参考信号的精确跟踪。

本发明的控制效果主要表现在：采用离散幂次吸引律，将齐次谐波和偶次谐波干扰补偿项“嵌入”到吸引律中，依此设计出并联型半周期重复控制器，该控制技术不仅能跟踪上给定的参考信号，而且可以实现对齐次谐波和偶次谐波的消除；同时，采用新型幂次吸引律，减少系统颤振的同时又保证系统的鲁棒性。

附图说明

图1是逆变器的齐次和偶次谐波干扰的对称特性图；

图2是基于吸引律的控制系统设计流程图；

图3是本发明的非线性幂次函数Φ(e_k,δ)、幂次函数fal(e_k,δ)和幂次型符号函数|e_k|^λsgn(e_k)的比较图；

图4是本发明提出的吸引律e_k+1＝(1-ρ)e_k-εΦ(e_k,δ)和基于幂次函数的吸引律e_k+1＝(1-ρ)e_k-εfal(e_k,δ)和幂次吸引律e_k+1＝(1-ρ)e_k-ε|e_k|^λsgn(e_k)的误差收敛速度比较图；

图5是并联型半周期重复控制系统方框图；

图6是当控制器参数ρ＝0.3,ε＝0.3,δ＝0.5,λ＝0.5时，在齐次谐波半周期重复控制器(15)(不考虑干扰补偿项)作用下的给定参考信号、输出信号和跟踪误差信号；

图7是当控制器参数ρ＝0.3,ε＝0.3,δ＝0.5,λ＝0.5时，在并联型半周期重复控制器(21)(不考虑干扰补偿项)作用下的给定参考信号、输出信号和跟踪误差信号；

图8是当控制器参数ρ＝0.3,ε＝0.3,δ＝0.5,λ＝0.5时，在并联型半周期重复控制器(21)(不考虑干扰补偿项)作用下的给定参考信号、输出信号、跟踪误差信号以及三个边界；

图9是当控制器参数ρ＝0.3,ε＝0.3,δ＝0.5,λ＝0.5时，在并联型半周期重复控制器(21)(考虑干扰补偿项)作用下的给定参考信号、输出信号、跟踪误差信号以及三个边界。

具体实施方式

下面结合附图对本发明具体实施方式做进一步描述。

参照图1-9，一种采用周期干扰差分抑制策略的并联型半周期重复控制方法，被控对象为逆变器，所述的逆变器控制系统数学模型为：

y_k+1+a₁y_k+a₂y_k-1＝b₁u_k+b₂u_k-1+w_k+1 (1)

其中，y_k+1,y_k,y_k-1为逆变器k+1，k，k-1时刻的输出信号，u_k,u_k-1表示逆变器k，k-1时刻的控制输入信号，a₁,a₂,b₁,b₂为系统参数；w_k+1＝m_k+1+n_k+1+Δw_k+1为k+1时刻的系统干扰信号，包括齐次谐波m_k+1、偶次谐波n_k+1以及其他间谐波和参数摄动Δw_k+1；

给定参考信号r_k是周期为N的正弦信号，满足

r_k＝Asin(2πk/N)，r_k＝r_k-N (2)

其中，r_k,r_k-N分别为k，k-N时刻的给定参考信号，A为给定参考信号r_k的幅值；

齐次谐波m_k满足半周期的负对称特性(见图1)：

m_k＝-m_k-N/2 (3)

其中，m_k,m_k-N/2分别为第k和k-N/2时刻的齐次谐波；偶次谐波n_k满足半周期的正对称特性(见图1)：

n_k＝n_k-N/2 (4)

其中，n_k,n_k-N/2分别为第k和k-N/2时刻的偶次谐波；根据齐次/偶次谐波的半周期对称，分别构造两种等效干扰。对于式(3)，齐次谐波的等效干扰d_k,m为

d_k,m＝w_k+w_k-N/2 (5)

其中，w_k,w_k-N/2分别为k和k-N/2时刻的干扰信号；d_k,m为第k时刻的齐次谐波等效干扰；对于式(4)，偶次谐波的等效干扰d_k,n为

d_k,n＝w_k-w_k-N/2 (6)

其中，d_k,n为第k时刻的偶次谐波等效干扰；

一种采用周期干扰差分抑制策略的并联型半周期重复控制方法的设计流程见图2。构造下述离散吸引律：

e_k+1＝(1-ρ)e_k-εΦ(e_k,δ) (7)

其中，Φ(e_k,δ)为非线性幂次函数，具有“大误差幂次、小误差放大”的特性，其具体形式如下：

δ为非线性幂次函数的边界层参数；非线性幂次函数、幂次函数|e_k|^λsgn(e_k)以及fal()函数的比较图(见图3)；e_k+1,e_k分为k+1，k时刻的跟踪误差，且e_k＝r_k-y_k；0＜ρ＜1,δ＞0,0＜λ＜1为用于调整吸引速度的参数；吸引律(7)、幂次吸引律e_k+1＝(1-ρ)e_k-ε|e_k|^λsgn(e_k)以及基于fal函数的吸引律e_k+1＝(1-ρ)e_k-εfal(e_k,δ)的误差收敛速度比较(见图4)。

为了提高系统对齐次和偶次谐波干扰的抑制能力，可将离散吸引律(7)修正为

e_k+1＝(1-ρ)e_k-εΦ(e_k,δ)+d_k,m-d_k+1,m (9)

e_k+1＝(1-ρ)e_k-εΦ(e_k,δ)+d_k,n-d_k+1,n (10)

其中，d_k+1,m和d_k+1,n可分别抑制齐次谐波和偶次谐波干扰信号；d_k,m和d_k,n分别用于补偿d_k+1,m和d_k+1,n，实现对分数次谐波和非周期干扰信号的抑制。

由跟踪误差e_k＝r_k-y_k和系统(1)知，

将w_k+1+w_k+1-N/2表达为

记等效干扰d_k,m＝w_k+w_k-N/2，

将式(13)代入式(9)，则齐次谐波半周期重复控制器的表达式为

为了区别齐次谐波半周期重复控制器和偶次谐波半周期重复控制器，依据齐次谐波的半周期对称性而获得控制器记为u_k,m，则相应的齐次谐波半周期重复控制器表达为

由跟踪误差e_k＝r_k-y_k和系统(1)知，

将w_k+1-w_k+1-N/2表达为

记等效干扰d_k,n＝w_k-w_k-N/2，

将式(18)代入式(10)，则偶次谐波半周期重复控制器的表达式为

为了区别齐次谐波半周期重复控制器和偶次谐波半周期重复控制器，依据偶次谐波的半周期对称性而获得控制器记为u_k,n，则相应的偶次谐波半周期重复控制器表达为

为了同时抑制齐次谐波和偶次谐波干扰，将式(15)和(20)以并联形式组合成并联型半周期重复控制器(见图5)，其具体形式如下：

u_k＝k₁u_k,m+k₂u_k,n (21)

其中，0≤k₁,k₂≤1，k₁+k₂＝1。将u_k作为逆变器的控制输入信号，可量测获得逆变器系统输出信号y_k，跟随参考信号r_k变化，且系统跟踪误差的动态行为由式(22)表征；

e_k+1＝(1-ρ)e_k-εΦ(e_k,δ)+k₁(d_k,m-d_k+1,m)+k₂(d_k,n-d_k+1,n) (22)

进一步，为表征吸引律的吸引过程，本发明给出单调减区域、绝对吸引层以及稳态误差带这3个指标的表达式；这3个指标可用于指导控制器参数整定，其中单调减区域、绝对吸引层、稳态误差带的定义如下：

1)单调减区域Δ_MDR

2)绝对吸引层Δ_AAL

|e_k+1|＜|e_k|，当|e_k|＞Δ_AAL (24)

3)稳态误差带Δ_SSE

|e_k+1|≤Δ_SSE，当|e_k|≤Δ_SSE (25)

这里，Δ_MDR为单调减区域边界，Δ_AAL为绝对吸引层边界，Δ_SSE为稳态误差带边界。

在并联型半周期重复控制器(21)作用下，且干扰补偿误差满足

|k₁(d_k,m-d_k+1,m)+k₂(d_k,n-d_k+1,n))|≤Δ (26)

时，各指标的表达式如下：

1)单调减区域Δ_MDR表示为：

Δ_MDR＝max{Δ_MDR1,Δ_MDR2} (27)

式中，Δ_MDR1，Δ_MDR2为正实数，由下式确定

其中，Δ为干扰补偿误差的上确界；

2)绝对吸引层Δ_AAL表示为：

Δ_AAL＝max{Δ_AAL1,Δ_AAL2} (29)

式中，Δ_AAL1，Δ_AAL2为正实数，且满足

3)稳态误差带Δ_SSE表示为：

Δ_SSE＝max{Δ_SSE1,Δ_SSE2,Δ}

(31)

式中，Δ_SSE1，Δ_SSE2为正实数，由下式确定

其中，当ξ＞δ时，

当ξ≤δ时，

再进一步，离散多周期重复控制器设计完成后，需要整定其中的控制器参数。其可调参数包括ρ，ε，δ，λ，k₁,k₂，其参数k₁,k₂的整定可根据齐次谐波和偶次谐波的幅值来进行，其它参数ρ，ε，δ，λ的整定可根据表征吸引律吸引过程的三个边界指标进行。

实施例

针对逆变器输出波形进行闭环控制。所采用的逆变器系统由给定正弦信号部分，并联型半周期重复控制器，PWM调制部分，逆变器主控电路以及采样电路构成。其中给定正弦信号、并联型半周期重复控制器以及PWM调制模块均由DSP控制板实现，其余部分均由逆变器硬件电路实现。整个逆变器控制系统由DSP给定需要输出的期望信号，经过PWM调制后驱动逆变器的功率开关管的高低脉冲信号，实现导通和关断。逆变器输出信号经过LC滤波电流还原成正弦信号，由采样电流采样所需的电压和电流等信号数据返回到DSP，然后经过并联型半周期重复控制器作用后修正输入信号，实现逆变器波形跟踪控制，降低逆变器输出波形的THD值。

下面给出逆变器的并联型半周期重复控制器的设计过程。

首先建立系统数学模型。将逆变器的主控电路、LC滤波电路以及采样电路作为对象进行建模，得到如下二阶差分方程模型

y_k+1+a₁y_k+a₂y_k-1＝b₁u_k+b₂u_k-1+w_k+1 (33)

其中，y_k+1,y_k,y_k-1分别表示k+1，k，k-1时刻的逆变器输出电压，u_k,u_k-1表示k和k-1时刻逆变器的控制量，w_k+1为逆变器系统的不确定特性，由外部干扰以及未建模特性等所组成。模型中的系统参数a₁,a₂,b₁,b₂通过机理建模所获得，其具体取值如下：

a₁＝-0.5775,a₂＝0.2804,b₁＝0.4102,b₂＝0.2589 (34)

实施例中，逆变器的给定参考信号r_k+1＝220sin(2πfkT_s)，单位为V，信号频率为50Hz，采样周期为T_s＝0.0001s，参考信号周期为0.02s。逆变器系统的扰动信号为

其中，第一项用于模拟逆变器齐次谐波干扰信号，第二项用于模拟逆变器偶次谐波干扰信号，第三项用于模拟逆变器间(分数次)谐波干扰信号，第四项为随机扰动信号。

取上述系统参数进行数值仿真，用于检验并联型半周期重复控制器在逆变器系统上的实施结果，并与半周期重复控制器的实施结果进行对比：

(1)采用齐次谐波半周期控制器(15)且不考虑干扰补偿项(令d_k,m＝0)；仿真时，逆变器系统干扰只考虑齐次谐波和偶次谐波干扰，则h₁＝5,h₂＝1,h₃＝0,h₄＝0；控制器参数选取为ρ＝0.3,ε＝0.3,δ＝0.5,λ＝0.5，仿真结果见图6。

(2)采用并联型半周期重复控制器(21)且不考虑干扰补偿项(令d_k,m＝d_k,n＝0)，其中k₁＝0.7,k₂＝0.3；仿真时，逆变器系统干扰只考虑齐次谐波和偶次谐波干扰，则h₁＝5,h₂＝1,h₃＝0,h₄＝0；控制器参数选取为ρ＝0.3,ε＝0.3,δ＝0.5,λ＝0.5，仿真结果见图7。

上述数值仿真结果验证了本发明专利给出的并联型半周期重复控制器在半个周期(0.01s)就起作用，实现了齐次谐波和偶次谐波干扰的抑制(见图7)。从图6可以看出，当采用常规的半周期重复控制器只能抑制齐次谐波干扰，甚至放大了偶次谐波干扰。

为了检验本发明给出的单调减区域Δ_MDR，绝对吸引层边界Δ_AAL，稳态误差带边界Δ_SSE的边界值以及干扰补偿措施的有效性，就ρ，ε，δ，λ的取值和是否嵌入干扰补偿项分别进行讨论。仿真时，考虑齐次谐波、偶次谐波、间谐波以及其他非周期干扰，置h₁＝5,h₂＝1,h₃＝1,h₄＝0.1。

(1)采用并联型半周期重复控制器(21)且不考虑干扰补偿项(令d_k,m＝d_k,n＝0)，控制器参数设置为ρ＝0.3,ε＝0.3,δ＝0.5,λ＝0.5。此时，|k₁(d_k,m-d_k+1,m)+k₂(d_k,n-d_k+1,n)|≤Δ＝0.746，相应的Δ_MDR＝1.6094，Δ_AAL＝Δ_SSE＝1.3324，仿真结果见图8。

(2)采用并联型半周期重复控制器(21)且考虑干扰补偿项，控制器参数设置为ρ＝0.3,ε＝0.3,δ＝0.5,λ＝0.5。此时，|k₁(d_k,m-d_k+1,m)+k₂(d_k,n-d_k+1,n)|≤Δ＝0.44，相应的Δ_MDR＝1.0724，Δ_SSE＝Δ_AAL＝0.6564，仿真结果见图9。

数值仿真结果表明，采用并联型半周期重复控制器能到达预期控制效果，实现了快速、有效地抑制逆变器系统的齐次谐波、偶次谐波、间谐波以及其他非周期干扰信号。从图8和图9可以看出，嵌入了干扰补偿项的控制器在抑制间谐波以及其他非周期干扰信号上有更大的优势，稳态误差带的界更小。

Claims (2)

1.一种采用周期干扰差分抑制策略的并联型半周期重复控制方法，其特征在于：所述方法包括以下步骤：

步骤1：给定参考信号r_k是周期为N的正弦信号，满足

r_k＝A sin(2πk/N)，r_k＝r_k-N (1)

其中，r_k,r_k-N分别为k，k-N时刻的给定参考信号，A为给定参考信号r_k的幅值；

步骤2：齐次谐波m_k满足半周期的负对称特性：

m_k＝-m_k-N/2 (2)

其中，m_k,m_k-N/2分别为第k和k-N/2时刻的齐次谐波；偶次谐波n_k满足半周期的正对称特性：

n_k＝n_k-N/2 (3)

其中，n_k,n_k-N/2分别为第k和k-N/2时刻的偶次谐波；根据齐次/偶次谐波的半周期对称，分别构造两种等效干扰；对于式(2)，齐次谐波的等效干扰d_k,m为

d_k,m＝w_k+w_k-N/2 (4)

其中，w_k,w_k-N/2分别为k和k-N/2时刻的干扰信号；d_k,m为第k时刻的齐次谐波等效干扰；对于式(3)，偶次谐波的等效干扰d_k,n为

d_k,n＝w_k-w_k-N/2 (5)

其中，d_k,n为第k时刻的偶次谐波等效干扰；

步骤3：构造下述离散吸引律：

e_k+1＝(1-ρ)e_k-εΦ(e_k,δ) (6)

其中，Φ(e_k,δ)为改进型非线性幂次函数，具有“大误差幂次、小误差放大”的特性，其具体形式如下：

δ为非线性幂次函数的边界层参数，e_k+1,e_k分为k+1，k时刻的跟踪误差，且e_k＝r_k-y_k；0＜ρ＜1,δ＞0,0＜λ＜1为用于调整吸引速度的参数；

为了提高系统对齐次和偶次谐波干扰的抑制能力，将周期干扰抑制补偿项嵌入到吸引律(6)中，形成如下理想误差动态

e_k+1＝(1-ρ)e_k-εΦ(e_k,δ)+k₁(d_k,m-d_k+1,m)+k₂(d_k,n-d_k+1,n) (8)

其中，0≤k₁,k₂≤1，k₁+k₂＝1；d_k+1,m和d_k+1,n可分别抑制齐次谐波和偶次谐波干扰信号；d_k,m和d_k,n分别用于补偿d_k+1,m和d_k+1,n，实现对分数次谐波和非周期干扰信号的抑制；

步骤4：依据吸引律(8)，并联型半周期重复控制器的表达式为

u_k＝k₁u_k,m+k₂u_k,n (9)

其中，u_k,m，u_k,n分别为第k时刻的齐次谐波和偶次谐波半周期重复控制器，具体形式如下：

其中，u_k,m，u_k-1,m，u_k-N/2,m，u_k-1-N/2,m分别为第k，k-1，k-N/2，k-1-N/2时刻的齐次谐波半周期重复控制器；u_k,n，u_k-1,n，u_k-N/2,n，u_k-1-N/2,n分别为第k，k-1，k-N/2，k-1-N/2时刻的偶次谐波半周期重复控制器；r_k+1分别为k+1时刻的给定参考信号；y_k,y_k-1,y_k-1-N/2,y_k-N/2,y_k+1-N/2分别为第k，k-1，k-1-N/2，k-N/2，k+1-N/2时刻的输出信号；式(10)和式(11)中的a₁,a₂,b₁,b₂为逆变器的系统参数，逆变器的数学模型如下：

y_k+1+a₁y_k+a₂y_k-1＝b₁u_k+b₂u_k-1+w_k+1 (12)

其中，y_k+1,y_k,y_k-1为逆变器k+1，k，k-1时刻的输出信号，u_k,u_k-1表示逆变器k，k-1时刻的控制输入信号，a₁,a₂,b₁,b₂为系统参数；w_k+1＝m_k+1+n_k+1+Δw_k+1为k+1时刻的系统干扰信号，包括齐次谐波m_k+1、偶次谐波n_k+1以及其他间谐波和参数摄动Δw_k+1；

步骤5：将u_k作为逆变器的控制输入信号，可量测获得逆变器系统输出信号y_k，跟随参考信号r_k变化，且系统跟踪误差的动态行为由式(8)表征。

2.如权利要求1所述的一种采用周期干扰差分抑制策略的并联型半周期重复控制方法，其特征在于：采用并联型半周期重复控制器(9)，系统跟踪误差的吸引过程由单调减区域、绝对吸引层以及稳态误差带这三个边界指标表征，这3个指标可用于指导控制器参数整定；

在并联型半周期重复控制器(9)作用下，且干扰补偿误差满足

|k₁(d_k,m-d_k+1,m)+k₂(d_k,n-d_k+1,n)|≤Δ (13)

时，各指标的表达式如下：

1)单调减区域Δ_MDR表示为：

Δ_MDR＝max{Δ_MDR1,Δ_MDR2} (14)

式中，Δ_MDR1，Δ_MDR2为正实数，由下式确定

其中，Δ为干扰补偿误差的上确界；

2)绝对吸引层Δ_AAL表示为：

Δ_AAL＝max{Δ_AAL1,Δ_AAL2} (16)

式中，Δ_AAL1，Δ_AAL2为正实数，且满足

3)稳态误差带Δ_SSE表示为：

Δ_SSE＝max{Δ_SSE1,Δ_SSE2,Δ} (18)

式中，Δ_SSE1，Δ_SSE2为正实数，由下式确定

其中，当ξ＞δ时，

当ξ≤δ时，

再进一步，参数k₁,k₂的整定可根据齐次谐波和偶次谐波的幅值来进行。

Images