CN109282804B - 一种单轴光纤陀螺寻北算法 - Google Patents

一种单轴光纤陀螺寻北算法 Download PDF

Info

Publication number
CN109282804B
CN109282804B CN201811015886.1A CN201811015886A CN109282804B CN 109282804 B CN109282804 B CN 109282804B CN 201811015886 A CN201811015886 A CN 201811015886A CN 109282804 B CN109282804 B CN 109282804B
Authority
CN
China
Prior art keywords
north
seeking
data
degrees
module
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
Application number
CN201811015886.1A
Other languages
English (en)
Other versions
CN109282804A (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 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 CN201811015886.1A priority Critical patent/CN109282804B/zh
Publication of CN109282804A publication Critical patent/CN109282804A/zh
Application granted granted Critical
Publication of CN109282804B publication Critical patent/CN109282804B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C19/00Gyroscopes; Turn-sensitive devices using vibrating masses; Turn-sensitive devices without moving masses; Measuring angular rate using gyroscopic effects
    • G01C19/58Turn-sensitive devices without moving masses
    • G01C19/64Gyrometers using the Sagnac effect, i.e. rotation-induced shifts between counter-rotating electromagnetic beams

Abstract

本发明属于惯性寻北定向测量领域,具体涉及一种单轴光纤陀螺寻北算法。具体包括如下步骤:采集四个位置的陀螺、加速度计数据;然后解算出粗寻北值;再解算出精寻北值。相应的寻北装置包括采集模块、控制模块、解算模块和通信模块,加速度计信号由A/D采样模块进行信号采集,而光纤陀螺输出则直接通过串口采集,采集的加速度计信号与光纤陀螺信号被送给控制模块,经过控制模块的初步处理再传送给解算模块,最终由解算模块运用寻北算法解算出寻北结果并传送给控制模块,控制模块再通过串口通信模块将结果发出,其他设备可通过串口通信模块与寻北设备进行通信,本发明精度更高,稳定行更好,应用前景广阔。

Description

一种单轴光纤陀螺寻北算法
技术领域
本发明属于惯性寻北定向测量领域,具体涉及一种单轴光纤陀螺寻北算法。
背景技术
光纤陀螺是基于Sagnac效应的新一代角速率传感器,寻北定向是光纤陀螺的重要应用之一。在军事领域,陀螺寻北仪为武器系统快速、机动、精确打击提供了保障;在民用领域,陀螺寻北仪为矿山开采、资源勘测、隧道施工等提供了方位基准。其应用前景越来越广阔。寻北技术通常情况下是采用静态寻北方案,即在某几个位置处采集陀螺数据,根据不同位置陀螺敏感的地球自转角速度的分量,列出方程求解北向夹角。但存在以下缺陷:现有的寻北仪一般采用机械陀螺作为寻北核心部件,但机械陀螺存在可靠性差,寿命短,容易受加速度影响,导致测试精度不够,稳定性不好;现有的寻北仪一般采用单陀螺和双加速度计,增大了系统的体积,也增加了制造成本;寻北仪的精度与寻北位置的准确性有很大关系,现有寻北仪系统采用光电编码器来确定转位系统的四个位置,保证寻北位置的准确性,但光电编码器存在体积大,价格高的缺陷;同时光电编码器在使用时需要安装在转轴轴心处,故提高了寻北仪装置整体结构设计的难度,最终导致寻北仪具有体积大,成本高的缺陷。
发明内容
本发明的目的在于提供一种精度高的单轴光纤陀螺寻北算法。
一种单轴光纤陀螺寻北算法,具体包括如下步骤:
步骤1、采集四个位置的陀螺、加速度计数据;
步骤2、根据步骤1数据解算出粗寻北值;
步骤3、解算出精寻北值。
所述一种单轴光纤陀螺寻北算法,步骤1具体包括如下步骤:
步骤1.1、寻北仪把初始位置作为第一个位置并在这个位置静止采集数据;
步骤1.2、采集一段时间数据后,寻北仪顺时针旋转90°到达第二个位置并在此位置采集同样时长的数据;
步骤1.3、第二个位置数据采集完成后,寻北仪再顺时针旋90°到达第三个位置并在此位置采集同样时长的数据;
步骤1.4、第三个位置数据采集完成后,寻北仪再顺时针旋90°到达第四个位置并在此位置采集同样时长的数据。
所述一种单轴光纤陀螺寻北算法,步骤2具体包括如下步骤:
步骤2.1、运用上面四个位置的数据,将陀螺仪测得的其敏感轴上的地球自转角速度的分量用如下公式表示:
Figure BDA0001786079420000021
其中
Figure BDA0001786079420000022
为载体坐标系下角速度三轴分量,ωie=15.0411°/h,为地球自转角速度,L为地理纬度,α,θ,γ分别为航向角,俯仰角和横滚角;
步骤2.2、取陀螺敏感轴方向为东向,分别用(0)、(1)表示起始位置和终止位置且这两个位置相距180°,则有如下关系式;
Figure BDA0001786079420000023
Figure BDA0001786079420000024
步骤2.3、求出俯仰角与横滚角,公式如下:
Figure BDA0001786079420000025
其中
Figure BDA0001786079420000026
为载体坐标系下加速度在X轴,Y轴的分量;
步骤2.4、将式(2)与式(3)作差整理得:
A=Bsinα+Ccosα (5)
其中
Figure BDA0001786079420000027
B=cosγ0+cosγ1 (7)
C=sinγ0sinθ0+sinγ1sinθ1 (8)
步骤2.5、求出航向角α,公式如下:
Figure BDA0001786079420000028
步骤2.6、二、三位置解得一个寻北值α1;二,四位置解得一个寻北值α2,通过下式将两个寻北值加权平均求粗寻北值:
α=cos2α1·α1+cos2α2·α2 (10)
所述一种单轴光纤陀螺寻北算法,步骤3具体包括如下步骤:
步骤3.1、根据步骤2解算出的粗寻北值使寻北仪精寻北的初始位置在地理西向,以此为精寻北第一个位置并采集一段时间数据,采集完成后使寻北仪顺时针旋转180°到达第二个位置,并在此位置上采集同样时长的数据;
步骤3.2、根据上面两个位置的数据解算出精寻北值,寻北解算公式与步骤2中一样,将粗寻北值与精寻北值相加便得到最终的寻北结果。
本发明的有益效果在于:
本发明提供的一种粗,精分段式寻北的方法,相对于传统四位置寻北算法,传统四位置寻北算法解算出的结果收敛性不如分段式寻北;另一方面,本发明精度更高,稳定行更好,应用前景广阔。
附图说明
图1为光纤陀螺寻北的基本原理图;
图2为初始位置为1时分段式寻北算法效果比较图;
图3为初始位置为2时分段式寻北算法效果比较图;
图4为初始位置为3时分段式寻北算法效果比较图;
图5为初始位置为4时分段式寻北算法效果比较图;
图6为全温条件下分段式寻北算法效果比较图;
图7为全温温度曲线图;
图8为单轴光纤陀螺寻北算法的装置框图。
具体实施方式
下面结合附图对本发明做进一步描述:
如附图1所示,为光纤陀螺寻北的基本原理图;如图所示地球自转角速率在地理北向上的分量为ωecosL,设光纤陀螺测得的角速率为ωg,则可以得到寻北值为:
Figure BDA0001786079420000031
图2,图3,图4,图5给出了这两种算法的寻北结果,这4幅图是在四个不同的初始位置。以图2的初始位置为基准,图3,图4,图5的初始位置分别是在图2的初始位置上顺时针旋转90°,180°,270°得到。每个位置下做6次寻北,可以看到:
在每个初始位置下,比起分段式的寻北方案,传统四位置寻北解算出的结果收敛性不如分段式寻北。图2中传统四位置寻北的寻北结果重复性(RMSE)为0.031°,分段式寻北的寻北结果重复性(RMSE)为0.015°。图3中传统四位置寻北的寻北结果重复性(RMSE)为0.049°,分段式寻北的寻北结果重复性(RMSE)为0.012°。图4中传统四位置寻北的寻北结果重复性(RMSE)为0.045°,分段式寻北的寻北结果重复性(RMSE)为0.016°。图5中传统四位置寻北的寻北结果重复性(RMSE)为0.042°,分段式寻北的寻北结果重复性(RMSE)为0.009°。四个初始位置下的传统四位置寻北均值分别为:117.275°,207.300°,297.320°,27.321°。而次四个初始位置下分段式寻北均值分别为:117.317°,207.326°,297.304°,27.297°。若将图3,4,5三个位置的寻北均值换算到图2位置处则分别有:四个初始位置下传统四位置寻北均值为117.275°,117.300°,117.320°,117.321°,其均值极差为0.046°。四个位置下分段式寻北均值为117.317°,117.326°,117.304°,117.297°,其均值极差为0.029°。故无论从重复性上还是均值极差上来看,分段式的寻北方式效果显著。
另外在全温条件下进行两种方法的对比,如图6所示。全温温度曲线如图7所示。传统四位置寻北的寻北结果重复性(RMSE)为0.219°,分段式寻北的寻北结果重复性(RMSE)为0.03°。故在全温条件下分段式的寻北方式效果显著。
结合图8,单轴光纤陀螺寻北装置结构组成包括采集模块、控制模块、解算模块和通信模块。光纤陀螺与加速度计的信号通过采集模块采集,采集到的信号传送给控制模块,经过控制模块对信号的初步处理再发送给解算模块,最终由解算模块运用寻北算法解算出寻北结果并传送给控制模块,控制模块再通过串口通信模块将结果发出。
各模块的构成及功能如下:
模块一:采样模块
采样模块用于接收光纤陀螺和加速度计的原始数据。
模块二:控制模块
控制模块接收到来自采集模块的数据,先进行初步的处理,再将完整的传感器数据发送给解算模块。同时它还负责控制通信与其他设备进行通信。
模块三:解算模块
解算模块就是将每个位置的传感器数据存储下来,最后套用寻北算法计算出寻北值。
模块四:通信模块
通信模块负责与其他设备进行通信,比如将寻北结果向外发送,或从外设接收一些指令。
一种单轴光纤陀螺寻北算法,具体包括如下步骤:
步骤1、采集四个位置的陀螺、加速度计数据;
步骤2、根据步骤1数据解算出粗寻北值;
步骤3、解算出精寻北值。
所述一种单轴光纤陀螺寻北算法,步骤1具体包括如下步骤:
步骤1.1、寻北仪把初始位置作为第一个位置并在这个位置静止采集数据;
步骤1.2、采集一段时间数据后,寻北仪顺时针旋转90°到达第二个位置并在此位置采集同样时长的数据;
步骤1.3、第二个位置数据采集完成后,寻北仪再顺时针旋90°到达第三个位置并在此位置采集同样时长的数据;
步骤1.4、第三个位置数据采集完成后,寻北仪再顺时针旋90°到达第四个位置并在此位置采集同样时长的数据。
所述一种单轴光纤陀螺寻北算法,步骤2具体包括如下步骤:
步骤2.1、运用上面四个位置的数据,将陀螺仪测得的其敏感轴上的地球自转角速度的分量用如下公式表示:
Figure BDA0001786079420000051
其中
Figure BDA0001786079420000052
为载体坐标系下角速度三轴分量,ωie=15.0411°/h,为地球自转角速度,L为地理纬度,α,θ,γ分别为航向角,俯仰角和横滚角;
步骤2.2、取陀螺敏感轴方向为东向,分别用(0)、(1)表示起始位置和终止位置且这两个位置相距180°,则有如下关系式;
Figure BDA0001786079420000053
Figure BDA0001786079420000054
步骤2.3、求出俯仰角与横滚角,公式如下:
Figure BDA0001786079420000055
其中
Figure BDA0001786079420000056
为载体坐标系下加速度在X轴,Y轴的分量;
步骤2.4、将式(2)与式(3)作差整理得:
A=Bsinα+Ccosα (5)
其中
Figure BDA0001786079420000061
B=cosγ0+cosγ1 (7)
C=sinγ0sinθ0+sinγ1sinθ1 (8)
步骤2.5、求出航向角α,公式如下:
Figure BDA0001786079420000062
步骤2.6、二、三位置解得一个寻北值α1;二,四位置解得一个寻北值α2,通过下式将两个寻北值加权平均求粗寻北值:
α=cos2α1·α1+cos2α2·α2 (10)
所述一种单轴光纤陀螺寻北算法,步骤3具体包括如下步骤:
步骤3.1、根据步骤2解算出的粗寻北值使寻北仪精寻北的初始位置在地理西向,以此为精寻北第一个位置并采集一段时间数据,采集完成后使寻北仪顺时针旋转180°到达第二个位置,并在此位置上采集同样时长的数据;
步骤3.2、根据上面两个位置的数据解算出精寻北值,寻北解算公式与步骤2中一样,将粗寻北值与精寻北值相加便得到最终的寻北结果。

Claims (2)

1.一种单轴光纤陀螺寻北算法,其特征在于,具体包括如下步骤:
步骤1、采集四个位置的陀螺、加速度计数据;
步骤2、根据步骤1数据解算出粗寻北值;
步骤3、解算出精寻北值;
所述步骤1具体包括如下步骤:
步骤1.1、寻北仪把初始位置作为第一个位置并在这个位置静止采集数据;
步骤1.2、采集一段时间数据后,寻北仪顺时针旋转90°到达第二个位置并在此位置采集同样时长的数据;
步骤1.3、第二个位置数据采集完成后,寻北仪再顺时针旋90°到达第三个位置并在此位置采集同样时长的数据;
步骤1.4、第三个位置数据采集完成后,寻北仪再顺时针旋90°到达第四个位置并在此位置采集同样时长的数据;
所述步骤2具体包括如下步骤:
步骤2.1、运用上面四个位置的数据,将陀螺仪测得的其敏感轴上的地球自转角速度的分量用如下公式表示:
Figure FDA0003235231240000011
其中
Figure FDA0003235231240000012
为载体坐标系下角速度三轴分量,ωie=15.0411°/h,为地球自转角速度,L为地理纬度,α,θ,γ分别为航向角,俯仰角和横滚角;
步骤2.2、取陀螺敏感轴方向为东向,分别用(0)、(1)表示起始位置和终止位置且这两个位置相距180°,则有如下关系式;
Figure FDA0003235231240000013
Figure FDA0003235231240000014
步骤2.3、求出俯仰角与横滚角,公式如下:
Figure FDA0003235231240000015
其中
Figure FDA0003235231240000016
为载体坐标系下加速度在X轴,Y轴的分量;
步骤2.4、将式(2)与式(3)作差整理得:
A=Bsinα+Ccosα (5)
其中
Figure FDA0003235231240000021
B=cosγ0+cosγ1 (7)
C=sinγ0sinθ0+sinγ1sinθ1 (8)
步骤2.5、求出航向角α,公式如下:
Figure FDA0003235231240000022
步骤2.6、二、三位置解得一个寻北值α1;二,四位置解得一个寻北值α2,通过下式将两个寻北值加权平均求粗寻北值:
α=cos2α1·α1+cos2α2·α2 (10)
2.根据权利要求1所述一种单轴光纤陀螺寻北算法,其特征在于,所述步骤3具体包括如下步骤:
步骤3.1、根据步骤2解算出的粗寻北值使寻北仪精寻北的初始位置在地理西向,以此为精寻北第一个位置并采集一段时间数据,采集完成后使寻北仪顺时针旋转180°到达第二个位置,并在此位置上采集同样时长的数据;
步骤3.2、根据上面两个位置的数据解算出精寻北值,寻北解算公式与步骤2中一样,将粗寻北值与精寻北值相加便得到最终的寻北结果。
CN201811015886.1A 2018-09-01 2018-09-01 一种单轴光纤陀螺寻北算法 Active CN109282804B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811015886.1A CN109282804B (zh) 2018-09-01 2018-09-01 一种单轴光纤陀螺寻北算法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811015886.1A CN109282804B (zh) 2018-09-01 2018-09-01 一种单轴光纤陀螺寻北算法

Publications (2)

Publication Number Publication Date
CN109282804A CN109282804A (zh) 2019-01-29
CN109282804B true CN109282804B (zh) 2021-11-09

Family

ID=65183473

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811015886.1A Active CN109282804B (zh) 2018-09-01 2018-09-01 一种单轴光纤陀螺寻北算法

Country Status (1)

Country Link
CN (1) CN109282804B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110926447B (zh) * 2019-12-16 2022-02-22 重庆华渝电气集团有限公司 一种具有自主导航功能的单轴光纤陀螺寻北方法及航姿导航方法
CN111765880B (zh) * 2020-07-06 2022-07-29 中国人民解放军96901部队22分队 一种基于单光纤陀螺的高精度四位置寻北方法
CN112729264B (zh) * 2020-12-24 2021-11-02 中南大学 一种任意四位置单陀螺寻北方法
CN113432605B (zh) * 2021-06-30 2023-03-31 河北科技大学 基于光纤陀螺的寻北方法及系统
WO2024003894A1 (en) * 2022-06-29 2024-01-04 Israel Aerospace Industries Ltd. Improved system, method and computer program product for north-finding
CN115164942B (zh) * 2022-09-07 2022-11-25 中国船舶重工集团公司第七0七研究所 一种捷联式陀螺寻北仪水平寻北精度自动化测试方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105659809B (zh) * 2005-06-02 2010-05-05 北京航天时代光电科技有限公司 基于倾斜补偿和粗精组合位置法的光纤陀螺寻北仪及寻北方法
CN102207386A (zh) * 2011-03-07 2011-10-05 东南大学 基于方位效应误差补偿的寻北方法
CN102840856A (zh) * 2011-06-24 2012-12-26 西安测绘研究所 一种动态旋转调制的陀螺仪寻北方法
CN103512570A (zh) * 2013-10-12 2014-01-15 北京信息科技大学 一种预置45°带解算失效判定的寻北方法
CN103822624A (zh) * 2014-03-12 2014-05-28 西安北方捷瑞光电科技有限公司 一种单轴激光陀螺寻北仪及其寻北方法
CN103983259A (zh) * 2014-05-29 2014-08-13 中国人民解放军第二炮兵工程大学 基于摆式陀螺最大进动速度粗寻北原理的全方位寻北方法
CN106705997A (zh) * 2016-12-01 2017-05-24 浙江大学 一种光纤陀螺寻北启动误差抑制方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105659809B (zh) * 2005-06-02 2010-05-05 北京航天时代光电科技有限公司 基于倾斜补偿和粗精组合位置法的光纤陀螺寻北仪及寻北方法
CN102207386A (zh) * 2011-03-07 2011-10-05 东南大学 基于方位效应误差补偿的寻北方法
CN102840856A (zh) * 2011-06-24 2012-12-26 西安测绘研究所 一种动态旋转调制的陀螺仪寻北方法
CN103512570A (zh) * 2013-10-12 2014-01-15 北京信息科技大学 一种预置45°带解算失效判定的寻北方法
CN103822624A (zh) * 2014-03-12 2014-05-28 西安北方捷瑞光电科技有限公司 一种单轴激光陀螺寻北仪及其寻北方法
CN103983259A (zh) * 2014-05-29 2014-08-13 中国人民解放军第二炮兵工程大学 基于摆式陀螺最大进动速度粗寻北原理的全方位寻北方法
CN106705997A (zh) * 2016-12-01 2017-05-24 浙江大学 一种光纤陀螺寻北启动误差抑制方法

Also Published As

Publication number Publication date
CN109282804A (zh) 2019-01-29

Similar Documents

Publication Publication Date Title
CN109282804B (zh) 一种单轴光纤陀螺寻北算法
CN107270893B (zh) 面向不动产测量的杆臂、时间不同步误差估计与补偿方法
CN101793523B (zh) 一种组合导航和光电探测一体化系统
CN103900611B (zh) 一种惯导天文高精度复合两位置对准及误差标定方法
CN102564452B (zh) 一种基于惯性导航系统的在线自主标定方法
CN107063254B (zh) 一种陀螺地磁组合的姿态解算方法
CN201159646Y (zh) 捷联式多位置陀螺罗盘
CN103697878B (zh) 一种单陀螺单加速度计旋转调制寻北方法
CN110926447B (zh) 一种具有自主导航功能的单轴光纤陀螺寻北方法及航姿导航方法
CN103743413A (zh) 倾斜状态下调制寻北仪安装误差在线估计与寻北误差补偿方法
CA2738665A1 (en) Method and apparatus for precision azimuth measurement
CN110672871B (zh) 一种基于组合惯导信息和光电转塔视频跟踪的运动目标测速测向方法
CN109459008A (zh) 一种小型中高精度光纤陀螺寻北装置及方法
CN109489661B (zh) 一种卫星初始入轨时陀螺组合常值漂移估计方法
CN102937450A (zh) 一种基于陀螺测量信息的相对姿态确定方法
CN112378399A (zh) 基于捷联惯导和数字全站仪的煤矿巷道掘进机器人精确定位定向方法
CN110514200B (zh) 一种惯性导航系统及高转速旋转体姿态测量方法
CN104296741B (zh) 采用距离平方和距离平方变化率的wsn/ahrs紧组合方法
CN111722295A (zh) 一种水下捷联式重力测量数据处理方法
CN111912402B (zh) 基于地磁信息辅助gps的高旋转载体的测姿方法及装置
CN103487053B (zh) 一种任意二位置捷联寻北方法
CN105486275A (zh) 一种基于9轴惯性测量单元的磁偏角计算方法
Lei et al. Initial alignment for SINS based on low-cost IMU
Liu et al. High precision fiber-optic gyroscope resolution test method based on low precision turntable
CN104897157A (zh) 基于足部航姿参考和肩部电子罗盘的个人导航系统及方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant