CN113504730A - 一种考虑执行器饱和的非线性飞行器鲁棒控制方法 - Google Patents

一种考虑执行器饱和的非线性飞行器鲁棒控制方法 Download PDF

Info

Publication number
CN113504730A
CN113504730A CN202110887313.3A CN202110887313A CN113504730A CN 113504730 A CN113504730 A CN 113504730A CN 202110887313 A CN202110887313 A CN 202110887313A CN 113504730 A CN113504730 A CN 113504730A
Authority
CN
China
Prior art keywords
state
actuator saturation
nonlinear
designed
considering
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.)
Pending
Application number
CN202110887313.3A
Other languages
English (en)
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.)
Northwestern Polytechnical University
AVIC Chengdu Aircraft Design and Research Institute
Original Assignee
Northwestern Polytechnical University
AVIC Chengdu Aircraft Design and Research Institute
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 Northwestern Polytechnical University, AVIC Chengdu Aircraft Design and Research Institute filed Critical Northwestern Polytechnical University
Priority to CN202110887313.3A priority Critical patent/CN113504730A/zh
Publication of CN113504730A publication Critical patent/CN113504730A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明涉及一种考虑执行器饱和的非线性飞行器鲁棒控制方法,属于飞行控制领域。包括以下步骤:将飞行器模型分解为考虑所有干扰和不确定的线性主系统和考虑执行器饱和的非线性辅系统;设计观测器用以估计辅系统的状态和主系统的输出;针对主系统和辅系统分别设计控制器,设计好主、辅系统的控制器后,综合二者得到原系统的控制器。本发明充分考虑了系统的非线性信息及不确定和干扰的影响,可以取得良好的轨迹跟踪效果;充分考虑了执行器饱和对系统控制性能的影响,可以通过辅系统的补偿实现在主系统中仍采用经典的线性控制方法。本方法简单有效,灵活性和可靠性较高。

Description

