CN101694390A - 基于光纤惯性测量系统的舰船升沉运动测量方法 - Google Patents

基于光纤惯性测量系统的舰船升沉运动测量方法 Download PDF

Info

Publication number
CN101694390A
CN101694390A CN200910073076A CN200910073076A CN101694390A CN 101694390 A CN101694390 A CN 101694390A CN 200910073076 A CN200910073076 A CN 200910073076A CN 200910073076 A CN200910073076 A CN 200910073076A CN 101694390 A CN101694390 A CN 101694390A
Authority
CN
China
Prior art keywords
sin
cos
gamma
beta
semi
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
Application number
CN200910073076A
Other languages
English (en)
Other versions
CN101694390B (zh
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 Ship Navigation Technology Co., Ltd.
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 CN2009100730766A priority Critical patent/CN101694390B/zh
Publication of CN101694390A publication Critical patent/CN101694390A/zh
Application granted granted Critical
Publication of CN101694390B publication Critical patent/CN101694390B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Navigation (AREA)

Abstract

本发明提供的是一种基于光纤惯性测量系统的舰船升沉运动测量方法。实时采集光纤陀螺捷联惯性系统的输出数据;得到载体坐标系b与地理坐标系t的关系矩阵;得到载体坐标系与半固定坐标系之间的关系矩阵;得到半固定坐标系下垂直轴方向的平动加速度;得到该坐标系下第N个采样点的速度量;对速度量进行滤波处理;得到半固定坐标系下垂直轴方向平动位移值。本发明利用现有的惯性测量系统,通过一定的量测方法,提供舰船升沉运动信息,从而提高舰船的可操纵性、完善导航信息内容。

Description

基于光纤惯性测量系统的舰船升沉运动测量方法
(一)技术领域
本发明涉及的是一种测量舰船升沉运动的方法,具体地说是在没有外界信息帮助下,利用系统的姿态信息与加速度计输出信息,自主测量舰船的升沉运动的即时位移地方法。
(二)背景技术
应用于舰船的惯性测量系统是由惯性测量组件和导航计算机组成,它实时地测量舰船的位置信息以及舰船的姿态信息。由于其在测量过程中不需要任何外部信息,也不向外辐射能量,自主性强,隐蔽性好,测量精度高且具有实时性,所以它被广泛地应用各类舰船的导航系统中,并发挥着重要作用,在国防建设中具有重要意义和巨大的经济效益。
惯性测量技术是一项发展较为成熟、精度较高、稳定性优良的自主导航技术,在军用和民用领域都发挥着越来越重要的作用。原有的惯性测量系统可以测得舰船的位置、速度、姿态信息,但并不能测出升沉运动的距离量。由于海况复杂,这种由海浪引起的运动已经成为舰船运动中重要的一部分,在船舶多种海上作业中均会产生影响。舰船升沉运动信息可用于舰载机的起落、武器发射、船舶靠港入岸等领域,还可以增强舰船的操纵性。在以往的方法中,为舰船操纵者提供舰船升沉运动信息多是通过海浪监测、船舶运动建模等手段完成。这种方法中,由于模型复杂,运算量大,其间不可避免的简化过程会带来较大的误差,并不能真实有效的反应舰船六自由度(横摇、纵摇、艏摇、升沉、纵荡、横荡)运动信息。如果惯性测量系统在提供舰船三个姿态信息基础上,还能够提供三个轴向平动(横荡、纵荡及升沉)中对舰船影响最大的升沉运动信息,那么就可以规避建模复杂问题,实现对舰船任意时刻的四自由度(同样方法可推广得到六自由度)信息同步、真实的测量,为操纵者提供准确及时地参考量。
(三)发明内容
本发明的目的在于提供一种能够提高舰船的可操纵性、完善导航信息内容的基于光纤惯性测量系统的舰船升沉运动测量方法。
本发明的目的是这样实现的:主要包括如下实施步骤:
步骤1、实时采集光纤陀螺捷联惯性系统充分预热后的光纤陀螺和各轴加速度计的输出数据;
步骤2、利用导航计算机中的已有方法,经过初始对准后,得到较为精确的舰船实时姿态信息(横摇角α,纵摇角β,俯仰角γ),进一步得到载体坐标系b与地理坐标系t的关系矩阵
C b t = cos β cos γ - sin β sin α sin γ - cos α sin γ sin β cos γ + cos β sin α sin γ cos β sin γ + sin β sin α cos γ cos α cos γ sin β sin γ - cos β sin α cos γ - sin β cos α sin α cos β cos α ;
步骤3、根据舰船行使操控时设定的航向,得到舰船的主航向角
Figure G2009100730766D0000022
,得到半固定坐标系D与地理坐标系t的关系矩阵Ct D,进而得到载体坐标系与半固定坐标系之间的关系矩阵
Figure G2009100730766D0000023
步骤4、结合步骤3中的关系矩阵Cb D,将由传感器测的载体坐标系下三个轴向的加速度通过坐标转化,得到半固定坐标系下垂直轴方向的平动加速度ai D
步骤5、将半固定坐标系下的垂直轴方向平动加速度进行一次积分,采样时间间隔为h,得到该坐标系下第N个采样点的速度量:
v D ( N ) = Σ i = 1 N a i D * h ( i = 1,2,3 , . . . , N ) ;
步骤6、选用高通滤波器
H ( j ω c ) = ( j ω c 1 + j ω c ) 3
根据海浪情况的不同设置截止频率,截止频率ωc要大于海浪频率的5倍,对速度量vD进行滤波处理,得到速度量
v=vD*H(jωc)
进而得到每个时刻的半固定坐标系下的垂直轴方向平动即升沉的速度vi
步骤7、将半固定坐标系垂直轴方向平动速度值进行积分,得到半固定坐标系下垂直轴方向平动位移值:
s i ( N ) = Σ i = 0 N v i * h ( i = 1,2,3 , . . . , N )
即舰船操作者需要的舰船的升沉位移值。
本发明技术具有以下优点:在不增加新的传感器,不需要外部信息的情况下,利用舰船中安装的中高精度光纤惯性系统输出的姿态信息和加速度信息,实时的真实的提供舰船升沉运动信息,增强原有惯性测量系统功能。
对本发明的有益效果通过如下多种动态条件下利用高精度光纤陀螺惯性系统针对垂直轴平动信息即升沉运动的测量试验结果加以验证。
试验条件:
(1)光纤陀螺捷联惯性系统器件精度陀螺常值漂移为0.01度/小时,加速度计的随机常值偏置为0.0001g。
(2)平台可以模拟规则海浪周期运动以及瞬时运动。
试验结果:单次试验时间视不同的运动模式而定。利用四自由度转台(可模拟升沉、横摇、纵摇、艏摇运动)系统模拟舰船升沉周期运动,以及瞬时平动。选取高精度光纤陀螺捷联惯性系统,将其安装与转台台面上,模拟舰船运动。分别进行振荡幅值为0.2米、振荡周期1.5秒的升沉测量,振荡幅值为0.6米、振荡周期15秒的升沉测量,以及瞬时0.2米平动的实验。实验所得结果证明对于周期运动的测量结果误差低于2%(|实测值-真实值|/真实值),瞬时运动的误差在5%以下。测量值的延迟时间较短,可以忽略不计。
(四)附图说明
图1基于光纤捷联惯性测量系统的舰船升沉运动测量技术流程图。
图2平台振荡幅值为0.2米、振荡周期1.5秒的升沉运动时本方法的测量结果。
图3平台振荡幅值为0.6米、振荡周期15秒的升沉运动时本方法的测量结果。
图4平台平动距离为0.2米,瞬时下降时本方法的测量结果。
(五)具体实施方式
下面结合附图举例对本发明做更详细地描述:
本实施方式中,将舰船捷联惯性系统安装在舰船中心位置,选取精度较高的光纤陀螺作为其主要元器件。结合图1,舰船升沉运动测量的具体实施步骤如下:
步骤1、将光纤陀螺捷联惯性系统充分预热后,实时采集其光纤陀螺和各轴加速度计输出数据。预热时间根据具体系统设定。
步骤2、利用导航计算机中的已有方法,经过初始对准后,得到较为精确的舰船实时姿态信息(俯仰角α、横滚角β、艏摇角γ,均为小角度),进一步得到,载体坐标系b与地理坐标系t的关系矩阵
C b t = cos β cos γ - sin β sin α sin γ - cos α sin γ sin β cos γ + cos β sin α sin γ cos β sin γ + sin β sin α cos γ cos α cos γ sin β sin γ - cos β sin α cos γ - sin β cos α sin α cos β cos α
步骤3、根据舰船行使操控时设定的航向,确定舰船的主航向角
Figure G2009100730766D0000042
(该角有别于航向角,是由舰船操纵者提供)。得到半固定坐标系D与地理坐标系t的关系矩阵Ct D。进而得到载体坐标系与半固定坐标系之间的关系矩阵
Figure G2009100730766D0000043
步骤4、结合步骤3中的方向余弦矩阵Cb D,将由传感器测的载体坐标系下三个轴向的加速度通过坐标转化,得到半固定坐标系下垂直轴方向的平动加速度ai D
步骤5、将半固定坐标系下的垂直轴方向平动加速度进行一次积分,采样时间间隔为h,得到该坐标系下第N个采样点的速度量:
v D ( N ) = Σ i = 1 N a i D * h ( i = 1,2,3 , . . . , N )
步骤6、选用高通滤波器
H ( j ω c ) = ( j ω c 1 + j ω c ) 3
根据海浪情况的不同设置截止频率,截止频率ωc要大于海浪频率的5倍,对速度量vD进行滤波处理,得到速度量
v=vD*H(jωc)
进而得到每个时刻的半固定坐标系下的垂直轴方向平动(升沉)的速度vi
步骤7、将半固定坐标系垂直轴方向平动速度值进行积分,得到半固定坐标系下垂直轴方向平动位移值:
s i ( N ) = Σ i = 0 N v i * h ( i = 1,2,3 , . . . , N )
即舰船操作者需要的舰船的升沉位移值。
本发明还可以包括如下特征:
1)所得到的舰船升沉信息(垂直轴向平动信息)与光纤捷联惯性系统传感器精度相关性较大,也就是说需要的系统精度尤其是其传感器达到中等精度以上水平。
2)如果要得到精确的舰船升沉运动测量值,需要在引入高通滤波器之后稳定5分钟左右。
3)如果舰船以周期振荡形式进行升沉往复运动,该测量系统会得准确的呈周期性的量测值。如果舰船沿半固定坐标系的垂直轴方向进行短距离平移运动,该测量系统会输出以起始o点为基准的瞬时升沉值,然后逐渐归零(如图4)。归零时间与选取的滤波器的截止频率ωc有关。

