CN103389090A - 一种惯性导航系统极区导航模式初始速度的测量方法 - Google Patents

一种惯性导航系统极区导航模式初始速度的测量方法 Download PDF

Info

Publication number
CN103389090A
CN103389090A CN2013103223223A CN201310322322A CN103389090A CN 103389090 A CN103389090 A CN 103389090A CN 2013103223223 A CN2013103223223 A CN 2013103223223A CN 201310322322 A CN201310322322 A CN 201310322322A CN 103389090 A CN103389090 A CN 103389090A
Authority
CN
China
Prior art keywords
inertial navigation
navigation system
coordinate system
theta
initial velocity
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.)
Pending
Application number
CN2013103223223A
Other languages
English (en)
Inventor
奔粤阳
鲍桂清
尹冬寒
赵维珩
阮双双
杨祥龙
叶攀
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN2013103223223A priority Critical patent/CN103389090A/zh
Publication of CN103389090A publication Critical patent/CN103389090A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Navigation (AREA)

Abstract

本发明公开了一种惯性导航系统极区导航模式初始速度的测量方法,包括以下步骤:采集惯性导航系统正常模式下输出信息;测量惯性导航系统与东经90°之间的经度差;计算辅助系数W;测量地理坐标系与横地理坐标系之间的转角;计算地理坐标系与横地理坐标系之间的方向余弦矩阵;测量惯性惯性导航系统极区模式的初始速度。本发明基于地球椭球模型求取正常模式与极区模式导航坐标系之间的转换矩阵,避免了简单以常规球体模拟地球时求取转换矩阵的误差,从而减小惯性导航系统模式转换的误差。本发明只需惯性导航系统正常模式下的输出信息就可以测量出极区模式所需的初始速度信息,不需要其他外部设备和信息,测量方法简单方便,有利于实际应用。

Description