一种考虑执行器饱和的非线性飞行器鲁棒控制方法
技术领域
本发明涉及一种飞行器控制方法,特别是涉及一种考虑执行器饱和的飞行器鲁棒跟踪控制方法,属于飞行控制领域。
背景技术
在航空航天领域,受物理特性限制,飞行器执行机构的作动范围有限,当控制器给出的控制信号超出执行器的最大作动能力就会出现执行器饱和的情况,比如在高超声速飞行器中,升降舵偏角的作动范围为δe∈[-15°,15°],油门阀开度的工作范围为Φ∈[0.05,1.2]。执行器饱和是一种强非线性,会影响系统稳定性和动态性能。因此,在控制器设计过程中必须考虑执行器饱和的影响,保障系统的稳定性和工作性能。
处理执行器饱和的一种简单方法是在不考虑饱和的情况下,先按照线性系统设计控制器,再引入补偿环节抵消饱和给系统带来的不利影响。但是这种间接处理方法的处理效果有限,更优的选择是采取直接方法,即考虑饱和非线性,使用非线性控制理论设计非线性控制器。这种方法可以取得更好的控制效果,然而却需要完全摒弃工程上经常应用的经典线性控制器。
发明内容
要解决的技术问题
为了避免现有技术的不足之处,本发明提出一种考虑执行器饱和的飞行器鲁棒跟踪控制方法,该方法能够充分利用工程上成熟的线性控制方法,同时克服饱和非线性对飞行的影响。
技术方案
一种考虑执行器饱和的非线性飞行器鲁棒控制方法,其特征在于步骤如下:
步骤1:将飞行器模型分解为考虑所有干扰和不确定的线性主系统和考虑执行器饱和的非线性辅系统;
步骤2:设计观测器用以估计辅系统的状态和主系统的输出;
步骤3:针对主系统和辅系统分别设计控制器,设计好主、辅系统的控制器后,综合二者得到原系统的控制器。
本发明进一步的技术方案为:步骤1具体为:
其中主系统为
Figure BDA0003194679470000021
yp=Cxp,xp(0)=x0
其中φ(xd)是一个与时间相关,与状态无关的函数,xp,up,yp为主系统的状态、输入和输出;
Figure BDA0003194679470000022
为系统矩阵,
Figure BDA0003194679470000023
为输入矩阵,
Figure BDA0003194679470000024
为输出矩阵,
Figure BDA0003194679470000025
为未知干扰;
定义辅系统的状态、输入和输出xs,us,ys
Figure BDA0003194679470000026
辅系统是原系统与主系统的差
Figure BDA0003194679470000027
ys=Cxs,xs(0)=0
根据加性状态分解,有如下关系
x=xp+xs,u=up+us,y=yp+ys
本发明进一步的技术方案为:步骤2具体为:
设计如下观测器:
Figure BDA0003194679470000028
Figure BDA0003194679470000029
用以估计辅系统的状态和主系统的输出,即xs,yp
本发明进一步的技术方案为:步骤3具体为:
考虑主系统,设计具有扰动抑制能力的主控制器为
up=H(yp-yd)
使得当t→∞时yp(t)-yd(t)→0,其中H(·)是一个线性函数;
考虑辅系统,设计具有处理系统非线性和执行器饱和的辅控制器为
us=L(xs)
使得当t→∞时,xs(t)→0,L(·)是待设计的函数;
最终设计的鲁棒跟踪控制器为
Figure BDA0003194679470000031
Figure BDA0003194679470000032
u=H(yp-yd)+L(xs)
那么,原飞行器系统的状态满足当t→∞时,y(t)-yd(t)→0。
有益效果
本发明提出的一种考虑执行器饱和的飞行器鲁棒控制方法,它是一种在原有线性控制基础上充分考虑了系统饱和非线性的控制方法,在有饱和出现的情况下通过将饱和分离出来单独进行处理,整体上可以获得比原有线性控制方法更好的飞行控制效果。本方法简单有效,灵活性和可靠性较高。
与现有技术相比有益效果为:
(1)本发明充分考虑了系统的非线性信息及不确定和干扰的影响,可以取得良好的轨迹跟踪效果。
(2)本发明充分考虑了执行器饱和对系统控制性能的影响,可以通过辅系统的补偿实现在主系统中仍采用经典的线性控制方法。
(3)本发明通过问题分解降低了控制器的设计难度,两个子系统控制器的设计相比原系统控制器的设计更加简单灵活。
附图说明
附图仅用于示出具体实施例的目的,而并不认为是对本发明的限制,在整个附图中,相同的参考符号表示相同的部件。
图1是本发明一种考虑执行器饱和的飞行器鲁棒跟踪控制方法的流程图
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图和实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。此外,下面描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
本发明针对的考虑执行器饱和的飞行器模型如下:
Figure BDA0003194679470000041
其中,
Figure BDA0003194679470000042
为飞行器系统状态,系统输出为
Figure BDA0003194679470000043
控制输入为
Figure BDA0003194679470000044
为系统矩阵,
Figure BDA0003194679470000045
为输入矩阵,
Figure BDA0003194679470000046
是和x相关的一个非线性函数向量,
Figure BDA0003194679470000047
为输出矩阵,
Figure BDA0003194679470000048
为未知干扰,期望输出表示为yd,饱和函数为sat(·)=[sat(u1),…,sat(ui),…,sat(um)]T,i=1,2,…,m,定义如下
Figure BDA0003194679470000049
针对系统(1),(A,C)可观,系统状态可测。
系统控制目标为:在任意初始条件x0下,当系统存在输入饱和及干扰和不确定d时,实现系统的鲁棒跟踪控制,即当t→∞时,y(t)-yd(t)→0。
本发明提出一种考虑执行器饱和的飞行器鲁棒控制方法,具体如图1所示,本控制方法包括以下三个部分:
(1)面向控制的模型建立
因为系统(1)为包含饱和非线性的飞行器模型,为了能够继续使用经典的线性控制方法,这里将饱和非线性分离出来进行单独处理。为了实现这一目标,首先对考虑执行器饱和的飞行器模型进行模型变换。
将系统(1)分解为考虑所有干扰和不确定的线性主系统和考虑执行器饱和的非线性辅系统。其中主系统为
Figure BDA0003194679470000051
其中φ(xd)是一个与时间相关,与状态无关的函数,xp,up,yp为主系统的状态,输入和输出。因此,主系统(3)是线性的。然后,定义辅系统的状态,输入和输出xs,us,ys
Figure BDA0003194679470000052
辅系统是原系统与主系统的差
Figure BDA0003194679470000053
根据加性状态分解,有如下关系
x=xp+xs,u=up+us,y=yp+ys.
至此,模型变换完成,主系统(3)和辅系统(4)为面向控制的模型。下面介绍总体的控制器设计框架。
(2)总体控制框架建立
主系统(3)是一个多输入多输出的线性系统,有着外部信号φ(yd)。针对主系统,指定控制任务为实现跟踪控制,使得当t→∞时,系统输出yp(t)-yd(t)→0。
辅系统(4)是一个精确的非线性系统,包含了原系统非线性及执行器饱和。当主系统跟踪上了目标,那么辅系统的平衡点为原点。针对辅系统,指定控制任务为镇定飞行器系统的非线性及饱和,使得当t→∞时,系统状态xs(t)→0。
由于主系统和辅系统是设计模型,而非真实存在的模型,因此需要设计观测器来获取这两个系统的状态值。设计如下观测器
Figure BDA0003194679470000061
用以估计辅系统的状态和主系统的输出,即xs,yp
基于面向控制的模型,原非线性系统的输出跟踪问题分解为线性主系统的输出跟踪问题和非线性辅系统的镇定控制问题。主系统实现轨迹跟踪控制,即控制目标为当t→∞时,yp(t)-yd(t)→0。主系统跟踪问题就变成了输出反馈跟踪问题。辅系统只需要解决存在执行器饱和情况下的镇定问题,控制目标为当t→∞时,xs(t)→0。镇定问题比原系统的跟踪问题简单,因为辅系统的状态已知且不需要解决跟踪问题。当yp(t)-yd(t)→0且xs(t)→0时,有y(t)-yd(t)→0。因此,总体的控制框架为:观测器提供辅系统状态和主系统输出的观测值,在此基础上,主系统控制器和辅系统控制器共同作用实现存在执行器饱和的飞行器的鲁棒飞行控制。下面介绍主系统及辅系统的控制器设计。
(3)控制器设计
针对主系统和辅系统分别设计控制器,设计好主、辅系统的控制器后,可以综合二者得到原系统的控制器。
考虑主系统,设计具有扰动抑制能力的主控制器为
up=H(yp-yd) (6)
使得当t→∞时yp(t)-yd(t)→0,其中H(·)是一个线性函数。
考虑辅系统,设计具有处理系统非线性和执行器饱和的辅控制器为
us=L(xs) (7)
使得当t→∞时,xs(t)→0,L(·)是待设计的函数。
只要两个子问题可以被很好的解决,那么原问题也就被解决了。最终设计的原系统的鲁棒跟踪控制器为
Figure BDA0003194679470000071
那么,原飞行器系统的状态满足当t→∞时,y(t)-yd(t)→0。
以高超声速飞行器为例,进行考虑执行器饱和的飞行器鲁棒跟踪控制器的设计。
马赫数为15,高度为110000ft的巡航段高超声速飞行器纵向通道Winged-Cone模型具体如下:
Figure BDA0003194679470000072
此模型包括5个状态变量x=[V,h,α,γ,q]T和两个控制输入u=[δe,β]T。其中,V是速度,γ是航迹角,h是高度,α是攻角,q是俯仰角速率,δe是升降舵舵偏角,β为节流阀开度。T,D,L,Myy分别代表推力,阻力,升力和俯仰转动力矩。m,Iyy,μ,r分别代表飞行器的质量,俯仰转动惯量,重力常数和地心的径向距离。
考虑控制输入的执行器饱和为β∈[0.05,1.2],δe∈[-15°,15°],设计非线性鲁棒跟踪控制器保证在执行器饱和影响下高超声速飞行器依然可以取得很好的跟踪性能。
鲁棒跟踪控制方法具体步骤如下:
(1)面向控制的模型建立
高超声速飞行器模型可以写成如下的形式
Figure BDA0003194679470000073
其中Φ(·,·)为表示模型(9)的一个综合函数,上式通过模型变换可以写成式(1)的形式。
考虑平衡点xtrim,utrim,基于小角度线性化并且考虑常规线性化忽略的非线性得到
Figure BDA0003194679470000081
其中
Figure BDA0003194679470000082
取输出y=[V,h]T
然后通过系统分解得到如下主系统和辅系统
主系统:
Figure BDA0003194679470000083
辅系统:
Figure BDA0003194679470000084
(2)总体控制框架建立
设计如下观测器用于观测两个子系统状态
Figure BDA0003194679470000085
主系统实现轨迹跟踪控制,即控制目标为当t→∞时,yp(t)-yd(t)→0。主系统跟踪问题就变成了输出反馈跟踪问题。辅系统只需要解决镇定问题,控制目标为当t→∞时,xs(t)→0。当yp(t)-yd(t)→0且xs(t)→0时有y(t)-yd(t)→0。因此,总体的控制框架为:观测器提供辅系统状态和主系统输出的观测值,再此基础上,主系统和辅系统共同作用实现存在执行器饱和的飞行器的鲁棒飞行控制。
(3)控制器设计
考虑鲁棒跟踪控制器的设计。为主系统(12)设计线性二次积分器。一般为了消除跟踪误差,控制器中会考虑积分项
Figure BDA0003194679470000086
则可得如下扩维系统
Figure BDA0003194679470000087
设计状态反馈控制器
up=-Kxxp-Keq (17)
其中Kx和Ke为反馈矩阵,则可保证系统的鲁棒跟踪目标。这里可采用线性二次调节器方法来确定反馈矩阵Kx和Ke
为辅系统(13)设计状态反馈控制器。设计控制输入为
us=lTxs (18)其中
Figure BDA0003194679470000091
是一个为控制器参数,则可保证辅系统的渐进稳定。
最终的控制器由主系统控制器、辅系统控制器和观测器组合而成
Figure BDA0003194679470000092
本发明未详细说明部分属于领域技术人员公知常识。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明公开的技术范围内,可轻易想到各种等效的修改或替换,这些修改或替换都应涵盖在本发明的保护范围之内。

