CN113029127A - 基于分布式循环结构的飞行器自主姿态估计方法 - Google Patents
基于分布式循环结构的飞行器自主姿态估计方法 Download PDFInfo
- Publication number
- CN113029127A CN113029127A CN202110261113.7A CN202110261113A CN113029127A CN 113029127 A CN113029127 A CN 113029127A CN 202110261113 A CN202110261113 A CN 202110261113A CN 113029127 A CN113029127 A CN 113029127A
- Authority
- CN
- China
- Prior art keywords
- attitude
- attitude estimation
- aircraft
- magnetometer
- equation
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000012937 correction Methods 0.000 claims abstract description 25
- 239000011159 matrix material Substances 0.000 claims abstract description 13
- 238000005259 measurement Methods 0.000 claims description 28
- 230000005358 geomagnetic field Effects 0.000 claims description 7
- 238000005096 rolling process Methods 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 2
- 238000004891 communication Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/005—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 with correlation of navigation data from several sources, e.g. map or contour matching
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/04—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by terrestrial means
- G01C21/08—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by terrestrial means involving use of the magnetic field of the earth
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/10—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
- G01C21/12—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
- G01C21/16—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
- G01C21/165—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Automation & Control Theory (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Navigation (AREA)
Abstract
本发明为了解决基于微机械陀螺/TAM姿态估计方法在线运行困难,精度较低的问题,提出了一种基于分布式循环结构的飞行器自主姿态估计方法。首先,建立自主姿态估计滤波器并将微机械陀螺和无校正误差的三轴磁强计数据代入,求得当前时刻飞行器的姿态;然后,利用已求解的当前时刻姿态矩阵建立TAM校正***;最后,将TAM校正***输出后无误差的数值重新输入姿态估计***,对下一时刻的姿态进行估计。
Description
技术领域
本发明属于飞行器姿态估计领域,涉及一种基于分布式循环结构的飞行器自主姿态估计方法。
背景技术
全球卫星定位***自应用以来取得了巨大的成就。但是随着反制技术的发展,其隐蔽性差、易受敌方电子战干扰的突出缺点也在多场局部战争中显露无疑。因此为了能进一步适应当前战争的形式,发展一种不需要数据链通信、可靠的自主导航方式已成为亟待解决的问题。微机械陀螺仪工作过程中不需要发射和接收信号,隐蔽性好、短期精度高、可提供全天候连续实时导航信息,已成为导航***的必选项。但同时具有积累误差,需要其他辅助信息及时修正。三轴磁强计(Three Axis Magnetometer,TAM)是目前常用的矢量磁测量传感器。与微机械陀螺仪组合后有更大的优势,惯保持了纯惯性导航设备全天候、连续导航的特点;被动测量,有很好的隐蔽性,避免通信时暴露和受到敌方诱骗、干扰;***精度高,能够为载体提供姿态、航向和位置等多维度高精度导航信息,在军事、安保等领域有很大的应用潜力。
但是,由于这两种传感器精度相对较低,如何获得较高的姿态精度是需要考虑的问题。TAM实际工作环境和线下测试时有很大区别,极端的温度梯度、机械压力、复杂的电磁环境等均对传感器的状态有较大影响,对姿态精度起着重要的作用。一些实际案例也证实了对TAM进行在线校正以确保获得一定精度的姿态确定结果来说是非常必要的。但对于在线校正问题,以上方法是将参数化成额外增加的状态。这样一个全阶滤波器的实现会有以下几个方面的困难:(1)状态向量的维数较大,会给实时操作带来计算负担;(2)状态包含快速变化项和缓变项,滤波器参数很难调节;(3)由于TAM误差校正参数较好的可观测性,所以量测更新能够更好的调节它,这将导致尺度因素和偏差对量测噪声较敏感,而姿态则反应较慢。可见,如何在姿态估计的同时更好地在线校正TAM误差参数是当前需要解决的一个关键问题。
针对现有微机械陀螺/TAM姿态估计***的不足,本发明提出一种基于分布式循环结构的飞行器自主姿态估计新方法。该方法将姿态估计和传感器校正分为两部分,通过这种分布式双循环的结构降低***维数,提高姿态估计精度。
发明内容
本发明的目的在于突破微机械陀螺以及TAM精度较低,TAM易受电磁环境干扰,在线估计状态维数较多的问题,提出一种基于分布式循环结构的飞行器自主姿态估计新方法。其中要解决的问题包括:
(1)状态向量的维数较大,会给实时操作带来计算负担;
(2)状态包含快速变化项和缓变项,滤波器参数很难调节;
(3)由于TAM误差校正参数较好的可观测性,所以量测更新能够更好的调节它,这将导致尺度因素和偏差对量测噪声较敏感,而姿态则反应较慢。
本发明所述的基于分布式循环结构的飞行器自主姿态估计方法,其特征在于包括以下技术措施:
步骤一:利用离散四元数动力学方程建立飞行器姿态估计***状态方程和三轴磁强计测量模型建立姿态估计***观测方程,从而完整构建飞行器姿态估计***,并求解该姿态估计***,从而得到当前时刻的姿态四元数、姿态转换矩阵、三轴姿态,利用当前时刻姿态矩阵将当地地磁场矢量理论值转化到载体坐标系下,将转化到载体坐标系下的当地地磁场矢量值与校正后的三轴磁强计测量值做差,得到地磁场残差向量;
步骤二:利用加入测量噪声的残差向量构建三轴磁强计校正***观测方程,以及由尺度因素误差和零偏两个误差参数作为***状态构建的三轴磁强计校正***状态方程,由此构建完整的三轴磁强计校正***,并利用UKF滤波器估计该***当前时刻的状态,也即误差参数;
步骤三:利用估计出的当前时刻的误差参数校正三轴磁强计测量值,并将校正好的测量值代入飞行器姿态估计***,重新估计三轴姿态。
对比现有技术,本发明所述的基于分布式循环结构的飞行器自主姿态估计方法,有益效果在于:
1)将整体***分为姿态估计***和传感器校正***,降低了每个滤波器的维数,减少了计算量和实时操作的负担;
2)将整体***分为姿态估计***和传感器校正***,使得快速变化项和缓变项能够分开,易于参数调节;
3)在姿态估计的同时,能够校正TAM误差参数,降低机载电磁环境对TAM测量的影响,提高***精度。
4)采用微机械陀螺和TAM构建起组合姿态估计***,在显著降低成本的同时,实现自主姿态估计,提高***隐蔽性、抗诱骗的能力。
附图说明
附图1是基于分布式循环结构的飞行器自主姿态估计方法步骤流程图;
附图2是分布式循环***结构图。
具体实施方式
为了解决基于微机械陀螺/TAM姿态估计方法在线运行困难,精度较低的问题,提出了一种基于分布式循环结构的飞行器自主姿态估计方法。首先,建立自主姿态估计滤波器并将微机械陀螺和无校正误差的三轴磁强计数据代入,求得当前时刻飞行器的姿态;然后,利用已求解的当前时刻姿态矩阵建立TAM校正***;最后,将TAM校正***输出后无误差的数值重新输入姿态估计***,对下一时刻的姿态进行估计。
以下结合说明书附图1对本发明做进一步的详细描述。参照说明书附图1,本发明的处理流程分以下步骤:
(1)当前时刻姿态四元数的一步预测
定义姿态四元数q为一四维向量:
q≡[q0 ρT]T,
其中ρ≡[q1 q2 q3]T为向量部分;q0为标量部分,有
由四元数q表示的姿态矩阵A(·)如下:
其中,I3×3为3×3的单位矩阵,<ρ×>为斜对称矩阵。
四元数微分方程为
本文采用光纤陀螺提取角速度ω,其输出模型可以表示为
ω=ω-b-ηa (2)
方程(1)、(2)、(3)可采用一阶龙格库塔法进行离散化。
(2)姿态估计***观测方程的构建
在载体坐标系下的磁场向量Bk和地理坐标系下磁场向量Hk的关系为
(3)基于滚动时域估计的飞行器姿态估计
通过陀螺和TAM建立姿态估计***模型之后,就需要考虑如何能根据传感器的数据得到飞行器当前最优姿态估计。考虑到四元数的归一化性质,本发明建立了了带有等式约束的滚动时域估计模型。它的基本思想为,只考虑当前时刻以前固定数量的数据,同时将历史数据对估计的影响用一种方法来近似描述,则每次参与优化的数据量不变,这种方法实现起来就像一个滚动的观测“窗口”。这个“窗”包含过去有限的量测值,“窗”的最后一个输出为当前状态估计。
窗口最后的一个点(i=L)代表现在时刻k,为方便起见,表示窗口中的某时刻时,直接采用L2,…,L表示。带等式约束的完整MHE问题可表示为
xi+1=fi+1,i(xi,ui)+wi,for i=1,…,L-1 (5)
yi+1=hi+1(xi+1)+vi+1,for i=1,…,L-1 (6)
li(xi)=0,for i=1,…,L (7)
其中,xi是待估计的状态。式(5)描述了在离散情况下的状态递推过程。在该递推方程中ui是已知输入,它的过程噪声为wi,过程噪声方差为Qi。方程(6)为量测方程,yi为量测值,vi、Ri分别为量测噪声和量测噪声方差。在时域内状态和方差是对历史信息的总结。式(8)中cnt,k为代价函数,它描述了到达代价。值得注意的是i=1时刻的量测值并没有包含在代价函数中。假定从量测中得到的所有信息已经包含在先验估计中。等式约束li由式(7)给出。MHE问题的输出为xL:时域“窗”中最后一个时刻的估计状态向量。
这就是一个带权重的最小二乘问题,在该问题中权重就是方差阵的逆:
采用二范数记法:
把方程(5)、方程(6)、方程(7)带入到代价函数(4)中,得到下式:
s.t.
li(xi)=0,for i=1,…,L
这样就将状态估计问题转化成了非线性规划(NPL)问题,可以用数值优化的方法求解,代入陀螺数据,以及经过校正的三轴磁强计的测量值B′m,k,能够解得当前时刻的四元数最优估计值、姿态矩阵以及飞行器三轴姿态。姿态估计过程以磁强计测量周期ΔT进行更新。
(4)三轴磁强计测量残差的求解
制约磁强计精度的各种误差因素会引起姿态估计误差,因此必须在线估计传感器的参数以排除传感器误差的影响。传感器校正过程是通过使用姿态估计阶段的残差来估计k时刻的尺度因素sk和零偏dk。误差参数定义为Θk:
Θk=[sx,k sy,k sz,k dx,k dy,k dz,k]T
令为滤波后的四元数最优估计值,Bm,k为磁强计在k时刻的测量值。假设离散区间内角速率保持恒定,该时刻的磁场三分量测量误差向量Δyk即为当地地磁场向量在机体坐标系下的理论值与经过误差校正的三轴磁强计测量值之差
(5)三轴磁强计误差参数的估计
将Δyk积累一定数量之后就可以用来校正传感器。令n为一个循环中积累的量测值数量,那么量测方程可以写为
zi=Δyi+ri,i=1,…,n (9)
其中,zi是在一个循环内积累的第i个时刻的误差,Δyi为地磁场残差向量,ri是相应的噪声,它对姿态估计误差也有影响。
在离散形式下,传感器校正***的状态方程为
其中,bi是一个过程噪声项,用来表示尺度因素和零偏的变化量。对于通过式(9)和式(10)建立的传感器校正***可以用采用UKF来完成滤波。
(6)三轴磁强计测量值的校正
得到校正后的三轴磁强计测量值B′k,m。将该值代入姿态估计***循环,可以估计飞行器姿态。
Claims (5)
1.基于分布式循环结构的飞行器自主姿态估计方法,其特征在于,包括以下步骤:
步骤一:利用离散四元数动力学方程建立飞行器姿态估计***状态方程和三轴磁强计测量模型建立姿态估计***观测方程,从而完整构建飞行器姿态估计***,并求解该姿态估计***,从而得到当前时刻的姿态四元数、姿态转换矩阵、三轴姿态,利用当前时刻姿态矩阵将当地地磁场矢量理论值转化到载体坐标系下,将转化到载体坐标系下的当地地磁场矢量值与校正后的三轴磁强计测量值做差,得到地磁场残差向量;
步骤二:利用加入测量噪声的残差向量构建三轴磁强计校正***观测方程,以及由尺度因素误差和零偏两个误差参数作为***状态构建的三轴磁强计校正***状态方程,由此构建完整的三轴磁强计校正***,并利用UKF滤波器估计该***当前时刻的状态,也即误差参数;
步骤三:利用估计出的当前时刻的误差参数校正三轴磁强计测量值,并将校正好的测量值代入飞行器姿态估计***,重新估计三轴姿态。
3.根据权利要求1所述的基于分布式循环结构的飞行器自主姿态估计方法,其特征在于,步骤一中飞行器姿态估计***求解方法为:
将离散的飞行器姿态估计***代入滚动时域估计模型,采用高斯-牛顿迭代法并结合MEMS速率陀螺和校正好的三轴磁强计测量值求解该滚动时域估计模型,从而得到当前时刻的姿态四元数、姿态转换矩阵、三轴姿态。
4.根据权利要求1所述的基于分布式循环结构的飞行器自主姿态估计方法,其特征在于,步骤二中三轴磁强计校正***观测方程构建方法为:
令n为一个循环中积累的量测值数量,那么将观测方程写为
zi=Δyi+ri,i=1,…,n
其中,zi是在一个循环内积累的第i个时刻的误差,Δyi为地磁场残差向量,ri是相应的噪声。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110261113.7A CN113029127B (zh) | 2021-03-10 | 2021-03-10 | 基于分布式循环结构的飞行器自主姿态估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110261113.7A CN113029127B (zh) | 2021-03-10 | 2021-03-10 | 基于分布式循环结构的飞行器自主姿态估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113029127A true CN113029127A (zh) | 2021-06-25 |
CN113029127B CN113029127B (zh) | 2023-05-09 |
Family
ID=76469236
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110261113.7A Active CN113029127B (zh) | 2021-03-10 | 2021-03-10 | 基于分布式循环结构的飞行器自主姿态估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113029127B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TW498169B (en) * | 2000-11-29 | 2002-08-11 | American Gnc Corp | Interruption-free hand-held positioning method and system thereof |
CN106342284B (zh) * | 2008-08-18 | 2011-11-23 | 西北工业大学 | 一种飞行载体姿态确定方法 |
US8515672B1 (en) * | 2012-03-01 | 2013-08-20 | Honeywell International Inc. | Systems and methods to incorporate master navigation system resets during transfer alignment |
CN107544042A (zh) * | 2016-06-29 | 2018-01-05 | 中国科学院电子学研究所 | 一种磁力计阵列校正方法 |
CN108919156A (zh) * | 2018-06-27 | 2018-11-30 | 中国人民解放军海军航空大学 | 基于噪声补偿的三轴磁强计线下校正方法 |
CN109001787A (zh) * | 2018-05-25 | 2018-12-14 | 北京大学深圳研究生院 | 一种姿态角解算与定位的方法及其融合传感器 |
CN111076722A (zh) * | 2019-11-18 | 2020-04-28 | 广州南方卫星导航仪器有限公司 | 基于自适应的四元数的姿态估计方法及装置 |
CN111578931A (zh) * | 2020-05-21 | 2020-08-25 | 中国人民解放军海军航空大学 | 基于在线滚动时域估计的高动态飞行器自主姿态估计方法 |
CN111739063A (zh) * | 2020-06-23 | 2020-10-02 | 郑州大学 | 一种基于多传感器融合的电力巡检机器人定位方法 |
-
2021
- 2021-03-10 CN CN202110261113.7A patent/CN113029127B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TW498169B (en) * | 2000-11-29 | 2002-08-11 | American Gnc Corp | Interruption-free hand-held positioning method and system thereof |
CN106342284B (zh) * | 2008-08-18 | 2011-11-23 | 西北工业大学 | 一种飞行载体姿态确定方法 |
US8515672B1 (en) * | 2012-03-01 | 2013-08-20 | Honeywell International Inc. | Systems and methods to incorporate master navigation system resets during transfer alignment |
CN107544042A (zh) * | 2016-06-29 | 2018-01-05 | 中国科学院电子学研究所 | 一种磁力计阵列校正方法 |
CN109001787A (zh) * | 2018-05-25 | 2018-12-14 | 北京大学深圳研究生院 | 一种姿态角解算与定位的方法及其融合传感器 |
CN108919156A (zh) * | 2018-06-27 | 2018-11-30 | 中国人民解放军海军航空大学 | 基于噪声补偿的三轴磁强计线下校正方法 |
CN111076722A (zh) * | 2019-11-18 | 2020-04-28 | 广州南方卫星导航仪器有限公司 | 基于自适应的四元数的姿态估计方法及装置 |
CN111578931A (zh) * | 2020-05-21 | 2020-08-25 | 中国人民解放军海军航空大学 | 基于在线滚动时域估计的高动态飞行器自主姿态估计方法 |
CN111739063A (zh) * | 2020-06-23 | 2020-10-02 | 郑州大学 | 一种基于多传感器融合的电力巡检机器人定位方法 |
Non-Patent Citations (1)
Title |
---|
赵国荣等: "基于滚动时域估计的飞行器姿态估计及三轴磁强计在线校正", 《物理学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN113029127B (zh) | 2023-05-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107270893B (zh) | 面向不动产测量的杆臂、时间不同步误差估计与补偿方法 | |
CN103630137B (zh) | 一种用于导航***的姿态及航向角的校正方法 | |
CN104698485B (zh) | 基于bd、gps及mems的组合导航***及导航方法 | |
CN105910602B (zh) | 一种组合导航方法 | |
CN111189442B (zh) | 基于cepf的无人机多源导航信息状态预测方法 | |
CN108761512A (zh) | 一种弹载bds/sins深组合自适应ckf滤波方法 | |
CN109443349A (zh) | 一种姿态航向测量***及其融合方法、存储介质 | |
CN112146655A (zh) | 一种BeiDou/SINS紧组合导航***弹性模型设计方法 | |
CN112325886A (zh) | 一种基于重力梯度仪和陀螺仪组合的航天器自主定姿*** | |
CN112461262A (zh) | 一种校正三轴磁强计误差的装置和方法 | |
CN115507849A (zh) | 一种基于ins/gnss组合导航辅助的磁传感器校正方法和*** | |
CN110849360A (zh) | 面向多机协同编队飞行的分布式相对导航方法 | |
CN107860382B (zh) | 一种在地磁异常情况下应用ahrs测量姿态的方法 | |
Zou et al. | A nonlinear transfer alignment of distributed POS based on adaptive second-order divided difference filter | |
CN115900770A (zh) | 一种机载环境下磁传感器的在线校正方法和*** | |
Soken et al. | In flight magnetometer calibration via unscented Kalman filter | |
CN112284388B (zh) | 一种无人机多源信息融合导航方法 | |
CN114459466A (zh) | 一种基于模糊控制的mems多传感器数据融合处理方法 | |
Gong et al. | Airborne earth observation positioning and orientation by SINS/GPS integration using CD RTS smoothing | |
CN112683265B (zh) | 一种基于快速iss集员滤波的mimu/gps组合导航方法 | |
CN113029127A (zh) | 基于分布式循环结构的飞行器自主姿态估计方法 | |
Candan et al. | Estimation of attitude using robust adaptive Kalman filter | |
CN114812546A (zh) | 遮蔽空间单兵导航位姿修正方法及装置 | |
Wu et al. | Real-time magnetic disturbance determination for micro air vehicles via gravity and global navigation satellite system measurements | |
CN117948986B (zh) | 一种极区因子图构建方法和极区组合导航方法 |
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 |