CN101187561A - 适合于光纤陀螺的载体姿态测量方法 - Google Patents
适合于光纤陀螺的载体姿态测量方法 Download PDFInfo
- Publication number
- CN101187561A CN101187561A CNA2007101448462A CN200710144846A CN101187561A CN 101187561 A CN101187561 A CN 101187561A CN A2007101448462 A CNA2007101448462 A CN A2007101448462A CN 200710144846 A CN200710144846 A CN 200710144846A CN 101187561 A CN101187561 A CN 101187561A
- Authority
- CN
- China
- Prior art keywords
- omega
- carrier
- optical fibre
- renewal cycle
- posture renewal
- 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
Images
Landscapes
- Gyroscopes (AREA)
Abstract
本发明提供的是一种适合于光纤陀螺的载体姿态测量方法。包括:通过外部设备确定载体的初始位置参数与初始速度值;光纤陀螺捷联惯性导航***进行初始对准,确定载体相对导航坐标系的初始姿态,得到姿态四元数的初始值;确定姿态更新周期H=tm-tm-1;采集光纤陀螺输出的载体相对于惯性坐标系的角速度计算旋转矢量的增量Δφ;通过旋转矢量与四元数的关系,得到姿态更新周期H内姿态更新四元数q(H);由姿态四元数更新方程更新姿态四元数;计算载体坐标系b系相对于导航坐标系n系的捷联矩阵T;求载体相对导航坐标系的姿态角等步骤。本发明解决了在载体高动态环境或是高频率振动环境中,圆锥效应对于载体姿态测量精度产生影响的问题。
Description
(一)技术领域
本发明提供的是一种测量方法,具体地说是一种使用光纤陀螺测量载体相对于导航坐标系姿态参数的方法。
(二)背景技术
在捷联惯性导航***中,惯性传感器(陀螺和加速度计)直接固联在载体上。它们敏感载体相对于惯性坐标系的姿态变化及速度变化信息,这些信息经过导航计算机的处理输出导航参数,从而完成导航制导任务。由于采用惯性传感器固联的方式,惯性传感器会直接敏感载体的振动和环境的干扰,对捷联惯性导航***的测量精度带来一定的影响。在捷联惯性导航***中,采用的是数学平台,主要依靠陀螺的输出信息来构建数学平台,以描述载体相对导航坐标系的姿态。而陀螺的输出为角增量或角速度,不能直接用于描述坐标系的相对关系。在常规的捷联惯性导航***的姿态测量方法中,总是认为载体坐标系和导航坐标系间的转换是通过一系列的转动来实现的。所以,为了测量载体的姿态,需要由相应的角速度或角增量来求得描述坐标系关系的四元数或旋转矢量。但是在高动态环境下,在姿态更新周期内的每次陀螺采样对应的载体坐标系的位置是不同的。另外,传统方法中是假设从载体坐标系到参考坐标系间的转动的次序是可以不予考虑的,这是基于无限小转动是矢量的原理来得到的。但实际的工程中,特别是在载体高动态机动情形或恶劣振动环境下,载体的转动常常是有限转动,而有限转动不是矢量,其转动次序不能交换。因而将姿态更新周期中的所有陀螺输出的角增量求和后再求解四元数的传统姿态测量方法将带来较大的误差,即圆锥误差。圆锥误差类似于陀螺漂移,给载体的姿态测量带来负面影响。补偿圆锥误差,高质量地测量载体的姿态就必须设计高性能的姿态测量方法。
在传统的姿态测量方法中,利用陀螺输出的角增量去估算出圆锥补偿项,从而估算出旋转向量的增量。该基本思路是基于在传统陀螺的输出中,角增量的获取是很方便的。即要想得到旋转向量的增量必须知道陀螺输出的角增量。而在目前广泛使用的光纤陀螺捷联惯性导航***中,光纤陀螺的输出信号为角速度信息,为此,通常的做法是采用分段线性化,认为在每一个采样周期内载体的角速度是常量。这种假设在采样周期足够小且载体的角速度变化较慢的时候是成立的,但对于高动态的场合就必须另加考虑了。已有的研究成果表明,在光纤陀螺输出角速度的情形下,对于传统的基于旋转矢量的姿态测量方法来说,性能不如常规的基于四元数的姿态测量方法,无法体现在应对圆锥效应的优越性。
(三)发明内容
本发明的目的在于提供一种能够解决在载体高动态环境或是高频率振动环境中,圆锥效应对于载体姿态测量精度产生影响的问题的适合于光纤陀螺的载体姿态测量方法。
本发明的目的是这样实现的:
步骤1、通过外部设备确定载体的初始位置参数与初始速度值;
步骤2、光纤陀螺捷联惯性导航***进行初始对准,确定载体相对导航坐标系的初始姿态,得到姿态四元数的初始值;
步骤3、确定姿态更新周期H=tm-tm-1,所述姿态更新周期H等于N倍的光纤陀螺采样周期h=tl-tl-1,所述N为大于零的整数;
步骤4、采集光纤陀螺输出的载体相对于惯性坐标系的角速度计算旋转矢量的增量Δφ;
步骤5、通过旋转矢量与四元数的关系,得到姿态更新周期H内姿态更新四元数q(H),
其中为旋转矢量增量的模,
步骤6、由姿态四元数更新方程更新姿态四元数
Q(tm)=Q(tm-1)q(H)
其中姿态更新四元数q(H)由步骤5得到,Q(tm)、Q(tm-1)分别表示载体在tm、tm-1时刻的姿态四元数;
步骤7、利用步骤6获得的tm时刻姿态四元数Q(tm)=[q0q1q2q3]T计算载体坐标系b系相对于导航坐标系n系的捷联矩阵T,
步骤8、利用步骤7获得的捷联矩阵T求得载体相对导航坐标系的姿态角,导航坐标系选取为游动方位坐标系,则载体的格网航向角ψ、纵摇角θ、横摇角γ都可由捷联矩阵T中的元素来表示,
θ0=arcsinT32
其中格网航向角ψ的定义域为(0°,360°);纵摇角θ的定义域为(-90°,90°);横摇角γ的定义域为(-90°,90°),
由姿态角的定义域得到它们的真值,
格网航向角的真值为
纵摇角的真值为
θ=θ0
横摇角的真值为
本发明还可以包括:
1、在步骤1中所述的载体的初始位置由GPS装置或者是外部的高精度组合导航设备提供,所述载体的初始速度值由DVL多普勒计程仪或者是外部的高精度组合导航设备提供。
2、在步骤2中所述的载体相对导航坐标系的初始姿态确定分为两种情形:当载体处于静基座时,采用基于古典控制理论的自对准方法;当载体处于动基座时,采用基于卡尔曼滤波理论的组合对准方法。
3、在步骤4中所述的旋转矢量增量Δφ的计算具体步骤如下:
步骤41:由光纤陀螺输出的输出计算Δφ中的惯性可测项,即角增量α;由姿态更新周期和光纤陀螺采样周期的关系H=Nh;每个姿态更新周期H内已知N+1个光纤陀螺采样角速度ω0,ω1,L ωN;则姿态更新周期H内的角增量α为
其中
步骤42:由光纤陀螺输出的输出计算Δφ中的非惯性可测项,即圆锥补偿项β;由姿态更新周期和光纤陀螺采样周期的关系H=Nh;每个姿态更新周期H内已知N+1个光纤陀螺采样角速度ω0,ω1,L ωN;
在姿态更新周期H内,圆锥补偿项β由光纤陀螺采样角速度叉乘项的线性组合来计算
其中KN-i为优化系数,在典型圆锥环境中确定;
步骤43:将步骤41中得到的姿态更新周期H内角增量α,与步骤41中得到的姿态更新周期H内圆锥补偿项β相加,得到姿态更新周期H内旋转矢量增量Δφ;
Δφ=α+β
设置N=3,陀螺在t时刻采样角速度ω0,并且每个姿态更新周期H,在t+H/3、t+2H/3、t+H时刻采样ω1、ω2、ω3;得到姿态更新周期H内旋转矢量增量Δφ为
其中
设置N=2,陀螺在t时刻采样角速度ω0,并且每个计算周期H内,在t+H/2、t+H时刻采样两次ω1、ω2;得到姿态更新周期H内旋转矢量增量Δφ为
其中
本发明与传统方法相比较所的优点主要体现在:
(1)针对光纤陀螺输出为角速度的情形,不同于传统测量方法以角增量作为输入,发明方法直接以角速度输入,以角速度的叉乘项来拟和圆锥补偿项。避免了使用光纤陀螺时,传统方法对于划船补偿项的测量偏差。
(2)相对于提高捷联惯性导航***姿态计算精度的另一种途径:选用高性能的光纤陀螺。发明的测量方法只需要作导航软件设计,不需要增加捷联惯性导航***的制造成本。
(四)附图说明
图1为本发明的适合光纤陀螺捷联惯性导航***的姿态测量方法流程图。
图2为典型圆锥环境中,发明的姿态测量方法;传统的姿态测量方法的实验曲线。典型圆锥环境定义为:载体的两个正交轴上存在同频率,相位差90度的角振动。光纤陀螺的采样频率为100Hz。
(五)具体实施方式
下面结合附图举例对本发明做更详细地描述:
结合图1,本发明的适合于光纤陀螺捷联惯性导航***的姿态测量方法包括如下步骤:
步骤1、通过外部设备确定载体的初始位置参数与初始速度值。
步骤2、光纤陀螺捷联惯性导航***进行初始对准,确定载体相对导航坐标系的初始姿态,得到姿态四元数的初始值。
步骤3、确定姿态更新周期H=tm-tm-1,所述姿态更新周期H等于N倍的光纤陀螺采样周期h=tl-tl-1。所述N为大于零的整数。
步骤4、采集光纤陀螺输出的载体相对于惯性坐标系的角速度计算旋转矢量的增量Δφ。
步骤5、由旋转矢量与四元数的关系,得到姿态更新周期H内姿态更新四元数q(H)。
其中为旋转矢量增量的模。
步骤6、由姿态四元数更新方程更新姿态四元数
Q(tm)=Q(tm-1)q(H) (2)
其中姿态更新四元数q(H)由步骤5求得。Q(tm)、Q(tm-1)分别表示载体在tm、tm-1时刻的姿态四元数。
步骤7、利用步骤6获得的tm时刻姿态四元数Q(tm)=[q0 q1 q2 q3]T计算载体坐标系b系相对于导航坐标系n系的捷联矩阵T。
步骤8、利用步骤7获得的捷联矩阵T求得载体相对导航坐标系的姿态角。导航坐标系选取为游动方位坐标系,则载体的格网航向角ψ、纵摇角θ、横摇角γ都可由捷联矩阵T中的元素来表示。
θ0=arcsinT32 (4)
其中格网航向角ψ的定义域为(0°,360°);纵摇角θ的定义域为(-90°,90°);横摇角γ的定义域为(-90°,90°)。
由姿态角的定义域可以确定它们的真值。格网航向角的真值为
纵摇角的真值为
θ=θ0
横摇角的真值为
在步骤1中所述的载体的初始位置由GPS装置或者是外部的高精度组合导航设备提供,所述载体的初始速度值由DVL多普勒计程仪或者是外部的高精度组合导航设备提供。
在步骤2中所述的载体相对导航坐标系的初始姿态确定分为两种情形:当载体处于静基座时,可以采用基于古典控制理论的自对准方法;当载体处于动基座时,可以采用基于卡尔曼滤波理论的组合对准方法。
在步骤4中所述的旋转矢量增量Δφ的计算具体步骤如下:
步骤41:由光纤陀螺输出的输出计算Δφ中的惯性可测项,即角增量α。由姿态更新周期和光纤陀螺采样周期的关系H=Nh。每个姿态更新周期H内已知N+1个光纤陀螺采样角速度ω0,ω1,LωN。则姿态更新周期H内的角增量α为
其中
步骤42:由光纤陀螺输出的输出计算Δφ中的非惯性可测项,即圆锥补偿项β。由姿态更新周期和光纤陀螺采样周期的关系H=Nh。每个姿态更新周期H内已知N+1个光纤陀螺采样角速度ω0,ω1,LωN。
在姿态更新周期H内,圆锥补偿项β可以由光纤陀螺采样角速度叉乘项的线性组合来计算
其中KN-i为优化系数,可以在典型圆锥环境中确定。
步骤43:将步骤41中得到的姿态更新周期H内角增量α,与步骤41中得到的姿态更新周期H内圆锥补偿项β相加,得到姿态更新周期H内旋转矢量增量Δφ。
Δφ=α+β (8)
设置N=3,陀螺在t时刻采样角速度ω0,并且每个姿态更新周期H,在t+H/3、t+2H/3、t+H时刻采样ω1、ω2、ω3。得到姿态更新周期H内旋转矢量增量Δφ为
其中
设置N=2,陀螺在t时刻采样角速度ω0,并且每个计算周期H内,在t+H/2、t+H时刻采样两次ω1、ω2。得到姿态更新周期H内旋转矢量增量Δφ为
其中
在姿态更新周期H内,姿态测量方法所使用的光纤陀螺采样数越多,姿态测量方法应对圆锥补偿的效果就越好,光纤陀螺捷联惯性导航***输出的姿态精度就越高。
Claims (5)
1.一种适合于光纤陀螺的载体姿态测量方法,其特征是:
步骤1、通过外部设备确定载体的初始位置参数与初始速度值;
步骤2、光纤陀螺捷联惯性导航***进行初始对准,确定载体相对导航坐标系的初始姿态,得到姿态四元数的初始值;
步骤3、确定姿态更新周期H=tm-tm-1,所述姿态更新周期H等于N倍的光纤陀螺采样周期h=tl-tl-1,所述N为大于零的整数;
步骤4、采集光纤陀螺输出的载体相对于惯性坐标系的角速度计算旋转矢量的增量Δφ;
步骤5、通过旋转矢量与四元数的关系,得到姿态更新周期H内姿态更新四元数q(H),
其中为旋转矢量增量的模,
步骤6、由姿态四元数更新方程更新姿态四元数
Q(tm)=Q(tm-1)q(H)
其中姿态更新四元数q(H)由步骤5得到,Q(tm)、Q(tm-1)分别表示载体在tm、tm-1时刻的姿态四元数;
步骤7、利用步骤6获得的tm时刻姿态四元数Q(tm)=[q0q1q2q3]T计算
载体坐标系b系相对于导航坐标系n系的捷联矩阵T,
步骤8、利用步骤7获得的捷联矩阵T求得载体相对导航坐标系的姿态角,导航坐标系选取为游动方位坐标系,则载体的格网航向角ψ、纵摇角θ、横摇角γ都可由捷联矩阵T中的元素来表示,
θ0=arcsinT32
其中格网航向角ψ的定义域为(0°,360°);纵摇角θ的定义域为(-90°,90°);横摇角γ的定义域为(-90°,90°),
由姿态角的定义域得到它们的真值,
格网航向角的真值为
纵摇角的真值为
θ=θ0
横摇角的真值为
2.根据权利要求1所述的适合于光纤陀螺的载体姿态测量方法,其特征是:在步骤1中所述的载体的初始位置由GPS装置或者是外部的高精度组合导航设备提供,所述载体的初始速度值由DVL多普勒计程仪或者是外部的高精度组合导航设备提供。
3.根据权利要求1或2所述的适合于光纤陀螺的载体姿态测量方法,其特征是:在步骤2中所述的载体相对导航坐标系的初始姿态确定分为两种情形:当载体处于静基座时,采用基于古典控制理论的自对准方法;当载体处于动基座时,采用基于卡尔曼滤波理论的组合对准方法。
4.根据权利要求1或2所述的适合于光纤陀螺的载体姿态测量方法,其特征是:在步骤4中所述的旋转矢量增量Δφ的计算具体步骤如下:
步骤41:由光纤陀螺输出的输出计算Δφ中的惯性可测项,即角增量α;由姿态更新周期和光纤陀螺采样周期的关系H=Nh;每个姿态更新周期H内已知N+1个光纤陀螺采样角速度ω0,ω1,LωN;则姿态更新周期H内的角增量α为
其中
步骤42:由光纤陀螺输出的输出计算Δφ中的非惯性可测项,即圆锥补偿项β;由姿态更新周期和光纤陀螺采样周期的关系H=Nh;每个姿态更新周期H内已知N+1个光纤陀螺采样角速度ω0,ω1,LωN;
在姿态更新周期H内,圆锥补偿项β由光纤陀螺采样角速度叉乘项的线性组合来计算
其中KN-1为优化系数,在典型圆锥环境中确定;
步骤43:将步骤41中得到的姿态更新周期H内角增量α,与步骤41中得到的姿态更新周期H内圆锥补偿项β相加,得到姿态更新周期片内旋转矢量增量Δφ;
Δφ=α+β
设置N=3,陀螺在t时刻采样角速度ω0,并且每个姿态更新周期H,在t+H/3、t+2H/3、t+H时刻采样ω1、ω2、ω3;得到姿态更新周期H内旋转矢量增量Δφ为
其中
设置N=2,陀螺在t时刻采样角速度ω0,并且每个计算周期H内,在t+H/2、t+H时刻采样两次ω1、ω2;得到姿态更新周期H内旋转矢量增量Δφ为
其中
5.根据权利要求3所述的适合于光纤陀螺的载体姿态测量方法,其特征是:在步骤4中所述的旋转矢量增量Δφ的计算具体步骤如下:
步骤41:由光纤陀螺输出的输出计算Δφ中的惯性可测项,即角增量α;由姿态更新周期和光纤陀螺采样周期的关系H=Nh;每个姿态更新周期H内已知N+1个光纤陀螺采样角速度ω0,ω1,L ωN;则姿态更新周期H内的角增量α为
其中
步骤42:由光纤陀螺输出的输出计算Δφ中的非惯性可测项,即圆锥补偿项β;由姿态更新周期和光纤陀螺采样周期的关系H=Nh;每个姿态更新周期H内已知N+1个光纤陀螺采样角速度ω0,ω1,L ωN;
在姿态更新周期H内,圆锥补偿项β由光纤陀螺采样角速度叉乘项的线性组合来计算
其中KN-1为优化系数,在典型圆锥环境中确定;
步骤43:将步骤41中得到的姿态更新周期H内角增量α,与步骤41中得到的姿态更新周期H内圆锥补偿项β相加,得到姿态更新周期H内旋转矢量增量Δφ;
Δφ=α+β
设置N=3,陀螺在t时刻采样角速度ω0,并且每个姿态更新周期H,在t+H/3、t+2H/3、t+H时刻采样ω1、ω2、ω3;得到姿态更新周期H内旋转矢量增量Δφ为
其中
设置N=2,陀螺在t时刻采样角速度ω0,并且每个计算周期H内,在t+H/2、t+H时刻采样两次ω1、ω2;得到姿态更新周期H内旋转矢量增量Δφ为
其中
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200710144846A CN100588906C (zh) | 2007-12-18 | 2007-12-18 | 适合于光纤陀螺的载体姿态测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200710144846A CN100588906C (zh) | 2007-12-18 | 2007-12-18 | 适合于光纤陀螺的载体姿态测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101187561A true CN101187561A (zh) | 2008-05-28 |
CN100588906C CN100588906C (zh) | 2010-02-10 |
Family
ID=39480052
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200710144846A Expired - Fee Related CN100588906C (zh) | 2007-12-18 | 2007-12-18 | 适合于光纤陀螺的载体姿态测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100588906C (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102768043A (zh) * | 2012-06-14 | 2012-11-07 | 辽宁工程技术大学 | 一种无外观测量的调制型捷联***组合姿态确定方法 |
CN102997932A (zh) * | 2011-09-15 | 2013-03-27 | 北京自动化控制设备研究所 | 一种消除高精度惯导***标定中转台抖动影响的方法 |
CN103148868A (zh) * | 2013-02-22 | 2013-06-12 | 哈尔滨工程大学 | 匀速直航下基于多普勒计程仪地理系测速误差估计的组合对准方法 |
CN103792843A (zh) * | 2014-01-24 | 2014-05-14 | 北京航天控制仪器研究所 | 一种惯性平台快速转动控制方法 |
CN103900608A (zh) * | 2014-03-21 | 2014-07-02 | 哈尔滨工程大学 | 一种基于四元数ckf的低精度惯导初始对准方法 |
CN103954301A (zh) * | 2014-05-12 | 2014-07-30 | 北京航天控制仪器研究所 | 惯性测量***基于火箭橇轨道坐标系的定位方法 |
CN104634344A (zh) * | 2013-11-11 | 2015-05-20 | 精工爱普生株式会社 | 姿势计算装置及姿势计算方法 |
CN103954301B (zh) * | 2014-05-12 | 2016-11-30 | 北京航天控制仪器研究所 | 惯性测量***基于火箭橇轨道坐标系的定位方法 |
CN106997053A (zh) * | 2017-02-24 | 2017-08-01 | 北京航天自动控制研究所 | 一种共支架基座三捷联惯组导航一致性试验方法 |
CN109297487A (zh) * | 2017-07-25 | 2019-02-01 | 北京信息科技大学 | 一种角速度输入条件下的姿态解耦方法 |
CN110345942A (zh) * | 2019-07-17 | 2019-10-18 | 哈尔滨工程大学 | 一种基于角速率输入的插值三子样圆锥误差补偿算法 |
CN111351483A (zh) * | 2020-03-31 | 2020-06-30 | 北京控制工程研究所 | 一种递归多子样大动态惯性导航方法 |
CN112857366A (zh) * | 2021-01-14 | 2021-05-28 | 东南大学 | 一种基于压缩结构的光纤捷联惯导***姿态解算方法 |
CN113720206A (zh) * | 2021-09-02 | 2021-11-30 | 重庆零壹空间科技集团有限公司 | 火箭地面瞄准方法、***、计算机设备和存储介质 |
-
2007
- 2007-12-18 CN CN200710144846A patent/CN100588906C/zh not_active Expired - Fee Related
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102997932B (zh) * | 2011-09-15 | 2015-08-05 | 北京自动化控制设备研究所 | 一种消除高精度惯导***标定中转台抖动影响的方法 |
CN102997932A (zh) * | 2011-09-15 | 2013-03-27 | 北京自动化控制设备研究所 | 一种消除高精度惯导***标定中转台抖动影响的方法 |
CN102768043A (zh) * | 2012-06-14 | 2012-11-07 | 辽宁工程技术大学 | 一种无外观测量的调制型捷联***组合姿态确定方法 |
CN102768043B (zh) * | 2012-06-14 | 2015-03-04 | 辽宁工程技术大学 | 一种无外观测量的调制型捷联***组合姿态确定方法 |
CN103148868A (zh) * | 2013-02-22 | 2013-06-12 | 哈尔滨工程大学 | 匀速直航下基于多普勒计程仪地理系测速误差估计的组合对准方法 |
CN103148868B (zh) * | 2013-02-22 | 2015-12-09 | 哈尔滨工程大学 | 匀速直航下基于多普勒计程仪地理系测速误差估计的组合对准方法 |
CN104634344A (zh) * | 2013-11-11 | 2015-05-20 | 精工爱普生株式会社 | 姿势计算装置及姿势计算方法 |
CN103792843A (zh) * | 2014-01-24 | 2014-05-14 | 北京航天控制仪器研究所 | 一种惯性平台快速转动控制方法 |
CN103900608A (zh) * | 2014-03-21 | 2014-07-02 | 哈尔滨工程大学 | 一种基于四元数ckf的低精度惯导初始对准方法 |
CN103900608B (zh) * | 2014-03-21 | 2016-08-17 | 哈尔滨工程大学 | 一种基于四元数ckf的低精度惯导初始对准方法 |
CN103954301A (zh) * | 2014-05-12 | 2014-07-30 | 北京航天控制仪器研究所 | 惯性测量***基于火箭橇轨道坐标系的定位方法 |
CN103954301B (zh) * | 2014-05-12 | 2016-11-30 | 北京航天控制仪器研究所 | 惯性测量***基于火箭橇轨道坐标系的定位方法 |
CN106997053A (zh) * | 2017-02-24 | 2017-08-01 | 北京航天自动控制研究所 | 一种共支架基座三捷联惯组导航一致性试验方法 |
CN106997053B (zh) * | 2017-02-24 | 2019-11-12 | 北京航天自动控制研究所 | 一种共支架基座三捷联惯组导航一致性试验方法 |
CN109297487A (zh) * | 2017-07-25 | 2019-02-01 | 北京信息科技大学 | 一种角速度输入条件下的姿态解耦方法 |
CN110345942A (zh) * | 2019-07-17 | 2019-10-18 | 哈尔滨工程大学 | 一种基于角速率输入的插值三子样圆锥误差补偿算法 |
CN111351483A (zh) * | 2020-03-31 | 2020-06-30 | 北京控制工程研究所 | 一种递归多子样大动态惯性导航方法 |
CN112857366A (zh) * | 2021-01-14 | 2021-05-28 | 东南大学 | 一种基于压缩结构的光纤捷联惯导***姿态解算方法 |
CN113720206A (zh) * | 2021-09-02 | 2021-11-30 | 重庆零壹空间科技集团有限公司 | 火箭地面瞄准方法、***、计算机设备和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN100588906C (zh) | 2010-02-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100588906C (zh) | 适合于光纤陀螺的载体姿态测量方法 | |
CN100547352C (zh) | 适合于光纤陀螺捷联惯性导航***的地速检测方法 | |
US10670424B2 (en) | Method for initial alignment of an inertial navigation apparatus | |
CN101173858B (zh) | 一种月面巡视探测器的三维定姿与局部定位方法 | |
CN105300379B (zh) | 一种基于加速度的卡尔曼滤波姿态估计方法及*** | |
CN107560613B (zh) | 基于九轴惯性传感器的机器人室内轨迹跟踪***及方法 | |
CN103776446B (zh) | 一种基于双mems-imu的行人自主导航解算算法 | |
CN101629826A (zh) | 基于单轴旋转的光纤陀螺捷联惯性导航***粗对准方法 | |
Anjum et al. | Sensor data fusion using unscented kalman filter for accurate localization of mobile robots | |
CN107063254B (zh) | 一种陀螺地磁组合的姿态解算方法 | |
CN103697878B (zh) | 一种单陀螺单加速度计旋转调制寻北方法 | |
CN109959374B (zh) | 一种行人惯性导航全时全程逆向平滑滤波方法 | |
CN102003967A (zh) | 基于罗经原理的旋转式船用捷联惯导方位对准方法 | |
CN102288177B (zh) | 一种基于角速率输出的捷联***速度解算方法 | |
CN109931952A (zh) | 未知纬度条件下捷联惯导直接解析式粗对准方法 | |
CN103955005B (zh) | 一种火箭橇轨道重力实时测量方法 | |
Wongwirat et al. | A position tracking experiment of mobile robot with inertial measurement unit (imu) | |
CN103292810B (zh) | 一种旋转式惯导***信号同步补偿方法 | |
CN109724579B (zh) | 一种陀螺罗经标定方法、装置、计算设备及存储介质 | |
CN109827572B (zh) | 一种检测车位置预测的方法及装置 | |
CN104101345A (zh) | 基于互补重构技术的多传感器姿态融合方法 | |
Zhe et al. | Adaptive complementary filtering algorithm for imu based on mems | |
JPH04231813A (ja) | 角度および角度特性曲線の測定方法 | |
TW498169B (en) | Interruption-free hand-held positioning method and system thereof | |
Chen et al. | Foot-mounted Dual-sensor Single-board Pedestrian Inertial Navigation System Based on Position and Velocity Constraints. |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100210 Termination date: 20171218 |