一种惯性导航系统极区导航模式初始速度的测量方法
技术领域
本发明属于惯性导航系统极区导航技术领域,尤其涉及一种惯性导航系统极区导航模式初始速度的测量方法。
背景技术
惯性导航系统对于保障船舶、飞机等在极区的航行与作业至关重要。然而目前我国常用的固定指北方位惯性导航系统,其机械编排无法在极区使用。为克服惯性导航系统上述机械编排无法在极区进行导航解算这一缺点,导航领域的研究人员设计出自由方位和游动方位机械编排。自由方位系统和游动方位系统在平台的方位指令上采取了措施,解决了极区航行时平台方位角施矩速率信号过大的问题。它们可解决惯性导航系统极区工作的部分问题,但由于极区经线圈的快速汇聚会导致系统位置矩阵的误差,进一步使得系统导航精度随纬度的升高而下降。同时,在地理极点附近,由于真北方向失去定义导致自由方位角和游动方位角失去定义,航向无法给出,即采用这两种机械编排惯性导航系统在地理极点附近无法导航。
解决这个问题可以采用横坐标系导航参考框架。运载体航行至极区时,转入极区导航模式,即采用横坐标系参考框架。横坐标参考框架是指惯性导航系统以横地理坐标系为导航坐标系,在横地理坐标系内解算导航参数,包括速度,位置等信息。由于依然采用的是通过加速度积分来求解速度的方式,因此在进入极区模式时需要横坐标系参考框架下速度的初始值。但是惯性导航系统存在的正常模式转入极区模式时速度信息的初始化问题。
发明内容
本发明实施例的目的在于提供一种惯性导航系统极区导航模式初始速度的测量方法,旨在解决惯性导航系统由正常模式转入极区模式时速度信息的初始化的问题。
本发明实施例是这样实现的,一种惯性导航系统极区导航模式初始速度的测量方法,该惯性导航系统极区导航模式初始速度的测量方法包括以下步骤:
采集惯性导航系统正常模式下输出信息;
测量惯性导航系统与东经90°之间的经度差;
计算辅助系数W;
测量地理坐标系与横地理坐标系之间的转角;
计算地理坐标系与横地理坐标系之间的方向余弦矩阵;
测量惯性惯性导航系统极区模式的初始速度。
进一步,该惯性导航系统极区导航模式初始速度的测量方法具体步骤为:
步骤一:采集惯性导航系统正常模式输出的经度信息λ,纬度信息L和速度信息v;
步骤二:测量惯性导航系统所在位置与东经90°经线的经度差δλ;
步骤三:利用步骤一中的纬度信息L计算辅助系数W;
步骤四:利用步骤三中的W测量地理坐标系与横地理坐标系之间的转角θ;
步骤五:利用步骤六中得到的地理坐标系与横地理坐标系之间的转角θ构造地理坐标系与横地理坐标系之间的方向余弦矩阵
Figure BDA00003582704600031
步骤六:由步骤五中得到的
Figure BDA00003582704600032
测量惯性导航系统极区模式的初始速度。
进一步,在步骤二中,经度差δλ表示为:δλ=λ-90。
进一步,在步骤三中:辅助系数W表示为:
Figure BDA00003582704600033
进一步,在步骤四中:地理坐标系与横地理坐标系之间的转角θ表示为:
θ = arctan ( W sin δλ sin L cos δλ ( 1 - e 2 ) ) .
进一步,在步骤五中:地理坐标系与横地理坐标系之间的方向余弦矩阵
Figure BDA00003582704600035
表示为:
C t t ‾ = cos θ sin θ 0 - sin cos θ 0 0 0 1 .
进一步,在步骤六中:惯性导航系统极区模式的初始速度表示为:
v x ‾ v y ‾ v z ‾ = cos θ sin θ 0 - sin cos θ 0 0 0 1 v x v y v z
其中
Figure BDA00003582704600043
Figure BDA00003582704600044
是惯性导航系统进入极区模式所需横东向、横北向和横天向速度,vx、vy和vz是惯性导航系统正常模式输出的东向、北向和天向速度。
本发明提供的惯性导航系统极区导航模式初始速度的测量方法,基于地球椭球模型求取正常模式与极区模式导航坐标系之间的转换矩阵,避免了简单以常规球体模拟地球时求取转换矩阵的误差,从而减小惯性导航系统模式转换的误差。本发明只需惯性导航系统正常模式下的输出信息就可以测量出极区模式所需的初始速度信息,不需要其他外部设备和信息,测量方法简单方便,有利于实际应用。此外,本发明填补了惯性导航系统极区导航模式初始速度测量方法的空白,解决了惯性导航系统由正常模式转入极区模式时速度信息的初始化问题。
附图说明
图1是本发明实施例提供的惯性导航系统极区导航模式初始速度的测量方法的流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
图1示出了本发明提供的惯性导航系统极区导航模式初始速度的测量方法流程。为了便于说明,仅仅示出了与本发明相关的部分。
本发明的惯性导航系统极区导航模式初始速度的测量方法,该惯性导航系统极区导航模式初始速度的测量方法包括以下步骤:
采集惯性导航系统正常模式下输出信息;
测量惯性导航系统与东经90°之间的经度差;
计算辅助系数W;
测量地理坐标系与横地理坐标系之间的转角;
计算地理坐标系与横地理坐标系之间的方向余弦矩阵;
测量惯性惯性导航系统极区模式的初始速度。
作为本发明实施例的一优化方案,该惯性导航系统极区导航模式初始速度的测量方法具体步骤为:
步骤一:采集惯性导航系统正常模式输出的经度信息λ,纬度信息L和速度信息v;
步骤二:测量惯性导航系统所在位置与东经90°经线的经度差δλ;
步骤三:利用步骤一中的纬度信息L计算辅助系数W;
步骤四:利用步骤三中的W测量地理坐标系与横地理坐标系之间的转角θ;
步骤五:利用步骤六中得到的地理坐标系与横地理坐标系之间的转角θ构造地理坐标系与横地理坐标系之间的方向余弦矩阵
步骤六:由步骤五中得到的
Figure BDA00003582704600062
测量惯性导航系统极区模式的初始速度。
作为本发明实施例的一优化方案,在步骤二中,经度差δλ表示为:δλ=|λ-90|。
作为本发明实施例的一优化方案,在步骤三中:辅助系数W表示为: W = 1 - e 2 sin 2 L .
作为本发明实施例的一优化方案,在步骤四中:地理坐标系与横地理坐标系之间的转角θ表示为:
θ = arctan ( W sin δλ sin L cos δλ ( 1 - e 2 ) ) .
作为本发明实施例的一优化方案,在步骤五中:地理坐标系与横地理坐标系之间的方向余弦矩阵
Figure BDA00003582704600065
表示为:
C t t ‾ = cos θ sin θ 0 - sin cos θ 0 0 0 1 .
作为本发明实施例的一优化方案,在步骤六中:惯性导航系统极区模式的初始速度表示为:
v x ‾ v y ‾ v z ‾ = cos θ sin θ 0 - sin cos θ 0 0 0 1 v x v y v z
其中
Figure BDA00003582704600073
是惯性导航系统进入极区模式所需横东向、横北向和横天向速度,vx、vy和vz是惯性导航系统正常模式输出的东向、北向和天向速度。
下面结合附图及具体实施例对本发明的应用原理作进一步描述。
如图1所示,本发明实施例的惯性导航系统极区导航模式初始速度的测量方法包括以下步骤:
S101:采集惯性导航系统正常模式下输出信息;
S102:测量惯性导航系统与东经90°之间的经度差;
S103:计算辅助系数W;
S104:测量地理坐标系与横地理坐标系之间的转角;
S105:计算地理坐标系与横地理坐标系之间的方向余弦矩阵;
S106:测量惯性惯性导航系统极区模式的初始速度。
本发明的具体步骤如下:
步骤一:采集惯性导航系统正常模式输出的经度信息λ,纬度信息L和速度信息v;
步骤二:测量惯性导航系统所在位置与东经90°经线的经度差δλ:
δλ=|λ-90|
步骤三:利用步骤一中的纬度信息L计算辅助系数W:
W = 1 - e 2 sin 2 L .
步骤四:利用步骤三中的W测量地理坐标系与横地理坐标系之间的转角θ:
θ = arctan ( W sin δλ sin L cos δλ ( 1 - e 2 ) ) .
步骤五:利用步骤六中得到的地理坐标系与横地理坐标系之间的转角θ构造地理坐标系与横地理坐标系之间的方向余弦矩阵
Figure BDA00003582704600083
C t t ‾ = cos θ sin θ 0 - sin cos θ 0 0 0 1 .
步骤六:由步骤五中得到的
Figure BDA00003582704600085
测量惯性导航系统极区模式的初始速度:
v x ‾ v y ‾ v z ‾ = cos θ sin θ 0 - sin cos θ 0 0 0 1 v x v y v z
其中
Figure BDA00003582704600087
是惯性导航系统进入极区模式所需横东向、横北向和横天向速度,vx、vy和vz是惯性导航系统正常模式输出的东向、北向和天向速度;
惯性导航系统正常模式输出参数:经度信息λ,纬度信息L,东向速度vx、北向速度vy和天向速度vz;惯性导航系统极区模式输出参数:横东向速度
Figure BDA00003582704600091
横北向速度
Figure BDA00003582704600092
和横天向速度
Figure BDA00003582704600093
本发明具体实施方式:
步骤一,惯性导航系统行至纬度为70度时开始由正常模式转换为极区模式;
步骤二,采集惯性导航系统正常模式输出的位置信息:经度λ和地理纬度L;
步骤三,测量惯性导航系统所在位置与东经90°经线的经度差δλ:
δλ=|λ-90|
步骤四,利用步骤二中的纬度信息L计算辅助系数W:
W = 1 - e 2 sin 2 L .
步骤五,利用步骤四中的W测量地理坐标系与横地理坐标系之间的转角θ:
θ = arctan ( W sin δλ sin L cos δλ ( 1 - e 2 ) ) .
步骤六,利用步骤五中得到的地理坐标系与横地理坐标系之间的转角θ构造地理坐标系与横地理坐标系之间的方向余弦矩阵
Figure BDA00003582704600101
C t t ‾ = cos θ sin θ 0 - sin cos θ 0 0 0 1 .
步骤七,由步骤六中得到的测量惯性导航系统极区模式的初始速度:
v x ‾ v y ‾ v z ‾ = cos θ sin θ 0 - sin cos θ 0 0 0 1 v x v y v z
其中
Figure BDA00003582704600105
Figure BDA00003582704600106
是惯性导航系统进入极区模式所需横东向、横北向和横天向速度,vx、vy和vz是惯性导航系统正常模式输出的东向、北向和天向速度。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (7)