Claims (2)

1.一种基于光纤惯性测量系统的舰船升沉运动测量方法,其特征是包括如下实施步骤:
步骤1、实时采集光纤陀螺捷联惯性系统充分预热后的光纤陀螺和各轴加速度计的输出数据;
步骤2、利用导航计算机中的已有方法,经过初始对准后,得到较为精确的舰船实时姿态信息即横摇角α、纵摇角β、俯仰角γ,进一步得到载体坐标系b与地理坐标系t的关系矩阵
C b t = cos β cos γ - sin β sin α sin γ - cos α sin γ sin β cos γ + cos β sin α sin γ cos β sin γ + sin β sin α cos γ cos α cos γ sin β sin γ - cos β sin α cos γ - sin β cos α sin α cos β cos α ;
步骤3、根据舰船行使操控时设定的航向,得到舰船的主航向角
Figure F2009100730766C0000012
得到半固定坐标系D与地理坐标系t的关系矩阵Ct D,进而得到载体坐标系与半固定坐标系之间的关系矩阵
Figure F2009100730766C0000013
步骤4、结合步骤3中的关系矩阵,将由传感器测的载体坐标系下三个轴向的加速度通过坐标转化,得到半固定坐标系下垂直轴方向的平动加速度ai D
步骤5、将半固定坐标系下的垂直轴方向平动加速度进行一次积分,采样时间间隔为h,得到该坐标系下第N个采样点的速度量:
v D ( N ) = Σ i = 1 N a i D * h ( i = 1,2,3 , . . . , N ) ;
步骤6、选用高通滤波器
H ( jω c ) = ( jω c 1 + jω c ) 3
根据海浪情况的不同设置截止频率ωc,对速度量vD进行滤波处理,得到速度量
v=vD*H(jωc)
进而得到每个时刻的半固定坐标系下的垂直轴方向平动即升沉的速度vi
步骤7、将半固定坐标系垂直轴方向平动速度值进行积分,得到半固定坐标系下垂直轴方向平动位移值:
s i ( N ) = Σ i = 0 N v i * h ( i = 1,2,3 , . . . , N )
即舰船操作者需要的舰船的升沉位移值。
2.根据权利要求1所述的基于光纤惯性测量系统的舰船升沉运动测量方法,其特征是所述根据海浪情况的不同设置截止频率是:截止频率ωc大于海浪频率的5倍。
CN2009100730766A 2009-10-20 2009-10-20 基于光纤惯性测量系统的舰船升沉运动测量方法 Active CN101694390B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009100730766A CN101694390B (zh) 2009-10-20 2009-10-20 基于光纤惯性测量系统的舰船升沉运动测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009100730766A CN101694390B (zh) 2009-10-20 2009-10-20 基于光纤惯性测量系统的舰船升沉运动测量方法

