CN109992004B - 一种lpv系统异步切换状态反馈控制器设计方法 - Google Patents

一种lpv系统异步切换状态反馈控制器设计方法 Download PDF

Info

Publication number
CN109992004B
CN109992004B CN201910382619.6A CN201910382619A CN109992004B CN 109992004 B CN109992004 B CN 109992004B CN 201910382619 A CN201910382619 A CN 201910382619A CN 109992004 B CN109992004 B CN 109992004B
Authority
CN
China
Prior art keywords
controller
subsystem
time
switching
interval
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
CN201910382619.6A
Other languages
English (en)
Other versions
CN109992004A (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.)
Harbin University of Science and Technology
Original Assignee
Harbin University of Science and Technology
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 Harbin University of Science and Technology filed Critical Harbin University of Science and Technology
Priority to CN201910382619.6A priority Critical patent/CN109992004B/zh
Publication of CN109992004A publication Critical patent/CN109992004A/zh
Application granted granted Critical
Publication of CN109992004B publication Critical patent/CN109992004B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course or altitude of land, water, air, or space vehicles, e.g. automatic pilot
    • G05D1/10Simultaneous control of position or course in three dimensions
    • G05D1/101Simultaneous control of position or course in three dimensions specially adapted for aircraft

Abstract

本发明公开了一种LPV系统异步切换状态反馈控制器设计方法,包括以下步骤:首先将飞行器的机翼掠角的变化速率设为时变参数,将攻角、俯仰速度设为状态状态变量,设升降舵偏转角为控制输入,抽象出飞行器纵向短周期的LPV模型;其次把时变参数划分为有限个子区间,把各子系统运行的区间划分为控制器与子系统匹配的区间和不匹配的区间;然后根据各控制器模型依赖的滞后时间、系统参数和Lyapunov稳定性条件推导出合适的驻留时间设计参数,利用模型依赖的平均驻留时间策略设计出各个子区间的控制器;最后,将设计的控制器转化成为可解的线性矩阵不等式形式,利用凸优化的思想求解线性矩阵不等式,获得控制器控制增益。通过使用本发明所提供的方法,可以设计鲁棒稳定的异步切换状态反馈控制器,使飞行器在飞行过程中的飞行控制始终具有良好的性能。

Description