1.一种惯性导航系统极区导航模式初始速度的测量方法,其特征在于,该惯性导航系统极区导航模式初始速度的测量方法包括以下步骤:
采集惯性导航系统正常模式下输出信息;
测量惯性导航系统与东经90°之间的经度差;
计算辅助系数W;
测量地理坐标系与横地理坐标系之间的转角;
计算地理坐标系与横地理坐标系之间的方向余弦矩阵;
测量惯性惯性导航系统极区模式的初始速度。
2.如权利要求1所述的惯性导航系统极区导航模式初始速度的测量方法,其特征在于,该惯性导航系统极区导航模式初始速度的测量方法具体步骤为:
步骤一:采集惯性导航系统正常模式输出的经度信息λ,纬度信息L和速度信息v;
步骤二:测量惯性导航系统所在位置与东经90°经线的经度差δλ;
步骤三:利用步骤一中的纬度信息L计算辅助系数W;
步骤四:利用步骤三中的W测量地理坐标系与横地理坐标系之间的转角θ;
步骤五:利用步骤六中得到的地理坐标系与横地理坐标系之间的转角θ构造地理坐标系与横地理坐标系之间的方向余弦矩阵
步骤六:由步骤五中得到的
Figure FDA00003582704500026
测量惯性导航系统极区模式的初始速度。
3.如权利要求2所述的惯性导航系统极区导航模式初始速度的测量方法,其特征在于,在步骤二中,经度差δλ表示为:δλ=|λ-90|。
4.如权利要求2所述的惯性导航系统极区导航模式初始速度的测量方法,其特征在于,在步骤三中:辅助系数W表示为:
W = 1 - e 2 s in 2 L .
5.如权利要求2所述的惯性导航系统极区导航模式初始速度的测量方法,其特征在于,在步骤四中:地理坐标系与横地理坐标系之间的转角θ表示为:
θ = arctan ( W sin δλ sin L cos δλ ( 1 - e 2 ) ) .
6.如权利要求2所述的惯性导航系统极区导航模式初始速度的测量方法,其特征在于,在步骤五中:地理坐标系与横地理坐标系之间的方向余弦矩阵
Figure FDA00003582704500024
表示为:
C t t ‾ = cos θ sin θ 0 - sin θ cos θ 0 0 0 1 .
7.如权利要求2所述的惯性导航系统极区导航模式初始速度的测量方法,其特征在于,在步骤六中:惯性导航系统极区模式的初始速度表示为:
v x ‾ v y ‾ v z ‾ = cos θ sin θ 0 - sin θ cos θ 0 0 0 1 v x v y v z
其中
Figure FDA00003582704500032
Figure FDA00003582704500033
是惯性导航系统进入极区模式所需横东向、横北向和横天向速度,vx、vy和vz是惯性导航系统正常模式输出的东向、北向和天向速度。
CN2013103223223A 2013-07-29 2013-07-29 一种惯性导航系统极区导航模式初始速度的测量方法 Pending CN103389090A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2013103223223A CN103389090A (zh) 2013-07-29 2013-07-29 一种惯性导航系统极区导航模式初始速度的测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2013103223223A CN103389090A (zh) 2013-07-29 2013-07-29 一种惯性导航系统极区导航模式初始速度的测量方法

