CN114625159B - 一种基于被控变量的结冰飞机控制方法 - Google Patents
一种基于被控变量的结冰飞机控制方法 Download PDFInfo
- Publication number
- CN114625159B CN114625159B CN202210069645.5A CN202210069645A CN114625159B CN 114625159 B CN114625159 B CN 114625159B CN 202210069645 A CN202210069645 A CN 202210069645A CN 114625159 B CN114625159 B CN 114625159B
- Authority
- CN
- China
- Prior art keywords
- aircraft
- sensor
- control
- acceleration
- angular acceleration
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 230000001133 acceleration Effects 0.000 claims abstract description 59
- 230000000694 effects Effects 0.000 claims abstract description 4
- 238000005259 measurement Methods 0.000 claims description 13
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 claims description 6
- 230000005484 gravity Effects 0.000 claims description 3
- 238000009434 installation Methods 0.000 claims description 2
- 230000008859 change Effects 0.000 abstract description 18
- 238000009825 accumulation Methods 0.000 abstract description 11
- 230000004069 differentiation Effects 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 7
- 230000004044 response Effects 0.000 description 5
- 238000004088 simulation Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 238000005096 rolling process Methods 0.000 description 2
- 230000000087 stabilizing effect Effects 0.000 description 2
- 241001481833 Coryphaena hippurus Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/10—Simultaneous control of position or course in three dimensions
- G05D1/101—Simultaneous control of position or course in three dimensions specially adapted for aircraft
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Vehicle Body Suspensions (AREA)
Abstract
本说明书实施例公开了一种基于被控变量的结冰飞机控制方法。本发明通过线加速度和角加速度反馈的增量动态逆方法设计结冰飞机的控制律,飞机气动参数的改变都反映在线加速度和角加速度的变化里,当机身和舵面出现积冰引起飞机气动特性变化时,飞机闭环飞行品质不受影响。此外由于微分器对高频噪声有放大效应,本文没有直接通过角速度微分得到角加速度,而是采用在机身配置低成本线加速度计的方式,利用各线加速度计之间差值获取飞机角加速度。
Description
技术领域
本说明书涉及飞行控制领域,尤其涉及一种基于被控变量的结冰飞机控制方法。
背景技术
飞机在含有大量过冷水滴的大气中飞行时,表面会聚集冰层,冰层主要集中在机翼、操纵面和进气道唇口等对气动特性影响显著的部位。结冰导致的飞机气动特性改变会使飞行品质下降,更为严重的可能导致飞机失控。
飞机结冰对飞机气动特性的影响主要体现在机身及平尾失速特性、飞行性能、稳定性和操纵性。实际飞行过程中结冰速度和积冰几何形状受飞行环境和气动部件外形综合影响,目前并不能做出很好的预测,尤其是积冰几何形状基本无法预测,很难精确判定积冰会对气动部件的外流场产生何种影响。
因此,需要提供一种适用于结冰状态的飞机的飞行控制方案。
发明内容
本发明的目的在于,提供一种适用于结冰状态的飞机的飞行控制方案。
为解决上述技术问题,本发明采用如下技术方案:
一种基于被控变量的结冰飞机控制方法,包括:
确定飞机在三个控制通道上的控制变量:
其中,LCV为横向控制变量,ps为稳定轴系下的滚转角速度,p为体轴系下的滚转角速度,r为体轴系下的偏航角速度,α为攻角;MCV为纵向控制变量,q为体轴系下的俯仰角速度,Q为动压,S为飞机参考面积,CLα为升力线斜率,m为飞机质量,Vco为俯仰角速度和法向加速度在纵向控制变量中比重相同时的速度;NCV为航向控制变量,β为侧滑角;
根据所述被控变量,结合增量动态逆原理,推导出对应上述被控变量的增量动态逆控制方法,实现通航飞机结冰状态的鲁棒控制,控制律中需要用到的角加速度信号通过安装于机身不同位置的线加速度计测量值解算得到。
本说明书实施例采用的上述至少一个技术方案能够达到以下有益效果:
与现有技术相比,本发明通过线加速度和角加速度反馈的增量动态逆方法设计结冰飞机的控制律,飞机气动参数的改变都反映在线加速度和角加速度的变化里,当机身和舵面出现积冰引起飞机气动特性变化时,飞机闭环飞行品质不受影响。此外由于微分器对高频噪声有放大效应,本文没有直接通过角速度微分得到角加速度,而是采用在机身顶点配置低成本线加速度计的方式,利用各线加速度计之间差值解算出飞机角加速度。
附图说明
图1为本申请实施例所涉及的飞机外表的结冰形状示意图;
图2为本申请实施例中所提供的传感器的位置示意图;
图3a为本申请实施例所提供的纵向上的控制响应逻辑示意图,图3b为本申请实施例所提供的在横向上的控制响应逻辑示意图;
图4a至图4c为本申请实施例所提供的指令响应的三通道的仿真结果的示意图;
图5a至图5c为本申请实施例所提供的结冰状态与常规状态下本申请的方案和常规方案的三通道的对比示意图。
具体实施方式
为使本申请的目的、技术方案和优点更加清楚,下面将结合本申请具体实施例及相应的附图对本申请技术方案进行清楚、完整地描述。显然,所描述的实施例仅是本申请一部分实施例,而不是全部的实施例。基于本说明书中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
飞行过程中结冰速度和积冰几何形状受飞行环境和气动部件外形综合影响,目前并不能做出很好的预测,尤其是积冰几何形状基本无法预测,图1列出几种常见形状,从图中积冰形状和积冰位置看,积冰主要发生在气动部件前缘,会破环气动部件的外流场,不同冰型对流线的破环程度有较大的差异,其中角状冰对流场的破坏尤其严重,其导致的气动参数改变量比较大,有可能诱发飞机“海豚跳”,但很难精确判定积冰会对气动部件的外流场产生何种影响。
气动部件积冰导致飞机静稳定导数、舵面铰链力矩和操纵导数改变,会改变飞机飞行品质,进一部会导致严重的飞行问安全问题,因此需要克服飞机积冰带来的影响。
基于此,本申请实施例提供一种基于被控变量的结冰飞机控制方法,具体包括:
步骤S201,确定飞机的角加速度
在本申请实施例中,需要使用状态变化率作为反馈信号,但是状态变化率往往不容易直接测量,以飞机动力学***为例,状态变化率中的线加速度可通过加速度计直接测量,直接敏感角加速度的传感器由于量程和带宽受限,因此都采用间接方法来测量角加速度。
常见的角加速度间接测量方法本质上都是基于微分运算的,而微分运算会放大测量噪声,降低角加速度信号的信噪比,因此本文采用基于线加速度计和陀螺仪相结合的方法来获得飞机在质心上角加速度。
安装在机身上某处的线加速度可以下式表示:
其中下标“cg”表示质心处的值,依据两个不同位置(i,j)线加速度相减可以得到角加速度满足如下关系:
写成标量形式为
其中
显然式(3)左边矩阵的秩为2,因此一对线加速度传感器不足以直接计算出三轴角加速度,必须采用三个或三个以上的线加速度传感器布置在机身上,才可以计算出唯一的三轴角加速度值。考虑到一般运输类飞机机身实际几何形状,轴向和展向有较大的空间来布置线加速度传感器,参照式(2),较大的相对距离可产生较大的线加速度差量,可降低线加速度计测量噪声对最终结果的影响。四个传感器的具体安装位置如图2中的所示,分别位于机翼两端、机头和机尾,其中①号和③号(即第一传感器和第三传感器)位于飞机对称面内,且z坐标相同,②号和④号(即第二传感器和第四传感器)的x坐标和z坐标。
参照式(2)和(3),再结合上面传感器安装位置,可以计算出飞机的三轴角加速度。对①号和③号传感器而言,式(2)和(3)变为:
对②号和④号传感器则有:
偏航角加速度既可以采用①号和③号传感器计算得到,也可以通过②号和④号传感器得到,若将两者综合,则可以消除角速度测量误差对角加速度测量值的影响:
为了让线加速度测量值能够计算出角加速度,因此线加速度传感器必须分散布置非质心位置,因此测得线加速度必须进行角速率修正,采用(1)式描述飞机质心处的线加速度acg。
S203,根据所述在质心处的线加速度acg,采用如下方式确定机在三个不同控制通道上的控制变量。
在设计飞机的飞行控制律时需要对纵向、航向和横向三个通道进行综合控制,具体控制变量为如下形式:
其中,LCV为横向控制变量,ps为稳定轴系下的滚转角速度,p为滚转角速度,r为偏航角速度,α为攻角;MCV为纵向控制变量,q为俯仰角速度,Q为动压,S为飞机参考面积,CLα为升力线斜率,m为飞机质量,Vco为俯仰角速度和法向加速度在纵向控制变量中比重相同时的速度,来源于C*准则,在升力线性的情况下,上述纵向控制变量与C*准则等价,取值一般为120m/s左右,低速时偏重于俯仰角速度控制,高速时偏重于迎角控制(实际等同于过载);NCV为航向控制通道,β为侧滑角。
S205,根据所述控制变量输入相应指令,生成舵面指令以实现飞机的飞行控制。
在本申请实施例中,对于飞机的纵向和横航向分开设计。如图3a和图3b所示,图3a为本申请实施例所提供的纵向控制律示意图,图3b为本申请实施例所提供的在横航向上的控制律示意图。
对纵向控制而言,状态取为x=[α,q],因此纵向控制中需要用到迎角变化率,显然迎角和侧滑角的变化率并不能直接测量,只能用微分方式得到,但是气流角的测量精度较低,同时由于频带受限导致数据采样率偏低,因此通过气流角测量值微分得到的迎角变化率不能直接用于上面的纵向控制律当中。本文采用角速度和线加速度测量信号代入到迎角动力学方程中,直接计算得到迎角变化率,飞机迎角的动力学方程可写为如下形式:
若引入不扣除重力加速度测量信号(原始信号),上式可以变为:
显然由上式可以通过线加速度计、气流角传感器和三轴陀螺仪的测量值直接得到,不用引入跟踪微分器。增量动态逆方法中用到的俯仰角加速度则是通过第3节中的方案得到。
对横航向控制而言,状态取为:
x=[p,r]
因此横航向增量动态逆控制方法需要用到滚转角加速度和偏航角加速度,同样可通过第3节中的方案得到,这里不再详述。稳定轴系下偏航角速度rs作为侧滑角控制的内环,稳定轴系下偏航角速度指令rsc可用下式表示:
rsc=-KNCV(NCV-β)+fβ (11)
其中对于fβ同样可以采用传感器的实测信号代替实际的气动模型:
参照控制品质需求和时标分离思想各控制参数取为:
KMCV=10rad/s;KLCV=10rad/s;KNCV=2rad/s。
与现有技术相比,本发明通过线加速度和角加速度反馈的增量动态逆方法设计结冰飞机的控制律,飞机气动参数的改变都反映在线加速度和角加速度的变化里,当机身和舵面出现积冰引起飞机气动特性变化时,飞机闭环飞行品质不受影响。此外由于微分器对高频噪声有放大效应,本文没有直接通过角速度微分得到角加速度,而是采用在机身顶点配置低成本线加速度计的方式,利用各线加速度计之间差值获取飞机质心处的角加速度。
进一步的,本申请实施例还提供了相应的如下验证:
试验对象为一较大展弦比的通航飞机,机翼面积为24m2,展长为19.8m,平均气动弦长为2m,机身长度为15m,质量为3000kg,转动惯量参数为:
Ixx=21562kg·m2;Iyy=14231kg·m2;Izz=34256kg·m2;Ixz=239kg·m2
伺服***采用一阶惯性***建模,带宽为ωa=50rad/s。分别对三通道进行跟踪仿真,飞行速度为60m/s,飞行高度1000m,指令形式为偶极方波,仿真时不考虑模型偏差和测量误差。
三通道的指令响应特性可用如下传递函数来描述:
指令响应的仿真结果如图4a至图4c所示,可见在不考虑模型偏差和测量误差的情况下,三通道的被控变量都能较好的跟踪上预期指令。其中横向控制在实现滚转指令跟踪的同时,也能较好抑制侧滑角,在航向上,能实现带侧滑飞行,不仅能抑制侧风扰动,同时能让飞机在着陆阶段对准跑道。
飞机结冰时的气动参数会发生变化。如表1所示。
表1飞机结冰前后气动参数
结冰后气动参数会发生较大的变化,以对内环姿态控制最主要的几个参数幅值变化为例:Cmα减小17.6%,减小13.5%,Cnβ减少了27%,Clβ减少10%,/>减小26%,/>减小33.3%。图5为三通道结冰状态与标称状态仿真结果,同时给出了常规控制方法与本申请的方法进行比较,可以看出,对于本申请的方法,由于结冰状态气动参数与标称状态有较大差异,各通道舵面偏转与标称状态不同,在结冰状态下则能更好的跟踪被控变量。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于装置、设备和介质类实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可,这里就不再一一赘述。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下,在权利要求书中记载的动作或步骤或模块可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外,在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中,多任务处理和并行处理也是可以的或者可能是有利的。
Claims (5)
1.一种基于被控变量的结冰飞机控制方法,所述方法包括:
确定飞机在三个控制通道上的控制变量:
其中,LCV为横向控制变量,ps为稳定轴系下的滚转角速度,p为体轴系下的滚转角速度,r为体轴系下的偏航角速度,α为攻角;MCV为纵向控制变量,q为体轴系下的俯仰角速度,Q为动压,S为飞机参考面积,CLα为升力线斜率,m为飞机质量,Vco为俯仰角速度和法向加速度在纵向控制变量中比重相同时的速度;NCV为航向控制变量,β为侧滑角;
根据所述控制变量输入相应指令,生成舵面指令以实现飞机的飞行控制。
2.如权利要求1所述的方法,需要确定飞机角加速度包括:
根据设置于非质心位置上的多个线加速度传感器的测量值ai,采用如下方式确定飞机角加速度
3.如权利要求2所述的方法,需要对上述多个线加速度传感器的安装位置进行布置,包括:
设置于机翼两端、机头和机尾共4个传感器,其中,传感器1和传感器3位于飞机对称面内,且z方向坐标相同,传感器2和传感器4的x方向坐标和z方向坐标相同。
4.如前面要求3所述的方法,具体确定通航飞机三轴角加速度的具体值,包括:
采用传感器1和传感器3所感测得到的线加速度分别确定和/>
或者,采用传感器2和传感器4感测得到的线加速度分别确定和/>
5.如权利要求4所述的方法,其中,还包括:
综合传感器1和传感器3所确定的结果和传感器2和传感器4所确定的结果,可消除角速度测量误差的影响:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210069645.5A CN114625159B (zh) | 2022-01-21 | 2022-01-21 | 一种基于被控变量的结冰飞机控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210069645.5A CN114625159B (zh) | 2022-01-21 | 2022-01-21 | 一种基于被控变量的结冰飞机控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114625159A CN114625159A (zh) | 2022-06-14 |
CN114625159B true CN114625159B (zh) | 2023-07-28 |
Family
ID=81898775
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210069645.5A Active CN114625159B (zh) | 2022-01-21 | 2022-01-21 | 一种基于被控变量的结冰飞机控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114625159B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110316358A (zh) * | 2019-03-29 | 2019-10-11 | 南京航空航天大学 | 基于动态逆的战斗机大迎角控制方法 |
CN113778120A (zh) * | 2021-10-27 | 2021-12-10 | 北京航空航天大学 | 一种多传感器融合的无人机复杂天候飞行控制方法 |
Family Cites Families (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA839624A (en) * | 1970-04-21 | C. Amos Homer | Forming inorganic filaments from melt emanating from rotor | |
US5170969A (en) * | 1988-11-23 | 1992-12-15 | The Boeing Company | Aircraft rudder command system |
US9443358B2 (en) * | 1995-06-07 | 2016-09-13 | Automotive Vehicular Sciences LLC | Vehicle software upgrade techniques |
RU2364548C2 (ru) * | 2007-05-17 | 2009-08-20 | Закрытое акционерное общество "Гражданские самолеты Сухого" | Система управления летательным аппаратом |
ITPN20120028A1 (it) * | 2012-06-01 | 2013-12-02 | Giovanni Amos Blandina | Dispositivo per evitare il fenomeno di aquaplaning in un veicolo |
FR3035978A1 (fr) * | 2015-05-04 | 2016-11-11 | Airbus Helicopters | Systeme de commande de giravion, giravion associe et methode de commande correspondante |
GB2554062A (en) * | 2016-08-22 | 2018-03-28 | Norwegian Univ Of Science And Technology | Icing control system |
CN106707759B (zh) * | 2017-02-17 | 2019-04-09 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种飞机Herbst机动控制方法 |
US10160537B2 (en) * | 2017-03-13 | 2018-12-25 | The Boeing Company | Apparatus for using aircraft active vibration control system as pilot cueing aid |
US11987368B2 (en) * | 2017-06-29 | 2024-05-21 | Bombarder Inc. | Aircraft operation during steep approach |
WO2019157455A1 (en) * | 2018-02-09 | 2019-08-15 | Skydio, Inc. | Aerial vehicle smart landing |
US11492923B2 (en) * | 2018-04-09 | 2022-11-08 | Gulfstream Aerospace Corporation | Ice shedding aircraft engine |
KR101982071B1 (ko) * | 2018-08-24 | 2019-05-24 | 주식회사 에그 (Egg) | 자율 결빙방지장치의 제어방법 |
CN114967737A (zh) * | 2019-07-12 | 2022-08-30 | 深圳市道通智能航空技术股份有限公司 | 一种飞行器控制方法及飞行器 |
CN112046761B (zh) * | 2020-08-04 | 2021-10-15 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种基于统计检验和滤波的飞机结冰在线探测方法 |
CN113625732B (zh) * | 2021-07-31 | 2024-07-05 | 西北工业大学 | 基于角加速度估计的增量反馈逆角速度控制律设计方法 |
CN113777931B (zh) * | 2021-11-09 | 2022-02-11 | 中国空气动力研究与发展中心计算空气动力研究所 | 结冰翼型气动模型构造方法、装置、设备及介质 |
-
2022
- 2022-01-21 CN CN202210069645.5A patent/CN114625159B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110316358A (zh) * | 2019-03-29 | 2019-10-11 | 南京航空航天大学 | 基于动态逆的战斗机大迎角控制方法 |
CN113778120A (zh) * | 2021-10-27 | 2021-12-10 | 北京航空航天大学 | 一种多传感器融合的无人机复杂天候飞行控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114625159A (zh) | 2022-06-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109614633B (zh) | 一种复合式旋翼飞行器非线性建模及线性化配平方法 | |
US8155794B2 (en) | System and method for control of aeroelasticity effects | |
Ryan et al. | Current and future research in active control of lightweight, flexible structures using the X-56 aircraft | |
US7757993B1 (en) | Method for reducing the turbulence and gust influences on the flying characteristics of aircraft, and a control device for this purpose | |
CN109703768B (zh) | 一种基于姿态/轨迹复合控制的软式空中加油对接方法 | |
CN106043695B (zh) | 一种油动多旋翼无人机定桨距变转速***及控制技术 | |
CN111290278B (zh) | 一种基于预测滑模的高超声速飞行器鲁棒姿态控制方法 | |
WO2020088094A1 (zh) | 一种旋翼机协调转弯控制方法及*** | |
CN107101636A (zh) | 一种使用卡尔曼滤波器辨识多旋翼动力学模型参数的方法 | |
CN111580540A (zh) | 基于Kalman滤波的变体飞行器T-S模糊控制方法 | |
Li et al. | Simulation method for wind tunnel based virtual flight testing | |
CN114625159B (zh) | 一种基于被控变量的结冰飞机控制方法 | |
CN109445283B (zh) | 一种用于欠驱动浮空器在平面上定点跟踪的控制方法 | |
CN113885545B (zh) | 考虑缆绳张力变化的系留无人机控制方法 | |
CN111859540B (zh) | 一种大气扰动中飞机颠簸响应的计算方法 | |
Kotikalpudi et al. | Real-time shape estimation for flexible unmanned air vehicle via Kalman filtering | |
CN116611175B (zh) | 一种大展弦比飞机体自由度颤振的预测方法 | |
RU2790358C1 (ru) | Способ определения аэродинамических характеристик летательного аппарата по результатам летного эксперимента | |
CN115828420B (zh) | 一种吸气式高超声速无人机动态可靠性分析方法 | |
Serokhvostov et al. | Mav unsteady characteristics in-flight measurement with the help of smartap autopilot | |
ZEIS, JR et al. | Angle of attack estimation using an inertial reference platform | |
Park et al. | Determination of flight attitude using kalman filter for an agricultural unmanned helicopter | |
Bennett et al. | The aircraft spin-a mathematical approach and comparison to flight test | |
CN116280182A (zh) | 民用飞机防滑刹车控制方法 | |
Nguyen et al. | Unsteady aerodynamic parameter estimation for multirotor helicopters |
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 |