一种LPV系统异步切换状态反馈控制器设计方法
技术领域
本发明涉及飞行器飞行过程中的飞行控制领域,具体是一种LPV系统异步切换状态反馈控制器设计方法。
背景技术
在飞行器的飞行控制中,需要改变飞行器的机翼扫掠角以适应任务的要求,例如亚音速目标,超音速目标,侦查和攻击配置,因此,一些参数将会随着机翼扫掠角的变化而变化,如平均气动弦,跨度和机翼面积等,为了应对这些复杂的任务,不同的扫掠角下设计不同的控制器,以达到更好的控制目标;
LPV系统刚好可以满足这种控制需求,将飞行器的扫掠角作为控制系统的时变参数,来进行系统控制器的设计,但由于时变参数的变化范围过大,只设计一个控制器往往无法满足控制需要,因此需要将时变参数划分为有限个互相不重叠的区间,在这些子区间上分别进行控制器的设计,当时变参数变化到相应的子区间时,调用合适的控制器来对系统进行控制,这就需要引入切换控制的思想,切换控制系统由一簇子系统和控制这些子系统的切换逻辑构成,当时变参数变化到相应的区间,所设计的控制率则切换到相应的控制器进行控制,然而,在实际情况下,由于信号传输、系统检测等原因会导致控制器的切换滞后于子系统,这就涉及到了异步切换的思想,以往的异步切换设计方法往往是时序的任意子系统之间的切换,这并不符合LPV系统中时变参数变化曲线光滑可微的事实,因此现有的控制器设计有的没有考虑到异步切换的情况,有的则没有考虑到LPV系统参数变化区间的可微性,是存在缺陷的。
发明内容
有鉴于此,本发明提供一种LPV系统异步切换状态反馈控制器设计方法,充分考虑到个子系统模型的参数,利用模型依赖的平均驻留时间策略设计出相应的控制器,考虑模型依赖的控制器滞后时间对异步切换系统的影响,综合系统参数、模型依赖的控制器滞后时间和Lyapunov稳定性条件推导出合适的驻留时间设计参数,且所设计的控制器能使异步切换系统全局一致指数稳定,从而设计出更符合实际情况的控制器;
本发明提供如下技术方案:一种LPV系统异步切换状态反馈控制器设计方法,该方法包括:
A、将飞行器的机翼掠角的变化速率设为时变参数,将攻角、俯仰速度设为状态状态变量,设升降舵偏转角为控制输入,抽象出飞行器纵向短周期的LPV模型;
B、把时变参数划分为有限个子区间,把各子系统运行的区间划分为控制器与子系统匹配的区间和不匹配的区间;
C、根据各控制器模型依赖的滞后时间、系统参数和Lyapunov稳定性条件推导出合适的驻留时间设计参数,利用模型依赖的平均驻留时间策略设计出各个子区间的控制器;
D、将设计的控制器转化成为可解的线性矩阵不等式形式,利用凸优化的思想求解线性矩阵不等式,获得控制器控制增益。
较佳的,所述步骤A中包括:
考虑飞行器飞行控制过程抽象为LPV模型:
Figure BDA0002052481920000021
式中x(t)∈Rn为状态空间,u(t)∈Rm为控制输入,w(t)为扰动输入,z(t)为期望输出,y(t)为测量输出,参数矩阵A(ρ),B(ρ),C(ρ),D(ρ),E(ρ),均为关于时变参数ρ的连续函数矩阵,ρ(t)是关于时间t的连续函数,ρ的值实时可测且有界,即
Figure BDA0002052481920000022
ρ(t)的变化率也被认为是有界的即
Figure BDA0002052481920000023
为方便表示,用ρ表示ρ(t);
则其相应的闭环系统为
Figure BDA0002052481920000024
式中
Figure BDA0002052481920000025
相应的控制器为:
u(t)=Kx(t) (3)
其中K为控制器的控制增益;使得所述由LPV系统与控制器K组成的闭环系统满足如下条件:
(1)闭环系统全局一致指数稳定;
(2)期望输出z和扰动w满足H∞性能指标:
Figure BDA0002052481920000031
较佳的,考虑步骤B包括:
当参数变化范围过大时,通过此范围所设计的控制器往往会具有很大的保守性,因此我们将参数的变化范围Ω划分为有限个子空间Ωi,i∈M,M={1,2,3,…,N},并在这些子空间上分别设计控制器,设划分的区间为:
Figure BDA0002052481920000032
Figure BDA0002052481920000033
Figure BDA0002052481920000034
则相应的子系统可以描述为:
Figure BDA0002052481920000035
相应的子系统的状态反馈控制器为:
u(t)=Ki(ρ(t))x(t),i=1,2,3,...N (8)
则系统(1)的闭环控制系统为:
Figure BDA0002052481920000036
式中
Figure BDA0002052481920000037
在实际情况下,由于信号传输、系统检测等原因会导致控制器的切换滞后于子系统,基于这类情况,将子系统运行的区间划分为子系统与控制器相匹配的区间和非匹配的区间,假设ρ∈Ωi时,即第i个子系统运行时间为t∈[tn,tn+1),也就是第σ(tn)=i个子系统运行,对于系统(1),考虑其控制输入为u(t)=Kσ(t)x(t),σ’(t)表示控制器的切换信号,Kσ(t)为所设计的控制器增益,为了与实际情况相符,假设控制器落后子系统的时间Δn≤τn=tn+1-tn,n=0,1,2…,因此控制器的切换时间序列可以表示为{(t00,σ’(t0)),(t11,σ’(t1)),…(tn-1n-1,σ’(tn-1)),…,n=0,1,2,…};
假设t=tn时,第σ(tn)=i个子系统运行,t=tn+1时第σ(tn+1)=j个子系统运行,则相应的控制器的动作时间分别为tnn,tn+1n+1,由于控制器的切换落后于子系统的切换,故在系统运行时间(t0,t)中存在控制器与子系统匹配的区间T+(t0,t)和不匹配的区间T-(t0,t),两种情况可以分别进行讨论,由系统(7)可以得到整个异步切换控制系统可以描述为:
Figure BDA0002052481920000041
Figure BDA0002052481920000042
Figure BDA0002052481920000043
较佳的,考虑步骤C包括:
对于切换连续线性系统(1),如果存在常数c>0,δ>0使得相应的系统解在任意初始条件下对所有的t≥t0满足
Figure BDA0002052481920000044
则称该系统在平衡态x=0时在此切换信号下是全局一致指数稳定的;
若使用切换率σ(t)控制切换系统,当t∈(t0,t)时,使用Nσ(t0,t)表示系统在此区间上的切换次数,若存在正实数N0和τa使得
Figure BDA0002052481920000045
则称τa为此切换率的平均驻留时间,N0表示振动幅度;
平均驻留时间强调的是相邻的两次切换之间的时间不得小于一个常数τa
若使用切换率σ(t)控制具有M个子系统的切换系统,当t∈(t0,t)时,使用Nσ,i(t0,t)表示第i个子系统在此区间上的切换次数,若存在相应的正实数N0,i和τi使得
Figure BDA0002052481920000046
则称τi为第i个子系统在此切换率下的模型依赖的驻留时间,N0,i表示第i个子系统的振动幅度,Ti(t0,t)表示在时间(t0,t)内第i个子系统运行的总时间;
在平均驻留时间策略中,每个子系统均采用一样的参数来决定各子系统的驻留时间,也就是说,平均驻留时间策略并非是模型依赖的,而模型依赖的驻留时间策略充分利用了所有的系统中的参数,因此这种策略具有更小的保守性,本专利采用模型依赖的驻留时间策略,各个子系统的驻留时间由其系统参数决定,用τi,i∈M表示第i个子系统的驻留时间;
若存在正实数α,γ使得下面两个条件成立,则称系统具有指数的H性能γ:
(1)当扰动为零时,系统指数稳定;
(2)在零初始条件下,扰动不为零系统满足
Figure BDA0002052481920000051
定理1、给定正常数αi>0和βi>0,μi>1,i∈M,存在一组正定实对称矩阵Pi∈Rn×n,和一个正实数γ,i,j∈M,i≠j,有
Figure BDA0002052481920000052
Figure BDA0002052481920000053
Pi≤μiPj (15)
Figure BDA0002052481920000054
Figure BDA0002052481920000055
则闭环系统(10)在切换率σ’(t)下全局一致指数稳定,且满足指数H性能指标γ,式中ζi∈(0,αi),i∈M,
Figure BDA0002052481920000056
考虑在控制过程中,各子系统的控制器滞后的时间往往与子系统本身相关,因此,在基于模型依赖的切换策略下,将Δi,αi,βi和μi综合考虑,进而可以得到各子系统合适的驻留时间满足系统全局一致指数稳定;
考虑闭环系统(33)的Lyapunov函数为:
Vi(t)=xT(t)Pix(t),i∈M (18)
当扰动为零时,由式(13)和式(14)分别可知,当t∈[tnn,tn+1)时有
Figure BDA0002052481920000057
对上式求解微分方程可得
Figure BDA0002052481920000061
当t∈[tn,tnn)时有
Figure BDA0002052481920000062
对上式求解微分方程可得
Figure BDA0002052481920000063
因此当t∈[tn,tn+1)时为一整个切换工作周期,故当t∈[tn,tn+1)时有
Figure BDA0002052481920000064
Figure BDA0002052481920000065
由式(23)可得
Figure BDA0002052481920000066
由式(17)可知
Figure BDA0002052481920000067
因此可知
Figure BDA0002052481920000068
由式(16)式(26)和条件1可知当t→∞时,Vσ(t)(t)收敛于零,因此当式(15)式(16)和式(17)成立时,闭环系统(10)全局一致渐近稳定,另外取
Figure BDA0002052481920000069
时,式(26)满足形式
Figure BDA0002052481920000071
可知系统(10)全局一致指数稳定;
若扰动不为零,则t∈[tnn,tn+1)时有
Figure BDA0002052481920000072
可推出
Figure BDA0002052481920000073
当t∈[tn,tnn)时有
Figure BDA0002052481920000074
可推出
Figure BDA0002052481920000075
式中Γ(t)=zT(t)z(t)-γ2wT(t)w(t);
假设t∈[tn,tn+1)时,第σ(tn)=i个子系统运行,t∈[tn+1,tn+2)时第σ(tn+1)=j个子系统运行i,j∈M,n=0,1,2…故由式(15)式(28)和式(30)可知
Figure BDA0002052481920000076
此处主要研究系统输出z(t)对扰动w(t)的抑制性能,故取x(t0)=0,将Γ(s)展开根据对应项的小于关系和式(31)可得
Figure BDA0002052481920000081
将积分限合并,并根据式(17)可得
Figure BDA0002052481920000082
当t∈[t0,t)时式(33)两边同乘
Figure BDA0002052481920000083
可得
Figure BDA0002052481920000084
对上式两边同乘
Figure BDA0002052481920000085
可得
Figure BDA0002052481920000086
由式(17)可知
Figure BDA0002052481920000087
可推出
Figure BDA0002052481920000088
化形可得
Figure BDA0002052481920000089
对上式从t=0到∞的积分可得
Figure BDA0002052481920000091
至此可知,所设计的异步切换控制器满足指数稳定,且满足指数H性能指标γ。
较佳的,考虑步骤D包括:
将得到的稳定性条件转化为可解的线性矩阵不等式形式:
定理2、给定正常数αi>0和βi>0,μi>1,i∈M,存在正定实对称矩阵Si∈Rn×n,Qi∈Rm,i∈M,和一个正实数γ,i,j∈M,i≠j,有
Figure BDA0002052481920000092
Figure BDA0002052481920000093
Sj≤μiSi (41)
Figure BDA0002052481920000094
Figure BDA0002052481920000095
则闭环系统(10)在切换率σ’(t)下全局一致指数稳定,且满足H性能指标γ,式中
Figure BDA0002052481920000096
控制器的增益可由下式得到
Figure BDA0002052481920000097
在定理1中令
Figure BDA0002052481920000098
KiSi=Qi,由式(13)和式(14)可知式(39)和式(40)成立,由式(15)可知
PiiPj≤0 (45)
由上式可推出
Figure BDA0002052481920000101
对上式进行变换可得
-Pi -1-IT(-μiPj)-1I≤0 (47)
上式中令
Figure BDA0002052481920000102
可得Sj≤μiSi,可知式(41)成立;
至此,得到所设计的切换LPV系统异步控制器的控制增益。
如上所见,本发明中的LPV系统异步切换状态反馈控制器设计方法,针对飞行器的飞行控制LPV系统模型,考虑到在异步切换中不同子系统的控制器的滞后时间也往往各不相同,为得到更好的系统性能,使用模型依赖的驻留时间策略对各子系统的控制器及驻留时间分别进行设计,当控制器与子系统不匹配时,利用各子系统满足模型依赖的平均驻留时间切换策略,得到整个系统的异步切换工作过程的稳定性条件,本发明中显然可见,在不匹配区间时其Lyapunov函数并非严格递减,这降低了对此区间的条件约束,但也意味着在非匹配区间不能保证系统严格稳定,故通过限制子系统与控制器的匹配时间与非匹配时间的比值来使整个系统达到全局一致指数稳定,而且这个控制器滞后子系统的时间也可以是模型依赖的,最后,将稳定性条件转换为可解的LMIs,并获得相应的控制增益;
本发明的有益效果为:与现有技术相比,本发明考虑在控制过程中,各子系统的控制器滞后的时间往往与子系统本身相关,因此,在基于模型依赖的切换策略下,将Δi,αi,βi和μi综合考虑,进而可以得到各子系统合适的驻留时间满足系统全局一致指数稳定。
附图说明
图1为LPV系统异步切换状态反馈控制器设计方法的流程示意图。
图2为时变参数的变化曲线。
图3为切换信号。
图4为飞行器系统在所设计的控制器控制下的状态变化曲线。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下参照附图并举实施例,对本发明进一步详细说明;
本发明中,提供了一种LPV系统异步切换状态反馈控制器设计方法,该方法适用于飞行器飞行控制LPV系统;
图1为本发明中的LPV系统异步切换状态反馈控制器设计方法的流程示意图,如图1所示,本发明实施例中的LPV系统异步切换状态反馈控制器设计方法,包括以下步骤:
步骤一,将飞行器的机翼掠角的变化速率设为时变参数,将攻角、俯仰速度设为状态状态变量,设升降舵偏转角为控制输入,抽象出飞行器纵向短周期的LPV模型;
在本发明的技术方案中,可以先考虑将系统抽象成一个LPV模型;
例如,较佳的,在本发明的具体实施例中,上述的步骤一包括:
考虑飞行器飞行控制过程抽象为LPV模型:
Figure BDA0002052481920000111
式中x(t)∈Rn为状态空间,u(t)∈Rm为控制输入,w(t)为扰动输入,z(t)为期望输出,y(t)为测量输出,参数矩阵A(ρ),B(ρ),C(ρ),D(ρ),E(ρ),均为关于时变参数ρ的连续函数矩阵,ρ(t)是关于时间t的连续函数,ρ的值实时可测且有界,即
Figure BDA0002052481920000112
ρ(t)的变化率也被认为是有界的即
Figure BDA0002052481920000113
为方便表示,用ρ表示ρ(t);
则其相应的闭环系统为
Figure BDA0002052481920000114
式中
Figure BDA0002052481920000115
相应的控制器为:
u(t)=Kx(t) (3)
其中K为控制器的控制增益;使得所述由LPV系统与控制器K组成的闭环系统满足如下条件:
(1)闭环系统全局一致指数稳定;
(2)期望输出z和扰动w满足H∞性能指标:
Figure BDA0002052481920000116
步骤二,把时变参数划分为有限个子区间,把各子系统运行的区间划分为控制器与子系统匹配的区间和不匹配的区间;
在本发明的技术方案中,将时变参数的变化范围划分为有限个互不重叠的子区间,并在这些子空间上分别进行控制器的设计;
例如,较佳的,在本发明的具体实施例中,先说明一些证明的数学准备:
当参数变化范围过大时,通过此范围所设计的控制器往往会具有很大的保守性,因此我们将参数的变化范围Ω划分为有限个子空间Ωi,i∈M,M={1,2,3,…,N},并在这些子空间上分别设计控制器,设划分的区间为:
Figure BDA0002052481920000121
Figure BDA0002052481920000122
Figure BDA0002052481920000123
则相应的子系统可以描述为:
Figure BDA0002052481920000124
相应的子系统的状态反馈控制器为:
u(t)=Ki(ρ(t))x(t),i=1,2,3,...N (8)
则系统(1)的闭环控制系统为:
Figure BDA0002052481920000125
式中
Figure BDA0002052481920000126
在实际情况下,由于信号传输、系统检测等原因会导致控制器的切换滞后于子系统,基于这类情况,将子系统运行的区间划分为子系统与控制器相匹配的区间和非匹配的区间,假设ρ∈Ωi时,即第i个子系统运行时间为t∈[tn,tn+1),也就是第σ(tn)=i个子系统运行,对于系统(1),考虑其控制输入为u(t)=Kσ(t)x(t),σ’(t)表示控制器的切换信号,Kσ(t)为所设计的控制器增益,为了与实际情况相符,假设控制器落后子系统的时间Δn≤τn=tn+1-tn,n=0,1,2…,因此控制器的切换时间序列可以表示为{(t00,σ’(t0)),(t11,σ’(t1)),…(tn-1n-1,σ’(tn-1)),…,n=0,1,2,…};
假设t=tn时,第σ(tn)=i个子系统运行,t=tn+1时第σ(tn+1)=j个子系统运行,则相应的控制器的动作时间分别为tnn,tn+1n+1,由于控制器的切换落后于子系统的切换,故在系统运行时间(t0,t)中存在控制器与子系统匹配的区间T+(t0,t)和不匹配的区间T-(t0,t),两种情况可以分别进行讨论,由系统(7)可以得到整个异步切换控制系统可以描述为:
Figure BDA0002052481920000131
Figure BDA0002052481920000132
Figure BDA0002052481920000133
步骤三,根据各控制器模型依赖的滞后时间、系统参数和Lyapunov稳定性条件推导出合适的驻留时间设计参数,利用模型依赖的平均驻留时间策略设计出各个子区间的控制器;
在本发明的技术方案中,考虑在控制过程中,各子系统的控制器滞后的时间往往与子系统本身相关,因此,在基于模型依赖的切换策略下,将Δi,αi,βi和μi综合考虑,进而可以得到各子系统合适的驻留时间满足系统全局一致指数稳定;
例如,较佳的,在本发明的具体实施例中,所述步骤三包括:
对于切换连续线性系统(1),如果存在常数c>0,δ>0使得相应的系统解在任意初始条件下对所有的t≥t0满足
Figure BDA0002052481920000134
则称该系统在平衡态x=0时在此切换信号下是全局一致指数稳定的;
若使用切换率σ(t)控制切换系统,当t∈(t0,t)时,使用Nσ(t0,t)表示系统在此区间上的切换次数,若存在正实数N0和τa使得
Figure BDA0002052481920000135
则称τa为此切换率的平均驻留时间,N0表示振动幅度;
平均驻留时间强调的是相邻的两次切换之间的时间不得小于一个常数τa
若使用切换率σ(t)控制具有M个子系统的切换系统,当t∈(t0,t)时,使用Nσ,i(t0,t)表示第i个子系统在此区间上的切换次数,若存在相应的正实数N0,i和τi使得
Figure BDA0002052481920000141
则称τi为第i个子系统在此切换率下的模型依赖的驻留时间,N0,i表示第i个子系统的振动幅度,Ti(t0,t)表示在时间(t0,t)内第i个子系统运行的总时间;
在平均驻留时间策略中,每个子系统均采用一样的参数来决定各子系统的驻留时间,也就是说,平均驻留时间策略并非是模型依赖的,而模型依赖的驻留时间策略充分利用了所有的系统中的参数,因此这种策略具有更小的保守性,本专利采用模型依赖的驻留时间策略,各个子系统的驻留时间由其系统参数决定,用τi,i∈M表示第i个子系统的驻留时间;
若存在正实数α,γ使得下面两个条件成立,则称系统具有指数的H性能γ:
(1)当扰动为零时,系统指数稳定;
(2)在零初始条件下,扰动不为零系统满足
Figure BDA0002052481920000142
定理1、给定正常数αi>0和βi>0,μi>1,i∈M,存在一组正定实对称矩阵Pi∈Rn×n,和一个正实数γ,i,j∈M,i≠j,有
Figure BDA0002052481920000143
Figure BDA0002052481920000144
Pi≤μiPj (15)
Figure BDA0002052481920000145
Figure BDA0002052481920000146
则闭环系统(10)在切换率σ’(t)下全局一致指数稳定,且满足指数H性能指标γ,式中ζi∈(0,αi),i∈M,
Figure BDA0002052481920000147
考虑在控制过程中,各子系统的控制器滞后的时间往往与子系统本身相关,因此,在基于模型依赖的切换策略下,将Δi,αi,βi和μi综合考虑,进而可以得到各子系统合适的驻留时间满足系统全局一致指数稳定;
考虑闭环系统(33)的Lyapunov函数为:
Vi(t)=xT(t)Pix(t),i∈M (18)
当扰动为零时,由式(13)和式(14)分别可知,当t∈[tnn,tn+1)时有
Figure BDA0002052481920000151
对上式求解微分方程可得
Figure BDA0002052481920000152
当t∈[tn,tnn)时有
Figure BDA0002052481920000153
对上式求解微分方程可得
Figure BDA0002052481920000154
因此当t∈[tn,tn+1)时为一整个切换工作周期,故当t∈[tn,tn+1)时有
Figure BDA0002052481920000155
Figure BDA0002052481920000156
由式(23)可得
Figure BDA0002052481920000157
由式(17)可知
Figure BDA0002052481920000158
因此可知
Figure BDA0002052481920000161
由式(16)式(26)和条件1可知当t→∞时,Vσ(t)(t)收敛于零,因此当式(15)式(16)和式(17)成立时,闭环系统(10)全局一致渐近稳定,另外取
Figure BDA0002052481920000162
时,式(26)满足形式
Figure BDA0002052481920000163
可知系统(10)全局一致指数稳定;
若扰动不为零,则t∈[tnn,tn+1)时有
Figure BDA0002052481920000164
可推出
Figure BDA0002052481920000165
当t∈[tn,tnn)时有
Figure BDA0002052481920000166
可推出
Figure BDA0002052481920000167
式中Γ(t)=zT(t)z(t)-γ2wT(t)w(t);
假设t∈[tn,tn+1)时,第σ(tn)=i个子系统运行,t∈[tn+1,tn+2)时第σ(tn+1)=j个子系统运行i,j∈M,n=0,1,2…故由式(15)式(28)和式(30)可知
Figure BDA0002052481920000171
此处主要研究系统输出z(t)对扰动w(t)的抑制性能,故取x(t0)=0,将Γ(s)展开根据对应项的小于关系和式(31)可得
Figure BDA0002052481920000172
将积分限合并,并根据式(17)可得
Figure BDA0002052481920000173
当t∈[t0,t)时式(33)两边同乘
Figure BDA0002052481920000174
可得
Figure BDA0002052481920000175
对上式两边同乘
Figure BDA0002052481920000176
可得
Figure BDA0002052481920000181
由式(17)可知
Figure BDA0002052481920000182
可推出
Figure BDA0002052481920000183
化形可得
Figure BDA0002052481920000184
对上式从t=0到∞的积分可得
Figure BDA0002052481920000185
至此可知,所设计的异步切换控制器满足指数稳定,且满足指数H性能指标γ。
步骤四,将设计的控制器转化成为可解的线性矩阵不等式形式,利用凸优化的思想求解线性矩阵不等式,获得控制器控制增益;
在本发明的技术方案中,可以使用多种具体的实施方式来实现上述的步骤四;以下将以其中的一种实现方式为例,对本发明的技术方案进行详细的介绍;
例如,较佳的,在本发明的具体实施例中,所述步骤四包括:
将得到的稳定性条件转化为可解的线性矩阵不等式形式:
定理2、给定正常数αi>0和βi>0,μi>1,i∈M,存在正定实对称矩阵Si∈Rn×n,Qi∈Rm,i∈M,和一个正实数γ,i,j∈M,i≠j,有
Figure BDA0002052481920000186
Figure BDA0002052481920000187
Sj≤μiSi (41)
Figure BDA0002052481920000191
Figure BDA0002052481920000192
则闭环系统(10)在切换率σ’(t)下全局一致指数稳定,且满足H性能指标γ,式中
Figure BDA0002052481920000193
控制器的增益可由下式得到
Figure BDA0002052481920000194
在定理1中令
Figure BDA0002052481920000195
KiSi=Qi,由式(13)和式(14)可知式(39)和式(40)成立,由式(15)可知
PiiPj≤0 (45)
由上式可推出
Figure BDA0002052481920000196
对上式进行变换可得
-Pi -1-IT(-μiPj)-1I≤0 (47)
上式中令
Figure BDA0002052481920000197
可得Sj≤μiSi,可知式(41)成立;
至此,得到所设计的切换LPV系统异步控制器的控制增益;
综上可知,在本发明中的LPV系统异步切换状态反馈控制器设计方法中,将飞行器的机翼掠角的变化速率设为时变参数,将攻角、俯仰速度设为状态状态变量,设升降舵偏转角为控制输入,抽象出飞行器纵向短周期的LPV模型;把时变参数划分为有限个子区间,把各子系统运行的区间划分为控制器与子系统匹配的区间和不匹配的区间;根据各控制器模型依赖的滞后时间、系统参数和Lyapunov稳定性条件推导出合适的驻留时间设计参数,利用模型依赖的平均驻留时间策略设计出各个子区间的控制器;将设计的控制器转化成为可解的线性矩阵不等式形式,利用凸优化的思想求解线性矩阵不等式,获得控制器控制增益;为进一步说明本发明的优越性,提供相关的仿真:在整个变化区间上是无法设计可行的LPV控制器的,在有序切换过程中,将参数的变化范围分为[0,1.6),[1.6,3)三个变化区间,分别进行控制器的设计,图2为时变参数的变化曲线,图3为切换信号,图4为状态响应曲线;
另外,由于本方法具有一定的通用性,所以通过使用该方法可以对所有可抽象为切换LPV模型的实际物理系统设计异步切换状态反馈控制器,以达到良好的控制效果;
对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神或基本特征的情况下,能够以其他的具体形式实现本发明;因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内,不应将权利要求中的任何附图标记视为限制所涉及的权利要求;
此外,应当理解,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。