Publications (1)

Publication Number Publication Date
CN103389090A true CN103389090A (zh) 2013-11-13

Family

ID=49533435

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2013103223223A Pending CN103389090A (zh) 2013-07-29 2013-07-29 一种惯性导航系统极区导航模式初始速度的测量方法

Country Status (1)

Country Link
CN (1) CN103389090A (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103940428A (zh) * 2014-03-21 2014-07-23 哈尔滨工程大学 一种惯性导航系统横速度的实时测量方法
CN103983277A (zh) * 2014-05-16 2014-08-13 哈尔滨工程大学 一种适用于极区的惯导系统综合校正方法
CN103983277B (zh) * 2014-05-16 2016-11-30 哈尔滨工程大学 一种适用于极区的惯导系统综合校正方法
CN117470233A (zh) * 2023-11-10 2024-01-30 中国人民解放军国防科技大学 一种优化的地球椭球模型下横坐标系组合导航方法
CN117470234A (zh) * 2023-11-10 2024-01-30 中国人民解放军国防科技大学 基于Psi角误差模型的舰船跨极区滤波切换方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5448486A (en) * 1993-04-29 1995-09-05 Honeywell Inc. Orthogonal polar coordinate system to accommodate polar navigation
CN103090867A (zh) * 2012-11-02 2013-05-08 哈尔滨工程大学 相对地心惯性系旋转的光纤陀螺捷联惯性导航系统误差抑制方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5448486A (en) * 1993-04-29 1995-09-05 Honeywell Inc. Orthogonal polar coordinate system to accommodate polar navigation
CN103090867A (zh) * 2012-11-02 2013-05-08 哈尔滨工程大学 相对地心惯性系旋转的光纤陀螺捷联惯性导航系统误差抑制方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘文超等: "游移方位INS/GPS组合导航在极地地区的应用", 《火力与指挥控制》 *
周琪等: "极区飞行格网惯性导航算法原理", 《西北工业大学学报》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103940428A (zh) * 2014-03-21 2014-07-23 哈尔滨工程大学 一种惯性导航系统横速度的实时测量方法
CN103983277A (zh) * 2014-05-16 2014-08-13 哈尔滨工程大学 一种适用于极区的惯导系统综合校正方法
CN103983277B (zh) * 2014-05-16 2016-11-30 哈尔滨工程大学 一种适用于极区的惯导系统综合校正方法
CN117470233A (zh) * 2023-11-10 2024-01-30 中国人民解放军国防科技大学 一种优化的地球椭球模型下横坐标系组合导航方法
CN117470234A (zh) * 2023-11-10 2024-01-30 中国人民解放军国防科技大学 基于Psi角误差模型的舰船跨极区滤波切换方法
CN117470234B (zh) * 2023-11-10 2024-04-26 中国人民解放军国防科技大学 基于Psi角误差模型的舰船跨极区滤波切换方法
CN117470233B (zh) * 2023-11-10 2024-06-21 中国人民解放军国防科技大学 一种优化的地球椭球模型下横坐标系组合导航方法

Similar Documents

Publication Publication Date Title
CN101858748B (zh) 高空长航无人机的多传感器容错自主导航方法
CN106767752A (zh) 一种基于偏振信息的组合导航方法
CN104049269B (zh) 一种基于激光测距和mems/gps组合导航系统的目标导航测绘方法
CN104374388A (zh) 一种基于偏振光传感器的航姿测定方法
CN103471616A (zh) 一种动基座sins大方位失准角条件下初始对准方法
CN104111445A (zh) 用于室内导航的超声波阵列辅助定位方法与系统
CN105021198A (zh) 一种基于多传感器综合导航的位置估计方法
CN103439727B (zh) 一种地面坐标的测量方法
CN103217699B (zh) 一种基于偏振信息的组合导航系统递推优化初始对准方法
CN109059964A (zh) 一种基于重力峰的惯性导航与重力测量双校准方法
CN104215242A (zh) 一种基于横向游移坐标系的极区惯性导航方法
CN103900566B (zh) 一种消除地球自转角速度对旋转调制型捷联惯导系统精度影响的方法
CN103852081A (zh) 用于大气数据/捷联惯导组合导航系统的真空速解算方法
CN105928515A (zh) 一种无人机导航系统
CN106093992A (zh) 一种基于cors的亚米级组合定位导航系统及导航方法
CN105180953A (zh) 基于气压传感的车辆导航系统的控制方法
CN103454662B (zh) 一种基于ckf的sins/北斗/dvl组合对准方法
CN103398725A (zh) 一种基于星敏感器的捷联惯导系统初始对准的方法
CN103017787A (zh) 适用于摇摆晃动基座的初始对准方法
CN107966145A (zh) 一种基于稀疏长基线紧组合的auv水下导航方法
CN102707080B (zh) 一种星敏感器模拟捷联惯导陀螺的方法
CN103575297B (zh) 基于卫星导航接收机的gnss和mimu组合导航航向角估计方法
CN103453907A (zh) 基于分层大气模型的行星进入段导航滤波方法
CN103712621A (zh) 偏振光及红外传感器辅助惯导系统定姿方法
CN106643726B (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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20131113