CN110673477A - 一种采用干扰差分抑制策略的离散时间幂次吸引控制器设计方法 - Google Patents

一种采用干扰差分抑制策略的离散时间幂次吸引控制器设计方法 Download PDF

Info

Publication number
CN110673477A
CN110673477A CN201910920607.4A CN201910920607A CN110673477A CN 110673477 A CN110673477 A CN 110673477A CN 201910920607 A CN201910920607 A CN 201910920607A CN 110673477 A CN110673477 A CN 110673477A
Authority
CN
China
Prior art keywords
interference
error
time
controller
reference signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201910920607.4A
Other languages
English (en)
Inventor
孙明轩
李旺
王晗
邹胜祥
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University of Technology ZJUT
Original Assignee
Zhejiang University of Technology ZJUT
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang University of Technology ZJUT filed Critical Zhejiang University of Technology ZJUT
Priority to CN201910920607.4A priority Critical patent/CN110673477A/zh
Publication of CN110673477A publication Critical patent/CN110673477A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Feedback Control In General (AREA)

Abstract

一种采用干扰差分抑制策略的离散时间幂次吸引控制器设计方法,给定模块产生参考信号;依据给定参考信号的具体形式,构造相应的干扰差分补偿反馈环节,其输出信号用于数字控制器中的干扰补偿;基于幂次吸引律构建理想误差动态,依据理想误差动态设计数字控制器,将当前的控制器计算获得的信号作为伺服对象的输入。具体的控制器参数整定可依据表征系统收敛性能的指标进行,提供了表征跟踪误差收敛过程的稳态误差带、绝对吸引层、单调减区域和系统跟踪误差首次进入稳态误差带的最大收敛步数的具体表达式。本发明提供的幂次吸引控制器设计方法,根据给定参考信号,采用相应干扰差分补偿措施,通过抑制干扰提高伺服系统跟踪精度。

Description