Claims (4)

1.一种考虑执行器饱和的非线性飞行器鲁棒控制方法,其特征在于步骤如下:
步骤1:将飞行器模型分解为考虑所有干扰和不确定的线性主系统和考虑执行器饱和的非线性辅系统;
步骤2:设计观测器用以估计辅系统的状态和主系统的输出;
步骤3:针对主系统和辅系统分别设计控制器,设计好主、辅系统的控制器后,综合二者得到原系统的控制器。
2.根据权利要求1所述一种考虑执行器饱和的非线性飞行器鲁棒控制方法,其特征在于步骤1具体为:
其中主系统为
Figure FDA0003194679460000011
yp=Cxp,xp(0)=x0
其中φ(xd)是一个与时间相关,与状态无关的函数,xp,up,yp为主系统的状态、输入和输出;
Figure FDA0003194679460000012
为系统矩阵,
Figure FDA0003194679460000013
为输入矩阵,
Figure FDA0003194679460000014
为输出矩阵,
Figure FDA0003194679460000015
为未知干扰;
定义辅系统的状态、输入和输出xs,us,ys
Figure FDA0003194679460000016
辅系统是原系统与主系统的差
Figure FDA0003194679460000017
ys=Cxs,xs(0)=0
根据加性状态分解,有如下关系
x=xp+xs,u=up+us,y=yp+ys
3.根据权利要求2所述一种考虑执行器饱和的非线性飞行器鲁棒控制方法,其特征在于步骤2具体为:
设计如下观测器:
Figure FDA0003194679460000021
Figure FDA0003194679460000022
用以估计辅系统的状态和主系统的输出,即xs,yp
4.根据权利要求3所述一种考虑执行器饱和的非线性飞行器鲁棒控制方法,其特征在于步骤3具体为:
考虑主系统,设计具有扰动抑制能力的主控制器为up=H(yp-yd)
使得当t→∞时yp(t)-yd(t)→0,其中H(·)是一个线性函数;
考虑辅系统,设计具有处理系统非线性和执行器饱和的辅控制器为us=L(xs)
使得当t→∞时,xs(t)→0,L(·)是待设计的函数;
最终设计的鲁棒跟踪控制器为
Figure FDA0003194679460000023
Figure FDA0003194679460000024
u=H(yp-yd)+L(xs)
那么,原飞行器系统的状态满足当t→∞时,y(t)-yd(t)→0。
CN202110887313.3A 2021-08-03 2021-08-03 一种考虑执行器饱和的非线性飞行器鲁棒控制方法 Pending CN113504730A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110887313.3A CN113504730A (zh) 2021-08-03 2021-08-03 一种考虑执行器饱和的非线性飞行器鲁棒控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110887313.3A CN113504730A (zh) 2021-08-03 2021-08-03 一种考虑执行器饱和的非线性飞行器鲁棒控制方法

