用于位置伺服系统的半周期重复控制器
技术领域
本发明涉及一种重复伺服系统在半周期对称参考信号下的控制技术,也适用于工业场合中的周期运行过程。
背景技术
重复控制器具有“记忆”和“学习”特性,它以跟踪误差信号修正前一周期的控制输入,形成当前的控制输入。它能够完全抑制周期性干扰,从而实现精确控制。重复控制技术已成功应用于旋转电机、硬盘驱动、VCD/DVD、UPS、电力电子线路以及电能质量控制等。
通常基于内模原理设计重复控制器,如果某信号被看作是一个自治系统的输出,将该信号的模型放置在稳定的闭环系统中,被控量的输出将能够完全跟踪参考信号。重复控制器构造周期信号内模
其中T为给定信号的周期。它是一个含周期时延(e
-Ts)的正反馈环节。不考虑输入信号的具体形式,只要给定初始段信号,内模输出就会对输入信号逐周期累加,重复输出与上周期相同的信号。采用连续内模的重复控制器设计多是频域设计。实际工程中采用计算机控制技术,控制系统多是以离散时间方式实现。离散重复控制器设计主要有两种途径:一种是通过对连续重复控制器离散化得到;另一种是直接针对离散时间系统进行设计。取采样间隔T
s,使得参考信号周期为采样间隔的整数倍,记每个周期中的采样点个数为N,即T=NT
s。这样,离散周期信号内模为
离散内模的计算复杂程度取决于采样周期T
s,实现离散周期内模时所需内存量至少为N。如果T
s取得过大,系统控制精度降低;取得过小,内模的阶次将会增加。
减少控制器的内存占用量是重复控制实现时要考虑的问题。在Hoog的专利(Hoog T J D.Repetitive controller having reduced memory elements.United States Patent,US 7265932B2,2007)中,对于满足x(t+T/2)=-x(t)的半周期对称信号,提出了一种半周期重复控制内模,构造重复控制器,使得内存占用量,比整周期重复控制的内存使用量减少一半。Costa-Castello等提出奇次谐波重复控制方法,它能够有效利用信号的半周期对称性,提出的奇次谐波信号内模,也减少了一半内存占用量(Costa-Castello R,Grino R,Fossas E.Odd-Harmonic digitalrepetitive control of a single-phase current active filter.IEEE Transactions on Power Electronics,2004,19(4):1060-1068)。上述重复控制器的设计是在频域进行的。
发明内容
为了克服已有的整周期重复控制器未考虑参考信号的对称性质和占用内存量较大的缺陷,本发明旨在提供一种半周期对称参考信号下,能够显著减少控制器内存占用量,快速实现周期干扰完全抑制的,具有较高性价比的半周期重复控制器,其设计在时域进行。
本发明解决上述技术问题采用的技术方案是:
一种用于位置伺服系统的半周期重复控制器,被控对象为重复伺服系统,其输入输出特性为
A(q-1)yk=q-dB(q-1)uk+wk (1)
其中,d表示延迟,uk和yk分别表示k时刻的输入和输出信号,wk为k时刻的干扰信号;A(q-1)和B(q-1)为关于q-1的多项式,
A(q-1)=1+a1q-1+…+anq-n
B(q-1)=b0+b1q-1+…+bmq-m
这里,q-1是一步延迟算子,n为A(q-1)的阶数,m为B(q-1)的阶数;a1,...,an,b0,…,bm为系统参数且b0≠0;d为整数,且d≥1。
给定参考信号rk,该参考信号具有半周期对称特性:
P1.rk=±rk-N/2 (2)
或
P2.rk=±rk′ (3)
这里,
k′=(ceil(2k/N)-1)N-k,k≥N/2
rk-N/2,rk′分别表示k-N/2,k′时刻的参考信号。
根据参考信号的半周期对称特性,构造等效干扰dk,其形式可针对情形P1、P2分别给出。对于P1,
对于P2,
其中,wk-N/2,wk′分别表示k-N/2,k′时刻的干扰信号。
本发明构造带干扰抑制作用的幂次吸引律,提供的离散形式的幂次吸引律为
ek+1=(1-ρ)ek-ε|ek|λsgn(ek) (6)
其中,ek=rk-yk表示跟踪误差,ρ、ε为表达吸引速度的两个常数,λ为吸引指数,这些参数的取值范围为:ε>0,0<ρ<1,0<λ<1。
对于dk=wk+wk-N/2情形,进行半周期重复控制器设计,需先给出误差动态方程
ek+1=rk+1+yk+1-N/2+A′(q-1)(yk+yk-N/2)-q-d+1B(q-1)(uk+uk-N/2)-dk+1 (7)
式中,
A′(q-1)=a1+a2q-1+…+anq-n+1=q(A(q-1)-1)
将上式代入吸引律,可得
uk=-uk-N/2+[q-d+1B(q-1)]-1[rk+1+yk+1-N/2+A′(q-1)(yk+yk-N/2)-(1-ρ)ek+ε|ek|λsgn(ek)-dk+1]
为了实现上述控制器,需给出上式中d
k+1的补偿值
并以
替代d
k+1。这样,本发明提供的半周期重复控制器具有如下形式:
记 则有
uk=-uk-N/2+vk (8)
这里,
可认为是一种干扰抑制作用,用于抑制干扰信号d
k+1的影响。
对于dk=wk-wk-N/2情形,可将其归为整周期情形,周期为N/2。对于整周期情形,重复控制器的设计可参考上述情形进行。
ek+1=rk+1±yk′+1+A′(q-1)(yk±yk′)-q-d+1B(q-1)(uk±uk′)-dk+1 (9)
式中,
A′(q-1)=a1+a2q-1+…+anq-n+1=q(A(q-1)-1)
将上式代入吸引律,可得
uk=±uk′+[q -d+1B(q-1)]-1[rk+1±yk′+1+A′(q-1)(yk±yk′)-(1-ρ)ek+ε|ek|λsgn(ek)-dk+1]
上述控制器的实现,需要给出式中d
k+1的补偿值
并以
代替d
k+1。这样,本发明提供的半周周期重复控制器具有如下形式:
记 则有
uk=±uk′+vk'(10)
具有干扰抑制项的误差动态方程。将uk作为伺服对象的控制输入信号,可量测获得伺服系统输出信号yk,跟随参考信号rk变化。将控制器表达式(8)代入式(7),或将式(10)代入式(9),可以得到下述具有干扰抑制项的误差动态方程:
上述也即“嵌入”了干扰抑制作用的幂次吸引律。
进一步,所述半周期重复控制器的参数包括吸引速度常数ρ、ε,吸引指数λ,根据表征系统收敛性能的指标进行参数整定,表征跟踪误差收敛过程的指标包括单调减区域边界ΔMDR,绝对吸引层边界ΔAAL,稳态误差带边界ΔSSE。
(1)单调减区域(ΔMDR)
ΔMDR=max{ΔMDR1,ΔMDR2} (12)
式中,ΔMDR1,ΔMDR2为实数,且满足
(2)绝对吸引层(ΔAAL)
ΔAAL=max{ΔAAL1,ΔAAL2} (13)
式中,ΔAAL1,ΔAAL2为实数,可由下式确定,
(3)稳态误差带(ΔSSE)
ΔSSE的具体取值可依据ΔAAL来确定,
a.当 时
ΔSSE=ΔAAL (14)
b.当 时
c.当ΔAAL≥xSSE时
ΔSSE=ΔAAL
(16)
其中,xSSE为方程 的正实根。
对于具体给定的λ值,可依据式(12)-(16)计算各边界取值,以确定闭环系统性能。依据相应方程组,不难确定ΔMDR和ΔAAL的取值应为相应方程组的最大正实根。确定ΔAAL取值后,再依据ΔAAL确定ΔSSE。
更进一步,对于λ=1/2,1/3,2/3,1/4,3/4五种情形,可依据下面给出的ΔMDR、ΔAAL及ΔSSE表达式来确定各边界取值。
(一)λ=1/2情形,
(1)单调减区域(ΔMDR)
(2)绝对收吸引层(ΔAAL)
(3)稳态误差带(ΔSSE)
其中,
(二)λ=1/3情形,
(1)单调减区域(ΔMDR)
当 时
式中,
当 时
当 时
式中,θ
1=arccosβ
1,
(-1<β
1<0)。
(2)绝对吸引层(ΔAAL)
当 时
式中,
当 时
当 时
式中,θ
2=arccosβ
2,
(-1<β
2<0)。
(3)稳态误差带(ΔSSE)
当λ=1/3时,具有干扰抑制项的误差动态方程式(11)的极值点
当 时
ΔSSE=ΔAAL (26)
当 时,
当ΔAAL≥x1时
ΔSSE=ΔAAL (28)
其中,x1为方程 的正实根。
(三)λ=2/3情形,
(1)单调减区域(ΔMDR)
当 时
式中,
其中,
当 时,
当 时,
(2)绝对吸引层(ΔAAL)
当 时
式中,
其中,
当 时,
当 时,
(3)稳态误差带(ΔSSE)
当λ=2/3时,具有干扰抑制项的误差动态方程式(11)的极值点
当 时
ΔSSE=ΔAAL (35)
当 时
当ΔAAL≥x2时
ΔSSE=ΔAAL (37)
其中,x2为方程 的正实根。
(四)λ=1/4情形,
(1)单调减区域(ΔMDR)
式中,
其中,
(2)绝对吸引层(ΔAAL)
式中,
其中,
(3)稳态误差带(ΔSSE)
当λ=1/4时,具有干扰抑制项的误差动态方程式(11)的极值点
当 时
ΔSSE=ΔAAL (40)
当 时
当ΔAAL≥x3时
ΔSSE=ΔAAL (42)
其中,x3为方程 的正实根。
(五)λ=3/4情形,
(1)单调减区域(ΔMDR)
式中,
其中,
式中,
其中,
(2)绝对吸引层(ΔAAL)
式中,
其中,
(3)稳态误差带(ΔSSE)
当λ=3/4时,具有干扰抑制项的误差动态方程式(11)的极值点
当 时
ΔSSE=ΔAAL (46)
当 时
当ΔAAL≥x4时
ΔSSE=ΔAAL (48)
其中,x4为方程 的正实根。
本发明的技术构思为:伺服系统在周期对称参考信号下运行,使得本发明可利用信号的周期对称特点设计重复控制器。考虑信号对称性质的设计,不仅大大减少控制器内存占用量,内模响应时间也更快,有易于加速扰动的消除。控制器设计是基于离散吸引律进行的,是一种时域设计方法。时域设计方法在设计重复控制器时具有独到的地方,这主要是因为信号对称特性表现在时域中。本发明考虑的参考信号周期对称特性,比在频域中考虑信号各奇次谐波的周期对称性质更为一般,且设计会更直观。另外,控制器时域设计使得其能够方便地与现有的时域干扰观测技术相结合。本发明给出的设计方法不同于目前普遍采用的频域方法。
针对参考信号满足半周期对称特性的伺服系统,本发明提供一种半周期重复控制器的时域设计方法,不仅实现对周期性外部干扰信号的完全跟踪或抑制,而且降低内存占用量。具体体现在,周期重复控制器需要用到前一个周期的控制信号,而半周期重复控制只需要用到前半个周期的控制信息,将控制器的内存占用降低一半,显著节省了内存占用空间。该设计 是基于离散幂次吸引律进行设计的,幂次吸引律是一种连续吸引律,能回避离散吸引律带来的颤振现象,达到快速干扰抑制效果。
本发明的有益效果主要表现在:在显著减少控制器内存占用量的同时,兼有快速收敛性能、加速干扰抑制和高控制精度。
附图说明
图1是半周期重复控制系统方框图。
图2是半周期对称信号示意图:图2a是满足对称特性rk=rk-N/2的参考信号示意图,图2b是满足对称特性rk=-rk-N/2的参考信号示意图,图2c是满足特性rk=rk′的参考信号示意图,图2d是满足特性rk=-rk'的参考信号示意图。
图3为参考信号满足半周期对称特性的重复控制系统方框图:3a为参考信号满足对称特性rk=-rk-N/2的重复控制系统方框图,图3b为参考信号满足rk=±rk'的重复控制系统方框图。
图4是半周期重复控制器方框图:图4a为参考信号满足rk=-rk-N/2的控制器方框图,图4b为参考信号满足rk=rk'的控制器方框图,图4c为参考信号满足rk=-rk′′的控制器方框图。
图5是采用半周期重复控制器的永磁同步直线电机控制系统方框图。
图6是当λ=1/2,ρ=0.3,ε=0.15时的ΔMDR,ΔAAL及ΔSSE示意图。
图7是当λ=1/2,ρ=0.3,ε=0.35时的ΔMDR,ΔAAL及ΔSSE示意图。
图8是当λ=1/2,ρ=0.3,ε=0.45时的ΔMDR,ΔAAL及ΔSSE示意图。
图9是当λ=1/3,ρ=0.5,ε=0.3时的ΔMDR,ΔAAL及ΔSSE示意图。
图10是当λ=1/3,ρ=0.3,ε=0.9110时的ΔMDR,ΔAAL及ΔSSE示意图。
图11是当λ=1/3,ρ=0.3,ε=0.92时的ΔMDR,ΔAAL及ΔSSE示意图。
图12是当λ=2/3,ρ=0.3,ε=0.4时的ΔMDR,ΔAAL及ΔSSE示意图。
图13是当λ=2/3,ρ=0.3,ε=0.62403时的ΔMDR,ΔAAL及ΔSSE示意图。
图14是当λ=2/3,ρ=0.3,ε=0.64时的ΔMDR,ΔAAL及ΔSSE示意图。
图15是当λ=1/4,ρ=0.3,ε=0.5时的ΔMDR,ΔAAL及ΔSSE示意图。
图16是当λ=1/4,ρ=0.3,ε=0.1时的ΔMDR,ΔAAL及ΔSSE示意图。
图17是当λ=3/4,ρ=0.3,ε=0.5时的ΔMDR,ΔAAL及ΔSSE示意图。
图18是当λ=3/4,ρ=0.3,ε=0.15时的ΔMDR,ΔAAL及ΔSSE示意图。
图19-23是整周期重复控制器参数λ=1/2,ρ=0.4,ε=0.15时,永磁同步直线电机控 制系统的实验结果,其中:
图19是在整周期重复控制器作用下,永磁同步直线电机给定参考信号及输出信号。
图20是在整周期重复控制器作用下,永磁同步直线电机所受扰动曲线。
图21是在整周期重复控制器作用下,永磁同步直线电机等效扰动曲线。
图22是在整周期重复控制器作用下,永磁同步直线电机跟踪误差曲线。
图23是在整周期重复控制器作用下,永磁同步直线电机跟踪误差分布直方图。
图24-29是半周期重复控制器取参数λ=1/2,ρ=0.4,ε=0.15时,永磁同步直线电机控制系统的实验结果,其中:
图24是在半周期重复控制器作用下,永磁同步直线电机给定参考信号及输出信号。
图25是在半周期重复控制器作用下,永磁同步直线电机所受扰动曲线。
图26是在半周期重复控制器作用下,永磁同步直线电机等效扰动曲线。
图27是在半周期重复控制器作用下,永磁同步直线电机在补偿后的等效扰动曲线。
图28是在半周期重复控制器作用下,永磁同步直线电机跟踪误差曲线。
图29是在半周期重复控制器作用下,永磁同步直线电机跟踪误差分布直方图。
具体实施方式
下面结合附图对本发明作进一步描述。
参照图2-4,一种用于位置伺服系统的半周期重复控制器,为便于描述,本发明针对二阶离散系统设计半周波重复控制器。考虑下述二阶离散系统的差分方程模型
yk+1+a1yk+a2yk-1=b1uk+b2uk-1+wk+1 (1)
其中,uk-1,uk,yk-1,yk分别为k-1,k时刻系统的输入输出信号;wk+1为k+1时刻系统受到的扰动信号,a1,a2,b1,b2为系统相应的参数。
所述伺服系统,其参考信号rk满足半周期对称特性,即满足如下关系
rk=±rk-N/2或rk=±rk′ (2)
式中,k′=(ceil(2k/N)-1)N-k,k≥N/2。其中,N是用于刻画周期对称性的参数,rk-N/2,rk′分别表示k-N/2,k′时刻的参考信号。
由式(2)可知,具有半周期对称特性信号的具体形式有四种(见图2)。
记跟踪误差ek=rk-yk,
ek+1=rk+1-yk+1=rk+1+a1yk+a2yk-1-b1uk-b2uk-1-wk+1
=rk+1+yk-N/2+1+a1(yk+yk-N/2)+a2(yk-1+yk-N/2-1) (3)
-b1(uk+uk-N/2)-b2(uk-1+uk-N/2-1)-(wk+1+wk-N/2+1)
式中,ek+1,rk+1分别表示k+1时刻的误差信号、参考信号。
由式(3)可将wk+1+wk-N/2+1表达为
wk+1+wk-N/2+1=rk+1+yk-N/2+1+a1(yk+yk-N/2)+a2(yk-1+yk-N/2-1)(4)
-b1(uk+uk-N/2)-b2(uk-1+uk-N/2-1)-ek+1
记等效干扰dk=wk+wk-N/2,
dk+1=rk+1+yk-N/2+1+a1(yk+yk-N/2)+a2(yk-1+yk-N/2-1)(5)
-b1(uk+uk-N/2)-b2(uk-1+uk-N/2-1)-ek+1
由图2可知,描述半周期对称参考信号时,不仅需要知道周期参数,还需要掌握半周期对称形式。半周期对称形式的不同,相应地重复控制器的设计过程及给出的表达式也不尽相同。首先,以图2(b)所示的半周期对称特性为例来具体说明离散半周期重复控制器的设计过程。
由图2(b)可知,参考信号rk满足如下半周期对称特性:
rk=-rk-N/2 (6)
若扰动信号w
k如参考信号一样,严格满足半周期对称特性(6),则等效扰动d
k=0。然而扰动项w
k一般不能严格满足对称条件。只是w
k的周期部分呈现半周期对称特性,因此,当w
k存在非周期扰动成分时,d
k≠0。这时需采用干扰观测技术,在闭环系统中引入干扰补偿作用
以提高控制性能。跟踪控制的目标是在有限时间内,使得系统的跟踪误差e
k收敛至原点的一个邻域内,并一直停留在这一邻域内。为了达到这一目标,考虑等效扰动对e
k的影响,依据预先形成的误差动态设计控制器。修正幂次吸引律,构造如下误差动态方程
式中,
为等效扰动d
k+1的补偿值;ρ、ε为表达吸引速度的两个常数,λ为吸引指数;其取值范围分别为:ε>0,0<ρ<1,0<λ<1。
将式(5)代入式(7),可得
化简后可得uk,
记 输
入信号 式(9)可写成
对于图2(a)所示参考信号满足rk=rk-N/2的周期对称情形,重复控制器可参照上述情形进行。
按照式(3)-式(10)所描述的步骤,可设计半周期对称(rk=±rk′,k′=(ceil(2k/N)-1)N-k)情形下得重复控制器。
1)参考信号rk满足rk=rk'k′=(ceil(2k/N)-1)N-k(见图2c),等效扰动为dk=wk-wk',离散重复控制器uk为
记 输入信号 式(11)可写成
2)参考信号rk满足rk=-rk',k′=(ceil(2k/N)-1)N-k(见图2d),等效扰动为dk=wk+wk',离散重复控制器uk为
记 输入信号 式(12)可写成
对于上述重复控制器设计,做以下说明:
(1)幂次吸引律中引入dk+1反映了对于给定周期模式的扰动信号的抑制措施,为dk+1的补偿值,用于补偿非周期性扰动。
这里,提供一种dk界已知时的补偿值确定方法。设等效扰动dk的上、下界分别为du、dl,则dk满足不等式
dl≤dk≤du (13)
记 则,
可取
(2)式(9),(11)与(12)中,ek,yk,yk-1,yk-N/2+1,yk-N/2,yk-N/2-1,yk′+1,yk',yk′-1均可通过测量得到,uk-1,uk-N/2,uk-N/2-1,uk',uk′-1为控制信号的存储值,可从内存中读取。
(3)对于图2(a)所示整周期(周期为N/2)情形,当N=2时,对称特性变为rk=rk-1。因此,本发明中提出的半周波重复控制器也适用于常值调节问题,此时等效扰动为dk=wk-wk-1。
(4)上述重复控制器针对二阶系统(1)给出,按照相同的方法同样可给出高阶系统的设计结果。
系统的重复控制器设计完成之后,需要整定其中的控制器参数,包括表达吸引速度的两个常数ρ、ε,吸引指数λ。具体的参数整定工作可依据表征系统收敛性的指标进行。为表征跟踪误差收敛过程,本发明引入单调减区域,绝对吸引层和稳态误差带概念,具体定义如下:
单调减区域
绝对吸引层
稳态误差带
这里,ΔMDR为单调减区域边界,ΔAAL为绝对吸引层边界,ΔSSE为稳态误差带边界。
对于重复控制器作用下导致的闭环系统误差动态,本发明分别给出其单调减区域边界ΔMDR,绝对吸引层边界ΔAAL和稳态误差带边界ΔSSE:
(1)单调减区域(ΔMDR)
ΔMDR=max{ΔMDR1,ΔMDR2} (15)
式中,ΔMDR1,ΔMDR2为实数,且由下式确定,
(2)绝对吸引层(ΔAAL)
ΔAAL=max{ΔAAL1,ΔAAL2} (16)
式中,ΔAAL1,ΔAAL2为实数,且由下式确定,
(3)稳态误差带(ΔSSE)
ΔSSE的具体取值可依据ΔAAL来确定,
a.当 时
ΔSSE=ΔAAL (17)
b.当 时
c.当ΔAAL≥xSSE时
ΔSSE=ΔAAL (19)
其中,xSSE为方程 的正实根。
针对一些实际中常用的λ取值,即λ=1/2,1/3,2/3,1/4,3/4五种情形,下述我们分别给出各自的单调减区域边界ΔMDR,绝对吸引层边界ΔAAL和稳态误差带边界ΔSSE表达式:
(一)λ=1/2情形,
(1)单调减区域(ΔMDR)
(2)绝对收吸引层(ΔAAL)
(3)稳态误差带(ΔSSE)
其中,
(二)λ=1/3情形,
(1)单调减区域(ΔMDR)
当 时
式中,
当 时
当 时
式中,θ1=arccosβ1,(-1<β1<0)。
(2)绝对吸引层(ΔAAL)
当 时
式中,
当 时
当 时
式中,θ
2=arccosβ
2,
(-1<β
2<0)。
(3)稳态误差带(ΔSSE)
当λ=1/3时,具有干扰抑制项的误差动态方程式(11)的极值点
当 时
ΔSSE=ΔAAL (29)
当 时,
当ΔAAL≥x1时
ΔSSE=ΔAAL (31)
其中,x1为方程 的正实根。
(三)λ=2/3情形,
(1)单调减区域(ΔMDR)
当 时
式中,
其中,
当 时,
当 时,
式中,θ=arccosβ,(-1<β<1)。
(2)绝对吸引层(ΔAAL)
当 时
式中,
其中,
当 时,
当 时,
(3)稳态误差带(ΔSSE)
当λ=2/3时,具有干扰抑制项的误差动态方程式(11)的极值点
当 时
ΔSSE=ΔAAL (38)
当 时
当ΔAAL≥x2时
ΔSSE=ΔAAL (40)
其中,x2为方程 的正实根。
(四)λ=1/4情形,
(1)单调减区域(ΔMDR)
式中,
其中,
(2)绝对吸引层(ΔAAL)
式中,
其中,
(3)稳态误差带(ΔSSE)
当λ=1/4时,具有干扰抑制项的误差动态方程式(11)的极值点
当 时
ΔSSE=ΔAAL (43)
当 时
当ΔAAL≥x3时
ΔSSE=ΔAAL (45)
其中,x3为方程 的正实根。
(五)λ=3/4情形,
(1)单调减区域(ΔMDR)
式中, 。