Claims (1)

1.一种LPV系统异步切换状态反馈控制器设计方法,其特征在于,该方法包括:
A、将飞行器的机翼扫掠角的变化速率设为时变参数,将攻角、俯仰速度设为状态状态变量,设升降舵偏转角为控制输入,抽象出飞行器纵向短周期的LPV模型;
B、把时变参数划分为有限个子区间,把各子系统运行的区间划分为控制器与子系统匹配的区间和不匹配的区间;
C、根据各控制器模型依赖的滞后时间、系统参数和Lyapunov稳定性条件推导出合适的驻留时间设计参数,利用模型依赖的平均驻留时间策略设计出各个子区间的控制器;
D、将设计的控制器转化成为可解的线性矩阵不等式形式,利用凸优化的思想求解线性矩阵不等式,获得控制器控制增益;
所述步骤A包括:
考虑飞行器飞行控制过程抽象为LPV模型:
Figure FDA0003559610750000011
式中x(t)∈Rn为状态空间,u(t)∈Rm为控制输入,w(t)为扰动输入,z(t)为期望输出,y(t)为测量输出,参数矩阵A(ρ),B(ρ),C(ρ),D(ρ),E(ρ),均为关于时变参数ρ的连续函数矩阵,ρ(t)是关于时间t的连续函数,ρ的值实时可测且有界,即
Figure FDA0003559610750000012
ρ(t)的变化率也被认为是有界的即
Figure FDA0003559610750000013
为方便表示,用ρ表示ρ(t);
则其相应的闭环系统为
Figure FDA0003559610750000014
式中
Figure FDA0003559610750000015
相应的控制器为:
u(t)=Kx(t) (3)
其中K为控制器的控制增益;使得所述由LPV系统与控制器K组成的闭环系统满足如下条件:
(1)闭环系统全局一致指数稳定;
(2)期望输出z和扰动w满足H∞性能指标:
Figure FDA0003559610750000016
所述步骤B包括:
当参数变化范围过大时,通过此范围所设计的控制器往往会具有很大的保守性,因此我们将参数的变化范围Ω划分为有限个子空间Ωi,i∈M,M={1,2,3,…,N},并在这些子空间上分别设计控制器,设划分的区间为:
Figure FDA0003559610750000021
Figure FDA0003559610750000022
Figure FDA0003559610750000023
则相应的子系统可以描述为:
Figure FDA0003559610750000024
相应的子系统的状态反馈控制器为:
u(t)=Ki(ρ(t))x(t),i=1,2,3,...N (8)
则系统(1)的闭环控制系统为:
Figure FDA0003559610750000025
式中
Figure FDA0003559610750000026
在实际情况下,由于信号传输、系统检测等原因会导致控制器的切换滞后于子系统,基于这类情况,将子系统运行的区间划分为子系统与控制器相匹配的区间和非匹配的区间,假设ρ∈Ωi时,即第i个子系统运行时间为t∈[tn,tn+1),也就是第σ(tn)=i个子系统运行,对于系统(1),考虑其控制输入为u(t)=Kσ’(t)x(t),σ’(t)表示控制器的切换信号,Kσ’(t)为所设计的控制器增益,为了与实际情况相符,假设控制器落后子系统的时间Δn≤τn=tn+1-tn,n=0,1,2…,因此控制器的切换时间序列可以表示为{(t00,σ’(t0)),(t11,σ’(t1)),…(tn-1n-1,σ’(tn-1)),…,n=0,1,2,…};
假设t=tn时,第σ(tn)=i个子系统运行,t=tn+1时第σ(tn+1)=j个子系统运行,则相应的控制器的动作时间分别为tnn,tn+1n+1,由于控制器的切换落后于子系统的切换,故在系统运行时间(t0,t)中存在控制器与子系统匹配的区间T+(t0,t)和不匹配的区间T-(t0,t),两种情况可以分别进行讨论,由系统(7)可以得到整个异步切换控制系统可以描述为:
Figure FDA0003559610750000031
Figure FDA0003559610750000032
Figure FDA0003559610750000033
所述步骤C包括:
对于切换连续线性系统(1),如果存在常数c>0,δ>0使得相应的系统解在任意初始条件下对所有的t≥t0满足
Figure FDA0003559610750000034
则称该系统在平衡态x=0时在此切换信号下是全局一致指数稳定的;
若使用切换率σ(t)控制切换系统,当t∈(t0,t)时,使用Nσ(t0,t)表示系统在此区间上的切换次数,若存在正实数N0和τa使得
Figure FDA0003559610750000035
则称τa为此切换率的平均驻留时间,N0表示振动幅度;
平均驻留时间强调的是相邻的两次切换之间的时间不得小于一个常数τa
若使用切换率σ(t)控制具有M个子系统的切换系统,当t∈(t0,t)时,使用Nσ,i(t0,t)表示第i个子系统在此区间上的切换次数,若存在相应的正实数N0,i和τi使得
Figure FDA0003559610750000036
则称τi为第i个子系统在此切换率下的模型依赖的驻留时间,N0,i表示第i个子系统的振动幅度,Ti(t0,t)表示在时间(t0,t)内第i个子系统运行的总时间;
在平均驻留时间策略中,每个子系统均采用一样的参数来决定各子系统的驻留时间,也就是说,平均驻留时间策略并非是模型依赖的,而模型依赖的驻留时间策略充分利用了所有的系统中的参数,因此这种策略具有更小的保守性,本专利采用模型依赖的驻留时间策略,各个子系统的驻留时间由其系统参数决定,用τi,i∈M表示第i个子系统的驻留时间;
若存在正实数α,γ使得下面两个条件成立,则称系统具有指数的H性能γ:
(1)当扰动为零时,系统指数稳定;
(2)在零初始条件下,扰动不为零系统满足
Figure FDA0003559610750000041
定理1、给定正常数αi>0和βi>0,μi>1,i∈M,存在一组正定实对称矩阵Pi∈Rn×n,和一个正实数γ,i,j∈M,i≠j,有
Figure FDA0003559610750000042
Figure FDA0003559610750000043
Pi≤μiPj (15)
Figure FDA0003559610750000044
Figure FDA0003559610750000045
则闭环系统(10)在切换率σ’(t)下全局一致指数稳定,且满足指数H性能指标γ,
式中ζi∈(0,αi),i∈M,
Figure FDA0003559610750000046
考虑在控制过程中,各子系统的控制器滞后的时间往往与子系统本身相关,因此,在基于模型依赖的切换策略下,将Δi,αi,βi和μi综合考虑,进而可以得到各子系统合适的驻留时间满足系统全局一致指数稳定;
考虑闭环系统(33)的Lyapunov函数为:
Vi(t)=xT(t)Pix(t),i∈M (18)
当扰动为零时,由式(13)和式(14)分别可知,当t∈[tnn,tn+1)时有
Figure FDA0003559610750000051
对上式求解微分方程可得
Figure FDA0003559610750000052
当t∈[tn,tnn)时有
Figure FDA0003559610750000053
对上式求解微分方程可得
Figure FDA0003559610750000054
因此当t∈[tn,tn+1)时为一整个切换工作周期,故当t∈[tn,tn+1)时有
Figure FDA0003559610750000055
Figure FDA0003559610750000056
由式(23)可得
Figure FDA0003559610750000057
由式(17)可知
Figure FDA0003559610750000058
因此可知
Figure FDA0003559610750000059
由式(16)式(26)和条件1可知当t→∞时,Vσ(t)(t)收敛于零,因此当式(15)式(16)和式(17)成立时,闭环系统(10)全局一致渐近稳定,另外取
Figure FDA0003559610750000061
Figure FDA0003559610750000062
时,式(26)满足形式
Figure FDA0003559610750000063
可知系统(10)全局一致指数稳定;
若扰动不为零,则t∈[tnn,tn+1)时有
Figure FDA0003559610750000064
可推出
Figure FDA0003559610750000065
当t∈[tn,tnn)时有
Figure FDA0003559610750000066
可推出
Figure FDA0003559610750000067
式中Γ(t)=zT(t)z(t)-γ2wT(t)w(t);
假设t∈[tn,tn+1)时,第σ(tn)=i个子系统运行,t∈[tn+1,tn+2)时第σ(tn+1)=j个子系统运行i,j∈M,n=0,1,2…故由式(15)式(28)和式(30)可知
Figure FDA0003559610750000068
此处主要研究系统输出z(t)对扰动w(t)的抑制性能,故取x(t0)=0,将Γ(s)展开根据对应项的小于关系和式(31)可得
Figure FDA0003559610750000071
将积分限合并,并根据式(17)可得
Figure FDA0003559610750000072
当t∈[t0,t)时式(33)两边同乘
Figure FDA0003559610750000073
可得
Figure FDA0003559610750000074
对上式两边同乘
Figure FDA0003559610750000075
可得
Figure FDA0003559610750000076
由式(17)可知
Figure FDA0003559610750000077
可推出
Figure FDA0003559610750000078
化形可得
Figure FDA0003559610750000079
对上式从t=0到∞的积分可得
Figure FDA0003559610750000081
至此可知,所设计的异步切换控制器满足指数稳定,且满足指数H性能指标γ;
所述步骤D包括:
将得到的稳定性条件转化为可解的线性矩阵不等式形式:
定理2、给定正常数αi>0和βi>0,μi>1,i∈M,存在正定实对称矩阵Si∈Rn×n,Qi∈Rm,i∈M,和一个正实数γ,i,j∈M,i≠j,有
Figure FDA0003559610750000082
Figure FDA0003559610750000083
Sj≤μiSi (41)
Figure FDA0003559610750000084
Figure FDA0003559610750000085
则闭环系统(10)在切换率σ’(t)下全局一致指数稳定,且满足H性能指标γ,
式中
Figure FDA0003559610750000086
控制器的增益可由下式得到
Figure FDA0003559610750000087
在定理1中令
Figure FDA0003559610750000088
KiSi=Qi,由式(13)和式(14)可知式(39)和式(40)成立,由式(15)可知
PiiPj≤0 (45)
由上式可推出
Figure FDA0003559610750000091
对上式进行变换可得
-Pi -1-IT(-μiPj)-1I≤0 (47)
上式中令
Figure FDA0003559610750000092
可得Sj≤μiSi,可知式(41)成立;
至此,得到所设计的切换LPV系统异步控制器的控制增益。
CN201910382619.6A 2019-05-08 2019-05-08 一种lpv系统异步切换状态反馈控制器设计方法 Active CN109992004B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910382619.6A CN109992004B (zh) 2019-05-08 2019-05-08 一种lpv系统异步切换状态反馈控制器设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910382619.6A CN109992004B (zh) 2019-05-08 2019-05-08 一种lpv系统异步切换状态反馈控制器设计方法

