CN103454921A - 飞行控制系统非线性跟踪控制器设计的正切线性化方法 - Google Patents

飞行控制系统非线性跟踪控制器设计的正切线性化方法 Download PDF

Info

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
Application number
CN2013103983541A
Other languages
English (en)
Other versions
CN103454921B (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.)
No 2 Artillery Engineering University Of Chinese Pla
Original Assignee
No 2 Artillery Engineering University Of Chinese Pla
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 No 2 Artillery Engineering University Of Chinese Pla filed Critical No 2 Artillery Engineering University Of Chinese Pla
Priority to CN201310398354.1A priority Critical patent/CN103454921B/zh
Publication of CN103454921A publication Critical patent/CN103454921A/zh
Application granted granted Critical
Publication of CN103454921B publication Critical patent/CN103454921B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Feedback Control In General (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明涉及一种飞行器非线性控制跟踪控制系统的正切线性化方法。飞行器非线性跟踪控制系统的正切线性化方法,针对飞行控制系统的非线性跟踪控制器类,其特征在于:包括:对飞行器非线性控制系统进行状态扩展;计算正切线性化系统的状态微分方程;设计时变线性系统反馈镇定增益;求解微分方程得到非线性飞行控制系统跟踪控制律四大步骤。同现有技术相比,提供了基于正切线性化的飞行器控制系统跟踪控制器设计方法。可以对复杂飞行控制系统进行完全等价的线性化建模以及对控制系统进行跟踪控制器设计,为控制系统实时实现提供一种有效分析工具,进一步为满足飞行控制系统设计方面的高要求提供有力的理论依据和技术支撑,具有很好的工程应用价值。

Description

飞行控制系统非线性跟踪控制器设计的正切线性化方法
技术领域
本发明属于自动控制技术领域,涉及一种飞行器非线性控制跟踪控制系统的正切线性化方法。
背景技术
现代飞行器的飞行高度和飞行马赫数跨度范围大、飞行环境复杂,在飞行过程中飞行器的气动特性和气热特性剧烈变化,很难保持飞行器在平衡状态条件下稳定飞行,从而导致飞行器偏离预期飞行轨迹。因此,有必要寻求基于非平衡状态的线性化和轨迹跟踪控制系统的方法。
在传统飞行控制建模中,建模方法采用的是传统的雅可比线性化方法。然而,可使用雅可比线性化方法的前提条件是系统状态是在基于平衡状态附近做小扰动。但是,在某些实际被控对象中,系统工作状态可能是偏离平衡状态且无法满足小扰动条件。
正切线性化控制方法主要思想是将一些复杂的非线性控制问题转化为一类线性时变系统的控制问题,而这类线性时变系统的控制问题可采用现有方法有效解决。这可能主要是由于在应用该方法要求系统的状态变量不能变化太快的限定条件,例如要求原始系统的状态变化率必须满足随时间指数衰减等。然而,在实际应用过程中,某些非线性系统可看作慢时变的,例如在巡航状态下做机动飞行条件下的飞行器模型。因此,这种正切线性化控制方法为飞行器非线性系统的稳定性分析和控制器设计提供了一种较有应用前景的途径。但是到目前为止,在现有文献中利用正切线性化方法进行控制器设计的结果很少。
发明内容
针对上述现有技术状况,本发明的目的是:针对飞行器控制系统的非线性模型,提供一种基于正切线性化控制方法,给出姿态跟踪控制器设计的方法步骤,为飞行控制系统设计提供一种新的技术途径。该方法不仅适用于平衡点附近的飞行器控制设计,对于非平衡点处的控制器设计仍然有效。此外,该技术无需大量繁琐的数值计算和状态变量的非线性转化,并且被控系统的形式不受限制。
现将本发明构思及技术解决方案叙述如下:
本发明一种飞行器非线性跟踪控制系统的正切线性化方法,针对飞行控制系统的非线性跟踪控制器类,其特征在于:包括以下步骤:
步骤1:对飞行器非线性控制系统进行状态扩展
步骤1.1:考虑如下一类具有向量可微分的飞行器非线性控制系统:
x · = f ( x , u ) - - - ( 1 )
其中,x(t)∈Rn和u(t)∈Rp分别表示系统的状态向量和输入向量,f为满足如下条件的连续可微函数
Figure BSA0000094705490000022
步骤1.2:对于非线性控制系统(1),考虑如下跟踪控制问题:给定初始条件x(0)=x0,保证系统状态x(t)随着t→∞渐近收敛于目标状态x(∞);
步骤1.3:对任意t≥0,如果存在一个控制律u(t)及实数α,β>0,使得系统(1)的对应解x(t)满足如下条件
| | x · ( t ) | | ≤ βe - αt - - - ( 3 )
存在一个x(∞)使得下式成立
| | x ( t ) - x ( ∞ ) | | ≤ β α e - αt · - - - ( 4 )
如果
Figure BSA0000094705490000025
满足指数衰减,那么x(t)-x(∞)亦为指数衰减;
步骤1.4:一般情况下,非线性系统(1)的终值状态目标x(∞)并不能随意给定,但对于飞行控制系统而言,如果被控飞行系统的一个或多个状态需要精确地被跟踪达到期望目标,则可以应用下述策略进行实现;假定xi(t)为系统状态的某个待指定的状态且其期望目标为xi(∞),可以增加一个新的状态xn+1(t)使得
Figure BSA0000094705490000029
其中,如果是对多个状态实现跟踪控制,可参照上述方法依次增加新的状态量。那么,仅需针对如下状态扩展系统进行控制器设计
x ^ · ( t ) = f ^ ( x ^ ( t ) , u ( t ) ) · - - - ( 5 )
其中,
Figure BSA0000094705490000027
为扩展系统的状态向量,
Figure BSA0000094705490000028
仍满足如式(2)下连续可微函数条件;
步骤2:、计算正切线性化系统的状态微分方程
步骤2.1:由条件(3)可知,
Figure BSA0000094705490000031
所需满足的指数衰减条件可以利用下述微分方程表示的“微分系统”稳定性来保证
x ^ · · ( t ) = A ( t ) x ^ · ( t ) + B ( t ) u · ( t ) - - - ( 6 )
其中,
A ( t ) = ∂ f ^ ∂ x ^ ( x ^ ( t ) , u ( t ) ) , B ( t ) = ∂ f ^ ∂ u ( x ^ ( t ) , u ( t ) ) - - - ( 7 )
并且,
Figure BSA0000094705490000034
为扩展状态微分方程(5)在给定控制律u(t)情形下的解;
步骤2.2:为求解上述期望非线性控制律u(t),引入下述“正切线性化系统”,其微分方程如下
ζ · ( t ) = A ( t ) ζ ( t ) + B ( t ) ξ ( t ) - - - ( 8 )
其中,
Figure BSA0000094705490000036
Figure BSA0000094705490000037
A(t)、B(t)分别为式(7)中所定义的表达式;
步骤2.3:于是,可以利用现有的线性时变控制方法来保证正切线性化系统(8)的稳定性,针对正切线性化系统(8),设计如下的状态反馈控制律
ξ(t)=K(t)ζ(t)         (9)
综合考虑系统(8)和反馈控制律(9),可以得到如下闭环系统
ζ · ( t ) = [ A ( t ) + B ( t ) K ( t ) ] ζ ( t ) - - - ( 10 )
步骤2.4:针对对于正切线性化系统(8),在任意初始状态ζ(0)=ζ0下,如果存在实数α,β>0使得扩展系统(5)中的非线性控制律u(t)满足如下条件
ξ ( t ) = u · ( t ) = K ( t ) f ^ ( x ^ ( t ) , u ( t ) ) - - - ( 11 )
则被控系统(5)的对应解
Figure BSA00000947054900000310
必定满足如下条件
| | ζ ( t ) | | = | | x ^ · ( t ) | | ≤ βe - αt - - - ( 12 )
也即是,
Figure BSA00000947054900000312
随时间指数衰减;
步骤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 · = f ( x , u ) - - - ( 1 )
其中,x(t)∈Rn和u(t)∈Rp分别表示系统的状态向量和输入向量,f为满足如下条件的连续可微函数
Figure FSA0000094705480000012
步骤1.2:对于非线性控制系统(1),考虑如下跟踪控制问题:给定初始条件x(0)=x0,保证系统状态x(t)随着t→∞渐近收敛于目标状态x(∞);
步骤1.3:对任意t≥0,如果存在一个控制律u(t)及实数α,β>0,使得系统(1)的对应解x(t)满足如下条件
| | x · ( t ) | | ≤ βe - αt - - - ( 3 )
存在一个x(∞)使得下式成立
| | x ( t ) - x ( ∞ ) | | ≤ β α e - αt · - - - ( 4 )
如果
Figure FSA0000094705480000015
满足指数衰减,那么x(t)-x(∞)亦为指数衰减;
步骤1.4:非线性系统(1)的终值状态目标x(∞)不能随意给定;对于飞行控制系统被控飞行系统的一个或多个状态需要精确地被跟踪达到期望目标,则可以应用下述策略进行实现;假定xi(t)为系统状态的某个待指定的状态且其期望目标为xi(∞),可以增加一个新的状态xn+1(t)使得其中,如果是对多个状态实现跟踪控制,可参照上述方法依次增加新的状态量;那么,仅需针对如下状态扩展系统进行控制器设计
x ^ · ( t ) = f ^ ( x ^ ( t ) , u ( t ) ) · - - - ( 5 )
其中,
Figure FSA0000094705480000023
为扩展系统的状态向量,
Figure FSA0000094705480000024
仍满足如下连续可微函数条件(2)。
3.根据权利要求1所述的飞行器非线性跟踪控制系统的正切线性化方法,其特征在于:步骤2中所述的“计算正切线性化系统的状态微分方程”的具体步骤如下:
步骤2.1:由条件(3)可知,
Figure FSA0000094705480000025
所需满足的指数衰减条件可以利用下述微分方程表示的“微分系统”稳定性来保证
x ^ · · ( t ) = A ( t ) x ^ · ( t ) + B ( t ) u · ( t ) - - - ( 6 )
其中,
A ( t ) = ∂ f ^ ∂ x ^ ( x ^ ( t ) , u ( t ) ) , B ( t ) = ∂ f ^ ∂ u ( x ^ ( t ) , u ( t ) ) - - - ( 7 )
并且,
Figure FSA0000094705480000029
为扩展状态微分方程(5)在给定控制律u(t)情形下的解;
步骤2.2:为求解上述期望非线性控制律u(t),引入下述“正切线性化系统”,其微分方程如下
ζ · ( t ) = A ( t ) ζ ( t ) + B ( t ) ξ ( t ) - - - ( 8 )
其中,
Figure FSA00000947054800000211
Figure FSA00000947054800000212
A(t)、B(t)分别为式(7)中所定义的表达式;
步骤2.3:于是,可以利用现有的线性时变控制方法来保证正切线性化系统(8)的稳定性,针对正切线性化系统(8),设计如下的状态反馈控制律
ξ(t)=K(t)ζ(t)       (9)
综合考虑系统(8)和反馈控制律(9),可以得到如下闭环系统
ζ · ( t ) = [ A ( t ) + B ( t ) K ( t ) ] ζ ( t ) - - - ( 10 )
步骤2.4:针对对于正切线性化系统(8),在任意初始状态ζ(0)=ζ0下,如果存在实数σ,β>0使得扩展系统式(5)中的非线性控制律u(t)满足如下条件
ξ ( t ) = u · ( t ) = K ( t ) f ^ ( x ^ ( t ) , u ( t ) ) - - - ( 11 )
则被控系统式(5)的对应解
Figure FSA0000094705480000032
必定满足如下条件
| | ζ ( t ) | | = | | x ^ · ( t ) | | ≤ βe - αt - - - ( 12 )
也即是,随时间指数衰减。
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)与原始的飞行器控制系统形成闭环系统,实现了飞行控制系统的正切线性化跟踪控制设计。
CN201310398354.1A 2013-08-30 2013-08-30 飞行控制系统非线性跟踪控制器设计的正切线性化方法 Active CN103454921B (zh)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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 浙江大学 一种空分节能过程的非线性预测控制系统及方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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
CN103454921A (zh) 飞行控制系统非线性跟踪控制器设计的正切线性化方法
CN103558857A (zh) 一种btt飞行器的分布式复合抗干扰姿态控制方法
CN105182743A (zh) 一种基于鲁棒h无穷的变增益解耦控制方法
Durmaz et al. Sliding mode control for non-linear systems with adaptive sliding surfaces
CN103713517A (zh) 一种飞行控制系统自适应调参方法
CN107037727A (zh) 一种无人直升机大包线自适应增益调度方法
CN104122793B (zh) 一种满足预设性能的导弹过载控制方法
CN103809446B (zh) 飞行器多回路模型簇颤振抑制复合频率鲁棒控制器设计方法
CN103199546A (zh) 动态无功补偿装置最优二次调节器及其设计方法
CN103558761B (zh) 一种具有控制器输入饱和的非线性化学反应循环不确定时滞系统的控制方法
CN103809434B (zh) 纵向飞行模型簇复合根轨迹多级pid控制器设计方法
Pakmehr et al. Adaptive control of uncertain systems with gain scheduled reference models and constrained control inputs
CN104155876A (zh) 一种pid控制器的分离实现方法
CN113467255A (zh) 可重复使用运载器自适应多变量固定时间预设控制方法
CN110032066B (zh) 分数阶非线性系统轨迹跟踪的自适应迭代学习控制方法
CN103809445B (zh) 飞行器多回路模型簇复合pid控制器设计方法
CN103809448B (zh) 飞行器多回路模型簇颤振抑制复合根轨迹补偿鲁棒控制器设计方法
Sun et al. Tracking control via robust dynamic surface control for hypersonic vehicles with input saturation and mismatched uncertainties
CN104298111A (zh) 一种用于窑炉热工参数控制系统的模糊控制方法
CN103823375B (zh) 纵向飞行模型簇复合根轨迹多级pid鲁棒控制器设计方法
CN103809456B (zh) 纵向飞行模型簇人机闭环pid鲁棒控制器设计方法
Di et al. High performance L 1 adaptive control design for longitudinal dynamics of fixed-wing UAV
CN103809457B (zh) 纵向飞行模型簇复合频率控制器设计方法
CN104197792B (zh) 一类多平衡点饱和切换系统的离散增益调度控制器设计方法及实现btt导弹控制的方法
CN103199545B (zh) 动态无功补偿装置最优二次高斯控制器及其设计方法

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