CN101178313A - 适合于光纤陀螺捷联惯性导航***的地速检测方法 - Google Patents
适合于光纤陀螺捷联惯性导航***的地速检测方法 Download PDFInfo
- Publication number
- CN101178313A CN101178313A CNA2007101447309A CN200710144730A CN101178313A CN 101178313 A CN101178313 A CN 101178313A CN A2007101447309 A CNA2007101447309 A CN A2007101447309A CN 200710144730 A CN200710144730 A CN 200710144730A CN 101178313 A CN101178313 A CN 101178313A
- Authority
- CN
- China
- Prior art keywords
- speed
- ground
- carrier
- increment
- update cycle
- 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
适合于光纤陀螺捷联惯性导航***的地速检测方法,涉及一种划船补偿方法。它解决了在高动态环境或是高频率的振动环境中,划船效应对地速检测精度产生影响的问题。本发明通过在地速更新周期内,分别由比力信号和角速度信号获得载体相对于惯性坐标系的速度测量增量和角增量测量值,进而获得速度旋转项;同时通过比力信号和角速度信号的叉乘项拟和得到划船补偿项;由速度旋转项和划船补偿项计算更新周期内载体相对于惯性坐标系的速度增量,然后将获得的速度增量投影到导航坐标系n系中,在导航坐标系n系中计算重力加速度和哥氏加速度引起的速度增量,最终通过惯性导航***的基本方程计算得到地速值。本发明适合于载体处于高频振动或是高机动的情形。
Description
技术领域
本发明涉及一种光纤陀螺捷联惯性导航***的地速检测方法。
背景技术
在捷联惯性导航***中,包括两类关键的计算。其一就是更新载体姿态的计算,另一个则是更新载体地速(导航坐标系n系相对地球坐标系的速度值)的计算。在姿态更新计算中,如果存在载体角速度向量旋转的情形,就会引入圆锥误差。在更新载体地速的计算中,情况就更为复杂,因为,要获取载体的地速信息,就需要在导航参考坐标系中对比力信号进行积分,这个积分过程通常分两步进行:首先是将固联在载体系上的加速度计敏感的比力信号用该时刻的姿态方向余弦矩阵或四元数转换到导航参考坐标系中,然后在导航参考坐标系中对比力信号进行积分。由于捷联惯性导航***中的加速度计是直接固联在载体上的,经历了载体的姿态角运动,所以,如果存在载体角速度向量旋转或比力向量旋转的情形,或是比力和角速度的大小不恒定的情形,都会引入不可交换性误差的常值漂移积分项,这就是速度划船效应。如果在地速检测过程中不仅比力信号的大小和方向发生了变化,而且姿态角速度的大小和方向也发生了变化,则此过程将引入划船效应误差,这使高质量地完成比力积分的两个步骤变得十分困难,要高质量地完成捷联惯性导航***的比力积分,就需要高精度的加速度计和高性能的地速检测方法,而高精度的加速度计需要高水平的制造工艺使其成本增加,要想在提高捷联惯性导航***的地速检测精度的前提下,还不增加***成本,就需要尽量消除划船误差对地速检测的影响,即加入有效的划船补偿方法。
对于传统的划船补偿方法,思路是利用角增量、速度增量拟合出地速检测过程中的划船误差,这是建立在传统的惯性传感器输出为角增量、速度增量的基础上,传统方法可以直接使用。但是对于目前广泛使用的光纤陀螺捷联***,光纤陀螺输出为角速度,石英加速度计输出为比力,使在传统方法中用到的角增量和速度增量不能直接获取的,传统方法的实用性出现问题。由角速度获得角增量,比力获得速度增量可以采取数值分析中的辛普森积分,其实质是对角速度和比力进行分段线性化,这种处理方式只有在采样周期足够短,并且角速度、比力矢量变化较慢时才可以使用。在高动态环境或是高频振动环境中,积分过程使传统划船补偿算法对于划船补偿项的计算出现偏差,即传统划船补偿方法的系数在角速度、比力输入时不再是最优。而目前广泛使用的二阶龙格库达(Lgkd)法可以直接利用比力信息进行地速检测,但对于划船误差没有补偿效果。在高动态环境中或是高频振动环境中,划船效应将严重影响地速检测的精度。
发明内容
为了解决在高动态环境或是高频率的振动环境中,划船效应对***地速检测精度产生影响的问题,本发明提供了一种适合于光纤陀螺捷联惯性导航***的地速检测方法。
适合于光纤陀螺捷联惯性导航***的地速检测方法,具体步骤如下:
步骤一、通过外部设备确定载体的位置参数和初始的地速值;
步骤二、捷联惯性导航***进行初始对准,确定载体相对于导航坐标系n系的初始姿态;
步骤三、确定捷联惯性导航***的地速更新周期H=tm-tm-1,所述地速更新周期H等于划船补偿周期;设置H为N倍的惯性传感器采样周期,所述N为大于0的整数;
步骤四、分别采集石英加速度计输出的比力信号f和光纤陀螺输出的角速度信号ω,计算获得载体相对于惯性坐标系的速度增量在载体坐标系b系上的投影ΔvSm b;
步骤五、由姿态更新算法可以得到地速更新周期H起始时刻的捷联姿态矩阵Cb n,通过Cb n将步骤四获得的载体相对于惯性坐标系的速度增量在载体坐标系b系上的投影ΔvSm b转换到导航坐标系n系上的投影得到ΔvSm n;
步骤六、在导航坐标系n系上,计算地速更新周期H内重力加速度和哥氏加速度引起的速度增量ΔvG/Corm n;
步骤七、根据惯性导航***的基本方程,由步骤五获得的载体相对于惯性坐标系的速度增量在导航坐标系n系上的投影ΔvSm n和步骤六获得的重力加速度与哥氏加速度引起的速度增量ΔvG/Corm n,计算得到地速更新周期H内导航坐标系n系相对地球坐标系的速度增量;
步骤八、将步骤七获得的地速更新周期H内导航坐标系n系相对地球坐标系的速度增量与地速更新周期H初始时刻的地速值相加,更新得到H终止时刻的地速值vm n。
所述步骤四中,计算获得载体相对惯性坐标系的速度的过程为:
步骤四一:在地速更新周期H内,对采集到的石英加速度计输出的比力信号f积分得到地速更新周期H内载体相对于惯性坐标系的速度测量增量vm;
步骤四二:在地速更新周期H内,对采集到的光纤陀螺输出的角速度信号ω积分得到地速更新周期H内载体相对于惯性坐标系的角增量测量值θm;
步骤四三:通过角增量测量θm和速度测量增量vm获得地速更新周期H内的速度旋转项Δvrotm;
步骤四四:用地速更新周期H内的N+1个光纤陀螺和石英加速度计采样值的叉乘项拟和划船补偿项Δvscul;
步骤四五:将步骤四一获得的速度测量增量vm、步骤四三获得的速度旋转项Δvrotm、步骤四四获得的划船补偿项Δvscul相加,获得地速更新周期H内载体相对于惯性坐标系的速度在载体坐标系b系上的投影ΔvSm b。
本发明的适合于光纤陀螺捷联惯性导航***的地速检测方法具有以下优点:
一、本发明的方法是以角速度、比力作为输入量,既是直接利用光纤陀螺和石英加速度计的输出,用角速度、比力表示地速的微分方程,以角速率和比力的叉乘项拟合划船误差项,进而消除划船效应对地速更新的影响。
二、本发明的方法避免了使用传统划船补偿方法在角速度、比力作为输入时,积分误差对于划船补偿项的影响,又对高频振动和高动态环境中的划船误差予以补偿。
三、本发明的方法在不增加***成本的情况下,有效地提高了捷联惯性导航***的地速输出精度。
四、本发明的方法采用每N个采样周期计算一次地速增量,相对于传统的采用二阶龙格库达法的地速检测方法,在提高精度的同时,其地速的更新频率降低为原来的1/N,即计算频率降低为原来的1/N,从而减少了导航计算机的计算量,节约了导航计算机的有限资源。
本发明的方法尤其适合于载体处于高频振动或是高机动的情形。
附图说明
图1为本发明的流程图。
图2为步骤四中通过比力信号和角速度信号获得载体相对于惯性坐标系的速度增量的流程图。
图3为典型划船环境中,分别采用本发明的方法、传统的划船补偿方法以及二阶龙格库达法获得的速度误差仿真曲线。典型划船环境定义为载体的两个正交轴存在同频且同相位的角振动和线振动。光纤陀螺和石英加速度计的采样频率为100Hz。
具体实施方式
本实施方式所述的适合于光纤陀螺捷联惯性导航***的地速检测方法,具体步骤如下:
步骤一、通过外部设备确定载体的位置参数和初始的地速值;
步骤二、捷联惯性导航***进行初始对准,确定载体相对于导航坐标系n系的初始姿态;
步骤三、确定地速更新周期H=tm-tm-1,所述地速更新周期H等于划船补偿周期;设置H为N倍的惯性传感器采样周期,所述N为大于0的整数;
步骤四、分别采集石英加速度计输出的比力信号f和光纤陀螺输出的角速度信号ω,计算获得载体相对于惯性坐标系的速度增量在载体坐标系b系上的投影ΔvSm b;
步骤五、由姿态更新算法可以得到地速更新周期H起始时刻的捷联姿态矩阵Cb n,通过Cb n将步骤四中获得的载体相对于惯性坐标系的速度增量在载体坐标系b系上的投影ΔvSm n转换到导航坐标系n系上的投影得到ΔvSm n;
步骤六、在导航坐标系n系上,计算地速更新周期H内重力加速度和哥氏加速度引起的速度增量ΔvG/Corm n;
步骤七、根据惯性导航***的基本方程,由步骤五获得的载体相对于惯性坐标系的速度增量在导航坐标系n系上的投影ΔvSm n和步骤六获得的重力加速度和哥氏加速度引起的速度增量ΔvG/Corm n,计算得到地速更新周期H内导航坐标系n系相对地球坐标系的速度增量;
步骤八、将步骤七获得的地速更新周期H内导航坐标系n系相对地球坐标系的速度增量与地速更新周期H初始时刻的地速值相加,更新得到H终止时刻的地速值vm n。
所述步骤四中的计算获得载体相对于惯性坐标系的速度的过程参见图2所示,具体步骤为:
步骤四二:在地速更新周期H内,对采集到的石英加速度计输出的比力信号f积分得到地速更新周期H内载体相对于惯性坐标系的速度测量增量vm;
步骤四三:在地速更新周期H内,对采集到的光纤陀螺输出的角速度信号ω积分得到地速更新周期H内载体相对于惯性坐标系的角增量测量值θm;
步骤四四:通过角增量测量θm和速度测量增量vm获得地速更新周期H内的速度旋转项Δvrotm;
步骤四五:地速更新周期H内,用N+1个光纤陀螺和石英加速度计采样值的叉乘项拟和划船补偿项;
步骤四六:将步骤四一获得的速度测量增量vm、步骤四三获得的速度旋转项Δvrotm、步骤四四获得的划船补偿项Δvscul相加,获得计算地速更新周期H内载体相对于惯性坐标系的速度及其在载体坐标系b系上的投影ΔvSm b。
在步骤一中所述的载体的位置参数由GPS装置或者是外部的高精度组合导航设备提供,所述载体的初始速度由DVL多普勒计程仪或者是外部的高精度组合导航设备提供。
在步骤四一中,通过对石英加速度计输出比力信号f进行积分得到速度测量增量vm:
在步骤四二中,通过对光纤陀螺输出的角速度信号ω进行积分得到测量角增量θm:
在步骤四三中所述的速度旋转项Δvrotm是速度旋转效应补偿量,由速度测量增量vm和测量角增量θm获得:
在步骤四四中所述的划船补偿项Δvsculm是动态积分项,是非惯性可测项,由采集到的石英加速度计输出比力信号f和光纤陀螺输出的角速度信号ω获得:
在步骤四五中所述的载体相对于惯性坐标系的速度增量在载体坐标系b系上的投影Δvsm b是由比力加速度引起,由步骤四一、步骤四三和步骤四四的结果获得:
在划船补偿方法的标准测试输入环境,即典型划船环境中,数值积分误差对载体相对于惯性***坐标的速度测量增量vm和速度旋转项Δvrotm的计算不产生积累误差,只会造成划船补偿项Δvsculm的漂移误差,为了避免漂移误差的产生,在划船补偿项Δvsculm的计算公式中使用角速度ω和比力f的叉乘项进行拟和计算。
例如:当H=2h的时候,采样3次光纤陀螺和石英加速度计输出值,划船补偿项Δvsculm的离散计算为
当H=3h的时候,采样4次光纤陀螺和石英加速度计输出值,Δvsculm的离散计算为
在地速更新周期H内,划船补偿方法所使用的惯性传感器采样数越多,划船补偿方法的输出精度就越高。
在步骤五中,根据所述捷联矩阵Cb n可由姿态更新获得,利用Cb n将步骤四中获得的载体相对于惯性坐标系的速度增量转换到导航坐标系n系中获得ΔvSm n为:
在步骤六中重力加速度和哥氏加速度引起的速度增量ΔvG/Corm n为
其中gP n为重力,ωen n为位置角速度,ωie n为地球转速。在地速更新周期H内纬度只有微小的变化,位置角速度和地球转速的变化量都是小量,所以把它们在地速更新周期H内的平均值代入进行离散数字计算,
其中
在步骤八中,选取地理坐标系作为导航坐标系n系,由惯性导航***的基本方程,地速变化率为:
导航计算机中,对地速进行离散递推计算
得到地速更新周期终止时刻的地速值,完成地速的更新。
其中地速的初始值vm-1 n由上一个地速更新周期计算给出;载体相对于惯性坐标系的速度在导航坐标系n系的投影ΔvSm n由步骤五获得。
Claims (7)
1.适合于光纤陀螺捷联惯性导航***的地速检测方法,具体步骤如下:
步骤一、通过外部设备确定载体的位置参数和初始的地速值;
步骤二、捷联惯性导航***进行初始对准,确定载体相对于导航坐标系n系的初始姿态;
步骤三、确定地速更新周期H=tm-tm-1,所述地速更新周期H等于划船补偿周期;设置H为N倍的惯性传感器采样周期,所述N为大于0的整数;
步骤四、分别采集石英加速度计输出的比力信号f和光纤陀螺输出的角速度信号ω,计算获得载体相对于惯性坐标系的速度增量在载体坐标系b系上的投影ΔvSm b;
步骤五、由姿态更新算法可以得到地速更新周期H起始时刻的捷联姿态矩阵Cb n,通过Cb n将步骤四获得的载体相对于惯性坐标系的速度增量在载体坐标系b系上的投影ΔvSm n转换到导航坐标系n系上的投影得到ΔvSm n;
步骤六、在导航坐标系n系上,计算地速更新周期H内重力加速度和哥氏加速度引起的速度增量ΔvG/Corm n;
步骤七、根据惯性导航***的基本方程,由步骤五获得的载体相对于惯性坐标系的速度增量在导航坐标系n系上的投影ΔvSm n和步骤六获得的重力加速度和哥氏加速度引起的速度增量ΔvG/Corm n,计算地速更新周期H内导航坐标系n系相对地球坐标系的速度增量;
步骤八、将步骤七获得的地速更新周期H内导航坐标系n系相对地球坐标系的速度增量与地速更新周期H初始时刻的地速值相加,更新得到地速更新周期H终止时刻的地速值vm n,
其特征在于,所述步骤四中的计算获得载体相对于惯性坐标系的速度的过程为:
步骤四一:在地速更新周期H内,对采集到的石英加速度计输出的比力信号f积分得到地速更新周期H内载体相对于惯性坐标系的速度测量增量vm;
步骤四二:在地速更新周期H内,对采集到的光纤陀螺输出的角速度信号ω积分得到地速更新周期H内载体相对于惯性坐标系的角增量测量值θm;;
步骤四三:通过角增量测量θm和速度测量增量vm获得地速更新周期H内的速度旋转项Δvrotm;
步骤四四:地速更新周期H内,用N+1个光纤陀螺和石英加速度计采样值的叉乘项拟和划船补偿项;
步骤四五:将步骤四一获得的速度测量增量vm、步骤四三获得的速度旋转项Δvrotm、步骤四四获得的划船补偿项Δvscul相加,获得计算地速更新周期H内载体相对于惯性坐标系的速度及其在载体坐标系b系上的投影ΔvSm b。
2.根据权利要求1所述的适合于光纤陀螺捷联惯性导航***的地速检测方法,其特征在于,在步骤一中所述的载体的位置参数由GPS装置或者是外部的高精度组合导航设备提供,所述载体的初始速度由DVL多普勒计程仪或者是外部的高精度组合导航设备提供。
3.根据权利要求l所述的适合于光纤陀螺捷联惯性导航***的地速检测方法,其特征在于,在步骤四三中所述的速度旋转项Δvrotm是速度旋转效应补偿量,由速度测量增量vm和测量角增量θm获得:
5.根据权利要求1所述的适合于光纤陀螺捷联惯性导航***的地速检测方法,其特征在于,在步骤四五中所述的载体相对于惯性坐标系的速度增量在载体坐标系b系上的投影Δvsm b是由比力加速度引起,由步骤四一、步骤四三和步骤四四的结果获得:
6.根据权利要求1所述的适合于光纤陀螺捷联惯性导航***的地速检测方法,其特征在于,在步骤六中重力加速度和哥氏加速度引起的速度增量ΔvG/Corm n为 其中gP n为重力,ωen n为位置角速度,ωie n为地球转速,在地速更新周期H内由离散数字计算获得
7.根据权利要求1所述的适合于光纤陀螺捷联惯性导航***的地速检测方法,其特征在于,在步骤八中,选取地理坐标系作为导航坐标系n系,由惯性导航***的基本方程,地速的变化率为: 导航计算机中,对地速进行离散递推,获得地速值vm n为
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007101447309A CN100547352C (zh) | 2007-12-03 | 2007-12-03 | 适合于光纤陀螺捷联惯性导航***的地速检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007101447309A CN100547352C (zh) | 2007-12-03 | 2007-12-03 | 适合于光纤陀螺捷联惯性导航***的地速检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101178313A true CN101178313A (zh) | 2008-05-14 |
CN100547352C CN100547352C (zh) | 2009-10-07 |
Family
ID=39404656
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2007101447309A Expired - Fee Related CN100547352C (zh) | 2007-12-03 | 2007-12-03 | 适合于光纤陀螺捷联惯性导航***的地速检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100547352C (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103175528A (zh) * | 2013-02-19 | 2013-06-26 | 哈尔滨工程大学 | 基于捷联惯导***的捷联罗经姿态测量方法 |
CN103292811A (zh) * | 2013-06-20 | 2013-09-11 | 哈尔滨工程大学 | 一种虚拟旋转地理坐标系下解算速度的捷联惯性导航方法 |
CN103323022A (zh) * | 2013-04-26 | 2013-09-25 | 哈尔滨工程大学 | 一种角增量速度增量捷联惯性导航***粗对准方法 |
CN104567873A (zh) * | 2015-01-15 | 2015-04-29 | 南京航空航天大学 | 针对高动态载体应用的高精度捷联惯导姿态角确定方法 |
CN104634344A (zh) * | 2013-11-11 | 2015-05-20 | 精工爱普生株式会社 | 姿势计算装置及姿势计算方法 |
CN104677356A (zh) * | 2015-03-16 | 2015-06-03 | 黄磊 | 一种基于角增量和比力输出的划桨速度计算方法 |
CN105511484A (zh) * | 2015-11-27 | 2016-04-20 | 深圳一电航空技术有限公司 | 控制无人机稳定飞行的方法及装置 |
CN110319833A (zh) * | 2019-07-09 | 2019-10-11 | 哈尔滨工程大学 | 一种无误差的光纤陀螺捷联惯导***速度更新方法 |
CN110763231A (zh) * | 2019-10-15 | 2020-02-07 | 哈尔滨工程大学 | 一种适用于光纤陀螺滤波信号的无误差的姿态更新方法 |
CN112363518A (zh) * | 2020-09-24 | 2021-02-12 | 北京航天自动控制研究所 | 一种减少火箭飞行软件前台运算时间的姿态角计算方法 |
CN113490855A (zh) * | 2019-03-29 | 2021-10-08 | 拓普康定位***公司 | 用于对石英振荡器的g灵敏度进行测量的方法和*** |
CN114008410A (zh) * | 2019-06-14 | 2022-02-01 | 赛峰电子与防务公司 | 用于监测惯性测量单元的性能的方法 |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102168978B (zh) * | 2010-12-17 | 2012-10-31 | 北京航空航天大学 | 一种船用惯性导航***摇摆基座开环对准方法 |
CN102155941B (zh) * | 2011-03-23 | 2013-05-01 | 重庆华渝电气仪表总厂 | 捷联垂直参考基准***上挠性陀螺仪动态启停控制*** |
-
2007
- 2007-12-03 CN CNB2007101447309A patent/CN100547352C/zh not_active Expired - Fee Related
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103175528A (zh) * | 2013-02-19 | 2013-06-26 | 哈尔滨工程大学 | 基于捷联惯导***的捷联罗经姿态测量方法 |
CN103175528B (zh) * | 2013-02-19 | 2015-04-08 | 哈尔滨工程大学 | 基于捷联惯导***的捷联罗经姿态测量方法 |
CN103323022A (zh) * | 2013-04-26 | 2013-09-25 | 哈尔滨工程大学 | 一种角增量速度增量捷联惯性导航***粗对准方法 |
CN103292811A (zh) * | 2013-06-20 | 2013-09-11 | 哈尔滨工程大学 | 一种虚拟旋转地理坐标系下解算速度的捷联惯性导航方法 |
CN104634344A (zh) * | 2013-11-11 | 2015-05-20 | 精工爱普生株式会社 | 姿势计算装置及姿势计算方法 |
CN104567873A (zh) * | 2015-01-15 | 2015-04-29 | 南京航空航天大学 | 针对高动态载体应用的高精度捷联惯导姿态角确定方法 |
CN104567873B (zh) * | 2015-01-15 | 2017-07-07 | 南京航空航天大学 | 针对高动态载体应用的高精度捷联惯导姿态角确定方法 |
CN104677356B (zh) * | 2015-03-16 | 2017-06-16 | 黄磊 | 一种基于角增量和比力输出的划桨速度计算方法 |
CN104677356A (zh) * | 2015-03-16 | 2015-06-03 | 黄磊 | 一种基于角增量和比力输出的划桨速度计算方法 |
CN105511484A (zh) * | 2015-11-27 | 2016-04-20 | 深圳一电航空技术有限公司 | 控制无人机稳定飞行的方法及装置 |
CN113490855A (zh) * | 2019-03-29 | 2021-10-08 | 拓普康定位***公司 | 用于对石英振荡器的g灵敏度进行测量的方法和*** |
CN113490855B (zh) * | 2019-03-29 | 2022-10-25 | 拓普康定位***公司 | 用于对石英振荡器的g灵敏度进行测量的方法和*** |
CN114008410A (zh) * | 2019-06-14 | 2022-02-01 | 赛峰电子与防务公司 | 用于监测惯性测量单元的性能的方法 |
CN114008410B (zh) * | 2019-06-14 | 2024-07-05 | 赛峰电子与防务公司 | 用于监测惯性测量单元的性能的方法 |
CN110319833A (zh) * | 2019-07-09 | 2019-10-11 | 哈尔滨工程大学 | 一种无误差的光纤陀螺捷联惯导***速度更新方法 |
CN110319833B (zh) * | 2019-07-09 | 2022-07-15 | 哈尔滨工程大学 | 一种无误差的光纤陀螺捷联惯导***速度更新方法 |
CN110763231A (zh) * | 2019-10-15 | 2020-02-07 | 哈尔滨工程大学 | 一种适用于光纤陀螺滤波信号的无误差的姿态更新方法 |
CN112363518A (zh) * | 2020-09-24 | 2021-02-12 | 北京航天自动控制研究所 | 一种减少火箭飞行软件前台运算时间的姿态角计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN100547352C (zh) | 2009-10-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100547352C (zh) | 适合于光纤陀螺捷联惯性导航***的地速检测方法 | |
CN105300379B (zh) | 一种基于加速度的卡尔曼滤波姿态估计方法及*** | |
CN102289306B (zh) | 姿态感知设备及其定位、鼠标指针的控制方法和装置 | |
CN103808331B (zh) | 一种mems三轴陀螺仪误差标定方法 | |
CN101290326B (zh) | 石英挠性加速度计测量组件的参数辨识标定方法 | |
CN101514900B (zh) | 一种单轴旋转的捷联惯导***初始对准方法 | |
CN101706284B (zh) | 提高船用光纤陀螺捷联惯导***定位精度的方法 | |
CN103090870B (zh) | 一种基于mems传感器的航天器姿态测量方法 | |
CN103900565B (zh) | 一种基于差分gps的惯导***姿态获取方法 | |
CN101949710B (zh) | Gnss辅助mems惯性传感器零偏的快速在线动态标定方法 | |
CN100588906C (zh) | 适合于光纤陀螺的载体姿态测量方法 | |
KR20170104623A (ko) | 관성 항법 장치의 초기 정렬 방법 | |
CN101701825A (zh) | 高精度激光陀螺单轴旋转惯性导航*** | |
CN103323625B (zh) | 一种mems-imu中加速度计动态环境下的误差标定补偿方法 | |
CN101246023A (zh) | 微机械陀螺惯性测量组件的闭环标定方法 | |
CN103278163A (zh) | 一种基于非线性模型的sins/dvl组合导航方法 | |
CN101571394A (zh) | 基于旋转机构的光纤捷联惯性导航***初始姿态确定方法 | |
CN101162147A (zh) | 大失准角下船用光纤陀螺捷联航姿***系泊精对准方法 | |
CN103674064B (zh) | 捷联惯性导航***的初始标定方法 | |
CN101706287A (zh) | 一种基于数字高通滤波的旋转捷联***现场标定方法 | |
CN101183004A (zh) | 一种在线实时消除光纤陀螺捷联惯导***振荡误差的方法 | |
CN101738203A (zh) | 挠性陀螺仪静态漂移零次和一次加速度相关项误差模型最优位置标定方法 | |
CN109931952A (zh) | 未知纬度条件下捷联惯导直接解析式粗对准方法 | |
Wang et al. | An improve hybrid calibration scheme for strapdown inertial navigation system | |
CN109084755B (zh) | 一种基于重力视速度与参数辨识的加速度计零偏估计方法 |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20091007 Termination date: 20111203 |