Publications (2)

Publication Number Publication Date
CN101694390A true CN101694390A (zh) 2010-04-14
CN101694390B CN101694390B (zh) 2012-08-22

Family

ID=42093376

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009100730766A Active CN101694390B (zh) 2009-10-20 2009-10-20 基于光纤惯性测量系统的舰船升沉运动测量方法

Country Status (1)

Country Link
CN (1) CN101694390B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105698789A (zh) * 2016-01-12 2016-06-22 西北工业大学 舰船升沉测量方法及其测量系统
CN106643728A (zh) * 2016-12-16 2017-05-10 哈尔滨工程大学 基于自适应频率估计的船舶升沉运动信息估计方法
KR20180016476A (ko) * 2015-06-18 2018-02-14 후아웨이 테크놀러지 컴퍼니 리미티드 고성능 이동 프런트홀을 위한 임베딩된 제어 신호에 의한 캐스케이드 파형 변조
CN109425339A (zh) * 2017-08-21 2019-03-05 哈尔滨工程大学 一种基于惯性技术的考虑杆臂效应的舰船升沉误差补偿方法
CN109631949A (zh) * 2018-12-29 2019-04-16 重庆华渝电气集团有限公司 一种光纤捷联惯导设备测试装置及测试方法
CN110319838A (zh) * 2019-07-09 2019-10-11 哈尔滨工程大学 一种自适应的运动姿态参考系统升沉测量方法
CN112629540A (zh) * 2020-12-16 2021-04-09 北京航天控制仪器研究所 一种基于载体姿态信息的升沉测量方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6697736B2 (en) * 2002-02-06 2004-02-24 American Gnc Corporation Positioning and navigation method and system thereof
CN100491915C (zh) * 2004-12-30 2009-05-27 中国科学院自动化研究所 一种自动舵航向控制系统及其方法

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20180016476A (ko) * 2015-06-18 2018-02-14 후아웨이 테크놀러지 컴퍼니 리미티드 고성능 이동 프런트홀을 위한 임베딩된 제어 신호에 의한 캐스케이드 파형 변조
KR101986964B1 (ko) 2015-06-18 2019-06-07 후아웨이 테크놀러지 컴퍼니 리미티드 고성능 이동 프런트홀을 위한 임베딩된 제어 신호에 의한 캐스케이드 파형 변조
CN105698789A (zh) * 2016-01-12 2016-06-22 西北工业大学 舰船升沉测量方法及其测量系统
CN106643728A (zh) * 2016-12-16 2017-05-10 哈尔滨工程大学 基于自适应频率估计的船舶升沉运动信息估计方法
CN109425339A (zh) * 2017-08-21 2019-03-05 哈尔滨工程大学 一种基于惯性技术的考虑杆臂效应的舰船升沉误差补偿方法
CN109631949A (zh) * 2018-12-29 2019-04-16 重庆华渝电气集团有限公司 一种光纤捷联惯导设备测试装置及测试方法
CN109631949B (zh) * 2018-12-29 2020-05-19 重庆华渝电气集团有限公司 一种光纤捷联惯导设备测试装置及测试方法
CN110319838A (zh) * 2019-07-09 2019-10-11 哈尔滨工程大学 一种自适应的运动姿态参考系统升沉测量方法
CN112629540A (zh) * 2020-12-16 2021-04-09 北京航天控制仪器研究所 一种基于载体姿态信息的升沉测量方法
CN112629540B (zh) * 2020-12-16 2024-02-09 北京航天控制仪器研究所 一种基于载体姿态信息的升沉测量方法

