CN104344836B - 一种基于姿态观测的冗余惯导系统光纤陀螺系统级标定方法 - Google Patents
一种基于姿态观测的冗余惯导系统光纤陀螺系统级标定方法 Download PDFInfo
- Publication number
- CN104344836B CN104344836B CN201410599074.1A CN201410599074A CN104344836B CN 104344836 B CN104344836 B CN 104344836B CN 201410599074 A CN201410599074 A CN 201410599074A CN 104344836 B CN104344836 B CN 104344836B
- Authority
- CN
- China
- Prior art keywords
- inertial navigation
- error
- optical fibre
- fibre gyro
- angle
- 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
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)
- Gyroscopes (AREA)
Abstract
一种基于姿态观测的冗余惯导系统光纤陀螺系统级标定方法,有八大步骤:一、将惯导安装转台上,确定载体的初始位置参数;二、确定光纤陀螺轴向与惯导本体坐标系安装关系;三、惯导预热,采集光纤陀螺输出数据进行精标定;四、使惯导位于固定位置东、北、天位置,采集光纤陀螺输出数据,标定光纤陀螺零偏误差并实时补偿;五、在东、北、天坐标系下绕X轴以角速度ω旋转一周,进行第一次标定参数修正;六、在东、北、天坐标系下绕Y轴以角速度ω旋转一周,进行第二次标定参数修正;七、在东、北、天坐标系下绕Z轴以角速度ω旋转一周,进行第三次标定参数修正;八、通过第三次标定参数修正,得到斜置光纤陀螺高精度的标定参数零偏、标度因数、失准角结果。
Description
技术领域
本发明属于惯性导航技术领域,特别是涉及一种基于姿态观测的冗余惯导系统光纤陀螺系统级标定方法。
背景技术
惯性技术是一项涉及多学科的综合技术,它是惯性导航和惯性制导技术、惯性仪表技术、惯性测量技术以及有关系统和装置技术的统称。惯性导航系统依靠自身的惯性敏感元件,不依赖任何外界信息测量导航参数,因此它不受天然的或人为的干扰,具有很好的隐蔽性,是一种完全自主式的导航系统。
惯性导航系统标定是通过比较系统中惯性器件的输出和已知参考输入,确定一组参数使惯导系统输出与输入相吻合的过程,惯导系统标定的理论基础是系统辨识和参数估计,其目的是确定惯性器件组合的数学误差模型或误差数学的模型参数。惯导系统使用之前必须进行标定,对器件零偏、标度因数、安装失准角等参数进行补偿。冗余惯导系统中器件安装方式与三轴惯导系统有较大差异,传统的标定方法在冗余惯导系统中实现起来特别繁琐,而且精度较低。因此,新的适用于冗余系统的标定方法已成为必然需求。
惯导系统常用标定方法主要有:分立式标定法、模观测标定方法、系统级标定方法等。
分立式标定方法也称为基于转台标定方法,需要转台为系统提供标准输入信息,对转台精度要求较高,同时,分立式标定过程依赖转台,一般只能在实验室进行。
模观测标定方法是指基于惯导系统输入加速度、角速度激励的模分别和加速度计比力测量、陀螺角速度测量的模相等的原理,以输入加速度、角速度的模作为观测,计算惯导系统参数的方法。目前模观测标定计算采用迭代算法,其收敛性严重依赖标定参数初值。
系统级标定方法主要基于导航解算误差的原理:惯导系统进入导航状态之后,其参数误差(惯性器件参数误差、初始对准姿态误差,初始位置误差等)经由导航解算会传递到导航结果(位置、速度、姿态等)中去,表现为导航误差,如能获取导航误差的全部或部分信息,就可能对惯导系统参数做出估计。系统级标定方法降低了对转台的精度要求,利用低精度转台就可以达到较高的标定精度,因此是现场标定的理想方法。
系统级标定方法相对于其他标定方法拥有较大优势,在现场标定和高精度标定的场合,系统级标定将占据重要地位。在冗余系统或者特定场合(如空间用单表陀螺等)中,光纤陀螺并不一定按照笛卡尔坐标系正交安装,而是采用特定的斜置安装方式以满足特定需求、提高系统可靠性和精度,传统意义上标定方法相对繁琐、精度低甚至不再适用,因此,斜置光纤陀螺高精度系统级标定方法拥有重大需求。
发明内容
本发明的目的在于提高冗余型光纤捷联惯性导航系统斜置光纤陀螺初始标定精度,提供了一种基于姿态观测的冗余惯导系统光纤陀螺系统级标定方法,它是一种适用于冗余光纤惯性导航系统斜置光纤陀螺系统级精标定方法。
本发明一种基于姿态观测的冗余惯导系统光纤陀螺系统级标定方法,该方法体步骤如下:
步骤一:将惯导系统安装在转台上,确定载体的初始位置参数,包括经度、纬度等;
步骤二:确定光纤陀螺轴向与惯导系统本体坐标系安装关系即安装角,计算安装矩阵;
步骤三:惯导系统预热,在已有光纤陀螺粗略标定参数(零偏、标度因数、失准角等)基础上(粗标定完成),准备采集光纤陀螺输出数据进行精标定;
光纤陀螺输出的数据为载体相对于惯性参考系的角速度
步骤四:使惯导系统位于某一固定位置(如东、北、天位置),采集光纤陀螺输出数据,以转台姿态角与惯导系统姿态角之间误差作为观测量,标定光纤陀螺零偏误差并实时补偿;
步骤五:在东、北、天坐标系下绕X轴以角速度ω(角速度ω应远大于地球转速,忽略地球转速影响)旋转一周,以转台姿态角与惯导系统姿态角之间的误差作为观测量,进行第一次标定参数修正;
步骤六:在东、北、天坐标系下绕Y轴以角速度ω(角速度ω应远大于地球转速,忽略地球转速影响)旋转一周,以转台姿态角与惯导系统姿态角之间的误差作为观测量,进行第二次标定参数修正;
步骤七:在东、北、天坐标系下绕Z轴以角速度ω(角速度ω应远大于地球转速,忽略地球转速影响)旋转一周,以转台姿态角与惯导系统姿态角之间的误差作为观测量,进行第三次标定参数修正;
步骤八:通过步骤七的第三次标定参数修正,得到斜置光纤陀螺高精度的标定参数(零偏、标度因数、失准角)结果。
其中,步骤四至步骤七中采用基于卡尔曼滤波技术的误差标定方法,利用转台姿态角与系统姿态角之间的误差作为观测量,通过卡尔曼滤波迭代,估计光纤陀螺零偏误差、标度因数误差及失准角,对光纤陀螺粗标定结果进行修正。具体步骤如下:
步骤一:建立光纤陀螺标定的系统状态方程和观测方程。
若以冗余系统中所有光纤陀螺标定参数构建卡尔曼滤波器,最终只能得到标定参数最小二乘解,而非真实解。此处采用任意三个编号为a,b,c的光纤陀螺和加速度计构建一套子惯导系统。
子惯导系统中光纤陀螺误差项系统级精标定卡尔曼滤波器的状态方程为:
其中15维状态矢量ψT为俯仰、横滚及航向姿态角误差;ΔBg T为光纤陀螺零偏残差矢量;ΔBg T=[ΔBga,ΔBgb,ΔBgc]T;为光纤陀螺标度因数误差残差矢量; 均为光纤陀螺安装失准角残差矢量: 表示系统状态的微分,状态方程中15阶方阵即状态矩阵,表示为如下形式:
其中Ag1和Ag2可表示为如下形式:
表示地球自转角速度在导航坐标系n系下的投影(下标的x,y及z表示沿导航坐标系的三个坐标轴)。为子惯导系统的配置矩阵
其中hi=[cos(αi)cos(βi)]·i+[sin(αi)cos(βi)]·j+[sin(βi)]·k,i=a,b,c,这里i,j和k表示轴Xb,Yb和Zb上的单位矢量,αi表示hi在Xb-Yb平面上的投影向量与Xb轴的夹角,βi表示hi与Xb-Yb平面的夹角(如图1),对角线矩阵 表示载体系(b系)相对惯性系(i系)角速度在b系下的投影(即是陀螺输出角速度值)。
其中pi=[sin(αi)cos(βi) -cos(αi)cos(βi) 0]。
其中qi=[-cos(αi)sin(βi) -sin(αi)sin(βi) cos(βi)](如图1、图2所示)。为转台所示捷联姿态矩阵。
假设Θg为系统噪声方差阵,式(1)中Wg为服从正态分布N(0,Θg)的系统噪声,满足如下条件:
以转台姿态角与惯导系统姿态角之间误差作为观测量构建卡尔曼滤波器的量测方程,形式如下:
Zg=FgXg+Vg (6)
上式中状态矢量Xg的定义与式(1)相同,观测量其中θx,θy,θz为系统解算姿态,为转台姿态(以姿态作为观测量)。量测矩阵Fg为15阶方阵,可表示为如下形式:
假设Rg为量测噪声方差阵,式(6)中Vg为服从正态分布N(0,Rg)的系统噪声,满足如下条件:
步骤二:对系统状态方程进行离散化。
对步骤一建立的系统状态变量进行估计,需要对系统状态方程进行离散化。离散化采用泰勒级数展开,则:
其中:Φ(k+1,k)为状态一步转移矩阵、I为十五阶单位阵、A(k)为状态转移矩阵,T为滤波周期。
系统模型噪声的方差为:
其中:Q(k)为离散系统噪声方差阵、Q连续系统噪声方程强度阵、A为状态转移矩阵。
步骤三:进行卡尔曼滤波状态估计。
对卡尔曼滤波器进行迭代,状态预测估计方程、方差预测方程、状态预测估计方程、方差迭代方程以及滤波增益方程可表示为如下形式:
Pk=(I-KkHk)Pk/k-1 (14)
最后估计得到冗余系统中斜置光纤陀螺零偏误差、标度因数误差以及失准角的标定结果,并对粗标定结果进行修正。
本发明的优点在于:
本发明所述方法在分立式标定(简称粗标定,粗标定之后光纤陀螺标定参数尚有残差)基础上,不增加传统标定程式,通过转动序列、滤波器设计的参数等的合理设计,对光纤陀螺零偏误差、标度因数误差及失准角进行进一步补偿,能够很大程度上提高斜置光纤陀螺标定精度。
附图说明
图1为理想传感器轴向与系统本体坐标系的安装关系示意图。
图2为实际传感器轴向与系统本体坐标系的安装关系示意图。
图3为系统级标定仿真平台示意图。
图4为六冗余RFINS的正十二面体安装方式示意图。
图5为系统级标定仿真过程中的子惯导系统示意图。
图6(a)为陀螺仪零偏误差随时间变化曲线示意图。
图6(b)为陀螺仪标度因数误差随时间变化曲线示意图。
图6(c)为陀螺仪失准角误差A随时间变化曲线示意图。
图6(d)为陀螺仪失准角误差B随时间变化曲线示意图。
图7为本发明流程框图。
图中符号说明如下:βi表示hi与Xb-Yb平面的夹角
OXbYbZb为系统本体坐标系;
Hi为理想传感器轴向,αi为Hi在XbOYb平面投影与Xb轴的夹角;βi表示hi与Xb-Yb平面的夹角;
Hi'为实际传感器轴向,δαi为Hi'与Hi夹角在XbOYb平面的投影,δβi为Hi'和Zb轴夹角与Hi和Zb轴夹角之差;
M1、M2、M3分别为XbOYb平面、XbOZb平面、YbOZb平面,α角为各平面内(M1、M2、M3)传感器轴向与最近系统本体坐标轴夹角,ABCDEF分别为六个传感器轴向;
ω表示系统绕坐标轴旋转角速度(角速度ω应远大于地球转速)。
具体实施方式
见图7,本发明一种基于姿态观测的冗余惯导系统光纤陀螺系统级标定方法,包括以下几个步骤:
步骤一:将光纤捷联惯组安装在转台上,确定载体的初始位置参数,包括经度、纬度等;
步骤二:确定光纤陀螺轴向与系统本体坐标系安装关系即安装角,计算安装矩阵;
步骤三:光纤捷联惯组预热,在已有光纤陀螺粗略标定参数(零偏、标度因数、失准角等)基础上(粗标定完成),准备采集光纤陀螺输出数据进行精标定;
光纤陀螺输出的数据为载体相对于惯性参考系的角速度
步骤四:使系统位于东、北、天位置,采集光纤陀螺输出数据,以转台姿态角与系统姿态角之间误差作为观测量,标定光纤陀螺零偏误差并实时补偿;
步骤五:在东、北、天坐标系下绕X轴以角速度ω(角速度ω应远大于地球转速,忽略地球转速影响)旋转一周,以转台姿态角与系统姿态角之间的误差作为观测量,进行第一次标定参数修正;
步骤六:在东、北、天坐标系下绕Y轴以角速度ω(角速度ω应远大于地球转速,忽略地球转速影响)旋转一周,以转台姿态角与系统姿态角之间的误差作为观测量,进行第二次标定参数修正;
步骤七:在东、北、天坐标系下绕Z轴以角速度ω(角速度ω应远大于地球转速,忽略地球转速影响)旋转一周,以转台姿态角与系统姿态角之间的误差作为观测量,进行第三次标定参数修正;
步骤八:通过步骤七的第三次标定参数修正,得到斜置光纤陀螺高精度的标定参数(零偏、标度因数、失准角)结果。
步骤四至步骤七中采用基于卡尔曼滤波技术的误差标定方法,利用转台姿态角与系统姿态角之间的误差作为观测量,通过卡尔曼滤波迭代,估计光纤陀螺零偏误差、标度因数误差及失准角,对光纤陀螺粗标定结果进行修正。具体步骤如下:
步骤一:建立光纤陀螺标定的系统状态方程和观测方程。
若以冗余系统中所有光纤陀螺标定参数构建卡尔曼滤波器,最终只能得到标定参数最小二乘解,而非真实解。此处采用任意三个编号为a,b,c的光纤陀螺和加速度计构建一套子惯导系统。
子惯导系统中光纤陀螺误差项系统级精标定卡尔曼滤波器的状态方程为:
其中15维状态矢量包括俯仰、横滚及航向姿态角误差ψT;光纤陀螺零偏残差矢量:ΔBg=[ΔBga,ΔBgb,ΔBgc]T,光纤陀螺标度因数误差残差矢量:光纤陀螺安装失准角残差矢量:
表示系统状态的微分,状态方程中15阶方阵(状态矩阵)可表示为如下形式:
其中Ag1和Ag2可表示为如下形式:
表示地球自转角速度在导航坐标系n系下的投影(下标的x,y及z表示沿导航坐标系的三个坐标轴)。为子惯导系统的配置矩阵
其中hi=[cos(αi)cos(βi)]·i+[sin(αi)cos(βi)]·j+[sin(βi)]·k(i=a,b,c),这里i,j和k表示轴Xb,Yb和Zb上的单位矢量,αi表示hi在Xb-Yb平面上的投影向量与Xb轴的夹角,βi表示hi与Xb-Yb平面的夹角(如图1),对角线矩阵 表示载体系(b系)相对惯性系(i系)角速度在b系下的投影(即是陀螺输出角速度值)。
其中pi=[sin(αi)cos(βi) -cos(αi)cos(βi) 0]。
其中qi=[-cos(αi)sin(βi) -sin(αi)sin(βi) cos(βi)](如图1、图2所示)。为转台所示捷联姿态矩阵。
假设Θg为系统噪声方差阵,式(1)中Wg为服从正态分布N(0,Θg)的系统噪声,满足如下条件:
以转台姿态角与惯导系统姿态角之间误差作为观测量构建卡尔曼滤波器的量测方程,形式如下:
Zg=FgXg+Vg (6)
上式中状态矢量Xg的定义与式(1)相同,观测量其中θx,θy,θz为系统解算姿态,为转台姿态(以姿态作为观测量)。量测矩阵Fg为15阶方阵,可表示为如下形式:
假设Rg为量测噪声方差阵,式(6)中Vg为服从正态分布N(0,Rg)的系统噪声,满足如下条件:
步骤二:对系统状态方程进行离散化。
对步骤一建立的系统状态变量进行估计,需要对系统状态方程进行离散化。离散化采用泰勒级数展开,则:
其中:Φ(k+1,k)为状态一步转移矩阵、I为十五阶单位阵、A(k)为状态转移矩阵,T为滤波周期。
系统模型噪声的方差为:
其中:Q(k)为离散系统噪声方差阵、Q连续系统噪声方程强度阵、A为状态转移矩阵。
步骤三:进行卡尔曼滤波状态估计。
对卡尔曼滤波器进行迭代,状态预测估计方程、方差预测方程、状态预测估计方程、方差迭代方程以及滤波增益方程可表示为如下形式:
Pk=(I-KkHk)Pk/k-1 (14)
最后估计得到冗余系统中斜置光纤陀螺零偏误差、标度因数误差以及失准角的标定结果,并对粗标定结果进行修正。
实施例仿真:
下面结合实施例及标定仿真对本发明做进一步说明,本发明实施例以某六冗余型光纤捷联惯导系统的标定为例。
(1)系统级标定仿真平台
系统级标定仿真平台是在分立式标定仿真平台的基础上加入光纤陀螺误差精标定模块,包括光纤陀螺误差精标定卡尔曼滤波器。标定仿真平台结构框图如图3所示。
(2)系统级标定仿真条件及结论
六冗余RFINS的正十二面体安装方式由图4所示,六冗余斜置RFINS系统结构中的ABC光纤陀螺与加速度计构成,如图5所示。其中AB轴位于面M2,C轴位于面M1,其与坐标轴的夹角为α=31°43'2.9”。
根据图5的配置方式,可得系统的安装矩阵为:
子惯导系统中光纤陀螺的零偏、标度及安装失准角等各项标定参数如下表所示,仿真过程中采样周期为20ms。
表1冗余光纤陀螺标定参数
光纤陀螺 | |
零偏 | 0.1°/h |
标度因数(ppm) | 100 |
失准角A(角分) | 1 |
失准角B(角分) | 1 |
噪声方差 | 0.005°/h |
在光纤陀螺误差系统级标定过程中,利用前述精标定卡尔曼滤波器以及转台转动序列。系统每次转动过程持续5s,转动完成后静止1min。下图描述了子惯导系统中光纤陀螺的零偏误差、标度因数误差和两类安装失准角误差的估计曲线。
从图6a-d可以看出,当以系统姿态与转台姿态的差值作为观测量时(系统速度置0),由于系统静止过程中光纤陀螺标度因数及安装失准角不被激发,因此在第一个位置可以较为准确地标定陀螺零偏误差;当系统分别绕本体坐标系转动一周后,光纤陀螺的标度因数误差和安装失准角误差可以准确标定;当系统转至第三个位置时A类安装失准角即可准确标定,是由于系统从绕天向转动时A类失准角未影响系统姿态误差,而到绕东向转动时A类失准角将作用于系统姿态误差,而此时标度因数误差和B类安装失准角误差已经到达非真值的稳态,所以A类安装失准角可以准确求得。
Claims (1)
1.一种基于姿态观测的冗余惯导系统光纤陀螺系统级标定方法,其特征在于:该方法体步骤如下:
步骤一:将惯导系统安装在转台上,确定载体的初始位置参数,包括经度、纬度;
步骤二:确定光纤陀螺轴向与惯导系统本体坐标系安装关系即安装角,计算安装矩阵;
步骤三:惯导系统预热,在已有光纤陀螺粗略标定参数零偏、标度因数、失准角基础上,准备采集光纤陀螺输出数据进行精标定;
步骤四:使惯导系统位于某一固定位置东、北、天位置,采集光纤陀螺输出数据,以转台姿态角与惯导系统姿态角之间误差作为观测量,标定光纤陀螺零偏误差并实时补偿;
步骤五:在东、北、天坐标系下绕X轴以角速度ω旋转一周,以转台姿态角与惯导系统姿态角之间的误差作为观测量,进行第一次标定参数修正;
步骤六:在东、北、天坐标系下绕Y轴以角速度ω旋转一周,以转台姿态角与惯导系统姿态角之间的误差作为观测量,进行第二次标定参数修正;
步骤七:在东、北、天坐标系下绕Z轴以角速度ω旋转一周,以转台姿态角与惯导系统姿态角之间的误差作为观测量,进行第三次标定参数修正;
步骤八:通过步骤七的第三次标定参数修正,得到斜置光纤陀螺高精度的标定参数零偏、标度因数、失准角结果;
其中,步骤四至步骤七中采用基于卡尔曼滤波技术的误差标定方法,利用转台姿态角与系统姿态角之间的误差作为观测量,通过卡尔曼滤波迭代,估计光纤陀螺零偏误差、标度因数误差及失准角,对光纤陀螺粗标定结果进行修正,具体步骤如下:
步骤一:建立光纤陀螺标定的系统状态方程和观测方程;
若以冗余系统中所有光纤陀螺标定参数构建卡尔曼滤波器,最终只能得到标定参数最小二乘解,而非真实解,此处采用任意三个编号为a,b,c的光纤陀螺和加速度计构建一套子惯导系统;
子惯导系统中光纤陀螺误差项系统级精标定卡尔曼滤波器的状态方程为:
其中15维状态矢量ψT为俯仰、横滚及航向姿态角误差;ΔBg T为光纤陀螺零偏残差矢量;ΔBg T=[ΔBga,ΔBgb,ΔBgc]T;为光纤陀螺标度因数误差残差矢量; 均为光纤陀螺安装失准角残差矢量: 表示系统状态的微分,状态方程中15阶方阵即状态矩阵,表示为如下形式:
其中Ag1和Ag2表示为如下形式:
表示地球自转角速度在导航坐标系n系下的投影,下标的x,y及z表示沿导航坐标系的三个坐标轴;为子惯导系统的配置矩阵
其中hi=[cos(αi)cos(βi)]·i+[sin(αi)cos(βi)]·j+[sin(βi)]·k,i=a,b,c,这里i,j和k表示轴Xb,Yb和Zb上的单位矢量,αi表示hi在Xb-Yb平面上的投影向量与Xb轴的夹角,βi表示hi与Xb-Yb平面的夹角,对角线矩阵 表示载体系即b系相对惯性系即i系角速度在b系下的投影,即是陀螺输出角速度值;
其中pi=[sin(αi)cos(βi) -cos(αi)cos(βi) 0];
其中qi=[-cos(αi)sin(βi) -sin(αi)sin(βi) cos(βi)],为转台所示捷联姿态矩阵;
假设Θg为系统噪声方差阵,式(1)中Wg为服从正态分布N(0,Θg)的系统噪声,满足如下条件:
以转台姿态角与惯导系统姿态角之间误差作为观测量构建卡尔曼滤波器的量测方程,形式如下:
Zg=FgXg+Vg (6)
上式中状态矢量Xg的定义与式(1)相同,观测量其中θx,θy,θz为系统解算姿态,为转台姿态;量测矩阵Fg为15阶方阵,表示为如下形式:
假设Rg为量测噪声方差阵,式(6)中Vg为服从正态分布N(0,Rg)的系统噪声,满足如下条件:
步骤二:对系统状态方程进行离散化;
对步骤一建立的系统状态变量进行估计,需要对系统状态方程进行离散化;离散化采用泰勒级数展开,则:
其中:Φ(k+1,k)为状态一步转移矩阵、I为十五阶单位阵、A(k)为状态转移矩阵,T为滤波周期;
系统模型噪声的方差为:
其中:Q(k)为离散系统噪声方差阵、Q连续系统噪声方程强度阵、A为状态转移矩阵;
步骤三:进行卡尔曼滤波状态估计;
对卡尔曼滤波器进行迭代,状态预测估计方程、方差预测方程、状态预测估计方程、方差迭代方程以及滤波增益方程表示为如下形式:
Pk=(I-KkHk)Pk/k-1 (14)
最后估计得到冗余系统中斜置光纤陀螺零偏误差、标度因数误差以及失准角的标定结果,并对粗标定结果进行修正。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410599074.1A CN104344836B (zh) | 2014-10-30 | 2014-10-30 | 一种基于姿态观测的冗余惯导系统光纤陀螺系统级标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410599074.1A CN104344836B (zh) | 2014-10-30 | 2014-10-30 | 一种基于姿态观测的冗余惯导系统光纤陀螺系统级标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104344836A CN104344836A (zh) | 2015-02-11 |
CN104344836B true CN104344836B (zh) | 2017-04-19 |
Family
ID=52500940
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410599074.1A Active CN104344836B (zh) | 2014-10-30 | 2014-10-30 | 一种基于姿态观测的冗余惯导系统光纤陀螺系统级标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104344836B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106839877B (zh) * | 2015-12-07 | 2018-08-24 | 上海机电工程研究所 | 一种光学多波段复合目标模拟系统的误差补偿方法 |
CN106884645A (zh) * | 2015-12-16 | 2017-06-23 | 航天科工惯性技术有限公司 | 陀螺测斜仪的标定方法 |
CN106767915A (zh) * | 2016-12-07 | 2017-05-31 | 上海航天控制技术研究所 | 一种运载火箭带冗余斜轴的光纤惯组标定方法 |
CN106969783B (zh) * | 2017-02-13 | 2021-03-26 | 哈尔滨工业大学 | 一种基于光纤陀螺惯性导航的单轴旋转快速标定技术 |
CN107655493B (zh) * | 2017-09-06 | 2021-04-06 | 东南大学 | 一种光纤陀螺sins六位置系统级标定方法 |
WO2019085526A1 (zh) * | 2017-11-03 | 2019-05-09 | 北京凌宇智控科技有限公司 | 面向三维空间的定位校正方法、组合定位方法及其设备 |
CN111561948B (zh) * | 2019-12-05 | 2023-07-28 | 北京计算机技术及应用研究所 | 四轴冗余捷联惯导的系统级标定方法 |
CN111314607B (zh) * | 2020-02-21 | 2021-07-13 | 昆山丘钛微电子科技有限公司 | 一种光学防抖的标定方法、装置、电子设备及介质 |
CN111486870B (zh) * | 2020-04-23 | 2021-01-05 | 中南大学 | 一种斜置捷联惯组系统级标定方法 |
CN111721288B (zh) * | 2020-06-19 | 2022-03-29 | 哈尔滨工业大学 | 一种mems器件零偏修正方法、装置及存储介质 |
CN115060258B (zh) * | 2022-08-18 | 2022-11-01 | 中国船舶重工集团公司第七0七研究所 | 一种基于谐振惯导系统的卡尔曼滤波精对准方法 |
CN116026330B (zh) * | 2023-03-30 | 2023-06-16 | 中国船舶集团有限公司第七〇七研究所 | 基于光纤陀螺数字信号的三轴旋转式框架施矩方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE1623557B2 (de) * | 1966-08-09 | 1977-09-08 | Sperry Rand Corp., New York, N.Y. (V.St.A.) | Redundantes kreiselsystem mit drei kreiselgeraeten fuer luftfahrzeuge |
US5479161A (en) * | 1994-03-25 | 1995-12-26 | Honeywell Inc. | Automatic calibration of redundant sensors |
JP2002090174A (ja) * | 2000-09-19 | 2002-03-27 | Toshiba Corp | 冗長航法装置及び冗長航法方法 |
CN102735232A (zh) * | 2012-06-26 | 2012-10-17 | 北京航天时代光电科技有限公司 | 光纤陀螺组合体惯性测量装置及其标定方法 |
CN102927994A (zh) * | 2012-10-23 | 2013-02-13 | 北京航空航天大学 | 一种斜置冗余捷联惯性导航系统的快速标定方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7106189B2 (en) * | 2004-04-29 | 2006-09-12 | Tracetech Incorporated | Tracking system and methods thereof |
-
2014
- 2014-10-30 CN CN201410599074.1A patent/CN104344836B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE1623557B2 (de) * | 1966-08-09 | 1977-09-08 | Sperry Rand Corp., New York, N.Y. (V.St.A.) | Redundantes kreiselsystem mit drei kreiselgeraeten fuer luftfahrzeuge |
US5479161A (en) * | 1994-03-25 | 1995-12-26 | Honeywell Inc. | Automatic calibration of redundant sensors |
JP2002090174A (ja) * | 2000-09-19 | 2002-03-27 | Toshiba Corp | 冗長航法装置及び冗長航法方法 |
CN102735232A (zh) * | 2012-06-26 | 2012-10-17 | 北京航天时代光电科技有限公司 | 光纤陀螺组合体惯性测量装置及其标定方法 |
CN102927994A (zh) * | 2012-10-23 | 2013-02-13 | 北京航空航天大学 | 一种斜置冗余捷联惯性导航系统的快速标定方法 |
Non-Patent Citations (5)
Title |
---|
一种余度捷联惯导系统标定结果的验证方法;芦佳振等;《系统工程与电子技术》;20081115(第11期);2245-2248 * |
多惯组冗余系统安装误差的空中标定技术;张金亮等;《传感技术学报》;20130331;第26卷(第3期);367-371 * |
捷联惯性导航系统多传感器冗余技术研究;刘莎红;《中国硕士学位论文全文数据库信息科技辑》;20130215(第2期);48-77 * |
斜置惯性测量单元的一体化标定技术;郭鹏飞等;《中国惯性技术学报》;20070615(第03期);377-381 * |
空间冗余陀螺仪在轨标定方法;杨华波等;《宇航学报》;20100131(第01期);104-110 * |
Also Published As
Publication number | Publication date |
---|---|
CN104344836A (zh) | 2015-02-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104344836B (zh) | 一种基于姿态观测的冗余惯导系统光纤陀螺系统级标定方法 | |
CN104344837B (zh) | 一种基于速度观测的冗余惯导系统加速度计系统级标定方法 | |
CN101706284B (zh) | 提高船用光纤陀螺捷联惯导系统定位精度的方法 | |
CN105371844B (zh) | 一种基于惯性/天文互助的惯性导航系统初始化方法 | |
CN104457754B (zh) | 一种基于sins/lbl紧组合的auv水下导航定位方法 | |
CN101788296B (zh) | 一种sins/cns深组合导航系统及其实现方法 | |
CN106767900A (zh) | 一种基于组合导航技术的船用光纤捷联惯导系统的在线标定方法 | |
CN104019828A (zh) | 高动态环境下惯性导航系统杆臂效应误差在线标定方法 | |
CN109974697A (zh) | 一种基于惯性系统的高精度测绘方法 | |
CN106885570A (zh) | 一种基于鲁棒sckf滤波的紧组合导航方法 | |
CN101261130B (zh) | 一种船用光纤捷联惯导系统传递对准精度评估方法 | |
CN105180968A (zh) | 一种imu/磁强计安装失准角在线滤波标定方法 | |
CN103852085B (zh) | 一种基于最小二乘拟合的光纤捷联惯导系统现场标定方法 | |
CN103575299A (zh) | 利用外观测信息的双轴旋转惯导系统对准及误差修正方法 | |
CN101571394A (zh) | 基于旋转机构的光纤捷联惯性导航系统初始姿态确定方法 | |
CN103852086B (zh) | 一种基于卡尔曼滤波的光纤捷联惯导系统现场标定方法 | |
CN103076025B (zh) | 一种基于双解算程序的光纤陀螺常值误差标定方法 | |
CN103792561B (zh) | 一种基于gnss通道差分的紧组合降维滤波方法 | |
CN104374401A (zh) | 一种捷联惯导初始对准中重力扰动的补偿方法 | |
CN103245359A (zh) | 一种惯性导航系统中惯性传感器固定误差实时标定方法 | |
CN106940193A (zh) | 一种基于Kalman滤波的船舶自适应摇摆标定方法 | |
CN102519485B (zh) | 一种引入陀螺信息的二位置捷联惯性导航系统初始对准方法 | |
CN101915579A (zh) | 一种基于ckf的sins大失准角初始对准新方法 | |
CN112595350B (zh) | 一种惯导系统自动标定方法及终端 | |
CN112146655B (zh) | 一种BeiDou/SINS紧组合导航系统弹性模型设计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |