CN103454921A - 飞行控制系统非线性跟踪控制器设计的正切线性化方法 - Google Patents
飞行控制系统非线性跟踪控制器设计的正切线性化方法 Download PDFInfo
- Publication number
- CN103454921A CN103454921A CN2013103983541A CN201310398354A CN103454921A CN 103454921 A CN103454921 A CN 103454921A CN 2013103983541 A CN2013103983541 A CN 2013103983541A CN 201310398354 A CN201310398354 A CN 201310398354A CN 103454921 A CN103454921 A CN 103454921A
- Authority
- CN
- China
- Prior art keywords
- control system
- state
- control
- nonlinear
- tangent
- 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.)
- Granted
Links
Landscapes
- Feedback Control In General (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明涉及一种飞行器非线性控制跟踪控制系统的正切线性化方法。飞行器非线性跟踪控制系统的正切线性化方法,针对飞行控制系统的非线性跟踪控制器类,其特征在于:包括:对飞行器非线性控制系统进行状态扩展;计算正切线性化系统的状态微分方程;设计时变线性系统反馈镇定增益;求解微分方程得到非线性飞行控制系统跟踪控制律四大步骤。同现有技术相比,提供了基于正切线性化的飞行器控制系统跟踪控制器设计方法。可以对复杂飞行控制系统进行完全等价的线性化建模以及对控制系统进行跟踪控制器设计,为控制系统实时实现提供一种有效分析工具,进一步为满足飞行控制系统设计方面的高要求提供有力的理论依据和技术支撑,具有很好的工程应用价值。
Description
技术领域
本发明属于自动控制技术领域,涉及一种飞行器非线性控制跟踪控制系统的正切线性化方法。
背景技术
现代飞行器的飞行高度和飞行马赫数跨度范围大、飞行环境复杂,在飞行过程中飞行器的气动特性和气热特性剧烈变化,很难保持飞行器在平衡状态条件下稳定飞行,从而导致飞行器偏离预期飞行轨迹。因此,有必要寻求基于非平衡状态的线性化和轨迹跟踪控制系统的方法。
在传统飞行控制建模中,建模方法采用的是传统的雅可比线性化方法。然而,可使用雅可比线性化方法的前提条件是系统状态是在基于平衡状态附近做小扰动。但是,在某些实际被控对象中,系统工作状态可能是偏离平衡状态且无法满足小扰动条件。
正切线性化控制方法主要思想是将一些复杂的非线性控制问题转化为一类线性时变系统的控制问题,而这类线性时变系统的控制问题可采用现有方法有效解决。这可能主要是由于在应用该方法要求系统的状态变量不能变化太快的限定条件,例如要求原始系统的状态变化率必须满足随时间指数衰减等。然而,在实际应用过程中,某些非线性系统可看作慢时变的,例如在巡航状态下做机动飞行条件下的飞行器模型。因此,这种正切线性化控制方法为飞行器非线性系统的稳定性分析和控制器设计提供了一种较有应用前景的途径。但是到目前为止,在现有文献中利用正切线性化方法进行控制器设计的结果很少。
发明内容
针对上述现有技术状况,本发明的目的是:针对飞行器控制系统的非线性模型,提供一种基于正切线性化控制方法,给出姿态跟踪控制器设计的方法步骤,为飞行控制系统设计提供一种新的技术途径。该方法不仅适用于平衡点附近的飞行器控制设计,对于非平衡点处的控制器设计仍然有效。此外,该技术无需大量繁琐的数值计算和状态变量的非线性转化,并且被控系统的形式不受限制。
现将本发明构思及技术解决方案叙述如下:
本发明一种飞行器非线性跟踪控制系统的正切线性化方法,针对飞行控制系统的非线性跟踪控制器类,其特征在于:包括以下步骤:
步骤1:对飞行器非线性控制系统进行状态扩展
步骤1.1:考虑如下一类具有向量可微分的飞行器非线性控制系统:
其中,x(t)∈Rn和u(t)∈Rp分别表示系统的状态向量和输入向量,f为满足如下条件的连续可微函数
步骤1.2:对于非线性控制系统(1),考虑如下跟踪控制问题:给定初始条件x(0)=x0,保证系统状态x(t)随着t→∞渐近收敛于目标状态x(∞);
步骤1.3:对任意t≥0,如果存在一个控制律u(t)及实数α,β>0,使得系统(1)的对应解x(t)满足如下条件
存在一个x(∞)使得下式成立
步骤1.4:一般情况下,非线性系统(1)的终值状态目标x(∞)并不能随意给定,但对于飞行控制系统而言,如果被控飞行系统的一个或多个状态需要精确地被跟踪达到期望目标,则可以应用下述策略进行实现;假定xi(t)为系统状态的某个待指定的状态且其期望目标为xi(∞),可以增加一个新的状态xn+1(t)使得其中,如果是对多个状态实现跟踪控制,可参照上述方法依次增加新的状态量。那么,仅需针对如下状态扩展系统进行控制器设计
步骤2:、计算正切线性化系统的状态微分方程
其中,
步骤2.2:为求解上述期望非线性控制律u(t),引入下述“正切线性化系统”,其微分方程如下
步骤2.3:于是,可以利用现有的线性时变控制方法来保证正切线性化系统(8)的稳定性,针对正切线性化系统(8),设计如下的状态反馈控制律
ξ(t)=K(t)ζ(t) (9)
综合考虑系统(8)和反馈控制律(9),可以得到如下闭环系统
步骤2.4:针对对于正切线性化系统(8),在任意初始状态ζ(0)=ζ0下,如果存在实数α,β>0使得扩展系统(5)中的非线性控制律u(t)满足如下条件
步骤3:设计时变线性系统反馈镇定增益
为求取如式(9)所示的状态反馈控制律,可针对实际控制系统情形,利用现有时变线性系统控制理论与方法进行求解;这里仅针对飞行器在巡航机动飞行过程中控制系统的特点,参考有关慢时变线性系统控制结论,可以给出如式(9)中状态反馈控制增益如下
K(t)=-R-1BT(t)P(t) (13)
其中,P(t)为下述时间冻结代数Riccati方程的解
(AT(t)+σI)P(t)+P(t)(A(t)+σI)-P(t)B(t)R-1BT(t)P(t)+Q=0 (14)
其中,Q,R均为合适维数的正定对称矩阵,σ>0是使得闭环化系统(10)的解ζ(t)满足式(12)中期望性能的可调自由参数。
步骤4:求解微分方程得到非线性飞行控制系统跟踪控制律
将式(13)中求得的反馈控制增益代入到式(11)中,通过给定的初始条件,联立系统状态方程(5),实时求解出非线性控制律u(t),然后将控制律u(t)与原始的飞行器控制系统形成闭环系统,实现了飞行控制系统的正切线性化跟踪控制设计。
本发明同现有技术相比,提供了基于正切线性化的飞行器控制系统跟踪控制器设计方法。不仅可以对复杂飞行控制系统进行完全等价的线性化建模,还可以对控制系统进行跟踪控制器设计,为控制系统实时实现提供一种有效分析工具,进一步为满足飞行控制系统设计方面的高要求提供有力的理论依据和技术支撑,具有很好的工程应用价值
具体实施方式
实施例
基于正切线性化方法的一类飞行控制系统非线性跟踪控制器设计方法,包括以下步骤:
步骤1:根据式(1)、(2)、(3)、(4)、(5)对飞行器非线性控制系统进行状态扩展
步骤2:根据式(6)、(7)、(8、)(9)、(10)、(11)、(12)计算正切线性化系统的状态微分方程;
步骤3:根据式(13)、(14)、(15)设计时变线性系统反馈镇定增益;
步骤4:将式(13)中求得的反馈控制增益代入到式(11)中,通过给定的初始条件,联立系统状态方程(5),实时求解出非线性控制律u(t),然后将控制律u(t)与原始的飞行器控制系统形成闭环系统,求解微分方程得到非线性飞行控制系统跟踪控制律;实现飞行控制系统的正切线性化跟踪控制设计。
Claims (5)
1.飞行器非线性跟踪控制系统的正切线性化方法,针对飞行控制系统的非线性跟踪控制器类,其特征在于:包括以下步骤:
步骤1:对飞行器非线性控制系统进行状态扩展;
步骤2:、计算正切线性化系统的状态微分方程;
步骤3:设计时变线性系统反馈镇定增益;
步骤4:求解微分方程得到非线性飞行控制系统跟踪控制律。
2.根据权利要求1所述的飞行器非线性跟踪控制系统的正切线性化方法,其特征在于:步骤1中所述的“对飞行器非线性控制系统进行状态扩展”的具体步骤如下:
步骤1.1:考虑如下一类具有向量可微分的飞行器非线性控制系统:
其中,x(t)∈Rn和u(t)∈Rp分别表示系统的状态向量和输入向量,f为满足如下条件的连续可微函数
步骤1.2:对于非线性控制系统(1),考虑如下跟踪控制问题:给定初始条件x(0)=x0,保证系统状态x(t)随着t→∞渐近收敛于目标状态x(∞);
步骤1.3:对任意t≥0,如果存在一个控制律u(t)及实数α,β>0,使得系统(1)的对应解x(t)满足如下条件
存在一个x(∞)使得下式成立
步骤1.4:非线性系统(1)的终值状态目标x(∞)不能随意给定;对于飞行控制系统被控飞行系统的一个或多个状态需要精确地被跟踪达到期望目标,则可以应用下述策略进行实现;假定xi(t)为系统状态的某个待指定的状态且其期望目标为xi(∞),可以增加一个新的状态xn+1(t)使得其中,如果是对多个状态实现跟踪控制,可参照上述方法依次增加新的状态量;那么,仅需针对如下状态扩展系统进行控制器设计
3.根据权利要求1所述的飞行器非线性跟踪控制系统的正切线性化方法,其特征在于:步骤2中所述的“计算正切线性化系统的状态微分方程”的具体步骤如下:
其中,
步骤2.2:为求解上述期望非线性控制律u(t),引入下述“正切线性化系统”,其微分方程如下
步骤2.3:于是,可以利用现有的线性时变控制方法来保证正切线性化系统(8)的稳定性,针对正切线性化系统(8),设计如下的状态反馈控制律
ξ(t)=K(t)ζ(t) (9)
综合考虑系统(8)和反馈控制律(9),可以得到如下闭环系统
步骤2.4:针对对于正切线性化系统(8),在任意初始状态ζ(0)=ζ0下,如果存在实数σ,β>0使得扩展系统式(5)中的非线性控制律u(t)满足如下条件
也即是,随时间指数衰减。
4.根据权利要求1所述的飞行器非线性跟踪控制系统的正切线性化方法,其特征在于:步骤3中所述的“设计时变线性系统反馈镇定增益”具体如下:
为求取如式(9)所示的状态反馈控制律,针对实际控制系统情形,利用现有时变线性系统控制理论与方法进行求解;这里仅针对飞行器在巡航机动飞行过程中控制系统的特点,参考有关慢时变线性系统控制结论,可以给出如式(9)中状态反馈控制增益如下
K(T)=-R-1BT(t)P(t) (13)
其中,P(t)为下述时间冻结代数Riccati方程的解
(AT(t)+σI)P(t)+P(t)(A(t)+σI)-P(t)B(t)R-1BT(t)P(t)+Q=0 (14)
其中,Q,R均为合适维数的正定对称矩阵,σ>0是使得闭环化系统(10)的解ζ(t)满足式(12)中期望性能的可调自由参数。
5.根据权利要求1所述的飞行器非线性跟踪控制系统的正切线性化方法,其特征在于:步骤4中所述的“求解微分方程得到非线性飞行控制系统跟踪控制律”具体如下:
将式(13)中求得的反馈控制增益代入到式(11)中,通过给定的初始条件,联立系统状态方程(5),实时求解出非线性控制律u(t),然后将控制律u(t)与原始的飞行器控制系统形成闭环系统,实现了飞行控制系统的正切线性化跟踪控制设计。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310398354.1A CN103454921B (zh) | 2013-08-30 | 2013-08-30 | 飞行控制系统非线性跟踪控制器设计的正切线性化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310398354.1A CN103454921B (zh) | 2013-08-30 | 2013-08-30 | 飞行控制系统非线性跟踪控制器设计的正切线性化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103454921A true CN103454921A (zh) | 2013-12-18 |
CN103454921B CN103454921B (zh) | 2017-04-12 |
Family
ID=49737417
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310398354.1A Active CN103454921B (zh) | 2013-08-30 | 2013-08-30 | 飞行控制系统非线性跟踪控制器设计的正切线性化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103454921B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104503239A (zh) * | 2014-12-17 | 2015-04-08 | 苏州大学 | 一种Lur’e型微分包含系统的控制器设计方法 |
CN105388756A (zh) * | 2014-08-27 | 2016-03-09 | 施耐德电气建筑有限公司 | 用于控制输入到建筑物内的能量的系统和方法 |
CN105511261A (zh) * | 2015-10-29 | 2016-04-20 | 重庆科技学院 | 基于非线性三自由度纵列式双旋翼直升机模型系统 |
CN107491080A (zh) * | 2017-07-12 | 2017-12-19 | 西北工业大学 | 一种基于非线性反馈和微分跟踪的飞行器姿态控制方法 |
CN110104218A (zh) * | 2019-04-19 | 2019-08-09 | 北京航天自动控制研究所 | 火箭发动机机架变形角的预偏角非线性补偿方法和装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030137332A1 (en) * | 2002-01-10 | 2003-07-24 | Masato Tanaka | Feedback control method and feedback control apparatus |
CN1612077A (zh) * | 2003-10-30 | 2005-05-04 | 株式会社东芝 | 参考模型跟踪控制系统和方法 |
CN1742241A (zh) * | 2002-12-09 | 2006-03-01 | 搭篷技术公司 | 具有可变动态特性的过程的自适应控制的一种系统和方法 |
CN101067743A (zh) * | 2007-06-13 | 2007-11-07 | 袁璞 | 变结构非线性模型预估控制器 |
CN101763037A (zh) * | 2009-12-17 | 2010-06-30 | 浙江大学 | 一种空分节能过程的非线性预测控制系统及方法 |
-
2013
- 2013-08-30 CN CN201310398354.1A patent/CN103454921B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030137332A1 (en) * | 2002-01-10 | 2003-07-24 | Masato Tanaka | Feedback control method and feedback control apparatus |
CN1742241A (zh) * | 2002-12-09 | 2006-03-01 | 搭篷技术公司 | 具有可变动态特性的过程的自适应控制的一种系统和方法 |
CN1612077A (zh) * | 2003-10-30 | 2005-05-04 | 株式会社东芝 | 参考模型跟踪控制系统和方法 |
CN101067743A (zh) * | 2007-06-13 | 2007-11-07 | 袁璞 | 变结构非线性模型预估控制器 |
CN101763037A (zh) * | 2009-12-17 | 2010-06-30 | 浙江大学 | 一种空分节能过程的非线性预测控制系统及方法 |
Non-Patent Citations (2)
Title |
---|
胡兆玮等: "使用控制模型的状态扩展模型研究", 《吉林师范大学学报自然科学版》, vol. 34, no. 3, 10 August 2013 (2013-08-10), pages 133 - 136 * |
蔡光斌等: "非线性系统控制器综合的一种新方法", 《中国科学技术大学学报》, vol. 42, no. 5, 31 May 2012 (2012-05-31), pages 383 - 390 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105388756A (zh) * | 2014-08-27 | 2016-03-09 | 施耐德电气建筑有限公司 | 用于控制输入到建筑物内的能量的系统和方法 |
CN105388756B (zh) * | 2014-08-27 | 2018-07-06 | 施耐德电气建筑有限公司 | 用于控制输入到建筑物内的能量的系统和方法 |
CN104503239A (zh) * | 2014-12-17 | 2015-04-08 | 苏州大学 | 一种Lur’e型微分包含系统的控制器设计方法 |
CN105511261A (zh) * | 2015-10-29 | 2016-04-20 | 重庆科技学院 | 基于非线性三自由度纵列式双旋翼直升机模型系统 |
CN107491080A (zh) * | 2017-07-12 | 2017-12-19 | 西北工业大学 | 一种基于非线性反馈和微分跟踪的飞行器姿态控制方法 |
CN107491080B (zh) * | 2017-07-12 | 2020-04-03 | 西北工业大学 | 一种基于非线性反馈和微分跟踪的飞行器姿态控制方法 |
CN110104218A (zh) * | 2019-04-19 | 2019-08-09 | 北京航天自动控制研究所 | 火箭发动机机架变形角的预偏角非线性补偿方法和装置 |
CN110104218B (zh) * | 2019-04-19 | 2021-04-30 | 北京航天自动控制研究所 | 火箭发动机机架变形角的预偏角非线性补偿方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN103454921B (zh) | 2017-04-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104698846B (zh) | 一种机械臂伺服系统的指定性能反演控制方法 | |
CN106406102B (zh) | 一种含干扰观测器的高超声速飞行器跟踪控制方法 | |
CN103454921A (zh) | 飞行控制系统非线性跟踪控制器设计的正切线性化方法 | |
CN105093931A (zh) | 一种航空发动机非线性系统控制器设计方法 | |
CN103558857A (zh) | 一种btt飞行器的分布式复合抗干扰姿态控制方法 | |
CN105182984A (zh) | 飞行器俯仰姿态的线性自抗扰控制器设计与参数整定 | |
Durmaz et al. | Sliding mode control for non-linear systems with adaptive sliding surfaces | |
CN107037727A (zh) | 一种无人直升机大包线自适应增益调度方法 | |
CN103248045A (zh) | 不确定电力系统等速滑模负荷频率控制器设计方法 | |
CN103809446B (zh) | 飞行器多回路模型簇颤振抑制复合频率鲁棒控制器设计方法 | |
CN103809434B (zh) | 纵向飞行模型簇复合根轨迹多级pid控制器设计方法 | |
CN104155876A (zh) | 一种pid控制器的分离实现方法 | |
CN104806302B (zh) | 一种基于非线性干扰观测器的汽轮发电机主汽门开度预测控制方法 | |
CN113467255A (zh) | 可重复使用运载器自适应多变量固定时间预设控制方法 | |
CN110032066B (zh) | 分数阶非线性系统轨迹跟踪的自适应迭代学习控制方法 | |
CN103809433B (zh) | 飞行器多回路模型簇复合根轨迹多级pid鲁棒控制器设计方法 | |
CN103777523B (zh) | 飞行器多回路模型簇复合pid鲁棒控制器设计方法 | |
CN103823376B (zh) | 纵向飞行模型簇复合pid控制器设计方法 | |
CN103809442A (zh) | 飞行器多回路模型簇复合频率鲁棒控制器设计方法 | |
CN104197792B (zh) | 一类多平衡点饱和切换系统的离散增益调度控制器设计方法及实现btt导弹控制的方法 | |
CN104298111A (zh) | 一种用于窑炉热工参数控制系统的模糊控制方法 | |
CN103809449B (zh) | 飞行器多回路模型簇颤振抑制复合pid鲁棒控制器设计方法 | |
CN103823367B (zh) | 纵向飞行模型簇颤振抑制复合频率鲁棒控制器设计方法 | |
CN103823375B (zh) | 纵向飞行模型簇复合根轨迹多级pid鲁棒控制器设计方法 | |
CN103809445B (zh) | 飞行器多回路模型簇复合pid控制器设计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |