CN102520726B - 大攻角飞行状态下的大气攻角及侧滑角估计方法 - Google Patents
大攻角飞行状态下的大气攻角及侧滑角估计方法 Download PDFInfo
- Publication number
- CN102520726B CN102520726B CN 201110426439 CN201110426439A CN102520726B CN 102520726 B CN102520726 B CN 102520726B CN 201110426439 CN201110426439 CN 201110426439 CN 201110426439 A CN201110426439 A CN 201110426439A CN 102520726 B CN102520726 B CN 102520726B
- Authority
- CN
- China
- Prior art keywords
- cos
- sin
- beta
- angle
- alpha
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 55
- 230000001133 acceleration Effects 0.000 claims abstract description 46
- 238000005259 measurement Methods 0.000 claims abstract description 46
- 230000008859 change Effects 0.000 claims abstract description 35
- 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 abstract description 12
- 239000000446 fuel Substances 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 36
- 238000004364 calculation method Methods 0.000 claims description 16
- 238000004422 calculation algorithm Methods 0.000 claims description 15
- 239000010763 heavy fuel oil Substances 0.000 claims description 11
- 238000001914 filtration Methods 0.000 claims description 7
- 230000006872 improvement Effects 0.000 claims description 6
- 239000000203 mixture Substances 0.000 claims description 6
- 238000005096 rolling process Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 230000005484 gravity Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 6
- 230000000694 effects Effects 0.000 description 10
- 238000010586 diagram Methods 0.000 description 9
- 239000000295 fuel oil Substances 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 239000000443 aerosol Substances 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000010729 system oil Substances 0.000 description 1
Images
Landscapes
- Navigation (AREA)
Abstract
本发明公开了一种大攻角飞行状态下的大气攻角及侧滑角估计方法,在不改变飞行器结构外形、不增添额外测量装置和硬件设备的前提下,充分利用现有机载捷联惯性导航***、GPS***、飞行控制***以及燃油测量***输出的参数,解算飞行器的角加速度、转动惯量、惯性积、力矩方程组的系数、总质量,基于飞行动力学模型,通过构建扩展卡尔曼滤波器,融飞行动力学模型求解和状态估计的过程为一体,实现当前时刻攻角、侧滑角和真空速的精确估计,进而解算出当前时刻的风速及风速变化率,并作为反馈信息输入到扩展卡尔曼滤波器中,完成下一时刻攻角及侧滑角的实时精确估计,通过递推求解的方式实现攻角和侧滑角的实时精确估计,提高了大气数据***的测量范围和测量精度以及对复杂飞行环境的适应性。
Description
技术领域
本发明涉及一种大攻角飞行状态下的大气攻角及侧滑角估计方法,属于大气数据、飞行器控制、惯性导航及GPS***等技术领域。
背景技术
攻角和侧滑角是飞行器飞行状态和飞行力学的两个重要参数,其精度和可靠性直接关系到飞行器控制、导航、发动机及火控等***的正常工作及性能发挥。现代高性能战机的主要战术性能要求包括超机动能力、超音速巡航、隐身性、高信息优势等。高机动飞行意味着飞行器具备大攻角飞行能力,大攻角飞行引起机体表面气流发生分离,飞行器处于非线性漩涡气流中,目前的攻角和侧滑角主要是通过安装在飞行器上的风标传感器、压差式传感器和零压差式传感器等来测量,这些类型的攻角和侧滑角传感器对气流扰动比较敏感,导致攻角和侧滑角传感器性能严重下降,难以满足大攻角飞行范围内攻角和侧滑角的精确测量,严重制约了高性能飞行器的性能发挥。另外,突出的传感器测量装置势必会降低飞行器的隐身性,因此解决高性能飞行器攻角和侧滑角的测量问题显得尤为重要。
目前国内外都已开展了大攻角下攻角和侧滑角的测量方法研究。1991年美国NASA报告刊登的一篇“Development of a Pneumatic High-Angle-of-Attack Flush Airdata Sensing(HI-FADS)System”论文,文中提出了一种依靠嵌入在飞行器前端不同位置上的压力传感器阵列来测量飞行器表面的压力分布,并由压力分布推算大气参数的嵌入式大气数据***,该方法能够实现大攻角飞行环境下的攻角和侧滑角的测量;2009年美国专利US7,564,539B2:“Optical Air Data Systems Methods”,提出了一种采用激光与大气分子和气溶胶分别发生瑞利散射和米氏散射效应,根据后向散射信号的强度和多普勒频移解算大气参数的光学式大气数据***,该方法在测量大气参数时不受飞行状态影响,满足大攻角飞行状态下攻角和侧滑角的测量;以上两种采用研究新型的嵌入式大气数据***和光学式大气数据***的方法技术难度大且代价较高;2006年AIAA刊登了会议“AIAA Guidance,Navigation,and Control Conference and Exhibit”上的一篇论文“Flight Testing of the X-45AJ-UCAS Computational Alpha-Beta System”,提供了一种利用惯性导航***、飞行控制***以及大气数据***进行攻角和侧滑角估计方法,以此避开了大攻角特性对飞行攻角、侧滑角测量的影响,但是该方法由于没有考虑飞行器周围实际风速的影响,估计的攻角和侧滑角存在***误差。
发明内容
本发明目的是针对现有大攻角大气参数测量技术存在的不足,充分利用机载捷联惯性导航***、GPS、飞行控制***以及燃油测量***的输出参数,基于飞行动力学模型,提出一种大攻角飞行状态下的大气攻角及侧滑角估计方法,该方法具有不受飞行环境和飞行状态制约,经济有效且易于工程实现的特点。
为实现上述目的,本发明采用如下技术方案:
一种大攻角飞行状态下的大气攻角及侧滑角估计方法,其特征在于:利用飞行器机载设备输出的参数,包括:
(1)捷联惯性导航***输出的姿态角、角速度以及线加速度;
(2)GPS***输出的地速;
(3)飞行控制***输出的空气动力、外合力矩、发动机推力;
(4)燃油测量***输出的剩余燃油质量;
由角速度通过微分方法求解出角加速度;
由角速度、角加速度和外合力矩,通过最小二乘估计方法求解出飞行器的转动惯量和惯性积,进而得到力矩方程组的系数;
由剩余燃油质量和已知的飞行器机体质量、机载设备质量、乘员质量以及武器质量相加,得到飞行器总质量;
根据飞行动力学模型,构建扩展卡尔曼滤波器,利用扩展卡尔曼滤波算法,实现大攻角飞行状态下当前时刻攻角、侧滑角和真空速的精确估计;
由当前时刻的姿态角、地速以及估计出的当前时刻的攻角、侧滑角、真空速,对当前时刻的风速进行解算,并结合上一时刻的风速,对当前时刻的风速变化率进行解算;
将得到的当前时刻的风速和风速变化率反馈给扩展卡尔曼滤波模算法模块,用于完成下一时刻的攻角及侧滑角的实时精确估计,通过递推求解的方式实现攻角和侧滑角的实时精确估计。
包括以下步骤:
(1)以周期ΔT读取捷联惯性导航***中的三个角速度、三个姿态角和三个线加速度信息,三个角速度信息分别为机体坐标系下x轴、y轴和z轴方向的横滚角速度p、俯仰角速度q、偏航角速度r;三个姿态角信息分别为俯仰角θ、横滚角φ、偏航角ψ;三个线加速度信息分别为机体坐标系下x轴方向的线加速度ax、y轴方向的线加速度ay、z轴方向的线加速度az,其中机体坐标系的x轴、y轴和z轴的指向分别为向前、向右、向下;
(2)以周期ΔT读取GPS***输出的地理坐标系下的地速,地理坐标系方向定义为北、东、地,此三个方向的地速分别为ug、vg、wg;
(3)以周期ΔT读取飞行控制***输出的飞行器所受空气动力、外合力矩、发动机推力,其中飞行器所受空气动力在气流坐标系下xa轴负方向、ya轴方向和za轴负方向的分量分别为阻力D、侧力Y、升力L,发动机推力在机体坐标系下x轴、y轴和z轴方向的分量分别为Tx、Ty、Tz,外合力矩在机体坐标系下x轴、y轴和z轴方向的分量分别为滚转力矩俯仰力矩M、偏航力矩N;
(4)以周期ΔT读取飞行器燃油测量***输出的剩余燃油质量mf;
(5)根据步骤(1)获取的机体坐标系下的横滚角速度p、俯仰角速度q、偏航角速度r,其中t时刻三个角速度信息分别为横滚角速度p(t)、俯仰角速度q(t)、偏航角速度r(t),t+ΔT时刻三个角速度信息分别为横滚角速度p(t+ΔT)、俯仰角速度q(t+ΔT)、偏航角速度r(t+ΔT),t时刻的机体坐标系下x轴、y轴和z轴方向的三个角加速度信息,横滚角加速度俯仰角加速度偏航角加速度按微分方法计算,即:
(6)将已知的飞行器角运动方程组 的三个方程式相加,得到方程EJ=U,其中 J=[Ix,Iy,Iz,Ixz]T′、T′表示转置,Ix、Iy、Iz为飞行器分别绕机体坐标系的x轴、y轴和z轴方向的转动惯量,飞行器对x轴、y轴的惯性积为Ixy,对x轴、z轴的惯性积为Ixz,对y轴、z轴的惯性积为Iyz,飞行器关于机体坐标系的对称面Oxz对称,惯性积Ixy=Iyz=0,Ixz为非零;从ta时刻开始到tb时刻结束的时间段内,利用步骤(1)连续获取n组横滚角速度p、俯仰角速度q、偏航角速度r,利用步骤(3)连续获取n组滚转力矩俯仰力矩M、偏航力矩N,利用步骤(5)连续求取n组横滚角加速度俯仰角加速度偏航角加速度根据得到的n组横滚角速度p、俯仰角速度q、偏航角速度r、横滚角加速度俯仰角加速度偏航角加速度滚转力矩俯仰力矩M、偏航力矩N和方程EJ=U,利用最小二乘估计方法,求出飞行器的转动惯量Ix、Iy、Iz和惯性积Ixz,其中n>4;
(7)根据步骤(6)得到的转动惯量Ix、Iy、Iz和惯性积Ixz,对已知的飞行器力矩方程组 中的9个系数c1、c2、c3、c4、c5、c6、c7、c8、c9进行求解,其中 其中
(8)将步骤(4)得到的剩余燃油质量mf与已知的飞行器机体质量、机载设备质量、乘员质量以及武器质量相加,得到飞行器总质量m;
(9)根据飞行动力学模型,选取飞行器的俯仰角、横滚角、横滚角速度、俯仰角速度、偏航角速度、攻角、侧滑角、真空速作为状态量,进而建立扩展卡尔曼滤波器状态方程;选取俯仰角、横滚角、横滚角速度、俯仰角速度、偏航角速度、三维线加速度以及真空速作为量测量,进而建立扩展卡尔曼滤波器量测方程;根据步骤(1)和步骤(2)获取的当前时刻即tk+1时刻的量测信息,步骤(3)获取的上一时刻即tk时刻的空气动力、外合力矩、发动机推力,步骤(7)求取的9个力矩方程组系数,步骤(8)得到的tk时刻的飞行器总质量,以及tk时刻的三维风速和三维风速变化率,利用扩展卡尔曼滤波方程得到tk+1时刻状态量的最优估计值,从而实现大攻角飞行状态下tk+1时刻的攻角及侧滑角的实时精确估计,其中tk>tb;
(10)根据步骤(1)得到tk+1时刻的俯仰角、横滚角、偏航角和步骤(2)得到的tk+1时刻的三维地速参数,结合步骤(9)估计得到的tk+1时刻的攻角、侧滑角和真空速,计算出tk+1时刻的三维风速,并结合tk时刻的三维风速,计算出tk+1时刻的三维风速变化率;
(11)将步骤(10)得到的tk+1时刻的三维风速和三维风速变化率反馈给扩展卡尔曼滤波模算法模块,用于完成步骤(9)中的下一时刻即tk+2时刻攻角及侧滑角的估计。
所述步骤(9)中卡尔曼滤波算法的具体步骤是:
(a)扩展卡尔曼滤波器状态方程的建立
根据飞行动力学模型,选取飞行器的俯仰角θ、横滚角φ、横滚角速度p、俯仰角速度q、偏航角速度r、攻角α、侧滑角β、真空速V作为状态量X,即X=[θ,φ,p,q,r,α,β,V]T′,进而建立扩展卡尔曼滤波器的状态方程:
状态方程的具体形式为:
其中,u(t)为控制输入量,设机体坐标系下沿x轴、y轴和z轴方向的三维风速分别为uw、vw、ww,三维风速uw、vw、ww的变化率分别为***噪声 分别由飞行器的外合力矩M、N、空气动力D、Y、L、发动机推力Tx、Ty、Tz、风速uw、vw、ww、风速变化率以及总质量m引起,g是重力加速度,***噪声系数阵G(t)定义为:
其中,
(b)扩展卡尔曼滤波器量测方程的建立
根据GPS得到的三维地速ug、vg、wg,利用解算出的V作为真空速量测量,并结合捷联惯性导航***得到的俯仰角θ、横滚角φ、横滚角速度p、俯仰角速度q、偏航角速度r、三维线加速度ax、ay、az组成总的量测量Z,即Z=[θ,φ,p,q,r,ax,ay,az,V]T′,进而建立扩展卡尔曼滤波器量测方程:
Z(t)=h[x(t),u(t),t]+v(t)
量测方程具体可表示为:
其中,量测噪声 由惯性器件的陀螺仪和线加速度计以及GPS计算的真空速的量测噪声引起;
(c)状态方程和量测方程转换为连续型线性干扰方程:
利用泰勒级数展开的方法将步骤(a)中的状态方程和步骤(b)中的量测方程转换为连续型线性干扰方程形式,即将f[X(t),u(t),t]和h[X(t),u(t),t]在状态最优估计值处进行泰勒展开,并略去二次以上项,得到连续型线性干扰方程:
式中,
(d)连续型线性干扰方程离散化:
取采样周期ΔT=tk+1-tk对步骤(c)中得到的连续型线性干扰方程进行离散化得到离散型线性干扰方程:
(e)扩展卡尔曼滤波方程:
根据步骤(3)获取的tk时刻的空气动力D、Y、L、外合力矩M、N、发动机推力Tx、Ty、Tz,步骤(7)求取的力矩方程组系数c1、c2、c3、c4、c5、c6、c7、c8、c9,步骤(8)计算的tk时刻的飞行器总质量m,以及tk时刻的三维风速uw、vw、ww和三维风速变化率参数,并结合tk时刻的状态最优估计值利用求取tk+1时刻状态量的一步预测值
根据步骤(2)获取tk+1时刻的三维地速ug、vg、wg,利用得到tk+1时刻真空速的量测量,并结合步骤(1)得到的tk+1时刻俯仰角θ、横滚角φ、横滚角速度p、俯仰角速度q、偏航角速度r、三维线加速度ax、ay、az参数组成tk+1时刻总的量测量Zk+1,结合求解得到的tk+1时刻状态量的一步预测值和滤波增益矩阵Kk+1,利用公式 进行tk+1时刻状态量改善值的求解;
扩展卡尔曼滤波器中,状态量中的俯仰角θ、横滚角φ、横滚角速度p、俯仰角速度q和偏航角速度r的初始值由***内部的捷联惯性导航***提供,状态量中攻角α、侧滑角β、真空速V的初始值由***外部的大气数据***提供,从而构成初始状态初始状态估值误差方差阵P0,***噪声方差阵初值Q0,量测噪声方差阵初值R0矩阵初值都由***外部直接输入,***噪声矩阵Qk由***统噪声方差阵初值Q0确定,量测噪声矩阵Rk+1由量测噪声方差阵初值R0确定,另外,三维风速uw、vw、ww的初始值利用***外部的风速仪提供,三维风速变化率的初始值均取为零。
所述步骤(10)中风速及风速变化率解算的具体方法是:
根据步骤(1)得到的tk+1时刻的俯仰角θ,横滚角φ和偏航角ψ,得到机体tk+1时刻的姿态转换矩阵Sθψφ:
结合步骤(2)获取的tk+1时刻的三维地速ug、vg、wg,得到tk+1时刻机体坐标系下沿x轴、y轴和z轴方向的三维地速ue,ve,we为:
结合步骤(9)估计得到的tk+1时刻的攻角α、侧滑角β、真空速V参数,利用下式将真空速V转换为机体坐标系下沿x轴、y轴和z轴方向的三维真空速ub、vb、wb,即
那么tk+1时刻机体坐标系下沿x轴、y轴和z轴方向的三维风速uw、vw、ww利用下式进行计算:
利用计算出的tk+1时刻的三维风速uw(tk+1)、vw(tk+1)、ww(tk+1),结合tk时刻的三维风速uw(tk)、vw(tk)、ww(tk),利用微分方法求取tk+1时刻三维风速变化率 即:
本发明的优点及显著效果:本发明在不改变飞行器结构外形、不增添额外测量装置和硬件设备的前提下,充分利用现有的机载惯性导航***、GPS、飞行控制***以及燃油测量***的输出参数,解算飞行器角加速度、转动惯量、惯性积、力矩方程组的系数和总质量,基于飞行动力学模型,通过构建扩展卡尔曼滤波器,融飞行动力学模型求解和状态估计的过程为一体,实现当前时刻攻角、侧滑角和真空速的精确估计,进而解算出当前时刻的风速及风速变化率,并作为反馈信息输入到扩展卡尔曼滤波器中,完成下一时刻攻角及侧滑角的实时精确估计,通过递推求解的方式实现攻角和侧滑角的实时精确估计。本方法不仅有利于飞行器隐身,同时提高了大气数据***的测量范围和测量精度以及对复杂飞行环境的适应性,是一种相对经济且易于工程实现的方法。对于提高我国新一代航空飞行器大攻角飞行状态下攻角和侧滑角的测量范围、改善飞行器飞行性能具有重要的现实意义。
附图说明
图1是本发明方法的原理流程示意图;
图2是力矩方程组系数计算模块的原理框图;
图3是扩展卡尔曼滤波算法模块的原理框图;
图4是风速及风速变化率计算模块的原理框图;
图5是扩展卡尔曼滤波算法模块与风速及风速变化率计算模块之间信息交互原理图;
图6是本发明攻角估计效果曲线示意图;
图7是本发明侧滑角估计效果曲线示意图;
图8是本发明攻角估计误差曲线示意图;
图9是本发明侧滑角估计误差曲线示意图。
具体实施方式
参看图1-5,本发明利用现有的机载惯性导航***、GPS、飞行控制***以及燃油测量***的输出参数,解算飞行器角加速度、转动惯量、惯性积、力矩方程组的系数和总质量,基于飞行动力学模型,通过构建扩展卡尔曼滤波器,融飞行动力学模型求解和状态估计的过程为一体,实现当前时刻攻角、侧滑角和真空速的精确估计,进而解算出当前时刻的风速及风速变化率,并作为反馈信息输入到扩展卡尔曼滤波器中,完成下一时刻攻角及侧滑角的实时精确估计,通过递推求解的方式实现攻角和侧滑角的实时精确估计。
为实现大攻角飞行状态下攻角和侧滑角的估计,需要完成的工作:
(1)捷联惯性导航***数据获取:
以周期ΔT读取捷联惯性导航***中的三个角速度、三个姿态角和三个线加速度信息,三个角速度信息分别为机体坐标系下x轴、y轴和z轴方向的横滚角速度p、俯仰角速度q、偏航角速度r;三个姿态角信息分别为俯仰角θ、横滚角φ、偏航角ψ;三个线加速度信息分别为机体坐标系下x轴方向的线加速度ax、y轴方向的线加速度ay、z轴方向的线加速度az,其中机体坐标系的x轴、y轴和z轴的指向分别为向前、向右、向下;
(2)GPS***数据获取:
以周期ΔT读取GPS***输出的地理坐标系下的地速,地理坐标系方向定义为北、东、地,此三个方向的地速分别为ug、vg、wg;
(3)飞行控制***数据获取:
以周期ΔT读取飞行控制***输出的飞行器所受空气动力、外合力矩、发动机推力,其中飞行器所受空气动力在气流坐标系下xa轴负方向、ya轴方向和za负方向的分量分别为阻力D、侧力Y、升力L,发动机推力在机体坐标系下x轴、y轴和z轴方向的分量分别为Tx、Ty、Tz,外合力矩在机体坐标系下x轴、y轴和z轴方向的分量分别为滚转力矩俯仰力矩M、偏航力矩N;
(4)飞行器剩余燃油质量获取:
以周期ΔT读取飞行器燃油测量***输出的剩余燃油质量mf;
(5)角加速度解算:
根据步骤(1)获取的机体坐标系下的横滚角速度p、俯仰角速度q、偏航角速度r,其中t时刻三个角速度信息分别为横滚角速度p(t)、俯仰角速度q(t)、偏航角速度r(t),t+ΔT时刻三个角速度信息分别为横滚角速度p(t+ΔT)、俯仰角速度q(t+ΔT)、偏航角速度r(t+ΔT),t时刻的机体坐标系下x轴、y轴和z轴方向的三个角加速度信息,横滚角加速度俯仰角加速度偏航角加速度按微分方法计算,即:
(6)飞行器转动惯量和惯性积解算:
将已知的飞行器角运动方程组 的三个方程式相加,得到方程EJ=U,其中 J=[Ix,Iy,Iz,Ixz]T′、T′表示转置,Ix、Iy、Iz为飞行器分别绕机体坐标系的x轴、y轴和z轴方向的转动惯量,飞行器对x轴、y轴的惯性积为Ixy,对x轴、z轴的惯性积为Ixz,对y轴、z轴的惯性积为Iyz,飞行器关于机体坐标系的对称面Oxz对称,惯性积Ixy=Iyz=0,Ixz为非零;从ta时刻开始到tb时刻结束的时间段内,利用步骤(1)连续获取n组横滚角速度p、俯仰角速度q、偏航角速度r,利用步骤(3)连续获取n组滚转力矩俯仰力矩M、偏航力矩N,利用步骤(5)连续求取n组横滚角加速度俯仰角加速度偏航角加速度根据得到的n组横滚角速度p、俯仰角速度q、偏航角速度r、横滚角加速度俯仰角加速度偏航角加速度滚转力矩俯仰力矩M、偏航力矩N和方程EJ=U,利用最小二乘估计方法,求出飞行器的转动惯量Ix、Iy、Iz和惯性积Ixz,其中n>4;
(7)力矩方程组系数解算:
根据步骤(6)得到的转动惯量Ix、Iy、Iz和惯性积Ixz,对已知的飞行器力矩方程组 中的9个系数c1、c2、c3、c4、c5、c6、c7、c8、c9进行求解,其中 其中
(8)飞行器总质量解算:
将步骤(4)得到的剩余燃油质量mf与已知的飞行器机体质量、机载设备质量、乘员质量以及武器质量相加,得到飞行器总质量m
(9)扩展卡尔曼滤波算法:
根据飞行动力学模型,选取飞行器的俯仰角、横滚角、横滚角速度、俯仰角速度、偏航角速度、攻角、侧滑角、真空速作为状态量,进而建立扩展卡尔曼滤波器状态方程;选取俯仰角、横滚角、横滚角速度、俯仰角速度、偏航角速度、三维线加速度以及真空速作为量测量,进而建立扩展卡尔曼滤波器量测方程;根据步骤(1)和步骤(2)获取的当前时刻即tk+1时刻的量测信息,步骤(3)获取的上一时刻即tk时刻的空气动力、外合力矩、发动机推力,步骤(7)求取的9个力矩方程组系数,步骤(8)得到的tk时刻的飞行器总质量,以及tk时刻的三维风速和三维风速变化率,利用扩展卡尔曼滤波方程得到tk+1时刻状态量的最优估计值,从而实现大攻角飞行状态下tk+1时刻的攻角及侧滑角的实时精确估计,其中tk>tb;
(a)扩展卡尔曼滤波器状态方程的建立
根据飞行动力学模型,选取飞行器的俯仰角θ、横滚角φ、横滚角速度p、俯仰角速度q、偏航角速度r、攻角α、侧滑角β、真空速V作为状态量X,即X=[θ,φ,p,q,r,α,β,V]T′,进而建立扩展卡尔曼滤波器的状态方程:
状态方程的具体形式为:
其中,u(t)为控制输入量,设机体坐标系下沿x轴、y轴和z轴方向的三维风速分别为uw、vw、ww,三维风速uw、vw、ww的变化率分别为***噪声 分别由飞行器的外合力矩M、N、空气动力D、Y、L、发动机推力Tx、Ty、Tz、风速uw、vw、ww、风速变化率以及总质量m引起,g是重力加速度,***噪声系数阵G(t)定义为:
其中,
(b)扩展卡尔曼滤波器量测方程的建立
根据GPS得到的三维地速ug、vg、wg,利用解算出的V作为真空速量测量,并结合捷联惯性导航***得到的俯仰角θ、横滚角φ、横滚角速度p、俯仰角速度q、偏航角速度r、三维线加速度ax、ay、az组成总的量测量Z,即Z=[θ,φ,p,q,r,ax,ay,az,V]T′,进而建立扩展卡尔曼滤波器量测方程:
Z(t)=h[X(t),u(t),t]+v(t)
量测方程具体可表示为:
其中,量测噪声 由惯性器件的陀螺仪和线加速度计以及GPS计算的真空速的量测噪声引起;
(c)状态方程和量测方程转换为连续型线性干扰方程:
利用泰勒级数展开的方法将步骤(a)中的状态方程和步骤(b)中的量测方程转换为连续型线性干扰方程形式,即将f[X(t),u(t),t]和h[X(t),u(t),t]在状态最优估计值处进行泰勒展开,并略去二次以上项,得到连续型线性干扰方程:
式中,
(d)连续型线性干扰方程离散化:
取采样周期ΔT=tk+1-tk对步骤(c)中得到的连续型线性干扰方程进行离散化得到离散型线性干扰方程:
(e)扩展卡尔曼滤波方程:
根据步骤(3)获取的tk时刻的空气动力D、Y、L、外合力矩M、N、发动机推力Tx、Ty、Tz,步骤(7)求取的力矩方程组系数c1、c2、c3、c4、c5、c6、c7、c8、c9,步骤(8)计算的tk时刻的飞行器总质量m,以及tk时刻的三维风速uw、vw、ww和三维风速变化率参数,并结合tk时刻的状态最优估计值利用求取tk+1时刻状态量的一步预测值
根据步骤(2)获取tk+1时刻的三维地速ug、vg、wg,利用得到tk+1时刻真空速的量测量,并结合步骤(1)得到的tk+1时刻俯仰角θ、横滚角φ、横滚角速度p、俯仰角速度q、偏航角速度r、三维线加速度ax、ay、az参数组成tk+1时刻总的量测量Zk+1,结合求解得到的tk+1时刻状态量的一步预测值和滤波增益矩阵Kk+1,利用公式 进行tk+1时刻状态量改善值的求解;
扩展卡尔曼滤波器中,状态量中的俯仰角θ、横滚角φ、横滚角速度p、俯仰角速度q和偏航角速度r的初始值由***内部的捷联惯性导航***提供,状态量中攻角α、侧滑角β、真空速V的初始值由***外部的大气数据***提供,从而构成初始状态初始状态估值误差方差阵P0,***噪声方差阵初值Q0,量测噪声方差阵初值R0矩阵初值都由***外部直接输入,***噪声矩阵Qk由***统噪声方差阵初值Q0确定,量测噪声矩阵Rk+1由量测噪声方差阵初值R0确定,另外,三维风速uw、vw、ww的初始值利用***外部的风速仪提供,三维风速变化率的初始值均取为零。
(10)风速及风速变化率解算:
根据步骤(1)得到tk+1时刻的俯仰角、横滚角、偏航角和步骤(2)得到的tk+1时刻的三维地速参数,结合步骤(9)估计得到的tk+1时刻的攻角、侧滑角和真空速,计算出tk+1时刻的三维风速,并结合tk时刻的三维风速,计算出tk+1时刻的三维风速变化率;
根据步骤(1)得到的tk+1时刻的俯仰角θ,横滚角φ和偏航角ψ,得到机体tk+1时刻的姿态转换矩阵Sθψφ:
结合步骤(2)获取的tk+1时刻的三维地速ug、vg、wg,得到tk+1时刻机体坐标系下沿x轴、y轴和z轴方向的三维地速ue,ve,we为:
结合步骤(9)估计得到的tk+1时刻的攻角α、侧滑角β、真空速V参数,利用下式将真空速V转换为机体坐标系下沿x轴、y轴和z轴方向的三维真空速ub、vb、wb,即
那么tk+1时刻机体坐标系下沿x轴、y轴和z轴方向的三维风速uw、vw、ww利用下式进行计算:
利用计算出的tk+1时刻的三维风速uw(tk+1)、vw(tk+1)、ww(tk+1),结合tk时刻的三维风速uw(tk)、vw(tk)、ww(tk),利用微分方法求取tk+1时刻三维风速变化率 即:
(11)风速和风速变化率反馈步骤:
将步骤(10)得到的tk+1时刻的三维风速和三维风速变化率反馈给扩展卡尔曼滤波模算法模块,用于完成步骤(9)中的下一时刻即tk+2时刻攻角及侧滑角的估计。
因此,通过递推求解的方式可以实现攻角和侧滑角的实时精确估计。
图1是本发明方法的基本原理框图,攻角/侧滑角计算***利用捷联惯性导航***、GPS***、飞行控制***和燃油测量***提供的参数,通过力矩方程组系数计算模块、风速及风速变化率计算模块、总质量计算模块以及扩展卡尔曼滤波算法模块,最终实现攻角和侧滑角的精确估计,其中各个模块之间的实线箭头代表了基本的逻辑连接关系。
图2是对图1中的力矩方程组系数计算模块的进一步细化,代表了力矩方程组系数计算的具体流程,更清晰的表明了信息的来源和传递过程。
图3是对图1中的扩展卡尔曼滤波算法模块的进一步细化,代表了攻角、侧滑角及真空速计算的具体流程,更清晰的表明了信息的来源和扩展卡尔曼滤波器的主要构成。
图4是对图1中的风速及风速变化率计算模块的进一步细化,更清晰的表明了当前时刻风速及风速变化率的计算流程和信息的来源。
图5是对图1中的扩展卡尔曼滤波算法模块和风速及风速变化率计算模块之间的信息交互的进一步细化,更清晰的表明了两个模块之间信息的在不同时刻相互传递和利用的过程,其中为状态量的初始值,t1为扩展卡尔曼滤波的开始时刻,是t1时刻的状态量的估计值,t2是t1的下一时刻,是t2时刻的状态量的估计值,是tk+1时刻的状态量的估计值,是tk+2时刻的状态量的估计值。
图6是本发明方法攻角估计效果曲线,横轴代表时间,纵轴表示攻角,飞行攻角变化范围-41°~75°,利用该发明估计的攻角与真实攻角曲线基本重合,吻合度较好,应用该发明不仅能够实现常规飞行攻角下时对攻角的精确估计,而且对大攻角飞行具有较好的跟踪效果。
图7是本发明方法侧滑角估计效果曲线,横轴代表时间,纵轴表示侧滑角,侧滑角变化范围为-12°~22°,利用该发明估计的侧滑角与真实侧滑角曲线基本重合,应用该发明能够对大攻角飞行状态下的侧滑角进行精确估计。
图8是本发明方法攻角估计误差曲线,横轴代表时间,纵轴代表攻角估计误差,应用该发明所估计的攻角误差均值为0.31°,估计攻角误差的标准差为0.20°,攻角估计最大误差控制在1.2°以内。
图9是本发明方法侧滑角估计误差曲线,横轴代表时间,纵轴代表侧滑角的估计误差,应用该发明所估计的侧滑角误差均值为0.05°,估计侧滑角误差的标准差为0.33°,侧滑角估计误差控制在1.2°以内。
具体实施例
本实施例基于飞行模拟软件X-Plane仿真的飞行数据,对大攻角飞行状态下的攻角和侧滑角估计方法进行实验,得出了有益的结论。
基于X-Plane***的仿真试验效果:从攻角和侧滑角估计效果曲线图6和图7及估计误差曲线图8和图9可以看出,本发明不仅能够在常规飞行时实现攻角和侧滑角的精确估计,而且能够对大攻角飞行具有较好的跟踪效果,其中攻角和侧滑角估计误差均值分别为0.31°和0.05°,估计误差的标准差分别为0.20°和0.33°,尽管在大攻角飞行阶段估计的攻角和侧滑角误差略有增大,但是所估计的攻角和侧滑角误差均能控制在1.2°以内。本发明具有很强的工程应用价值。
Claims (4)
1.一种大攻角飞行状态下的大气攻角及侧滑角估计方法,其特征在于:利用飞行器机载设备输出的参数,包括:
①捷联惯性导航***输出的姿态角、角速度以及线加速度;
②GPS***输出的地速;
③飞行控制***输出的空气动力、外合力矩、发动机推力;
④燃油测量***输出的剩余燃油质量;
由角速度通过微分方法求解出角加速度;
由角速度、角加速度和外合力矩,通过最小二乘估计方法求解出飞行器的转动惯量和惯性积,进而得到力矩方程组的系数;
由剩余燃油质量和已知的飞行器机体质量、机载设备质量、乘员质量以及武器质量相加,得到飞行器总质量;
根据飞行动力学模型,构建扩展卡尔曼滤波器,利用扩展卡尔曼滤波算法,实现大攻角飞行状态下当前时刻攻角、侧滑角和真空速的精确估计;
由当前时刻的姿态角、地速以及估计出的当前时刻的攻角、侧滑角、真空速,对当前时刻的风速进行解算,并结合上一时刻的风速,对当前时刻的风速变化率进行解算;
将得到的当前时刻的风速和风速变化率反馈给扩展卡尔曼滤波模算法模块,用于完成下一时刻的攻角及侧滑角的实时精确估计,通过递推求解的方式实现攻角和侧滑角的实时精确估计。
2.根据权利要求1所述的大攻角飞行状态下的大气攻角及侧滑角估计方法,包括以下步骤:
(1)以周期ΔT读取捷联惯性导航***中的三个角速度、三个姿态角和三个线加速度信息,三个角速度信息分别为机体坐标系下x轴、y轴和z轴方向的横滚角速度p、俯仰角速度q、偏航角速度r;三个姿态角信息分别为俯仰角θ、横滚角φ、偏航角ψ;三个线加速度信息分别为机体坐标系下x轴方向的线加速度ax、y轴方向的线加速度ay、z轴方向的线加速度az,其中机体坐标系的x轴、y轴和z轴的指向分别为向前、向右、向下;
(2)以周期ΔT读取GPS***输出的地理坐标系下的地速,地理坐标系方向定义为北、东、地,此三个方向的地速分别为ug、vg、wg;
(3)以周期ΔT读取飞行控制***输出的飞行器所受空气动力、外合力矩、发动机推力,其中飞行器所受空气动力在气流坐标系下xa轴负方向、ya轴方向和za轴负方向的分量分别为阻力D、侧力Y、升力L,发动机推力在机体坐标系下x轴、y轴和z轴方向的分量分别为Tx、Ty、Tz,外合力矩在机体坐标系下x轴、y轴和z轴方向的分量分别为滚转力矩俯仰力矩M、偏航力矩N;
(4)以周期ΔT读取飞行器燃油测量***输出的剩余燃油质量mf;
(5)根据步骤(1)获取的机体坐标系下的横滚角速度p、俯仰角速度q、偏航角速度r,其中t时刻三个角速度信息分别为横滚角速度p(t)、俯仰角速度q(t)、偏航角速度r(t),t+ΔT时刻三个角速度信息分别为横滚角速度p(t+ΔT)、俯仰角速度q(t+ΔT)、偏航角速度r(t+ΔT),t时刻的机体坐标系下x轴、y轴和z轴方向的三个角加速度信息,即横滚角加速度俯仰角加速度偏航角加速度按微分方法计算,即:
(6)将已知的飞行器角运动方程组 的三个方程式相加,得到方程EJ=U,其中 J=[Ix,Iy,Iz,Ixz]T'、T'表示转置,Ix、Iy、Iz为飞行器分别绕机体坐标系的x轴、y轴和z轴方向的转动惯量,飞行器对x轴、y轴的惯性积为Ixy,对x轴、z轴的惯性积为Ixz,对y轴、z轴的惯性积为Iyz,飞行器关于机体坐标系的对称面Oxz对称,惯性积Ixy=Iyz=0,Ixz为非零;从ta时刻开始到tb时刻结束的时间段内,利用步骤(1)连续获取n组横滚角速度p、俯仰角速度q、偏航角速度r,利用步骤(3)连续获取n组滚转力矩俯仰力矩M、偏航力矩N,利用步骤(5)连续求取n组横滚角加速度俯仰角加速度偏航角加速度根据得到的n组横滚角速度p、俯仰角速度q、偏航角速度r、横滚角加速度俯仰角加速度偏航角加速度滚转力矩俯仰力矩M、偏航力矩N和方程EJ=U,利用最小二乘估计方法,求出飞行器的转动惯量Ix、Iy、Iz和惯性积Ixz,其中n>4;
(7)根据步骤(6)得到的转动惯量Ix、Iy、Iz和惯性积Ixz,对已知的飞行器力矩方程组 中的9个系数c1、c2、c3、c4、c5、c6、c7、c8、c9进行求解,其中 其中
(8)将步骤(4)得到的剩余燃油质量mf与已知的飞行器机体质量、机载设备质量、乘员质量以及武器质量相加,得到飞行器总质量m;
(9)根据飞行动力学模型,选取飞行器的俯仰角、横滚角、横滚角速度、俯仰角速度、偏航角速度、攻角、侧滑角、真空速作为状态量,进而建立扩展卡尔曼滤波器状态方程;选取俯仰角、横滚角、横滚角速度、俯仰角速度、偏航角速度、三维线加速度以及真空速作为量测量,进而建立扩展卡尔曼滤波器量测方程;根据步骤(1)和步骤(2)获取的当前时刻即tk+1时刻的量测信息,步骤(3)获取的上一时刻即tk时刻的空气动力、外合力矩、发动机推力,步骤(7)求取的9个力矩方程组系数,步骤(8)得到的tk时刻的飞行器总质量,以及tk时刻的三维风速和三维风速变化率,利用扩展卡尔曼滤波方程得到tk+1时刻状态量的最优估计值,从而实现大攻角飞行状态下tk+1时刻的攻角及侧滑角的实时精确估计,其中tk>tb;
(10)根据步骤(1)得到tk+1时刻的俯仰角、横滚角、偏航角和步骤(2)得到的tk+1时刻的三维地速参数,结合步骤(9)估计得到的tk+1时刻的攻角、侧滑角和真空速,计算出tk+1时刻的三维风速,并结合tk时刻的三维风速,计算出tk+1时刻的三维风速变化率;
(11)将步骤(10)得到的tk+1时刻的三维风速和三维风速变化率反馈给扩展卡尔曼滤波模算法模块,用于完成步骤(9)中的下一时刻即tk+2时刻攻角及侧滑角的估计。
3.根据权利要求2所述的大攻角飞行状态下的大气攻角及侧滑角估计方法,其特征在于:所述步骤(9)中扩展卡尔曼滤波算法的具体步骤是:
(a)扩展卡尔曼滤波器状态方程的建立
根据飞行动力学模型,选取飞行器的俯仰角θ、横滚角φ、横滚角速度p、俯仰角速度q、偏航角速度r、攻角α、侧滑角β、真空速V作为状态量X,即X=[θ,φ,p,q,r,α,β,V]T',进而建立扩展卡尔曼滤波器的状态方程:
状态方程的具体形式为:
其中,u(t)为控制输入量,设机体坐标系下沿x轴、y轴和z轴方向的三维风速分别为uw、vw、ww,三维风速uw、vw、ww的变化率分别为***噪声 分别由飞行器的外合力矩M、N、空气动力D、Y、L、发动机推力Tx、Ty、Tz、风速uw、vw、ww、风速变化率以及总质量m引起,g是重力加速度,***噪声系数阵G(t)定义为:
其中,
(b)扩展卡尔曼滤波器量测方程的建立
根据GPS得到的三维地速ug、vg、wg,利用解算出的V作为真空速量测量,并结合捷联惯性导航***得到的俯仰角θ、横滚角φ、横滚角速度p、俯仰角速度q、偏航角速度r、三维线加速度ax、ay、az组成总的量测量Z,即Z=[θ,φ,p,q,r,ax,ay,az,V]T',进而建立扩展卡尔曼滤波器量测方程:
Z(t)=h[X(t),u(t),t]+v(t)
量测方程具体可表示为:
(c)状态方程和量测方程转换为连续型线性干扰方程:
利用泰勒级数展开的方法将步骤(a)中的状态方程和步骤(b)中的量测方程转换为连续型线性干扰方程形式,即将f[X(t),u(t),t]和h[X(t),u(t),t]在状态最优估计值处进行泰勒展开,并略去二次以上项,得到连续型线性干扰方程:
式中,
(d)连续型线性干扰方程离散化:
取采样周期ΔT=tk+1-tk对步骤(c)中得到的连续型线性干扰方程进行离散化得到离散型线性干扰方程:
(e)扩展卡尔曼滤波方程:
根据步骤(3)获取的tk时刻的空气动力D、Y、L、外合力矩M、N、发动机推力Tx、Ty、Tz,步骤(7)求取的力矩方程组系数c1、c2、c3、c4、c5、c6、c7、c8、c9,步骤(8)计算的tk时刻的飞行器总质量m,以及tk时刻的三维风速uw、vw、ww和三维风速变化率参数,并结合tk时刻的状态最优估计值利用求取tk+1时刻状态量的一步预测值
根据步骤(2)获取tk+1时刻的三维地速ug、vg、wg,利用得到tk+1时刻真空速的量测量,并结合步骤(1)得到的tk+1时刻俯仰角θ、横滚角φ、横滚角速度p、俯仰角速度q、偏航角速度r、三维线加速度ax、ay、az参数组成tk+1时刻总的量测量Zk+1,结合求解得到的tk+1时刻状态量的一步预测值和滤波增益矩阵Kk+1,利用公式 进行tk+1时刻状态量改善值的求解;
4.根据权利要求3所述的大攻角飞行状态下的大气攻角及侧滑角估计方法,其特征在于:所述步骤(10)中风速和风速变化率解算的具体方法是:
根据步骤(1)得到的tk+1时刻的俯仰角θ,横滚角φ和偏航角ψ,得到机体tk+1时刻的姿态转换矩阵Sθψφ:
结合步骤(2)获取的tk+1时刻的三维地速ug、vg、wg,得到tk+1时刻机体坐标系下沿x轴、y轴和z轴方向的三维地速ue,ve,we为:
结合步骤(9)估计得到的tk+1时刻的攻角α、侧滑角β、真空速V参数,利用下式将真空速V转换为机体坐标系下沿x轴、y轴和z轴方向的三维真空速ub、vb、wb,即
那么tk+1时刻机体坐标系下沿x轴、y轴和z轴方向的三维风速uw、vw、ww利用下式进行计算:
利用计算出的tk+1时刻的三维风速uw(tk+1)、vw(tk+1)、ww(tk+1),结合tk时刻的三维风速uw(tk)、vw(tk)、ww(tk),利用微分方法求取tk+1时刻三维风速变化率 即:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110426439 CN102520726B (zh) | 2011-12-19 | 2011-12-19 | 大攻角飞行状态下的大气攻角及侧滑角估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110426439 CN102520726B (zh) | 2011-12-19 | 2011-12-19 | 大攻角飞行状态下的大气攻角及侧滑角估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102520726A CN102520726A (zh) | 2012-06-27 |
CN102520726B true CN102520726B (zh) | 2013-07-03 |
Family
ID=46291678
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110426439 Expired - Fee Related CN102520726B (zh) | 2011-12-19 | 2011-12-19 | 大攻角飞行状态下的大气攻角及侧滑角估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102520726B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2784884C1 (ru) * | 2022-05-19 | 2022-11-30 | Федеральное государственное казенное военное образовательное учреждение высшего образования "Военная академия Ракетных войск стратегического назначения имени Петра Великого" МО РФ | Способ автоматического управления продольным движением беспилотного летательного аппарата при наличии ветрового возмущения |
Families Citing this family (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102809377B (zh) * | 2012-08-15 | 2015-08-12 | 南京航空航天大学 | 飞行器惯性/气动模型组合导航方法 |
CN103363993B (zh) * | 2013-07-06 | 2016-04-20 | 西北工业大学 | 一种基于无迹卡尔曼滤波的飞机角速率信号重构方法 |
CN103994748B (zh) * | 2014-05-27 | 2016-03-02 | 中国航天空气动力技术研究院 | 一种采用飞行和风洞试验数据估计无人机配平迎角的方法 |
US9880021B2 (en) * | 2014-10-20 | 2018-01-30 | Honeywell International Inc. | Systems and methods for attitude fault detection in one or more inertial measurement units |
CN105628086A (zh) * | 2014-10-29 | 2016-06-01 | 北京临近空间飞行器***工程研究所 | 一种基于锥面压力分布的超声速飞行来流参数解算方法 |
CN105136422B (zh) * | 2015-09-10 | 2017-10-13 | 中国航天空气动力技术研究院 | 风洞试验中修正飞行器模型侧滑弹性角的方法 |
CN105278536B (zh) * | 2015-11-05 | 2018-05-25 | 北京航天科颐技术有限公司 | 一种无人机姿态保持*** |
CN105912776B (zh) * | 2016-04-07 | 2019-03-08 | 中国人民解放军空军空降兵学院 | 一种基于多刚体模型的跳伞员坠落运动仿真方法 |
CN106840573B (zh) * | 2016-12-19 | 2019-04-30 | 中国航天空气动力技术研究院 | 一种嵌入式大气数据传感***标定方法 |
CN106682361A (zh) * | 2017-01-13 | 2017-05-17 | 沈阳航空航天大学 | 一种基于gps仿真的无人机飞行轨迹模拟***及方法 |
CN108475066B (zh) * | 2017-04-21 | 2021-02-19 | 深圳市大疆创新科技有限公司 | 无人飞行器姿态计算方法、飞行控制器及无人飞行器 |
CN107202664B (zh) * | 2017-05-24 | 2019-12-24 | 南京航空航天大学 | 一种用于嵌入式大气数据***的大气参数解算方法 |
CN110006425A (zh) * | 2019-04-11 | 2019-07-12 | 南京航空航天大学 | 基于载体动力学模型辅助的高动态角速度估计方法 |
CN110414110B (zh) * | 2019-07-19 | 2023-01-06 | 中仿智能科技(上海)股份有限公司 | 一种用于飞行失速状态下的飞机受力仿真方法 |
CN110427047B (zh) * | 2019-07-26 | 2022-06-03 | 深圳市道通智能航空技术股份有限公司 | 风速测算方法、风速估算器及无人机 |
CN110736854B (zh) * | 2019-09-29 | 2024-07-23 | 中航通飞华南飞机工业有限公司 | 一种基于机身两侧迎角传感器的飞行迎角获取方法 |
CN111122899B (zh) * | 2019-12-11 | 2020-11-17 | 南京航空航天大学 | 一种用于大气扰动中飞行的迎角侧滑角估计方法 |
CN111412887B (zh) * | 2020-03-31 | 2021-12-10 | 北京空天技术研究所 | 一种基于卡尔曼滤波的攻角、侧滑角辨识方法 |
CN114636842B (zh) * | 2022-05-17 | 2022-08-26 | 成都信息工程大学 | 一种高超声速飞行器的大气数据估计方法及装置 |
CN117251942B (zh) * | 2023-11-17 | 2024-03-08 | 成都凯天电子股份有限公司 | 一种估算飞行器空速、攻角和侧滑角的方法及*** |
CN117589190B (zh) * | 2024-01-18 | 2024-03-26 | 西北工业大学 | 基于惯导/飞控的大气参数解算方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102005058081B9 (de) * | 2005-12-06 | 2009-01-29 | Airbus Deutschland Gmbh | Verfahren zur Rekonstruktion von Böen und Strukturlasten bei Flugzeugen, insbesondere Verkehrsflugzeugen |
CN101413800B (zh) * | 2008-01-18 | 2010-09-29 | 南京航空航天大学 | 导航/稳瞄一体化***的导航、稳瞄方法 |
US8763955B1 (en) * | 2009-08-31 | 2014-07-01 | The Boeing Company | Method and apparatus for controlling a refueling drogue |
-
2011
- 2011-12-19 CN CN 201110426439 patent/CN102520726B/zh not_active Expired - Fee Related
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2784884C1 (ru) * | 2022-05-19 | 2022-11-30 | Федеральное государственное казенное военное образовательное учреждение высшего образования "Военная академия Ракетных войск стратегического назначения имени Петра Великого" МО РФ | Способ автоматического управления продольным движением беспилотного летательного аппарата при наличии ветрового возмущения |
Also Published As
Publication number | Publication date |
---|---|
CN102520726A (zh) | 2012-06-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102520726B (zh) | 大攻角飞行状态下的大气攻角及侧滑角估计方法 | |
Johansen et al. | On estimation of wind velocity, angle-of-attack and sideslip angle of small UAVs using standard sensors | |
CN108152529A (zh) | 一种基于飞行参数计算风速及风向的方法 | |
CN102809377B (zh) | 飞行器惯性/气动模型组合导航方法 | |
CN111766397B (zh) | 一种基于惯性/卫星/大气组合的气象风测量方法 | |
CN106527491A (zh) | 一种固定翼无人机控制***及横侧向飞行轨迹控制方法 | |
CN105005099B (zh) | 一种基于捷联惯导与飞行控制***的大气参数解算方法 | |
CN102425980B (zh) | 利用加速度计实现过载驾驶仪的控制方法 | |
Ramprasadh et al. | Multistage-fusion algorithm for estimation of aerodynamic angles in mini aerial vehicle | |
CN102607639A (zh) | 基于bp神经网络的大攻角飞行状态下大气数据测量方法 | |
CN103453907B (zh) | 基于分层大气模型的行星进入段导航滤波方法 | |
CN113111597A (zh) | 一种基于飞行数据的大气数据和扰动风估计方法 | |
CN111964688B (zh) | 结合无人机动力学模型和mems传感器的姿态估计方法 | |
CN106441372A (zh) | 一种基于偏振与重力信息的静基座粗对准方法 | |
CN101122637A (zh) | 一种sar运动补偿用sins/gps组合导航自适应降维滤波方法 | |
CN106249744A (zh) | 一种基于二级互补滤波的小型旋翼飞行器高度控制方法 | |
CN109541963B (zh) | 一种基于侧滑角信息的无人机测风建模方法 | |
CN104808673B (zh) | 一种基于卡尔曼滤波的四旋翼飞行器高度估计方法 | |
CN103033197B (zh) | 一种mems陀螺零位漂移的校正方法 | |
Lu et al. | Air data estimation by fusing navigation system and flight control system | |
CN112947522A (zh) | 一种基于有限时间观测器的硬式空中加油姿态控制方法 | |
CN114674345B (zh) | 一种惯导/相机/激光测速仪在线联合标定方法 | |
Emran et al. | A cascaded approach for quadrotor's attitude estimation | |
Witte et al. | Fundamental Turbulence Measurement with Unmanned Aerial Vehicles | |
Cao et al. | Research of attitude estimation of UAV based on information fusion of complementary filter |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130703 Termination date: 20151219 |
|
EXPY | Termination of patent right or utility model |