Publications (1)

Publication Number Publication Date
CN113504730A true CN113504730A (zh) 2021-10-15

Family

ID=78015443

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110887313.3A Pending CN113504730A (zh) 2021-08-03 2021-08-03 一种考虑执行器饱和的非线性飞行器鲁棒控制方法

Country Status (1)

Country Link
CN (1) CN113504730A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115051600A (zh) * 2022-07-18 2022-09-13 湖南科技大学 一种无刷直流电机伺服系统跟踪控制方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107203138A (zh) * 2017-06-27 2017-09-26 金陵科技学院 一种输入输出饱和的飞行器鲁棒控制方法
CN109062043A (zh) * 2018-08-01 2018-12-21 西北工业大学 考虑网络传输以及执行器饱和的航天器自抗扰控制方法
CN109143866A (zh) * 2018-09-25 2019-01-04 浙江工业大学 一种考虑执行器受限问题的刚性飞行器自适应固定时间姿态跟踪控制方法
CN113110543A (zh) * 2021-04-19 2021-07-13 西北工业大学 一种非线性非最小相位飞行器的鲁棒飞行控制方法
CN113126497A (zh) * 2021-04-14 2021-07-16 西北工业大学 考虑输入饱和的飞行器鲁棒跟踪控制方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107203138A (zh) * 2017-06-27 2017-09-26 金陵科技学院 一种输入输出饱和的飞行器鲁棒控制方法
CN109062043A (zh) * 2018-08-01 2018-12-21 西北工业大学 考虑网络传输以及执行器饱和的航天器自抗扰控制方法
CN109143866A (zh) * 2018-09-25 2019-01-04 浙江工业大学 一种考虑执行器受限问题的刚性飞行器自适应固定时间姿态跟踪控制方法
CN113126497A (zh) * 2021-04-14 2021-07-16 西北工业大学 考虑输入饱和的飞行器鲁棒跟踪控制方法
CN113110543A (zh) * 2021-04-19 2021-07-13 西北工业大学 一种非线性非最小相位飞行器的鲁棒飞行控制方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
REN JINRUI,ET AL.: "Docking control for probe-drogue refueling: An additive-state-decomposition-based output feedback iterative learning control method", 《CHINESE JOURNAL OF AERONAUTICS》 *
许斌 等: "基于时标分解的弹性高超声速飞行器智能控制", 《航空学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115051600A (zh) * 2022-07-18 2022-09-13 湖南科技大学 一种无刷直流电机伺服系统跟踪控制方法
CN115051600B (zh) * 2022-07-18 2024-07-12 湖南科技大学 一种无刷直流电机伺服系统跟踪控制方法

