CN110320795B - 一种采用自抗扰控制结构实现任意线性控制器的方法 - Google Patents
一种采用自抗扰控制结构实现任意线性控制器的方法 Download PDFInfo
- Publication number
- CN110320795B CN110320795B CN201910703107.5A CN201910703107A CN110320795B CN 110320795 B CN110320795 B CN 110320795B CN 201910703107 A CN201910703107 A CN 201910703107A CN 110320795 B CN110320795 B CN 110320795B
- Authority
- CN
- China
- Prior art keywords
- controller
- control
- linear
- ladrc
- state
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 20
- 230000009466 transformation Effects 0.000 claims abstract description 15
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000012938 design process Methods 0.000 claims description 6
- 238000013459 approach Methods 0.000 claims description 4
- 230000010354 integration Effects 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 238000011426 transformation method Methods 0.000 abstract description 2
- 238000013461 design Methods 0.000 description 5
- 230000036039 immunity Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B11/00—Automatic controllers
- G05B11/01—Automatic controllers electric
- G05B11/36—Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
- G05B11/42—Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential for obtaining a characteristic which is both proportional and time-dependent, e.g. P. I., P. I. D.
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
本发明提供一种采用自抗扰控制结构实现任意线性控制器的方法,所述方法在已知严格正则控制器K的前提下,通过引入线性变换,在状态空间下对控制器K与传统LADRC控制器进行描述,得到二者参数间的非对称Riccati方程,通过对非对称Riccati方程求解,获得与控制器K具有相同的抗扰及鲁棒性能的可调整的LADRC控制器,对可调整的LADRC控制器进行在线调整以获得理想控制性能。本发明在满足控制要求的情况下,采用自抗扰控制结构实现了任意线性控制器的控制,控制系统简单可靠;以线性系统理论为基础,在控制理论的状态空间描述下,通过线性变换的方法得出了相应的LADRC的实现,更为工程设计人员所接受。
Description
技术领域
本发明属于控制理论与应用技术领域,具体涉及一种采用自抗扰控制结构实现任意线性控制器的方法。
背景技术
自抗扰控制(Active Disturbance Rejection Control,简称ADRC)是中科院系统所韩京清研究员提出的一种通用的控制技术,它吸收并发扬了经典PID控制按误差进行调节的精髓,具有控制超调小、响应速度快、精度高、抗干扰能力强等优点而获得越来越广泛地关注。ADRC最本质的功能是扰动的估计和补偿能力,它主要由四个部分组成:跟踪微分器(tracking differentiator,TD)、扩张状态观测器(extended state observer,ESO)和非线性反馈(nonlinearstate error feedback,NLSEF),其核心思想是将系统内部的不确定性(定常或时变,线性或非线性)和外部不确定性(外部扰动)一起作为“总扰动”,通过构造扩张状态观测器对“总扰动”进行估计并实时补偿,从而获得较强的抗扰动能力。该方法不需要直接测量外扰,也不需要事先知道扰动作用规律,在了解了被控对象的相对阶次、输入输出通道个数、信号延迟时间等易获得且物理量清晰的特征量的基础上,就能够进行相关的控制器设计。ADRC技术是一种非线性控制方法,其结构还显得较为复杂,并且需要整定的参数较多,理论证明也较为困难,因此在实际应用中受到限制。
为解决上述问题,高志强教授等提出了线性自抗扰控制(linear ADRC,LADRC),利用所估计的误差及其各阶导数进行线性状态反馈,并且把ESO的设计和状态反馈的设计转化为观测器带宽和控制器带宽两个参数的选取,大大简化了ADRC的结构及参数整定,使得ADRC的实际应用成为可能。LADRC利用ESO对被控系统的内扰和外扰进行估计,然后采用状态反馈进行控制,不需要知道被控对象精确模型,确定对象相对阶后,控制结构固定,且只需整定2个或3个参数,是替代传统PID控制的潜在技术,其研究具有重大的理论和实际意义。目前对于自抗扰控制结构适用于哪些对象还未知,还缺乏一种采用自抗扰控制结构实现任意线性控制器的方法。
发明内容
为了解决上述问题,本发明提供一种采用自抗扰控制结构实现任意线性控制器的方法,所述方法在已知严格正则控制器K的前提下,通过引入线性变换,在状态空间下对控制器K与传统LADRC控制器进行描述,得到二者参数间的非对称Riccati方程,通过对非对称Riccati方程求解,获得与控制器K具有相同的抗扰及鲁棒性能的可调整的LADRC控制器,对可调整的LADRC控制器进行在线调整以获得理想控制性能;
进一步地,所述方法包括:
S1:给出传统LADRC控制器结构及设计过程;
S2:通过引入线性变换,在状态空间描述下得到已知严格正则控制器K与S1中的传统LADRC控制器的相互关系,给出二者参数间的非对称Riccati方程;
S3:通过对S2中非对称Riccati方程进行求解得到LADRC的各参数,获取与控制器K具有相同的抗扰及鲁棒性能的可调整的LADRC控制器;
进一步地,,所述S1具体包括:
假设对象模型为:
假设对象模型为:
传统LADRC控制器采用串联积分器模型,式(1)中,b0为对象模型高频增益的估计,y(t)为系统输出,u(t)为系统输入,f(t)为系统总扰动,d为外部干扰,t为时间;
线性自抗扰控制利用线性扩张状态观测器来估计被控对象的状态以及未知的广义扰动,设计线性扩张状态观测器如下:
对式(1)模型进行状态空间模型转换,得到:
且
Lo为线性扩张状态观测器增益:
Lo=[β1 β2…βp βp+1]T (6);
进一步地,基于假设对象模型所获得的被控对象的状态以及估计的广义扰动,采用如下线性状态反馈控制率控制:
而状态反馈增益Ko定义为:
Ko=[k1 k2…kp 1]/b0 (9);
线性自抗扰控制器具由如下状态空间实现:
通过以上设计过程可知,一个线性自抗扰控制需要设计如下三组参数:b0、ESO的观测器增益Lo以及p阶积分对象的控制增益Ko;其中b0取值应尽量逼近系统的实际高频增益b;通常b0越小,控制作用越强,但稳定裕度越小;
进一步地,所述S1中传统LADRC同时实现了对控制系统的跟踪与抗扰,其中与抗扰相关的控制器为
进一步地,所述S2具体包括:
对于被控对象,已知其严格正则控制器为K,状态空间实现为:
其中,Ak为p×p维,Bk为p×1维,Ck为1×p维,将线性自抗扰控制器(11),采用线性变换z=Pxk,则
则
由(14)可得
PBkCeP+PAk-AeP+BeCk=0 (15);
二者之间的线性变换阵P的求解对应为非对称代数Riccati方程的解;
进一步地,所述S3具体包括:
对于非对称代数Riccati方程的求解,通过构造相应的Hamiltonian矩阵来进行求解,其表达形式如下
此Hamiltonian矩阵为2p×2p维,求得其特征值为
选取其中的p个特征值对应的特征向量构成相应的变换阵为
其中V1为p×p维,V2为p×p维,则相应的非对称代数Riccati方程的解为P=V2V1 -1;
根据式(11)得到
最终获得到与控制器K相应的LADRC实现;
本发明的有益效果如下:
1):本发明在满足控制要求的情况下,采用自抗扰控制结构实现了任意线性控制器的控制,控制系统简单可靠;
2):本发明以线性系统理论为基础,在控制理论的状态空间描述下,通过线性变换的方法得出了相应的LADRC的实现,从而更为工程设计人员所接受;
3):采用本发明整定的LADRC控制器可以获得与已知给定控制器相当的控制性能,为很多工程上LADRC的整定提供了初始的调整参数,且便于根据对象实际运行状态及控制性能的要求来进行调整,模型的适应能力较强,可操作性好。
附图说明
图1为本发明所述二阶LADRC控制器结构图;
图2为本发明实施例1中所述LADRC与给定控制器K的频率响应曲线图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细描述。应当理解,此处所描述的具体实施例仅仅用于解释本发明,并不用于限定本发明。相反,本发明涵盖任何由权利要求定义的在本发明的精髓和范围上做的替代、修改、等效方法以及方案。进一步,为了使公众对本发明有更好的了解,在下文对本发明的细节描述中,详尽描述了一些特定的细节部分。对本领域技术人员来说没有这些细节部分的描述也可以完全理解本发明。
下面结合附图和具体实施例对本发明作进一步说明,但不作为对本发明的限定。下面为本发明的举出最佳实施例:
如图1-图2所示,本发明提供一种采用自抗扰控制结构实现任意线性控制器的方法,其中图1为二阶LADRC控制器结构图,同时也是本发明中进行参数整定的控制器结构图。
所述方法采用传统LADRC控制结构,根据已知控制器,采用状态空间下线性变换来得到LADRC的各参数。
所述线性自抗扰控制器的设计步骤为:
S1:给出传统LADRC控制器结构及设计过程;
S2:通过引入线性变换,在状态空间描述下得到已知严格正则控制器与S1中的LADRC控制器的相互关系,给出二者参数间的Riccati方程;
S3:通过对S2的非对称Riccati方程进行求解得到LADRC的各参数。
所述步骤S1具体包括:
假设对象模型为:
传统线性自抗扰控制采用串联积分器模型,式中,b0为对象模型高频增益的估计,y(t)为系统输出,u(t)为系统输入,f(t)为系统总扰动。
线性自抗扰控制利用线性扩张状态观测器来估计被控对象的状态以及未知的广义扰动。
其中
且
为所选取状态且Lo为观测器增益:
Lo=[β1 β2…βp βp+1]T (6);
在广义扰动f(t)有界的假设下,当Ae-LoCe渐进稳定时,z1(t),z2(t),…zp(t)趋近于输出并且zp+1(t)趋近于f(t)。因此可以利用zp+1(t)进行控制,使f(t)得到更快地抑制。
基于以上获得的被控对象的状态以及估计的广义扰动,采用如下线性状态反馈控制率控制
而状态反馈增益Ko定义为
Ko=[k1 k2…kp 1]/b0 (9);
综合上述,线性自抗扰控制器具有如下状态空间实现
通过以上设计过程可知,一个线性自抗扰控制需要设计以下三组参数:b0、ESO的观测器增益Lo以及p阶积分对象的控制增益Ko。其中b0取值应尽量逼近系统的实际高频增益b。通常b0越小,控制作用越强,但稳定裕度越小。
可以看到,LADRC同时实现了对控制系统的跟踪与抗扰。在工业过程控制中,我们更关心的是抗扰性能,而LADRC中与抗扰相关的控制器为
所述步骤S2具体包括:
对于对象(1),已知其严格正则控制器为K,状态空间实现为:
其中,Ak为p×p维,Bk为p×1维,Ck为1×p维。由线性自抗扰控制器(11),采用线性变换z=Pxk,则
则
由(13)可得
PBkCeP+PAk-AeP+BeCk=0 (14)。
由此可以看到,二者之间的线性变换阵P的求解为非对称代数Riccati方程的解。
所述步骤S3具体包括:
对于非对称代数Riccati方程的求解,可以通过构造相应的
Hamiltonian矩阵来进行求解,其表达形式如下
此Hamiltonian矩阵为2p×2p维,可以求得其特征值为
选取其中的p个特征值对应的特征向量构成相应的变换阵为
其中V1为p×p维,V2为p×p维,则相应的非对称代数Riccati方程的解为P=V2V1 -1。
根据式(11)可以得到
由此得到了与该控制器K相应的LADRC的实现。
综上所述,针对一般工业过程控制对象,任何已知严格正则控制器K,利用该发明所提供的方法都可以得到其相应的LADRC实现。经过以上步骤得到的LADRC可以获得与K相同的抗扰及鲁棒性能,利用得到的LADRC控制器参数进行控制,在需要时可以进行必要的在线调整以获得理想的控制性能。
实施例1:
给定控制器为:
通过S2和S3所述步骤,取b0=1得到LADRC控制器参数为:
K0=[-0.0711 -0.4564 0.6226 0.0385],
Lo=[1.7129 3.8064 2.9587 2.0458].
通过图2,给出了LADRC与给定控制器K的频率响应曲线,可以看到二者完全重合。因此通过图1所示结构对实际被控对象进行控制,可以得到与K相同抗扰性能。
以上所述的实施例,只是本发明较优选的具体实施方式的一种,本领域的技术人员在本发明技术方案范围内进行的通常变化和替换都应包含在本发明的保护范围内。
Claims (1)
1.一种采用自抗扰控制结构实现任意线性控制器的方法,其特征在于,所述方法针对一般工业过程控制对象,在任何已知严格正则控制器K的前提下,通过引入线性变换,在状态空间下对控制器K与传统LADRC控制器进行描述,得到二者参数间的非对称Riccati方程,通过对非对称Riccati方程求解,获得与控制器K具有相同的抗扰及鲁棒性能的可调整的LADRC控制器,对可调整的LADRC控制器进行在线调整以获得理想控制性能及相应的LADRC实现,所述方法包括:
S1:给出传统LADRC控制器结构及设计过程;
S2:通过引入线性变换,在状态空间描述下得到已知严格正则控制器K与S1中的传统LADRC控制器的相互关系,给出二者参数间的非对称Riccati方程;
S3:通过对S2中非对称Riccati方程进行求解得到LADRC的各参数,获取与控制器K具有相同的抗扰及鲁棒性能的可调整的LADRC控制器,所述S1具体包括:
假设对象模型为:
传统LADRC控制器采用串联积分器模型,式(1)中,b0为对象模型高频增益的估计,y(t)为系统输出,u(t)为系统输入,f(t)为系统总扰动,d为外部干扰,t为时间;
线性自抗扰控制利用线性扩张状态观测器来估计被控对象的状态以及未知的广义扰动,设计线性扩张状态观测器如下:
对式(1)模型进行状态空间模型转换,得到:
且
Lo为线性扩张状态观测器增益:
Lo=[β1 β2 … βp βp+1]T (6);
在广义扰动f(t)有界的假设下,当Ae-LoCe渐进稳定时,z1(t),z2(t),…zp(t)趋近于输出并且zp+1(t)趋近于f(t);因此通过zp+1(t)进行控制,使f(t)得到更快地抑制,基于假设对象模型所获得的被控对象的状态以及估计的广义扰动,采用如下线性状态反馈控制率控制:
状态反馈增益Ko定义为:
Ko=[k1 k2 … kp 1]/b0 (9);
线性自抗扰控制器具由如下状态空间实现:
通过以上设计过程可知,一个线性自抗扰控制需要设计如下三组参数:b0、ESO的观测器增益Lo以及p阶积分对象的控制增益Ko;其中b0取值应尽量逼近系统的实际高频增益b;通常b0越小,控制作用越强,但稳定裕度越小,所述S1中传统LADRC同时实现了对控制系统的跟踪与抗扰,其中与抗扰相关的控制器为:
所述S2具体包括:
对于被控对象,已知其严格正则控制器K,状态空间实现为:
其中,Ak为p×p维,Bk为p×1维,Ck为1×p维,将线性自抗扰控制器(11)采用线性变换z=Pxk,则
则
由(14)可得
PBkCeP+PAk-AeP+BeCk=0 (15);
二者之间的线性变换阵P的求解对应为非对称代数Riccati方程的解,所述S3具体包括:
对于非对称代数Riccati方程的求解,通过构造相应的Hamiltonian矩阵来进行求解,其表达形式如下
此Hamiltonian矩阵为2p×2p维,求得其特征值为
选取其中的p个特征值对应的特征向量构成相应的变换阵为
根据式(11)得到
最终获得到与控制器K相应的LADRC实现。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910703107.5A CN110320795B (zh) | 2019-07-31 | 2019-07-31 | 一种采用自抗扰控制结构实现任意线性控制器的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910703107.5A CN110320795B (zh) | 2019-07-31 | 2019-07-31 | 一种采用自抗扰控制结构实现任意线性控制器的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110320795A CN110320795A (zh) | 2019-10-11 |
CN110320795B true CN110320795B (zh) | 2021-06-01 |
Family
ID=68125044
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910703107.5A Expired - Fee Related CN110320795B (zh) | 2019-07-31 | 2019-07-31 | 一种采用自抗扰控制结构实现任意线性控制器的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110320795B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110703606A (zh) * | 2019-10-30 | 2020-01-17 | 曾喆昭 | 一种自耦pid协同控制理论新方法 |
CN111352345B (zh) * | 2020-04-08 | 2022-05-17 | 河南工业大学 | 一种基于极点配置的线性自抗扰控制参数整定方法 |
CN113110252B (zh) * | 2021-05-19 | 2022-03-18 | 广东电网有限责任公司 | 一种双参数新型基础控制器的控制方法 |
CN117895460B (zh) * | 2024-03-14 | 2024-07-16 | 国网四川省电力公司电力科学研究院 | 微电网储能变流器线性自抗扰控制参数整定方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104914725A (zh) * | 2015-06-12 | 2015-09-16 | 中国人民解放军63680部队 | 单入单出自抗扰控制系统鲁棒稳定性分析方法 |
CN106773648A (zh) * | 2016-12-19 | 2017-05-31 | 华侨大学 | 一种自抗扰控制的鲁棒保性能设计与参数整定方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8041436B2 (en) * | 2002-04-18 | 2011-10-18 | Cleveland State University | Scaling and parameterizing a controller |
-
2019
- 2019-07-31 CN CN201910703107.5A patent/CN110320795B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104914725A (zh) * | 2015-06-12 | 2015-09-16 | 中国人民解放军63680部队 | 单入单出自抗扰控制系统鲁棒稳定性分析方法 |
CN106773648A (zh) * | 2016-12-19 | 2017-05-31 | 华侨大学 | 一种自抗扰控制的鲁棒保性能设计与参数整定方法 |
Non-Patent Citations (2)
Title |
---|
Application and the parameter tuning of ADRC based on CPSO;LIU C F, ZANG B;《Proceedings of the 24th Chinese Control and Decision Conference》;20121231;第3277-3281页 * |
基于高阶控制器设计的线性自抗扰控制参数调整;傅彩芬,谭文;《控制理论与应用》;20170228;第34卷(第2期);第265-272页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110320795A (zh) | 2019-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110320795B (zh) | 一种采用自抗扰控制结构实现任意线性控制器的方法 | |
CN106802660B (zh) | 一种复合强抗扰姿态控制方法 | |
CN108303885B (zh) | 一种基于干扰观测器的电机位置伺服系统自适应控制方法 | |
CN108416126B (zh) | 一种dc/dc变换器稳定性分析方法及系统 | |
Hu et al. | Performance-oriented adaptive robust control of a class of nonlinear systems preceded by unknown dead zone with comparative experimental results | |
CN109459934A (zh) | 一种基于pid控制器整定降阶自抗扰控制器参数的方法 | |
CN107121932B (zh) | 电机伺服系统误差符号积分鲁棒自适应控制方法 | |
CN104242769A (zh) | 基于连续终端滑模技术的永磁同步电机速度复合控制方法 | |
CN104199294B (zh) | 电机伺服系统双神经网络摩擦补偿和有限时间协同控制方法 | |
Zhao et al. | Reinforcement learning based optimal control of linear singularly perturbed systems | |
Ran et al. | Self-tuning weighted measurement fusion Kalman filtering algorithm | |
CN108919652B (zh) | 一种自适应抗扰整形控制方法与系统 | |
CN109856970B (zh) | 一种具有网络诱导有界时滞和数据丢失的有限时镇定方法 | |
CN112769364B (zh) | 一种直流电机伺服系统的快速自适应抗扰控制方法 | |
CN108803325B (zh) | 永磁同步电机伺服系统鲁棒有限时间控制方法 | |
CN109240085B (zh) | 非高斯系统动态数据校正与系统控制性能优化方法 | |
CN105549385B (zh) | 多变量时滞非最小相位非方系统的解耦内模控制器、控制系统和控制方法 | |
Gao et al. | Neural network based dynamic surface integral nonsingular fast terminal sliding mode control for manipulators with disturbance rejection | |
CN111880470B (zh) | 压电驱动微定位平台的无抖振滑模控制方法 | |
Prakash et al. | ROBUST MODEL REFERENCE ADAPTIVE PI CONTROL. | |
Fridman et al. | High-order sliding-mode observation and identification for linear systems with unknown inputs | |
CN110095985B (zh) | 一种观测器设计方法和抗干扰控制系统 | |
CN116442223A (zh) | 一种机械手系统轨迹跟踪的非线性动态控制器设计方法 | |
CN113641193B (zh) | 一种非最小相位系统精确跟踪控制方法 | |
CN110209055B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210601 |
|
CF01 | Termination of patent right due to non-payment of annual fee |