CN109828469B - 相位优化的扩张状态观测器及抗扰控制系统 - Google Patents

相位优化的扩张状态观测器及抗扰控制系统 Download PDF

Info

Publication number
CN109828469B
CN109828469B CN201910183791.9A CN201910183791A CN109828469B CN 109828469 B CN109828469 B CN 109828469B CN 201910183791 A CN201910183791 A CN 201910183791A CN 109828469 B CN109828469 B CN 109828469B
Authority
CN
China
Prior art keywords
extended state
state observer
phase
output
optimized
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.)
Active
Application number
CN201910183791.9A
Other languages
English (en)
Other versions
CN109828469A (zh
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.)
Beijing Technology and Business University
Original Assignee
Beijing Technology and Business University
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 Beijing Technology and Business University filed Critical Beijing Technology and Business University
Priority to CN201910183791.9A priority Critical patent/CN109828469B/zh
Publication of CN109828469A publication Critical patent/CN109828469A/zh
Application granted granted Critical
Publication of CN109828469B publication Critical patent/CN109828469B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Feedback Control In General (AREA)

Abstract

本发明公开了一种相位优化的抗扰控制系统,包括:控制器、受控对象和相位优化的扩张状态观测器;其中,相位优化的扩张状态观测器包括扩张状态观测器和相位优化器;相位优化器包括乘法器和加法器;控制器的输出端接至扩张状态观测器的第一输入端,受控对象的输出接扩张状态观测器的第二输入端;扩张状态观测器的z3量输出端接加法器的第一输入端,z3的变化率输出端接乘法器,乘法器的输出端接加法器的第二输入端;加法器的输出端接控制器的第一反馈端,扩张状态观测器的z1、z2量输出端接控制器的第二反馈端。本发明的相位优化的扩张状态观测器及抗扰控制系统实现了对时变扰动的有效估计。

Description