Also Published As

Publication number Publication date
CN101694390B (zh) 2012-08-22

Similar Documents

Publication Publication Date Title
CN101694390B (zh) 基于光纤惯性测量系统的舰船升沉运动测量方法
CN103575299B (zh) 利用外观测信息的双轴旋转惯导系统对准及误差修正方法
CN103017755B (zh) 一种水下导航姿态测量方法
CN101706284B (zh) 提高船用光纤陀螺捷联惯导系统定位精度的方法
CN101793521B (zh) 基于光纤陀螺惯性测量系统的舰船横纵荡信息测量方法
CN101975872B (zh) 石英挠性加速度计组件零位偏置的标定方法
CN103852085B (zh) 一种基于最小二乘拟合的光纤捷联惯导系统现场标定方法
CN101424534B (zh) 惯性/重力组合导航半实物模拟装置
CN101162147A (zh) 大失准角下船用光纤陀螺捷联航姿系统系泊精对准方法
CN101246023A (zh) 微机械陀螺惯性测量组件的闭环标定方法
CN106767917A (zh) 一种斜置冗余惯导系统标定误差模型建模方法
CN103234560B (zh) 一种捷联惯导系统零位标定的方法
CN111537002B (zh) 一种激光捷联惯组安装误差的标定方法和定向方法
CN101178313A (zh) 适合于光纤陀螺捷联惯性导航系统的地速检测方法
CN102706361A (zh) 一种高精度多惯导系统姿态精度评定方法
CN105136166B (zh) 一种指定惯导位置精度的捷联惯导系统误差模型仿真方法
CN106767671B (zh) 基于三维电子罗盘的地质结构面产状计算方法
CN101183004A (zh) 一种在线实时消除光纤陀螺捷联惯导系统振荡误差的方法
CN103454662B (zh) 一种基于ckf的sins/北斗/dvl组合对准方法
CN102445176A (zh) 高速列车运行姿态参数测量系统
CN104482941A (zh) 一种船用光学惯导长航时保精度导航的系统级补偿方法
CN106767925A (zh) 带双轴转位机构的惯导系统三位置参数辨识对准方法
CN108362269B (zh) 一种基于航标体的波浪测量系统及测量方法
CN103575276A (zh) 一种双轴旋转惯性导航系统初始对准模型降阶方法
CN106123917B (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
C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20160920

Address after: 15 Heilongjiang, Nangang Province, Nantong street, building No. 258, building, ship, floor, No. 150001

Patentee after: Science Park Development Co., Ltd. of Harbin Engineering University

Patentee after: Zhao Yuxin

Address before: 150001 Heilongjiang, Nangang District, Nantong street, building No. 145, Harbin Engineering University, Intellectual Property Office

Patentee before: Harbin Engineering Univ.

C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20161021

Address after: 15 Heilongjiang, Nangang Province, Nantong street, building No. 258, building, ship, floor, No. 150001

Patentee after: Science Park Development Co., Ltd. of Harbin Engineering University

Patentee after: Harbin poly flame investment enterprise (limited partnership)

Address before: 15 Heilongjiang, Nangang Province, Nantong street, building No. 258, building, ship, floor, No. 150001

Patentee before: Science Park Development Co., Ltd. of Harbin Engineering University

Patentee before: Zhao Yuxin

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20170315

Address after: 150078 Harbin hi tech Industrial Development Zone Yingbin Road, the focus of the Russian park on the ground floor of the building 2D, No., East unit, level 2, level 22

Patentee after: Harbin Ship Navigation Technology Co., Ltd.

Address before: 15 Heilongjiang, Nangang Province, Nantong street, building No. 258, building, ship, floor, No. 150001

Patentee before: Science Park Development Co., Ltd. of Harbin Engineering University

Patentee before: Harbin poly flame investment enterprise (limited partnership)