Publications (2)

Publication Number Publication Date
CN109992004A CN109992004A (zh) 2019-07-09
CN109992004B true CN109992004B (zh) 2022-04-22

Family

ID=67136257

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910382619.6A Active CN109992004B (zh) 2019-05-08 2019-05-08 一种lpv系统异步切换状态反馈控制器设计方法

Country Status (1)

Country Link
CN (1) CN109992004B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111142550B (zh) * 2020-01-09 2021-07-27 上海交通大学 民用飞机辅助驾驶控制方法、系统及飞行品质评估方法
CN113848823A (zh) * 2021-08-17 2021-12-28 河南财经政法大学 一种动态智能生产系统的控制方法

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1567107A (zh) * 2003-06-09 2005-01-19 石油大学(北京) 复杂滞后过程的先进控制方法及其系统
CN102955428A (zh) * 2012-11-14 2013-03-06 华侨大学 基于lpv模型的满足设定点跟踪与扰动抑制性能的pi控制方法
WO2013187975A1 (en) * 2012-06-15 2013-12-19 Abb Research Ltd. Parallel computation of dynamic state estimation for power system
CN104505830A (zh) * 2015-01-14 2015-04-08 华北电力大学 时滞电力系统稳定性分析方法和装置
CN104698842A (zh) * 2015-02-06 2015-06-10 浙江大学 一种基于内点法的lpv模型非线性预测控制方法
CN104991566A (zh) * 2015-07-07 2015-10-21 北京航天自动控制研究所 一种用于高超声速飞行器的参数不确定性lpv系统建模方法
CN105094114A (zh) * 2015-06-23 2015-11-25 北京航天自动控制研究所 一种确定lpv控制器的稳定性的方法
CN105242540A (zh) * 2015-10-27 2016-01-13 西安石油大学 基于平均驻留时间的网络控制系统动态切换控制方法
CN106094511A (zh) * 2016-06-18 2016-11-09 哈尔滨理工大学 一种时滞lpv系统的鲁棒h∞状态反馈控制器的设计方法
EP3109725A1 (en) * 2015-06-12 2016-12-28 Innovative Solutions and Support, Inc. Upgraded flight management system for autopilot control and method of providing the same
CN106371313A (zh) * 2016-09-22 2017-02-01 哈尔滨理工大学 一种时滞lpv系统有记忆h∞状态反馈控制器设计方法
CN106843261A (zh) * 2016-10-25 2017-06-13 南京航空航天大学 一种变体飞行器过渡段的张量积插值建模与控制方法
CN107561939A (zh) * 2017-09-01 2018-01-09 哈尔滨理工大学 一种lpv系统输出反馈控制器
CN108628161A (zh) * 2017-03-23 2018-10-09 哈尔滨理工大学 一种切换lpv系统的设计方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110106345A1 (en) * 2009-11-03 2011-05-05 Takacs Robert S Low visibility landing system
CN105179164B (zh) * 2015-06-25 2018-11-09 江苏科技大学 基于t-s模糊模型的风能转换系统滑模控制方法及装置
CN109507878A (zh) * 2017-09-15 2019-03-22 中国科学院沈阳自动化研究所 基于平均驻留时间切换的状态反馈切换控制器设计方法
CN108762284B (zh) * 2018-05-17 2020-07-28 北京航空航天大学 一种基于lpv技术的航天器姿态跟踪控制方法与装置

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1567107A (zh) * 2003-06-09 2005-01-19 石油大学(北京) 复杂滞后过程的先进控制方法及其系统
WO2013187975A1 (en) * 2012-06-15 2013-12-19 Abb Research Ltd. Parallel computation of dynamic state estimation for power system
CN102955428A (zh) * 2012-11-14 2013-03-06 华侨大学 基于lpv模型的满足设定点跟踪与扰动抑制性能的pi控制方法
CN104505830A (zh) * 2015-01-14 2015-04-08 华北电力大学 时滞电力系统稳定性分析方法和装置
CN104698842A (zh) * 2015-02-06 2015-06-10 浙江大学 一种基于内点法的lpv模型非线性预测控制方法
EP3109725A1 (en) * 2015-06-12 2016-12-28 Innovative Solutions and Support, Inc. Upgraded flight management system for autopilot control and method of providing the same
CN105094114A (zh) * 2015-06-23 2015-11-25 北京航天自动控制研究所 一种确定lpv控制器的稳定性的方法
CN104991566A (zh) * 2015-07-07 2015-10-21 北京航天自动控制研究所 一种用于高超声速飞行器的参数不确定性lpv系统建模方法
CN105242540A (zh) * 2015-10-27 2016-01-13 西安石油大学 基于平均驻留时间的网络控制系统动态切换控制方法
CN106094511A (zh) * 2016-06-18 2016-11-09 哈尔滨理工大学 一种时滞lpv系统的鲁棒h∞状态反馈控制器的设计方法
CN106371313A (zh) * 2016-09-22 2017-02-01 哈尔滨理工大学 一种时滞lpv系统有记忆h∞状态反馈控制器设计方法
CN106843261A (zh) * 2016-10-25 2017-06-13 南京航空航天大学 一种变体飞行器过渡段的张量积插值建模与控制方法
CN108628161A (zh) * 2017-03-23 2018-10-09 哈尔滨理工大学 一种切换lpv系统的设计方法
CN107561939A (zh) * 2017-09-01 2018-01-09 哈尔滨理工大学 一种lpv系统输出反馈控制器

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
H∞ control for asynchronously switched linear parameter-varying systems with mode-dependent average dwell time;Qiugang Lu et al.;《IET Control Theory & Applications》;20130321;第7卷(第5期);第673-683页 *
LPV系统切换控制策略研究;郝现志;《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》;20180315(第3期);第1-41页 *
基于模型依赖平均驻留时间的线性切换系统有限时间H∞控制;王通 等;《控制与决策》;20150731;第30卷(第7期);第1189-1194页 *