相位优化的扩张状态观测器及抗扰控制系统
技术领域
本发明涉及先进控制技术领域,具体来讲涉及自抗扰控制中相位优化的扩张状态观测器及其设计方法以及相位优化的抗扰控制系统,解决了时变扰动估计准确性的问题。
背景技术
对未知、不确定性扰动的鲁棒性是控制器的重要性能指标之一,实际系统中,广泛存在着未知、不确定性扰动,并且这些扰动极难建模。目前,设计扰动估计器估计未知扰动,或者设计扩张状态观测器估计总扰动(包括内部参数变化、未建模动态以及外部未知、不确定扰动),通过控制律补偿所估扰动的扰动估计补偿方法和自抗扰控制方法非常有效。然而,对于时变扰动,现有扰动估计器的估计准确度较低,无法准确估计扰动,意味着无法保证闭环系统性能。
对于扩张状态观测器而言,要提高估计的准确度,必然提高观测器带宽。随之而来的问题是引入更多的量测噪声。特别地,扩张状态观测器仅能实现对常值扰动的稳态无差估计,对于时变扰动,其估计是有偏差的。
可见,基于扩张状态观测器的自抗扰控制性能极大地受限于扩张状态观测器的估计性能,为提高自抗扰控制或者基于扩张状态观测器的控制技术的性能,提高扩张状态观测器的估计准确度非常必要。
本发明为解决扩张状态观测器估计时变扰动偏差大的问题,提出了一种简单而有效的相位优化律,进而在几乎不增加计算成本的同时,设计出了相位优化的扩张状态观测器及抗扰控制系统,实现了对时变扰动的有效估计。
发明内容
为实现本发明之目的,采用以下技术方案予以实现:
一种相位优化的抗扰控制系统,包括:控制器、受控对象和相位优化的扩张状态观测器;其中,相位优化的扩张状态观测器包括扩张状态观测器和相位优化器;相位优化器包括乘法器和加法器;控制器的输出端接至扩张状态观测器的第一输入端,受控对象的输出接扩张状态观测器的第二输入端;扩张状态观测器的z3量输出端接加法器的第一输入端,z3的变化率输出端接乘法器,乘法器的输出端接加法器的第二输入端;加法器的输出端接控制器的第一反馈端,扩张状态观测器的z1、z2量输出端接控制器的第二反馈端。
所述的抗扰控制系统,其中:扩张状态观测器表示如下:
Figure BDA0001992161220000021
其中z1、z2、z3为扩张状态观测器的输出量,y是系统输出,z1是y的估计值,z2是y的估计值的变化率,z3是系统总干扰的估计值;d是外部扰动,u是控制输入,b0是控制增益,f是系统总干扰,ωo是观测器带宽。
所述的抗扰控制系统,其中:所述相位优化的扩张状态观测器设计为:
Figure BDA0001992161220000031
其中,z1,z2,z3PO为相位优化的扩张状态观测器的输出,z1,z2与扩张状态观测器的输出z1,z2相同,z3PO为扩张状态观测器的输出z3经相位优化律后所得的总扰动估计值。
所述的抗扰控制系统,其中控制器表示为:
Figure BDA0001992161220000032
其中u是控制输入,u0是控制器输出,b0是可调参数。
一种相位优化的扩张状态观测器的设计方法,包括如下步骤:
(1)设计扩张状态观测器:
Figure BDA0001992161220000033
其中z1、z2、z3为扩张状态观测器的输出量,y是输出,z1估计y,z2估计y的变化率,z3估计总扰动;d是外部扰动,u是控制输入,b0是控制增益,f是系统总干扰,ωo是观测器带宽;
(2)设计相位优化律
Figure BDA0001992161220000034
其中z3PO表示相位优化律的输出,
Figure BDA0001992161220000035
表示扩张状态观测器输出的总扰动估计值的变化率,c表示增益。
所述的扩张状态观测器的设计方法,其中::
根据式(3)得到z3的传递函数
Figure BDA0001992161220000036
其中s是拉普拉斯算子;
根据式(4)、(5),有z3PO的传递函数
Figure BDA0001992161220000041
令估计误差
Figure BDA0001992161220000042
取c=ωo/3可得:
Figure BDA0001992161220000043
Figure BDA0001992161220000044
y,u均取幅值为K的斜坡信号y(s)=K/s2,u(s)=K/s2,可得e3的稳态误差:
Figure BDA0001992161220000045
所述的扩张状态观测器的设计方法,其中:
根据式(7)、(8),用jω代替s,可得误差相位为
Figure BDA0001992161220000046
其中j为虚数单位,ω为角频率;
相位差为:
Figure BDA0001992161220000047
相位优化的扩张状态观测器设计为
Figure BDA0001992161220000051
其中,z1,z2,z3PO为相位优化的扩张状态观测器的输出,z1,z2与扩张状态观测器(3)的输出z1,z2相同,z3PO为扩张状态观测器(3)的输出z3经相位优化律后所得的总扰动估计值。
一种相位优化的抗扰控制系统的设计方法,包括如上一所述的相位优化的扩张状态观测器的设计方法,其特征在于还包括:
(3)设计设计控制器u0
Figure BDA0001992161220000052
其中,ωc为控制带宽,r是设定值。
基于相位优化的扩张状态观测器的控制律设计为
Figure BDA0001992161220000053
其中,z3PO为相位优化的扩张状态观测器的总扰动估计值。
附图说明
图1为相位优化的抗扰控制系统示意图;
图2为相位优化的抗扰控制系统设计流程示意图;
图3为控制和估计效果示意图。
具体实施方式
下面结合附图对本发明的具体实施方式进行详细说明:
本发明的相位优化的扩张状态观测器的设计方法包括:
(1)设计扩张状态观测器
许多物理系统,如机器人、航天器和机械系统,通常用二阶微分方程来描述,二阶非线性系统可描述为
Figure BDA0001992161220000061
其中y是系统输出,d是外部扰动,u是控制输入,b0是控制增益,f是系统总干扰,包含系统内部不确定性和外部扰动,将其扩展为系统的状态变量
Figure BDA0001992161220000062
则系统(1)状态方程为:
Figure BDA0001992161220000063
式中x1,x2,x3为系统状态变量。
建立扩张状态观测器,用式(3)表示:
Figure BDA0001992161220000064
选取合适观测器带宽ωo,扩张状态观测器即能实现对系统(2)中各状态的估计,即z1是y的估计值,z2是y的估计值的变化率,z3是系统总干扰的估计值,其中z1、z2、z3为扩张状态观测器的输出量(状态变量)。
(2)设计相位优化律
时变扰动下,扩张状态观测器无法给出准确估计结果。因此,需要设计一种快速、高效的相位优化的扩张状态观测器。本发明提出了一种相位优化律
Figure BDA0001992161220000065
其中z3PO表示相位优化律的输出,它是对时变扰动的估计。
Figure BDA0001992161220000066
表示扩张状态观测器输出的总扰动估计值的变化率,c表示增益。
根据式(3)可得z3的传递函数
Figure BDA0001992161220000071
其中s是拉普拉斯算子;
根据式(4)、(5),有z3PO的传递函数
Figure BDA0001992161220000072
令估计误差
Figure BDA0001992161220000073
取c=ωo/3可得
Figure BDA0001992161220000074
Figure BDA0001992161220000075
考虑分析的典型性,y,u均取幅值为K的斜坡信号y(s)=K/s2,u(s)=K/s2,可得e3的稳态误差:
Figure BDA0001992161220000076
可见,扩张状态观测器对于斜坡扰动估计有偏差,而使用相位优化律后,对斜坡扰动的稳态估计偏差为零。
根据式(7)、(8),用jω代替s,可得误差相位为
Figure BDA0001992161220000077
其中j为虚数单位,ω为角频率;
相位差为:
Figure BDA0001992161220000078
由于
Figure BDA0001992161220000079
可得
Figure BDA00019921612200000710
也就是说,z3PO的相位超前于z3
于是,相位优化的扩张状态观测器设计为
Figure BDA0001992161220000081
其中,z1,z2,z3PO为相位优化的扩张状态观测器的输出,z1,z2与扩张状态观测器(3)的输出z1,z2相同,z3PO为扩张状态观测器(3)的输出z3经相位优化律后所得的总扰动估计值。
进一步,为形成相位优化的抗扰控制系统,还包括以下步骤:
(3)设计控制器
设计控制器(比例微分控制器)u0
Figure BDA0001992161220000082
其中,ωc为控制带宽,r是设定值。
基于相位优化的扩张状态观测器的控制律设计为
Figure BDA0001992161220000083
其中,z3PO为相位优化的扩张状态观测器的总扰动估计值,u0为控制器的输出,b0为可调参数。当z3PO
Figure BDA0001992161220000084
估计准确时,系统(1)可简化为双积分器串联。
Figure BDA0001992161220000085
因采用相位优化律,使得只对常值扰动有很好估计效果的扩张状态观测器,对于时变扰动也有很好的估计效果。
如图1所示,相位优化的抗扰控制系统包括:控制器(Controller)、受控对象(Plant)和相位优化的扩张状态观测器(POESO);其中,相位优化的扩张状态观测器包括扩张状态观测器(ESO)和相位优化器(POL);相位优化器包括乘法器和加法器;控制器的输出端接至扩张状态观测器的第一输入端,受控对象的输出接扩张状态观测器的第二输入端;扩张状态观测器的z3量输出端接加法器的第一输入端,z3的变化率输出端接乘法器,乘法器的输出端接加法器的第二输入端;加法器的输出端接控制器的第一反馈端,扩张状态观测器的z1、z2量输出端接控制器的第二反馈端。
本发明的相位优化的抗扰控制系统的工作过程如下:将设定值r输入给控制器;控制器的输出u受外部扰动d影响,施加到受控对象中(如图1所示),控制器的输出u输入给扩张状态观测器(ESO);受控对象的输出y输入给扩张状态观测器(ESO);扩张状态观测器(ESO)的输出量z3的变化率经乘法器,乘以1/c后,输入给加法器并与扩张状态观测器(ESO)的输出z3量相加;加法器的输出量输入给控制器;扩张状态观测器(ESO)的输出量z1和z2输入给控制器。
扩张状态观测器的估计结果经相位优化律后,得到优化的总扰动估计值z3PO,它是相位优化的扩张状态观测器的一个输出。控制器根据设定值、系统状态估计值、总扰动估计值得到控制量。控制量施加给系统得到系统输出。将控制量与系统输出作为扩张状态观测器的输入,这样整个系统形成闭环。
相位优化的抗扰控制系统设计流程如图2所示:
步骤1:设计扩张状态观测器
按上述式(3),设计扩张状态观测器。其中,观测器带宽设为ωo=15。为测试不同观测器对不同类型干扰的估计效果,开始时不施加干扰。t=20s时施加阶跃干扰;t=40s时施加斜坡干扰;t=60s时施加正弦干扰,干扰如下式所示
Figure BDA0001992161220000091
步骤2:设计相位优化的相位优化的扩张状态观测器
按式(4)的相位优化律,优化扩张状态观测器对总扰动的估计值z3。其中,增益系数c=ωo/3=5。综合式(3)、(4),可得相位优化的扩张状态观测器,可用上述式12表示。
步骤3:设计控制律
根据第2步设计的相位优化的扩张状态观测器,得到对总扰动的优化估计值。由(13)设计控制器。其中,系统期望的系统输出为r=1,控制器带宽设为ωc=1。
步骤4:形成相位优化自抗扰控制
根据第3步得到的控制量,分别施加到被控对象与扩张观测器观测器中,形成闭环,可由式14表示。
仿真中,对比不使用观测器、使用扩张状态观测器和相位优化的扩张状态观测器的三种控制方法,即PD控制、基于扩张状态观测器的自抗扰控制(active disturbancerejection control,ADRC)和基于相位优化的扩张状态观测器的相位优化自抗扰控制(phase optimized active disturbance rejection control,POADRC)。采样时间为0.001秒,被控对象设置为两积分器串联,对比结果如图3所示。系统输出y如图3(a)所示,从中可见:首先,PD控制器可在没有干扰时,稳态无差地控制系统。然而,对系统施加干扰之后,可以看到它不能抵抗任何外扰。其次,ADRC可在没有干扰和存在常值扰动的情况下,稳态无差地控制系统。然而,可以看到它不能抵抗时变扰动。第三,POADRC可抵抗常值和时变扰动,稳态无差地控制系统。施加的扰动和估计值如图3(b)所示。可以观察到,相位优化的扩张状态观测器可准确估计时变扰动,而扩张状态观测器只能准确估计常值扰动。这是ADRC和POADRC控制性能差异的关键因素。

Claims (3)

1.一种相位优化的抗扰控制系统,包括:控制器、受控对象和相位优化的扩张状态观测器;其特征在于:相位优化的扩张状态观测器包括扩张状态观测器和相位优化器;相位优化器包括乘法器和加法器;控制器的输出端接至扩张状态观测器的第一输入端,受控对象的输出接扩张状态观测器的第二输入端;扩张状态观测器的z3量输出端接加法器的第一输入端,z3的变化率输出端接乘法器,乘法器的输出端接加法器的第二输入端;加法器的输出端接控制器的第一反馈端,扩张状态观测器的z1、z2量输出端接控制器的第二反馈端;相位优化的扩张状态观测器的相位优化律为:
Figure FDA0003228686220000011
其中z1是y的估计值,z2是y的估计值的变化率,z3是系统总干扰的估计值,z3PO表示相位优化律的输出,它是对时变扰动的估计,
Figure FDA0003228686220000012
表示扩张状态观测器输出的总扰动估计值的变化率,c表示增益;
基于相位优化的扩张状态观测器的控制律设计为:
Figure FDA0003228686220000013
其中,u0为控制器的输出,b0为可调参数。
2.一种如权利要求1所述的相位优化的扩张状态观测器的设计方法,其特征在于包括如下步骤:
(1)设计扩张状态观测器,用式(3)表示:
Figure FDA0003228686220000014
其中u是控制输入,y是系统输出,
选取合适观测器带宽ωo,扩张状态观测器即能实现对系统中各状态的估计,即z1是y的估计值,z2是y的估计值的变化率,z3是系统总干扰的估计值,其中z1、z2、z3为扩张状态观测器的输出量;
(2)设计相位优化律:
Figure FDA0003228686220000021
其中,
Figure FDA0003228686220000022
表示扩张状态观测器输出的总扰动估计值的变化率,c表示增益;
所述的扩张状态观测器的设计方法,其中:
根据式(3)得到z3的传递函数
Figure FDA0003228686220000023
其中s是拉普拉斯算子;
根据式(4)、(5),有z3PO的传递函数:
Figure FDA0003228686220000024
令估计误差
Figure FDA0003228686220000025
取c=ωo/3可得:
Figure FDA0003228686220000026
Figure FDA0003228686220000027
f是系统总干扰,d是外部扰动,y,u均取幅值为K的斜坡信号y(s)=K/s2,u(s)=K/s2,可得e3的稳态误差:
Figure FDA0003228686220000028
所述的扩张状态观测器的设计方法,其中:
根据式(7)、(8),用jω代替s,可得误差相位为:
Figure FDA0003228686220000031
其中j为虚数单位,ω为角频率;
相位差为:
Figure FDA0003228686220000032
相位优化的扩张状态观测器设计为:
Figure FDA0003228686220000033
其中,z1,z2与扩张状态观测器的输出z1,z2相同。
3.一种相位优化的抗扰控制系统的设计方法,包括根据权利要求2所述的相位优化的扩张状态观测器的设计方法,其特征在于还包括:
(3)设计控制器u0
Figure FDA0003228686220000034
其中,ωc为控制带宽,r是设定值,
基于相位优化的扩张状态观测器的控制律设计为:
Figure FDA0003228686220000035
其中,
u0为控制器的输出,b0为可调参数,当z3PO
Figure FDA0003228686220000041
估计准确时,相位优化的抗扰控制系统简化为双积分器串联:
Figure FDA0003228686220000042
CN201910183791.9A 2019-03-12 2019-03-12 相位优化的扩张状态观测器及抗扰控制系统 Active CN109828469B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910183791.9A CN109828469B (zh) 2019-03-12 2019-03-12 相位优化的扩张状态观测器及抗扰控制系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910183791.9A CN109828469B (zh) 2019-03-12 2019-03-12 相位优化的扩张状态观测器及抗扰控制系统

Publications (2)

Publication Number Publication Date
CN109828469A CN109828469A (zh) 2019-05-31
CN109828469B true CN109828469B (zh) 2021-11-02

Family

ID=66869159

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910183791.9A Active CN109828469B (zh) 2019-03-12 2019-03-12 相位优化的扩张状态观测器及抗扰控制系统

Country Status (1)

Country Link
CN (1) CN109828469B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110687800B (zh) * 2019-11-19 2022-05-13 大连海事大学 一种数据驱动自适应抗干扰控制器及其估计方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104142627A (zh) * 2014-07-03 2014-11-12 浙江工业大学 一种采用自抗扰控制技术的网络化无刷直流电机时延补偿和控制方法
CN106773672A (zh) * 2016-11-24 2017-05-31 滨州学院 提高观测精度的新型三阶线性扩张状态观测器构造方法
CN106873611A (zh) * 2017-01-22 2017-06-20 北京航空航天大学 一种多通道线性自抗扰控制器的设计方法
CN107272421A (zh) * 2017-08-10 2017-10-20 广东工业大学 一种自抗扰控制系统及方法
CN108919652A (zh) * 2018-10-10 2018-11-30 北京工商大学 一种自适应抗扰整形控制方法与系统
CN108983610A (zh) * 2018-08-03 2018-12-11 华侨大学 一种鲁棒自适应抗扰控制方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104142627A (zh) * 2014-07-03 2014-11-12 浙江工业大学 一种采用自抗扰控制技术的网络化无刷直流电机时延补偿和控制方法
CN106773672A (zh) * 2016-11-24 2017-05-31 滨州学院 提高观测精度的新型三阶线性扩张状态观测器构造方法
CN106873611A (zh) * 2017-01-22 2017-06-20 北京航空航天大学 一种多通道线性自抗扰控制器的设计方法
CN107272421A (zh) * 2017-08-10 2017-10-20 广东工业大学 一种自抗扰控制系统及方法
CN108983610A (zh) * 2018-08-03 2018-12-11 华侨大学 一种鲁棒自适应抗扰控制方法
CN108919652A (zh) * 2018-10-10 2018-11-30 北京工商大学 一种自适应抗扰整形控制方法与系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Compound Active Disturbance Rejection Control for Resonance Damping and Tracking of Nanopositioning;WEI Wei 等;《Proceedings of the 33rd Chinese Control Conference》;20140730;第5906-5909页 *
时滞系统的自抗扰控制综述;王丽君 等;《控制理论与应用》;20131231;第30卷(第12期);第1521-1533页 *