一种采用干扰差分抑制策略的离散时间幂次吸引控制器设计 方法
技术领域
本发明涉及一种采用干扰差分抑制策略的离散时间幂次吸引控制器设计方法,该方法适用于位置伺服系统,也适用于其他含有周期运行过程的工业场合。
背景技术
吸引律方法有别于离散滑模控制的趋近律方法。两者的主要区别表现在:吸引律方法将跟踪误差取代切换函数、原点取代切换面;趋近律方法要求有限时间达到切换面,而吸引律方法要求有限时间达到原点;吸引律方法设计的闭环系统仍具有关于参数漂移和外部干扰的鲁棒性能,只是滑模控制注重滑模运动的不变性,而吸引律方法追求系统稳态的不变性。
吸引律方法直接采用跟踪误差信号,控制器设计更为直接、简洁。吸引律反映了不考虑扰动时期望的系统误差动态特性;在存在干扰的情形下,直接依据吸引律设计的控制器无法实现。可将干扰抑制措施“嵌入”吸引律,构建具有扰动抑制作用的理想误差动态。依据构造的理想误差动态方程设计数字控制器,闭环系统动态过程由理想误差动态所决定,且具有理想误差动态所表征的期望跟踪性能。
以吸引律方法设计离散控制器时,刻画跟踪误差瞬态和稳态行为的指标可由理想误差动态给出,具体有下述四个指标:稳态误差带、绝对吸引层、单调递减区域以及首次收敛到稳态误差带所需的最大步数。四个指标的具体取值依赖于控制器参数和等效干扰信号的界,控制器参数和等效干扰信号的界不同,四个指标的取值也不同。一旦给定理想误差动态形式,可预先给出四个指标的具体表达式,用于控制器参数整定。
发明内容
本发明提出一种采用干扰差分抑制策略的离散时间幂次吸引控制器设计方法。为使得闭环系统具有预先设定的期望误差跟踪性能,采用干扰差分抑制策略,将其嵌入到幂次吸引律中,构建具有扰动抑制能力的理想误差动态,从而设计出控制器,使得闭环系统具有理想误差动态所刻画的特性,以提高控制性能,使得电机伺服系统实现高速、高精度跟踪。本发明具体给出稳态误差带、绝对吸引层、单调减区域以及首次收敛到稳态误差带所需的最大步数四个指标的具体表达式,皆可用于指导控制器参数整定。
本发明解决上述技术问题采用的技术方案是:
一种采用干扰差分抑制策略的离散时间幂次吸引控制器设计方法,包括以下步骤:
步骤1.给定参考信号rk为时间变量k的多项式,M表示该多项式的最高幂次,三种参考信号如下:
1)方波信号,M=0
Figure BDA0002217435960000021
2)三角波信号,M=1
Figure BDA0002217435960000022
3)S曲线,M=3
Figure BDA0002217435960000023
其中,A为幅值,N为参考信号在一个周期中的采样次数;
步骤2.构造理想误差动态,给出幂次吸引律,
ek+1=(1-ρ)ek-ε|ek|αsgn(ek) (4)
其中,ek=rk-yk表示k时刻,ek+1表示k+1时刻的跟踪误差,ρ、ε和分别表示吸引速度的常数,α为吸引指数,且相应的取值范围为:0<ρ<1,ε>0,0<α<1。
将干扰抑制措施嵌入上述吸引律中,构造如下理想误差动态
其中,dk+1为k+1时刻的等效干扰,
Figure BDA0002217435960000025
用于补偿等效干扰,为等效干扰补偿误差。
步骤3.干扰差分补偿策略
取等效干扰补偿作用为针对具体的参考信号定义等效干扰,其形式为
M=0时的等效干扰信号
dk+1=wk+1 (6)
M=1时的等效干扰信号
dk+1=wk+1-wk (7)
M=3时的等效干扰信号
dk+1=[(wk+1-wk)-(wk-wk-1)]-[(wk-wk-1)-(wk-1-wk-2)] (8)
其中,wk+1为k+1时刻的干扰。
定义干扰差分步数为L,L表示干扰补偿误差包含L个相继时刻的干扰,如式(6),两步干扰差分dk+1-dk=wk+1-wk,包含wk+1与wk两个时刻干扰,为有效抑制干扰,在选择等效干扰时,应满足如下条件
其中,
Figure BDA0002217435960000033
为不小于·的最小整数;
步骤4.控制器设计
依据理想误差动态(5)和等效干扰信号dk+1,给出如下控制器的表达式:
1)对于方波参考信号式
Figure BDA0002217435960000034
其中,F(q-1)=(q-d+1B(q-1)-b1)。
2)对于三角波参考信号式
Figure BDA0002217435960000035
其中,
Figure BDA0002217435960000036
3)对于S曲线参考信号式
Figure BDA0002217435960000041
其中
Figure BDA0002217435960000042
式(10)至式(12)中,A(q-1)、B(q-1)为伺服系统
A(q-1)yk=q-dB(q-1)uk+wk (13)
关于q-1的参数多项式:
A(q-1)=1+a1q-1+a2q-2+……+anq-n
B(q-1)=b1q-1+b2q-2+……+bmq-m
其中,uk与yk分别为伺服系统k时刻的输入及输出信号,d表示延迟因子,q-1为一步延迟算子,m、n分别为A(q-1)、B(q-1)的阶数,b1≠0,1≤m≤n。
进一步,所述方法还包括:
步骤5.性能分析
给出稳态误差带、绝对吸引层、单调减区域以及跟踪误差首次进入稳态误差带所需最多步数四个指标的具体表达式,用于刻画系统跟踪性能及指导控制器参数整定,其中,稳态误差带、绝对吸引层、单调减区域以及跟踪误差首次进入稳态误差带所需最多步数定义如下:
单调减区域ΔMDR:当ek大0于此边界时,ek同号递减,即满足如下条件:
Figure BDA0002217435960000043
绝对吸引层ΔAAL:当系统跟踪误差的绝对值|ek|大于此界时,其|ek|单调递减,即满足如条件:
Figure BDA0002217435960000044
稳态误差带ΔSSE:当系统误差一旦收敛进入该边界,那么误差就会稳定在此区域内,即满足如下条件:
Figure BDA0002217435960000059
最大收敛步数
Figure BDA0002217435960000051
跟踪误差最多经过
Figure BDA0002217435960000052
步进入稳态误差带。
等效干扰补偿误差满足
Figure BDA0002217435960000053
时,各指标的表达式如下:
(1)单调减区域ΔMDR
ΔMDR=max{ΔMDR1MDR2} (17)
Figure BDA0002217435960000054
其中,ΔMDR1和ΔMDR2均为实数,且由式(18)确定。
(2)绝对吸引层ΔAAL
ΔAAL=max{ΔAAL1AAL2} (19)
Figure BDA0002217435960000055
其中,ΔAAL1和ΔAAL2均为实数,且由式(20)确定。
(3)稳态误差带ΔSSE
ΔSSE=max{ΔSSE1SSE2} (21)
Figure BDA0002217435960000056
其中,ΔSSE1和ΔSSE2均为实数,且由式(22)确定。
针对于
Figure BDA0002217435960000057
各指标的表达式如下:
单调减区域ΔMDR
ΔMDR=max{ΔMDR1MDR2} (23)
1)当
Figure BDA0002217435960000058
Figure BDA0002217435960000061
2)当
Figure BDA0002217435960000064
3)当
Figure BDA0002217435960000066
其中
绝对吸引层ΔAAL
ΔAAL=max{ΔAAL1AAL2} (30)
1)当
Figure BDA0002217435960000071
Figure BDA0002217435960000072
Figure BDA0002217435960000073
2)当
Figure BDA0002217435960000076
3)当
其中
Figure BDA0002217435960000081
稳态误差带ΔSSE
1)当
Figure BDA0002217435960000082
或者ΔAAL≥δSSE
ΔSSE=ΔAAL (37)
2)
Figure BDA0002217435960000083
Figure BDA0002217435960000084
其中δSSE是方程的正实根;
收敛步数
Figure BDA00022174359600000810
Figure BDA0002217435960000087
其中,e0为跟踪误差初始值,表示不小于·的最小整数。
本发明技术构思为:提供一种采用干扰差分抑制策略的离散时间幂次吸引控制器设计方法。根据给定参考信号定义等效干扰,将干扰抑制措施嵌入幂次吸引律中,形成具有干扰抑制作用的理想误差动态。依据理想误差动态设计数字控制器,以实现对给定参考信号的精确跟踪。
本发明的有益效果主要表现在:根据给定参考信号,采用相应干扰差分补偿措施,通过抑制干扰提高跟踪精度。
附图说明
图1为伺服系统框图。
图2—图4为rk=10sin(2πfkTs)deg,f=1Hz,Ts=0.005s,Δ=0.3时,采用控制器式(33)的数值仿真结果,其中,图2是ρ=0.3、ε=0.1时的跟踪误差信号ek;图3是ρ=0.5、ε=0.15时的跟踪误差信号ek;图4是干扰误差补偿
图5—图7为参考信号rk如式(1),A=20,ρ=0.3、ε=0.1时,采用控制器式(33)的数值仿真结果,其中,图5是参考信号rk与输出信号yk;图6是跟踪误差信号ek;图7是干扰误差补偿
Figure BDA0002217435960000091
图8—图10为参考信号rk如式(2),A=20,ρ=0.3、ε=0.1时,采用控制器式(33)的数值仿真,其中,图8是参考信号rk与输出信号yk;图9是跟踪误差信号ek;图10是干扰误差补偿
Figure BDA0002217435960000092
图11—图13为参考信号rk如式(2),A=20,ρ=0.3、ε=0.1时,采用控制器式(34)的数值仿真结果,其中,图11是参考信号rk与输出信号yk;图12是跟踪误差信号ek;图13是干扰误差补偿
Figure BDA0002217435960000093
图14—图16为参考信号rk如式(3),A=20,ρ=0.3、ε=0.1时,采用控制器式(33)的数值仿真结果,其中,图14是参考信号rk与输出信号yk;图15是跟踪误差信号ek;图16是干扰误差补偿
Figure BDA0002217435960000094
图17—图19为参考信号rk如式(3),A=20,ρ=0.3、ε=0.1时,采用控制器式(34)的数值仿真结果,其中,图17是参考信号rk与输出信号yk;图18是跟踪误差信号ek;图19是干扰误差补偿
图20—图22为参考信号rk如式(3),A=20,ρ=0.3、ε=0.1时,采用控制器式(35)的数值仿真结果,其中,图20是参考信号rk与输出信号yk;图21是跟踪误差信号ek;图22是干扰误差补偿
Figure BDA0002217435960000096
图23—图25为参考信号rk如式(1),A=15deg,Ts=5ms,ρ=0.6、ε=0.3时,采用控制器式(33)的实验结果,其中,图23是参考信号rk与输出信号yk、图24是跟踪误差信号ek、图25是跟踪误差信号ek的直方图。
图26—图28为参考信号rk如式(2),A=110deg,Ts=5ms,ρ=0.6、ε=0.3时,采用控制器式(33)的实验结果,其中,图26是参考信号rk与输出信号yk、图27是跟踪误差信号ek、图28是跟踪误差信号ek的直方图。
图29—图31为参考信号rk如式(2),A=110deg,Ts=5ms,ρ=0.6、ε=0.3时,采用控制器式(34)的实验结果,其中,图29是参考信号rk与输出信号yk、图30是跟踪误差信号ek、图31是跟踪误差信号ek的直方图。
图32—图34为参考信号rk如式(3),A=60deg,Ts=5ms,ρ=0.6、ε=0.3时,采用控制器式(33)的实验结果,其中,图32是参考信号rk与输出信号yk、图33是跟踪误差信号ek、图34是跟踪误差信号ek的直方图。
图35—图37为参考信号rk如式(3),A=60deg,Ts=5ms,ρ=0.6、ε=0.3时,采用控制器式(34)的实验结果,其中,图35是参考信号rk与输出信号yk、图36是跟踪误差信号ek、图37是跟踪误差信号ek的直方图。
图38—图40为参考信号rk如式(3),A=60deg,Ts=5ms,ρ=0.6、ε=0.3时,采用控制器式(35)的实验结果,其中,图38是参考信号rk与输出信号yk、图39是跟踪误差信号ek、图40是跟踪误差信号ek的直方图。
具体实施方式
下面结合附图对本发明具体实施方式做进一步描述。
参照图1,一种采用干扰差分抑制策略的离散时间控制器幂次吸引律设计方法,包括以下步骤:
参考信号rk为时间变量k的多项式,M表示该多项式的最高幂次;三种参考信号如下:
1)方波信号,M=0
Figure BDA0002217435960000101
2)三角波信号,M=1
Figure BDA0002217435960000102
3)S曲线,M=3
Figure BDA0002217435960000103
其中,A为幅值,N为参考信号在一个周期中的采样次数;
步骤2.构造理想误差动态
针对幂次吸引律
ek+1=(1-ρ)ek-ε|ek|αsgn(ek) (4)
其中,ek=rk-yk表示k时刻,ek+1表示k+1时刻的跟踪误差,ρ、ε和分别表示吸引速度的常数,α为吸引指数,且相应的取值范围为:0<ρ<1,ε>0;
将干扰抑制措施嵌入上述吸引律中,构造如下理想误差动态
其中,dk+1为k+1时刻的等效干扰,
Figure BDA0002217435960000112
用于补偿等效干扰,
Figure BDA0002217435960000113
为等效干扰补偿误差;
步骤3.干扰差分补偿策略
取等效干扰补偿作用为
Figure BDA0002217435960000114
针对具体的参考信号定义等效干扰,其形式为
M=0时的等效干扰信号
dk+1=wk+1 (6)
M=1时的等效干扰信号
dk+1=wk+1-wk (7)
M=3时的等效干扰信号
dk+1=[(wk+1-wk)-(wk-wk-1)]-[(wk-wk-1)-(wk-1-wk-2)] (8)
其中,wk+1为k+1时刻的干扰;
定义干扰差分步数为L,L表示干扰补偿误差包含L个相继时刻的干扰,如式(6),两步干扰差分dk+1-dk=wk+1-wk,包含wk+1与wk两个时刻干扰,为有效抑制干扰,在选择等效干扰时,应满足如下条件
Figure BDA0002217435960000115
其中,
Figure BDA0002217435960000116
为不小于·的最小整数;
步骤4.控制器设计
依据理想误差动态(5)和等效干扰信号dk+1,给出如下控制器的表达式:
1)对于方波参考信号式
Figure BDA0002217435960000121
其中,F(q-1)=(q-d+1B(q-1)-b1)。
2)对于三角波参考信号式
Figure BDA0002217435960000122
其中,
Figure BDA0002217435960000123
3)对于S曲线参考信号式
Figure BDA0002217435960000124
其中
Figure BDA0002217435960000125
式(10)至式(12)中,A(q-1)、B(q-1)为伺服系统
A(q-1)yk=q-dB(q-1)uk+wk (13)
关于q-1的参数多项式:
A(q-1)=1+a1q-1+a2q-2+……+anq-n
B(q-1)=b1q-1+b2q-2+……+bmq-m
其中,uk与yk分别为伺服系统k时刻的输入及输出信号,d表示延迟因子,q-1为一步延迟算子,m、n分别为A(q-1)、B(q-1)的阶数,b1≠0,1≤m≤n。
进一步,所述方法还包括:
步骤5.性能分析
给出稳态误差带、绝对吸引层、单调减区域以及跟踪误差首次进入稳态误差带所需最多步数四个指标的具体表达式,用于刻画系统跟踪性能及指导控制器参数整定,其中,稳态误差带、绝对吸引层、单调减区域以及跟踪误差首次进入稳态误差带所需最多步数定义如下:
单调减区域ΔMDR:当ek大0于此边界时,ek同号递减,即满足如下条件:
Figure BDA0002217435960000131
绝对吸引层ΔAAL:当系统跟踪误差的绝对值|ek|大于此界时,其|ek|单调递减,即满足如条件:
Figure BDA0002217435960000132
稳态误差带ΔSSE:当系统误差一旦收敛进入该边界,那么误差就会稳定在此区域内,即满足如下条件:
Figure BDA0002217435960000133
最大收敛步数
Figure BDA0002217435960000134
跟踪误差最多经过
Figure BDA0002217435960000135
步进入稳态误差带。
针对于
Figure BDA0002217435960000136
时,各指标的表达式如下:
单调减区域ΔMDR
ΔMDR=max{ΔMDR1MDR2} (17)
1)当
Figure BDA0002217435960000137
Figure BDA0002217435960000139
2)当
Figure BDA00022174359600001310
Figure BDA00022174359600001311
Figure BDA0002217435960000141
3)当
Figure BDA0002217435960000142
Figure BDA0002217435960000143
其中
Figure BDA0002217435960000144
Figure BDA0002217435960000145
绝对吸引层ΔAAL
ΔAAL=max{ΔAAL1AAL2} (24)
1)当
Figure BDA0002217435960000146
Figure BDA0002217435960000147
Figure BDA0002217435960000148
2)当
Figure BDA0002217435960000151
Figure BDA0002217435960000152
Figure BDA0002217435960000153
3)当
Figure BDA0002217435960000154
Figure BDA0002217435960000155
其中
Figure BDA0002217435960000156
Figure BDA0002217435960000157
稳态误差带ΔSSE
1)当
Figure BDA0002217435960000158
或者ΔAAL≥δSSE
ΔSSE=ΔAAL (30)
2)
Figure BDA0002217435960000159
Figure BDA00022174359600001510
其中δSSE是方程
Figure BDA00022174359600001511
的正实根。
收敛步数
Figure BDA00022174359600001512
其中,e0为跟踪误差初始值,
Figure BDA0002217435960000162
表示不小于·的最小整数。
依据式(16)-(31)计算出各边界取值,以确定闭环系统的跟踪性能。
永磁同步电机装置执行位置精确跟踪任务,设计数字控制器用于位置环控制,其中电流环与速度环控制器由ELMO驱动器提供;位置环控制器由DSP开发板TMS320F2812提供。
通过最小二乘辨识方法,给出如下永磁同步伺服系统数学模型
yk+1-1.8949yk+0.8949yk-1=1.7908uk-0.5704uk-1+wk+1
对于M=0,由式(10)可得
Figure BDA0002217435960000163
对于M=1,由式(11)可得
Figure BDA0002217435960000164
对于M=3,由式(12)可得
Figure BDA0002217435960000165
通过数值仿真和实验结果验证本发明提供离散伺服系统中干扰差分补偿措施的有效性。
仿真分为两部分,第一部分验证三个边界值,第二部分验证干扰差分补偿措施的干扰抑制效果。
(1)给定参考信rk=10sin(2πfkTs)deg,频率f=1Hz,采样周期Ts=5ms,干扰为wk=0.15|mod(k,20)-10|+0.15|mod(k+5,20)-10|。在控制器式(33)作用下,参数ρ和ε不同取值,系统的性能指标也会不同,见图2—图3。
(i)当控制器参数为Δ=0.3,ρ=0.3,ε=0.1时(参见图2),性能指标为
ΔAAL=ΔSSE=ΔMDR=0.7035
(ii)当控制器参数为Δ=0.3,ρ=0.5,ε=0.15时(参见图3),性能指标为
ΔMDR=0.8889,ΔAAL=ΔSSE=0.3823
通过仿真,结果表明稳态误差带ΔSSE、绝对吸引层ΔAAL和单调减区域ΔMDR的准确性。
(2)参考信号分别为方波信号式(1)、三角波信号式(2)和S曲线式(3),幅值A分别为20、20、20。为验证干扰差分步数L在满式(9)时,设计的控制器能够实现对相应参考信号的精确跟踪,扰动信号选为wk=0.1rk,采样周期Ts=5ms,控制器参数ρ=0.3,ε=0.1。
1)参考信号rk为式(1),采用控制器式(33),仿真结果见图5—图7,图中,ΔSSE=0deg。
2)参考信号rk为式(2),采用控制器为式(33),仿真结果见图8—图10,图中,ΔSSE=0.02deg。
3)参考信号rk为式(2),采用控制器为式(34),仿真结果见图11—图13,图中,ΔSSE=0deg。
4)参考信号rk为式(3),采用控制器为式(33),仿真结果见图14—图16,图中,ΔSSE=0.03deg。
5)参考信号rk为式(3),采用控制器为式(34),仿真结果见图17—图19,图中,ΔSSE=0.002deg。
6)参考信号rk为式(3),采用控制器为式(35),仿真结果见图20—图22,图中,ΔSSE=0deg。
通过仿真(2),表明当M与L满足式(9)时,采用干扰差分补偿设计的控制器能够实现对给定参考信号的精确跟踪,当M与L不满足式(9)时,L越接近
Figure BDA0002217435960000171
跟踪效果越好,且控制过程无抖振。
在位置伺服装置上验证本发明提供的控制方法,图1为位置伺服系统框图。参考信号分别为方波信号式(1)、三角波信号式(2)和S曲线式(3)。分别采用控制器式(33)、式(34)和式(35)验证干扰差分补偿技术的效果。实验中采样周期Ts=5ms,控制器参数ρ=0.6,ε=0.5,参考信号式(1)、式(2)和式(3)的幅值A分别为15deg、110deg和60deg。实验结果如下:
1)参考信号rk如式(1),采用控制器(33),实验结果见图23—图25。系统跟踪误差为0.1deg。
2)参考信号rk如式(2),采用控制器(33),实验结果见图26—图28。系统跟踪误差为0.05deg。
3)参考信号rk如式(2),采用控制器(34),实验结果见图29—图31。系统跟踪误差为0.05deg。
4)参考信号rk如式(3),采用控制器(33),实验结果见图32—图34。系统跟踪误差为0.2deg。
5)参考信号rk如式(3),采用控制器(34),实验结果见图35—图37。系统跟踪误差为0.15deg。
6)参考信号rk如式(3),采用控制器(35),实验结果见图38—图40。系统跟踪误差为0.1deg。
实验结果表明,当M与L满足式(9)时,采用干扰差分补偿设计的控制器能够实现对给定参考信号的精确跟踪,当M与L不满足式(9)时,L越接近
Figure BDA0002217435960000181
系统跟踪性能越好。

Claims (2)

1.一种采用干扰差分抑制策略的离散时间幂次吸引控制器设计方法,其特征在于,所述方法包括以下步骤:
步骤1.给定参考信号rk为时间变量k的多项式,M表示该多项式的最高幂次,三种参考信号如下:
1)方波信号,M=0
2)三角波信号,M=1
Figure FDA0002217435950000012
3)S曲线,M=3
Figure FDA0002217435950000013
其中,A为幅值,N为参考信号在一个周期中的采样次数;
步骤2.构造理想误差动态
针对幂次吸引律
ek+1=(1-ρ)ek-ε|ek|αsgn(ek) (4)
其中,ek=rk-yk表示k时刻,ek+1表示k+1时刻的跟踪误差,ρ、ε和分别表示吸引速度的常数,α为吸引指数,且相应的取值范围为:0<ρ<1,ε>0;
将干扰抑制措施嵌入该吸引律,构造理想误差动态
Figure FDA0002217435950000021
其中,dk+1为k+1时刻的等效干扰,
Figure FDA0002217435950000022
用于补偿等效干扰,
Figure FDA0002217435950000023
为等效干扰补偿误差;
步骤3.干扰差分补偿策略
取等效干扰补偿作用为
Figure FDA0002217435950000024
针对具体的参考信号定义等效干扰,其形式为
M=0时的等效干扰信号
dk+1=wk+1 (6)
M=1时的等效干扰信号
dk+1=wk+1-wk (7)
M=3时的等效干扰信号
dk+1=[(wk+1-wk)-(wk-wk-1)]-[(wk-wk-1)-(wk-1-wk-2)] (8)
其中,wk+1为k+1时刻的干扰;
定义干扰差分步数为L,L表示干扰补偿误差包含L个相继时刻的干扰,如式(6),两步干扰差分dk+1-dk=wk+1-wk,包含wk+1与wk两个时刻干扰,为有效抑制干扰,在选择等效干扰时,应满足如下条件
Figure FDA0002217435950000025
其中,
Figure FDA0002217435950000026
为不小于·的最小整数;
步骤4.控制器设计
依据理想误差动态(5)和等效干扰信号dk+1,给出如下控制器的表达式:
1)对于方波参考信号式(1),
Figure FDA0002217435950000027
其中,F(q-1)=(q-d+1B(q-1)-b1);
2)对于三角波参考信号式(2),
Figure FDA0002217435950000028
其中,
Figure FDA0002217435950000031
3)对于S曲线参考信号式(3),
Figure FDA0002217435950000032
其中
Figure FDA0002217435950000033
式(10)至式(12)中,A(q-1)、B(q-1)为伺服系统
A(q-1)yk=q-dB(q-1)uk+wk (13)
关于q-1的参数多项式:
A(q-1)=1+a1q-1+a2q-2+……+anq-n
B(q-1)=b1q-1+b2q-2+……+bmq-m
其中,uk与yk分别为伺服系统k时刻的输入及输出信号,d表示延迟因子,q-1为一步延迟算子,m、n分别为A(q-1)、B(q-1)的阶数,b1≠0,1≤m≤n。
2.如权利要求1所述的一种采用干扰差分抑制策略的离散时间幂次吸引控制器设计方法,其特征在于,所述方法还包括:
步骤5.性能分析
给出稳态误差带、绝对吸引层、单调减区域以及跟踪误差首次进入稳态误差带所需最多步数四个指标的具体表达式,用于刻画系统跟踪性能及指导控制器参数整定,其中,稳态误差带、绝对吸引层、单调减区域以及跟踪误差首次进入稳态误差带所需最多步数定义如下:
1)稳态误差带ΔSSE
2)绝对吸引层ΔAAL
Figure FDA0002217435950000035
3)单调减区域ΔMDR
Figure FDA0002217435950000036
4)最大收敛步数跟踪误差最多经过
Figure FDA0002217435950000042
步进入稳态误差带;
等效干扰补偿误差满足时,针对于
Figure FDA0002217435950000044
时,各指标的表达式如下:
单调减区域ΔMDR
ΔMDR=max{ΔMDR1MDR2} (14)
1)当
Figure FDA0002217435950000046
Figure FDA0002217435950000047
2)当
Figure FDA0002217435950000048
Figure FDA0002217435950000049
Figure FDA00022174359500000410
3)当
Figure FDA00022174359500000411
Figure FDA00022174359500000412
其中
Figure FDA0002217435950000051
Figure FDA0002217435950000052
2)绝对吸引层ΔAAL
ΔAAL=max{ΔAAL1AAL2} (21)
1)当
Figure FDA0002217435950000054
Figure FDA0002217435950000055
2)当
Figure FDA0002217435950000056
Figure FDA0002217435950000057
Figure FDA0002217435950000058
3)当
Figure FDA0002217435950000059
Figure FDA0002217435950000061
其中
Figure FDA0002217435950000062
Figure FDA0002217435950000063
3)稳态误差带ΔSSE
1)当
Figure FDA0002217435950000064
或者ΔAAL≥δSSE
ΔSSE=ΔAAL (28)
2)
Figure FDA0002217435950000065
Figure FDA0002217435950000066
其中δSSE是方程
Figure FDA0002217435950000067
的正实根;
收敛步数
Figure FDA0002217435950000069
其中,e0为跟踪误差初始值,
Figure FDA00022174359500000610
表示不小于·的最小整数。
CN201910920607.4A 2019-09-27 2019-09-27 一种采用干扰差分抑制策略的离散时间幂次吸引控制器设计方法 Pending CN110673477A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910920607.4A CN110673477A (zh) 2019-09-27 2019-09-27 一种采用干扰差分抑制策略的离散时间幂次吸引控制器设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910920607.4A CN110673477A (zh) 2019-09-27 2019-09-27 一种采用干扰差分抑制策略的离散时间幂次吸引控制器设计方法

Publications (1)

Publication Number Publication Date
CN110673477A true CN110673477A (zh) 2020-01-10

Family

ID=69079510

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910920607.4A Pending CN110673477A (zh) 2019-09-27 2019-09-27 一种采用干扰差分抑制策略的离散时间幂次吸引控制器设计方法

Country Status (1)

Country Link
CN (1) CN110673477A (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108828958A (zh) * 2018-08-21 2018-11-16 浙江工业大学 一种采用干扰差分抑制策略的离散时间控制器无切换吸引律设计方法
CN110134010A (zh) * 2019-04-17 2019-08-16 浙江工业大学 一种采用等效扰动补偿伺服系统的幂次吸引重复控制方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108828958A (zh) * 2018-08-21 2018-11-16 浙江工业大学 一种采用干扰差分抑制策略的离散时间控制器无切换吸引律设计方法
CN110134010A (zh) * 2019-04-17 2019-08-16 浙江工业大学 一种采用等效扰动补偿伺服系统的幂次吸引重复控制方法

Similar Documents

Publication Publication Date Title
CN109188908B (zh) 基于指数型无切换吸引律的数字控制器设计方法
CN107544244B (zh) 基于椭圆吸引律和等效扰动扩张状态补偿的用于电机伺服系统的离散重复控制方法
CN107797448B (zh) 采用扰动扩张补偿的电机位置离散重复控制方法
CN109450307B (zh) 基于类正态分布吸引律和采用扰动扩张状态补偿的用于电机伺服系统的离散重复控制方法
CN107544245B (zh) 采用扰动扩张补偿的一阶惯性吸引律的用于电机伺服系统的离散重复控制方法
CN110134010B (zh) 一种采用等效扰动补偿伺服系统的幂次吸引重复控制方法
Hou et al. Discrete-time terminal sliding-mode tracking control with alleviated chattering
CN107121932B (zh) 电机伺服系统误差符号积分鲁棒自适应控制方法
CN107102634B (zh) 一种基于转台伺服系统的参数估计与跟踪控制方法
CN111546346B (zh) 一种柔性关节扰动观测方法、力矩控制方法和设备
CN108958041B (zh) 一种基于双曲正割吸引律的离散双周期重复控制方法
CN108828958B (zh) 一种采用干扰差分抑制策略的离散时间控制器无切换吸引律设计方法
CN110134014B (zh) 周期伺服系统幂次吸引重复控制的等效扰动补偿方法
CN105387859A (zh) Mems传感器组合温度漂移误差补偿方法
Durmaz et al. Sliding mode control for non-linear systems with adaptive sliding surfaces
CN110658719A (zh) 一种采用等效扰动补偿的伺服系统无切换吸引重复控制方法
CN110716430A (zh) 一种采用等效扰动补偿的伺服系统快速吸引重复控制方法
CN110389528A (zh) 基于扰动观测的数据驱动mems陀螺仪驱动控制方法
CN109100938B (zh) 一种采用干扰差分补偿的离散时间控制器吸引律设计方法
Trivedi et al. Roll stabilization: a higher-order sliding-mode approach
Cao et al. Robust flight control design using sensor-based backstepping control for unmanned aerial vehicles
CN111752158B (zh) 一种有限时间收敛的二阶滑模控制方法
CN110673477A (zh) 一种采用干扰差分抑制策略的离散时间幂次吸引控制器设计方法
CN104991445B (zh) 一种全局稳定的电机伺服系统自适应输出反馈鲁棒控制方法
Dorwarth et al. Nonlinear model order reduction for high Q MEMS gyroscopes

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20200110