Similar Documents

Publication Publication Date Title
CN109541941B (zh) 一种针对垂直起降运载器主动段飞行的自适应增广抗扰容错方法
Hu et al. Adaptive backstepping control for air-breathing hypersonic vehicles with input nonlinearities
CN109460050B (zh) 一种针对变体无人机的复合分层抗干扰控制方法
CN111290421A (zh) 一种考虑输入饱和的高超声速飞行器姿态控制方法
CN108427289B (zh) 一种基于非线性函数的高超声速飞行器跟踪控制方法
Jin et al. Robust adaptive neural network-based compensation control of a class of quadrotor aircrafts
Li et al. Predefined-time asymptotic tracking control for hypersonic flight vehicles with input quantization and faults
CN111367182A (zh) 考虑输入受限的高超声速飞行器抗干扰反步控制方法
Van Oort et al. Full-envelope modular adaptive control of a fighter aircraft using orthogonal least squares
Ding et al. Global smooth sliding mode controller for flexible air-breathing hypersonic vehicle with actuator faults
CN111158398A (zh) 考虑攻角约束的高超音速飞行器的自适应控制方法
Cordeiro et al. Robustness of incremental backstepping flight controllers: The boeing 747 case study
CN113504730A (zh) 一种考虑执行器饱和的非线性飞行器鲁棒控制方法
CN107943097B (zh) 飞行器的控制方法、装置和飞行器
CN113110543B (zh) 一种非线性非最小相位飞行器的鲁棒飞行控制方法
CN116360255A (zh) 一种非线性参数化高超声速飞行器的自适应调节控制方法
Zhao et al. Global adaptive neural backstepping control of a flexible hypersonic vehicle with disturbance estimation
Moin et al. State space model of an aircraft using Simulink
Lewis et al. Limited authority adaptive control architectures with dynamic inversion or explicit model following
Schumacher et al. PI control of a tailless fighter aircraft with dynamic inversion and neural networks
CN114637318A (zh) 一种高超声速飞行器滑模控制方法
Sobolic et al. Aerodynamic-free adaptive control of the NASA generic transport model
CN111061283A (zh) 一种基于特征模型的吸气式高超声速飞行器高度控制方法
CN113110581B (zh) 一种基于主辅系统结合的非线性飞行器位置保持控制方法
Lu et al. An ESO-based attitude control method for non-affine hypersonic flight vehicles

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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20211015