Also Published As

Publication number Publication date
CN109992004A (zh) 2019-07-09

Similar Documents

Publication Publication Date Title
CN109992004B (zh) 一种lpv系统异步切换状态反馈控制器设计方法
CN109033585B (zh) 不确定网络控制系统的pid控制器设计方法
Shi et al. Quantized learning control for flexible air-breathing hypersonic vehicle with limited actuator bandwidth and prescribed performance
Khalid et al. Conventional PID, adaptive PID, and sliding mode controllers design for aircraft pitch control
Jiao et al. Adaptive mode switching of hypersonic morphing aircraft based on type-2 TSK fuzzy sliding mode control
CN110298144B (zh) 基于交替事件触发的切换网络化飞控系统的输出调节方法
Wu et al. Adaptive neural flight control for an aircraft with time-varying distributed delays
CN111007724A (zh) 一种基于区间ii型模糊神经网络的高超音速飞行器指定性能量化跟踪控制方法
Pashilkar et al. Adaptive back-stepping neural controller for reconfigurable flight control systems
CN108459611B (zh) 一种近空间飞行器的姿态跟踪控制方法
Zhang et al. Fuzzy disturbance observer-based dynamic sliding mode control for hypersonic morphing vehicles
Tran et al. Back-stepping based flight path angle control algorithm for longitudinal dynamics
Tran et al. Integrator-backstepping control design for nonlinear flight system dynamics
CN110109357A (zh) 针对非标准型非线性航空器的半全局自适应控制方法
CN116360255A (zh) 一种非线性参数化高超声速飞行器的自适应调节控制方法
CN113110543B (zh) 一种非线性非最小相位飞行器的鲁棒飞行控制方法
Mobarez et al. Robust PID Flight Controller for Ultrastick-25e UAV
Tan et al. Robust control of air-breathing hypersonic vehicles with adaptive projection-based parameter estimation
Di et al. High performance L 1 adaptive control design for longitudinal dynamics of fixed-wing UAV
Ali et al. Robust level flight control design for scaled Yak-54 unmanned aerial vehicle using single sliding surface
CN106126940B (zh) 机器人分数阶pid控制器稳定性的形式化验证方法
Shen et al. Missile autopilot design based on robust LPV control
Yun et al. Smooth adaptive fixed time convergent controller design for BTT missiles with uncertainties
CN112034867A (zh) 一种无人机控制算法
Yu et al. Missile robust gain scheduling autopilot design using full block multipliers

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