Also Published As

Publication number Publication date
CN109828469A (zh) 2019-05-31

Similar Documents

Publication Publication Date Title
CN107168071B (zh) 一种基于干扰观测器的非线性系统自抗扰控制方法
Liu et al. Non-asymptotic fractional order differentiator for a class of fractional order linear systems
WO2019165801A1 (zh) 一种扰动感知控制方法
CN107272421B (zh) 一种自抗扰控制系统及方法
CN109870899A (zh) 一种基于扩张状态观测器的光电跟踪系统控制方法
CN111752262B (zh) 一种执行器故障观测器与容错控制器集成设计方法
CN105043348A (zh) 基于卡尔曼滤波的加速度计陀螺仪水平角度测量方法
CN108919652B (zh) 一种自适应抗扰整形控制方法与系统
Sun et al. Disturbance attenuation via double-domination approach for feedforward nonlinear system with unknown output function
CN109828469B (zh) 相位优化的扩张状态观测器及抗扰控制系统
CN106383443B (zh) 抗干扰控制方法及系统
CN104199993A (zh) 一种针对间歇过程的二维动态卡尔曼滤波器设计方法
CN108829946A (zh) 一种基于动态补偿技术的推力计算方法
CN108107728B (zh) 一种基于干扰补偿的电液位置伺服系统控制方法
CN110442026B (zh) 基于误差修正的扩张状态观测器、抗扰控制系统及设计方法
CN110569561B (zh) 一种分数阶pid控制器微分积分阶次估计方法
CN109507870A (zh) 结构自适应的分数阶比例积分或比例微分控制器设计方法
CN109426143B (zh) 负载转矩估算方法、系统、机电控制系统、方法及电机
CN110824922B (zh) 基于六阶b-样条小波神经网络的史密斯预估补偿方法
CN111682819B (zh) 基于改进ladrc的异步电机电流内环解耦控制方法
CN110632847A (zh) 一种自耦pd协同控制理论新方法
CN113671833B (zh) 不确定重复控制系统协调设计方法及系统
CN107421543B (zh) 一种基于状态扩维的隐函数量测模型滤波方法
Guo et al. Active Disturbance Rejection Control of Valve‐Controlled Cylinder Servo Systems Based on MATLAB‐AMESim Cosimulation
CN115236990A (zh) 一种广义估计器、广义抗扰控制器及设计方法

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
GR01 Patent grant
GR01 Patent grant