CN107677292A - 基于重力场模型的垂线偏差补偿方法 - Google Patents
基于重力场模型的垂线偏差补偿方法 Download PDFInfo
- Publication number
- CN107677292A CN107677292A CN201710894464.5A CN201710894464A CN107677292A CN 107677292 A CN107677292 A CN 107677292A CN 201710894464 A CN201710894464 A CN 201710894464A CN 107677292 A CN107677292 A CN 107677292A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- mtd
- msubsup
- inertial navigation
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
- G01C25/005—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass initial alignment, calibration or starting-up of inertial devices
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Navigation (AREA)
Abstract
本发明针对现有惯性导航系统精度受到垂线偏差影响而降低的问题，公开了一种基于重力场模型的垂线偏差补偿方法，属于惯性导航技术领域。本发明通过所述标定、计算正常重力矢量、由重力场模型计算垂线偏差、根据垂线偏差计算补偿后的重力矢量、初始对准和导航解算等步骤，补偿了垂线偏差对惯性导航系统的影响，提高了惯性导航定位精度，本发明步骤简单、易于实现，可在不增加额外的硬件成本的条件下减弱垂线偏差对惯性导航影响，提升惯性导航系统精度。
Description
技术领域：
本发明涉及重力辅助惯性导航技术，特别是涉及一种基于重力场模型的垂线偏差补偿方法，属于惯性导航领域。
背景技术：
惯性导航系统(Inertial Navigation System,INS)具有高精度、全天候、完全自主的优点，广泛地应用于航空、航天、舰船及地面车辆导航。垂线偏差所引起的水平重力扰动是惯性导航系统的主要误差源之一，对惯性导航系统的姿态测量精度、速度测量精度及定位精度有较大影响，特别是在长时间的惯性导航中表现尤为明显，垂线偏差引起惯性导航误差的原因在于定义导航坐标系时，假设的重力垂线方向(参考椭球的法线方向)与真实的重力垂线方向不一致，两者方向的差异正是垂线偏差。针对垂线偏差进行补偿以提高惯性导航系统精度是当前惯性导航技术领域的研究热点，垂线偏差补偿方法主要分为两类：基于垂线偏差实时测量数据的补偿方法和基于重力场模型的垂线偏差补偿方法。基于垂线偏差实时测量数据的补偿方法的原理是通过重力梯度仪或重力矢量仪实时测量垂线偏差，并将其补偿到惯性导航解算中以提高惯性导航系统精度，但重力梯度仪、重力矢量仪的使用受成本、获取渠道等诸多条件限制，难以广泛推广应用。基于重力场模型的垂线偏差补偿方法的原理是利用重力场模型计算垂线偏差并补偿到惯性导航解算中以提高惯性导航系统精度，目前已有多个适用于全球的重力场模型公布，为此类方法的推广应用奠定了基础。惯性导航是由初始对准与导航解算两个步骤组成，国内外学者在开展垂线偏差补偿方法研究时，主要关注了垂线偏差对导航解算的影响，而没有考虑其对初始对准的影响，但垂线偏差将影响初始对准所建立的导航坐标系，如果惯性导航解算中的各个向量未投影到统一的导航坐标系就会引起惯性导航误差。因此，进行垂线偏差补偿时，必须要从初始对准阶段就考虑垂线偏差的影响，并在初始对准与导航解算两个阶段都进行补偿，才能有效地补偿垂线偏差影响，提高惯性导航系统精度。
发明内容：
针对以上问题，本发明通过标定、计算正常重力矢量、由重力场模型计算垂线偏差、根据垂线偏差计算补偿后的重力矢量、初始对准和导航解算等步骤，补偿了垂线偏差对惯性导航系统的影响，提高了惯性导航定位精度。
一种基于重力场模型的垂线偏差补偿方法，其步骤为：
步骤S1：惯性导航系统工作前，通过高精度三轴转台对惯性导航系统进行标定，精确地估计出加速度计零偏，并在加速度计输出数据中扣除零偏；
步骤S2：定义导航系n和n′，其中n表示以参考椭球法线为垂向建立的导航坐标系，n′表示以真实重力为垂向建立的导航坐标系；使惯性导航系统处于静止状态，根据惯性导航系统的初始纬度L0、初始经度λ0和初始高度h0，使用正常重力模型计算初始位置处的标量正常重力值，以此构建未经过垂线偏差补偿的重力矢量向量上标n表示该向量投影在导航系n；
步骤S3：根据惯性导航系统的初始纬度L0、初始经度λ0和初始高度h0，使用EGM2008重力场球谐函数模型计算垂线偏差的北向分量ξ和东向分量η；
步骤S4：将步骤S2所述未经垂线偏差补偿的正常重力矢量绕旋转轴转动角度后，得到垂线偏差补偿后的重力矢量通过以下步骤实现：
步骤S41:根据步骤S3所述ξ和η，计算转动轴其中
步骤S42：计算转动角度
步骤S43：由转动轴和转动角度构建垂线偏差补偿旋转四元数Q＝[q0 q1 q2q3]T，其中，
步骤S44：由垂线偏差补偿旋转四元数Q＝[q0 q1 q2 q3]T，得到垂线偏差补偿姿态旋转矩阵
步骤S45：由垂线偏差补偿姿态旋转矩阵得到经过垂线偏差补偿后的重力矢量
步骤S5：将步骤S4所述经过垂线偏差补偿后的重力矢量用于解析粗对准，得到姿态矩阵的粗略估计值b表示载体坐标系；
步骤S6：建立五状态卡尔曼滤波器进行精对准，滤波器状态包括姿态误差和水平速度误差，并将步骤S4所述经过垂线偏差补偿后的重力矢量用于精对准中的惯性导航解算，得到姿态矩阵的精确估计值惯性导航初始对准完成；
步骤S7：惯性导航系统进入动态导航阶段，使用步骤S6所述的进行惯性导航姿态更新，同时使用步骤S4所述的经过垂线偏差补偿后的重力矢量用于惯性导航速度更新。
附图说明：
图1为导航坐标系n和n′示意图；
图2为由垂线偏差获得补偿后的重力矢量示意图；
图3为本发明所述基于重力场模型的垂线偏差补偿方法流程图；
图4为船载试验中，使用本发明补偿垂线偏差效果。
具体实施方式：
下面结合附图通过实施例对本发明的技术方案进行清楚、完整地描述。
步骤S1：惯性导航系统工作前，通过高精度三轴转台对惯性导航系统进行标定，精确地估计出加速度计零偏，并在加速度计输出数据中扣除零偏；
步骤S2：定义载体坐标系b为满足右手定律的正交坐标系，其原点在惯性导航系统质心，其轴向分别沿安装有惯性导航系统的运载体的俯仰轴、横滚轴及航向轴；如图1，定义导航坐标系n为满足右手定律的正交坐标系，其原点在载体所处位置，其xn轴指向东，yn轴指向北，zn轴指向天且与载体所处位置的参考椭球法线共线(参考椭球为WGS-84椭球模型)；如图1，定义导航坐标系n′为满足右手定律的正交坐标系，其原点在载体所处位置，其xn′轴指向东，yn′轴指向椭球北极点，zn′轴指向天且与载体所处位置的真实重力共线；使惯性导航系统处于静止状态，惯性导航系统的初始纬度L0、初始经度λ0和初始高度h0由卫星导航接收机或已知地标点确定，通过WGS-84坐标系下的正常重力模型计算初始位置处的标量正常重力值γ，以此构建导航坐标系n下未经过垂线偏差补偿的重力矢量其中，
偏心率e＝0.0818191908426；
γ0表示初始位置处未经过高度改正的正常重力值；
步骤S3：根据惯性导航系统的初始纬度L0、初始经度λ0和初始高度h0，使用EGM2008重力场球谐函数模型计算垂线偏差的北向分量ξ和东向分量η；
步骤S4：如图2，将步骤S2所述未经垂线偏差补偿的正常重力矢量绕旋转轴转动角度后，得到垂线偏差补偿后的重力矢量通过以下步骤实现：
步骤S41:根据步骤S3所述ξ和η，计算转动轴其中
步骤S42：计算转动角度
步骤S43：由转动轴和转动角度构建垂线偏差补偿旋转四元数Q＝[q0 q1 q2q3]T，其中：
步骤S44：由垂线偏差补偿旋转四元数Q＝[q0 q1 q2 q3]T，得到垂线偏差补偿姿态旋转矩阵
步骤S45：由垂线偏差补偿姿态旋转矩阵得到经过垂线偏差补偿后的重力矢量
步骤S5：将步骤S4所述经过垂线偏差补偿后的重力矢量用于解析粗对准，得到姿态矩阵的粗略估计值解析粗对准公式为：
其中，表示惯性导航系统中的加速度计测量的比力向量在所述载体坐标系b的投影；表示惯性导航系统中的陀螺测量的角速度向量在所述载体坐标系b的投影；表示初始位置的地球自转角速度向量在所述导航坐标系n的投影，Ω为地球自转角速率，在WGS-84坐标系下Ω＝7.292115×10-5rad/s：
步骤S6：建立五状态卡尔曼滤波器进行精对准，滤波器状态包括姿态误差和水平速度误差，并将步骤S4所述经过垂线偏差补偿后的重力矢量用于精对准中的惯性导航解算，得到姿态矩阵的精确估计值滤波器状态分别是：俯仰角误差θ、横滚角误差γ、航向角误差φ、东向速度误差δvE和北向速度误差δvN，所述的卡尔曼滤波器状态量如下：
所述卡尔曼滤波器状态方程如下
所述卡尔曼滤波器状态矩阵F如下
其中，vE、vN和vU分别是惯性导航系统计算的东向速度、北向速度和天向速度，L、h分别是惯性导航系统计算的纬度和高度，RN、RE分别是惯性导航系统计算的子午圈曲率半径、卯酉圈曲率半径，fE、fN及fU分别是惯性导航系统测量的比力向量在导航坐标系n投影的三个分量；状态方程噪声为w，包含陀螺测量噪声和加速度计测量噪声，状态方程噪声向量w的方差矩阵var(w)为：
δwGyro表示陀螺测量噪声标准差，δwAcc表示加速度计测量噪声标准差，δwGyro和δwAcc应根据所述惯性导航系统使用的陀螺与加速度计噪声水平设定，本实施例中设定为δwGyro＝0.01°/h和δwAcc＝5×10-4m/s2，所述卡尔曼滤波器观测量z及观测方程如下：
z＝[0-vE 0-vN]T＝[-vE -vN]T
其中，υ为观测噪声，其标准差可设置为10-1～10-2(m/s2)，本实施例中设置为10-2(m/s2)，H为观测矩阵，东向速度、北向速度观测均为0；
步骤S7：惯性导航系统进入动态导航阶段，使用步骤S6所述的进行惯性导航姿态更新，同时使用步骤S4所述的经过垂线偏差补偿后的重力矢量用于惯性导航速度更新。
将本发明所述方法用于船载高精度激光陀螺捷联式惯性导航系统的垂线偏差补偿，双频GPS(Global Position System，全球定位系统)接收机数据的精密单点定位结果为惯性导航计算提供初始位置和速度信息，并作为惯性导航系统位置的参考值，6小时的惯性导航定位误差如图4所示，其中实线为未使用垂线偏差补偿的惯性导航定位误差，虚线为使用本发明补偿后的定位误差，可以看到使用本发明所述方法补偿垂线偏差后，惯性导航系统精度有较为明显的提高，惯性导航6小时，未经垂线偏差补偿的导航误差为3450米，使用本发明所述方法补偿垂线偏差后，导航误差为2950米，惯性导航定位精度提升14.5％。
以上仅是本发明的优选实施方式，本发明的保护范围并不仅局限于上述实施例，凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理前提下的若干改进和润饰，应视为本发明的保护范围。
Claims (1)
1.基于重力场模型的垂线偏差补偿方法，其特征在于包含以下步骤：
步骤S1：惯性导航系统工作前，通过高精度三轴转台对惯性导航系统进行标定，精确地估计出加速度计零偏，并在加速度计输出数据中扣除零偏；
步骤S2：定义导航系n和n′，其中n表示以参考椭球法线为垂向建立的导航坐标系，n′表示以真实重力为垂向建立的导航坐标系；使惯性导航系统处于静止状态，根据惯性导航系统的初始纬度L0、初始经度λ0和初始高度h0，使用正常重力模型计算初始位置处的标量正常重力值，以此构建未经过垂线偏差补偿的重力矢量向量上标n表示该向量投影在导航系n；
步骤S3：根据惯性导航系统的初始纬度L0、初始经度λ0和初始高度h0，使用EGM2008重力场球谐函数模型计算垂线偏差的北向分量ξ和东向分量η；
步骤S4：将步骤S2所述未经垂线偏差补偿的正常重力矢量绕旋转轴转动角度θ后，得到垂线偏差补偿后的重力矢量通过以下步骤实现：
步骤S41:根据步骤S3所述ξ和η，计算转动轴其中
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mi>&eta;</mi>
<mo>/</mo>
<msqrt>
<mrow>
<msup>
<mi>&xi;</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>&eta;</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<mi>&xi;</mi>
<mo>/</mo>
<msqrt>
<mrow>
<msup>
<mi>&xi;</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>&eta;</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
步骤S42：计算转动角度
步骤S43：由转动轴和转动角度θ构建垂线偏差补偿旋转四元数Q＝[q0 q1 q2 q3]T，其中，
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>q</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mfrac>
<mi>&theta;</mi>
<mn>2</mn>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mfrac>
<mi>&theta;</mi>
<mn>2</mn>
</mfrac>
<mo>&CenterDot;</mo>
<msub>
<mi>u</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>q</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mfrac>
<mi>&theta;</mi>
<mn>2</mn>
</mfrac>
<mo>&CenterDot;</mo>
<msub>
<mi>u</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>q</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mfrac>
<mi>&theta;</mi>
<mn>2</mn>
</mfrac>
<mo>&CenterDot;</mo>
<msub>
<mi>u</mi>
<mn>3</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
步骤S44：由垂线偏差补偿旋转四元数Q＝[q0 q1 q2 q3]T，得到垂线偏差补偿姿态旋转矩阵
<mrow>
<msubsup>
<mi>C</mi>
<msup>
<mi>n</mi>
<mo>&prime;</mo>
</msup>
<mi>n</mi>
</msubsup>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>q</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>q</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>q</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>q</mi>
<mn>3</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>q</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>q</mi>
<mn>0</mn>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>q</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>q</mi>
<mn>3</mn>
</msub>
<mo>+</mo>
<msub>
<mi>q</mi>
<mn>0</mn>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>q</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>q</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>q</mi>
<mn>0</mn>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>q</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>q</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>q</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>q</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>q</mi>
<mn>3</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>q</mi>
<mn>2</mn>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>q</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<msub>
<mi>q</mi>
<mn>0</mn>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>q</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<msub>
<mi>q</mi>
<mn>0</mn>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>q</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>q</mi>
<mn>2</mn>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>q</mi>
<mn>3</mn>
</msub>
<mo>+</mo>
<msub>
<mi>q</mi>
<mn>0</mn>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>q</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>q</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>q</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>q</mi>
<mn>3</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
步骤S45：由垂线偏差补偿姿态旋转矩阵得到经过垂线偏差补偿后的重力矢量
<mrow>
<msup>
<mover>
<mi>g</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>n</mi>
</msup>
<mo>=</mo>
<msubsup>
<mi>C</mi>
<msup>
<mi>n</mi>
<mo>&prime;</mo>
</msup>
<mi>n</mi>
</msubsup>
<msup>
<mover>
<mi>&gamma;</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>n</mi>
</msup>
</mrow>
步骤S5：将步骤S4所述经过垂线偏差补偿后的重力矢量用于解析粗对准，得到姿态矩阵的粗略估计值b表示载体坐标系；
步骤S6：建立五状态卡尔曼滤波器进行精对准，滤波器状态包括姿态误差和水平速度误差，并将步骤S4所述经过垂线偏差补偿后的重力矢量用于精对准中的惯性导航解算，得到姿态矩阵的精确估计值惯性导航初始对准完成；
步骤S7：惯性导航系统进入动态导航阶段，使用步骤S6所述的进行惯性导航姿态更新，同时使用步骤S4所述的经过垂线偏差补偿后的重力矢量用于惯性导航速度更新。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710894464.5A CN107677292B (zh) | 2017-09-28 | 2017-09-28 | 基于重力场模型的垂线偏差补偿方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710894464.5A CN107677292B (zh) | 2017-09-28 | 2017-09-28 | 基于重力场模型的垂线偏差补偿方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107677292A true CN107677292A (zh) | 2018-02-09 |
| CN107677292B CN107677292B (zh) | 2019-11-15 |
Family
ID=61138862
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710894464.5A Active CN107677292B (zh) | 2017-09-28 | 2017-09-28 | 基于重力场模型的垂线偏差补偿方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107677292B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109059915A (zh) * | 2018-09-27 | 2018-12-21 | 清华大学 | 重力补偿方法、系统和装置 |
| CN110850382A (zh) * | 2019-12-18 | 2020-02-28 | 中国人民解放军61540部队 | 一种评估干涉雷达高度计测量精度的方法及系统 |
| CN111829553A (zh) * | 2020-06-18 | 2020-10-27 | 中国船舶重工集团公司第七0七研究所 | 一种基于pc-104的高精度惯导系统扰动重力补偿方法 |
| CN112798014A (zh) * | 2020-10-16 | 2021-05-14 | 中国人民解放军火箭军工程大学 | 一种基于重力场球谐模型补偿垂线偏差的惯导自对准方法 |
| CN114137625A (zh) * | 2021-11-24 | 2022-03-04 | 中国船舶重工集团公司第七0七研究所 | 一种基于两套惯导互观测的海洋垂线偏差测量方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102052922A (zh) * | 2010-11-19 | 2011-05-11 | 中国人民解放军海军工程大学 | 实际重力场影响惯性导航系统扰动重力补偿方法 |
| WO2013059989A1 (zh) * | 2011-10-25 | 2013-05-02 | 国防科学技术大学 | 一种惯性导航系统的运动对准方法 |
| CN104359496A (zh) * | 2014-11-26 | 2015-02-18 | 中国人民解放军国防科学技术大学 | 基于垂线偏差补偿的高精度姿态修正方法 |
| CN104374401A (zh) * | 2014-10-15 | 2015-02-25 | 哈尔滨工程大学 | 一种捷联惯导初始对准中重力扰动的补偿方法 |
| US9417091B2 (en) * | 2013-05-13 | 2016-08-16 | The Johns Hopkins University | System and method for determining and correcting field sensors errors |
-
2017
- 2017-09-28 CN CN201710894464.5A patent/CN107677292B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102052922A (zh) * | 2010-11-19 | 2011-05-11 | 中国人民解放军海军工程大学 | 实际重力场影响惯性导航系统扰动重力补偿方法 |
| WO2013059989A1 (zh) * | 2011-10-25 | 2013-05-02 | 国防科学技术大学 | 一种惯性导航系统的运动对准方法 |
| US9417091B2 (en) * | 2013-05-13 | 2016-08-16 | The Johns Hopkins University | System and method for determining and correcting field sensors errors |
| CN104374401A (zh) * | 2014-10-15 | 2015-02-25 | 哈尔滨工程大学 | 一种捷联惯导初始对准中重力扰动的补偿方法 |
| CN104359496A (zh) * | 2014-11-26 | 2015-02-18 | 中国人民解放军国防科学技术大学 | 基于垂线偏差补偿的高精度姿态修正方法 |
Non-Patent Citations (5)
| Title |
|---|
| RUONAN WU 等: "Gravity compensation using EGM2008 for high-precision long-term inertial navigation systems", 《SENSORS》 * |
| X. MU 等: "Direct compensation and kalman filter compensation of gravity vertical deviation based on grid database", 《2017 24TH SAINT PETERSBURG INTERNATIONAL CONFERENCE ON INTEGRATED NAVIGATION SYSTEMS (ICINS)》 * |
| 何昆鹏 等: "重力扰动对高精度惯性导航初始对准的影响与补偿", 《导航与控制》 * |
| 王皓: "高精度长航时惯导系统重力扰动补偿方法研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
| 赵忠 等: "高精度惯性导航系统垂线偏差影响与补偿", 《中国惯性技术学报》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109059915A (zh) * | 2018-09-27 | 2018-12-21 | 清华大学 | 重力补偿方法、系统和装置 |
| CN109059915B (zh) * | 2018-09-27 | 2020-12-01 | 清华大学 | 重力补偿方法、系统和装置 |
| CN110850382A (zh) * | 2019-12-18 | 2020-02-28 | 中国人民解放军61540部队 | 一种评估干涉雷达高度计测量精度的方法及系统 |
| CN111829553A (zh) * | 2020-06-18 | 2020-10-27 | 中国船舶重工集团公司第七0七研究所 | 一种基于pc-104的高精度惯导系统扰动重力补偿方法 |
| CN112798014A (zh) * | 2020-10-16 | 2021-05-14 | 中国人民解放军火箭军工程大学 | 一种基于重力场球谐模型补偿垂线偏差的惯导自对准方法 |
| CN114137625A (zh) * | 2021-11-24 | 2022-03-04 | 中国船舶重工集团公司第七0七研究所 | 一种基于两套惯导互观测的海洋垂线偏差测量方法 |
| CN114137625B (zh) * | 2021-11-24 | 2023-04-28 | 中国船舶重工集团公司第七0七研究所 | 一种基于两套惯导互观测的海洋垂线偏差测量方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107677292B (zh) | 2019-11-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108051866B (zh) | 基于捷联惯性/gps组合辅助水平角运动隔离的重力测量方法 | |
| CN107270893B (zh) | 面向不动产测量的杆臂、时间不同步误差估计与补偿方法 | |
| CN101949703B (zh) | 一种捷联惯性/卫星组合导航滤波方法 | |
| CN107677292B (zh) | 基于重力场模型的垂线偏差补偿方法 | |
| CN103090867B (zh) | 相对地心惯性系旋转的光纤陀螺捷联惯性导航系统误差抑制方法 | |
| CN103245359B (zh) | 一种惯性导航系统中惯性传感器固定误差实时标定方法 | |
| CN103900611B (zh) | 一种惯导天文高精度复合两位置对准及误差标定方法 | |
| CN102393201B (zh) | 航空遥感用位置和姿态测量系统（pos）动态杆臂补偿方法 | |
| CN101571394A (zh) | 基于旋转机构的光纤捷联惯性导航系统初始姿态确定方法 | |
| CN107655476A (zh) | 基于多信息融合补偿的行人高精度足部导航算法 | |
| CN105371844B (zh) | 一种基于惯性/天文互助的惯性导航系统初始化方法 | |
| CN103743413B (zh) | 倾斜状态下调制寻北仪安装误差在线估计与寻北误差补偿方法 | |
| CN103792561B (zh) | 一种基于gnss通道差分的紧组合降维滤波方法 | |
| CN103900608B (zh) | 一种基于四元数ckf的低精度惯导初始对准方法 | |
| CN105091907B (zh) | Sins/dvl组合中dvl方位安装误差估计方法 | |
| CN112097763A (zh) | 一种基于mems imu/磁力计/dvl组合的水下运载体组合导航方法 | |
| CN102169184A (zh) | 组合导航系统中测量双天线gps安装失准角的方法和装置 | |
| CN101514900A (zh) | 一种单轴旋转的捷联惯导系统初始对准方法 | |
| CN106289246A (zh) | 一种基于位置和姿态测量系统的柔性杆臂测量方法 | |
| CN113203418A (zh) | 基于序贯卡尔曼滤波的gnssins视觉融合定位方法及系统 | |
| CN102706366A (zh) | 一种基于地球自转角速率约束的sins初始对准方法 | |
| CN106405670A (zh) | 一种适用于捷联式海洋重力仪的重力异常数据处理方法 | |
| CN1995916A (zh) | 一种基于星敏感器标定的深综合组合导航方法 | |
| CN111947653A (zh) | 一种月面巡视探测器双模式惯性/视觉/天文导航方法 | |
| CN101881619A (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 |