CN103411623A - 速率陀螺校准方法 - Google Patents

速率陀螺校准方法 Download PDF

Info

Publication number
CN103411623A
CN103411623A CN2013102897322A CN201310289732A CN103411623A CN 103411623 A CN103411623 A CN 103411623A CN 2013102897322 A CN2013102897322 A CN 2013102897322A CN 201310289732 A CN201310289732 A CN 201310289732A CN 103411623 A CN103411623 A CN 103411623A
Authority
CN
China
Prior art keywords
axis
omega
rate
prime
carrier
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
CN2013102897322A
Other languages
English (en)
Other versions
CN103411623B (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.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical 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 Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN201310289732.2A priority Critical patent/CN103411623B/zh
Publication of CN103411623A publication Critical patent/CN103411623A/zh
Application granted granted Critical
Publication of CN103411623B publication Critical patent/CN103411623B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Gyroscopes (AREA)

Abstract

本发明公开了一种速率陀螺校准方法,用于解决现有速率陀螺误差标定过程复杂的技术问题。技术方案是将载体分别绕X轴、Y轴、Z轴旋转,由磁传感器的输出数据变化计算载体的旋转角速率,以此作为标定速率陀螺的基准角速率,分别校准速率陀螺的每个轴的标度因数和零偏误差。该方法无需速率转台,实现条件简单,适用于在线快速标定。经测试,陀螺静态精度能够达到十分钟漂移在2.5°以内,动态条件下,三个轴在输入角速率小于25°/s时,误差均能保持在±2°以内。

Description

速率陀螺校准方法
技术领域
本发明涉及一种陀螺校准方法,特别是涉及一种速率陀螺校准方法。
背景技术
速率陀螺的校准主要是在速率转台上采用位置标定法或速率标定法来完成。文献“微小型捷联惯性测量单元标定及补偿方法,李建利,房建成,盛蔚等,宇航学报,2008,03:947-951”公开了一种微小型捷联惯性测量单元标定及补偿方法。该方法根据微小型惯性测量单元角速度及加速度通道误差数学模型设计“六方位正反速率”试验,利用三轴速率转台设置微型捷联惯性测量单元的X、Y、Z轴分别与地球坐标系的天、地重合等六方位,在每个方位输入不同的正负转速。将六方位每个速率下的正反旋转采集数据分别求平均,根据6个方位的系统误差方程,标定MEMS陀螺的零偏、标度因数、安装误差和与比力有关项等误差系数。该方法虽然误差辨识精度较高,但是依赖于三轴速率转台的速率试验,而且要采集不同正反速率下转台转动数周的六方位数据,试验过程复杂,标定时间冗长,求解线性方程计算量也比较大,不利于在线标定。
发明内容
为了克服现有速率陀螺误差标定过程复杂的不足,本发明提供一种速率陀螺校准方法。该方法将载体分别绕X轴、Y轴、Z轴旋转,由磁传感器的输出数据变化计算载体的旋转角速率,以此作为标定速率陀螺的基准角速率,分别校准速率陀螺的每个轴的标度因数和零偏误差。该方法无需速率转台,实现条件简单,适用于在线快速标定。
本发明解决其技术问题所采用的技术方案是:一种速率陀螺校准方法,其特点是包括以下步骤:
步骤一、建立速率陀螺的误差模型为:
ω xm ω ym ω zm = N x 0 0 0 N y 0 0 0 N z ω x ω y ω z + B ωx B ωy B ωz - - - ( 1 )
式中,ωxm、ωym、ωzm分别为MEMS陀螺X、Y、Z轴输出角速度的测量值;ωx、ωy、ωz分别为陀螺X、Y、Z轴的实际输入角速度值;Nx、Ny、Nz分别为陀螺X、Y、Z轴的标度因数;Bωx、Bωy、Bωz分别为陀螺X、Y、Z轴的零偏误差。
步骤二、将安装有三轴磁传感器和三轴速率陀螺的航向姿态测量系统固定在水平的三轴无磁转台上,航姿测量系统的Z轴指天,X轴和Y轴与转台内框中轴线重合。
步骤三、固定航向角
Figure BDA00003492146700029
倾斜角γ=0°,载体绕X轴逆时针缓慢旋转。定义hx、hy、hz分别为T时刻三轴磁传感器输出的地磁场矢量在X、Y、Z轴上的分量,h'x、h'y、h'z为T+Ts时刻磁传感器的三轴输出数据。根据坐标旋转公式,求得Ts时间内载体绕X轴的旋转角度Δθ为:
Δθ = arctan ( h y ′ h z - h y h z ′ h y h y ′ + h z h z ′ ) - - - ( 2 )
选择磁传感器的数据更新频率为旋转时间差Ts,载体绕X轴旋转的基准角速率ωx为:
ω x = Δθ T s - - - ( 3 )
步骤四、固定航向角
Figure BDA00003492146700023
俯仰角θ=0°,载体绕Y轴逆时针缓慢旋转,根据坐标旋转公式,求得Ts时间内载体绕Y轴的旋转角度Δγ为:
Δγ = arctan ( - h x ′ h z - h x h z ′ h x h x ′ + h z h z ′ ) - - - ( 4 )
载体绕Y轴旋转的基准角速率ωy为:
ω y = Δγ T s - - - ( 5 )
步骤五、固定俯仰角θ=0°,倾斜角γ=0°,载体绕Z轴顺时针缓慢旋转,根据坐标旋转公式,求得Ts时间内绕Z轴的旋转角度
Figure BDA00003492146700026
为:
Figure BDA00003492146700027
载体绕Z轴旋转的基准角速率ωz为:
Figure BDA00003492146700028
步骤六、根据步骤三至步骤五得到的三个基准速率与陀螺输出的三轴角速率的测量值拟合(1)式,得到速率陀螺X、Y、Z轴的标度因数和零偏误差。
本发明的有益效果是:由于该方法将载体分别绕X轴、Y轴、Z轴旋转,由磁传感器的输出数据变化计算载体的旋转角速率,以此作为标定速率陀螺的基准角速率,分别校准速率陀螺的每个轴的标度因数和零偏误差。该方法无需速率转台,实现条件简单,适用于在线快速标定。经测试,陀螺静态精度能够达到十分钟漂移在2.5°以内,动态条件下,三个轴在输入角速率小于25°/s时,误差均能保持在±2°以内。
下面结合具体实施方式对本发明作详细说明。
具体实施方式
本发明速率陀螺校准方法具体步骤如下:
以室温条件下,由三轴MEMS加速度计、三轴MEMS角速率陀螺和三轴微型磁通门传感器构成的微型航向姿态测量系统在3SK-150型无磁转台上进行的实验为例对本发明进行说明。方法步骤如下:
(1)建立速率陀螺的误差模型为:
ω xm ω ym ω zm = N x 0 0 0 N y 0 0 0 N z ω x ω y ω z + B ωx B ωy B ωz - - - ( 1 )
其中,ωxm、ωym、ωzm分别为MEMS陀螺X、Y、Z轴输出角速度的测量值;ωx、ωy、ωz分别为X、Y、Z轴的实际输入角速度值;Nx、Ny、Nz分别为X、Y、Z轴的标度因数;Bωx、Bωy、Bωz分别为陀螺X、Y、Z轴的零偏误差。
(2)将航向姿态系统固定在水平放置的转台上,系统的Z轴指天,X轴和Y轴与转台内框中轴线重合。
(3)固定航向角
Figure BDA00003492146700032
倾斜角γ=0°,载体绕X轴逆时针缓慢旋转。定义hx、hy、hz分别为T时刻三轴磁传感器输出的地磁场矢量在X、Y、Z轴上的分量,h'x、h'y、h'z为T+Ts时刻磁传感器输出的地磁场矢量在X、Y、Z轴上的分量。根据坐标旋转公式,求得Ts时间内载体绕X轴的旋转角度Δθ为:
Δθ = arctan ( h y ′ h z - h y h z ′ h y h y ′ + h z h z ′ ) - - - ( 2 )
因陀螺数据更新率大于磁传感器数据更新率,所以选择磁传感器的数据更新频率为旋转时间差Ts,载体绕X轴旋转的基准角速率ωx为:
ω x = Δθ T s - - - ( 3 )
(4)固定航向角
Figure BDA00003492146700042
俯仰角θ=0°,载体绕Y轴逆时针缓慢旋转,根据坐标旋转公式,求得Ts时间内载体绕Y轴的旋转角度Δγ为:
Δγ = arctan ( - h x ′ h z - h x h z ′ h x h x ′ + h z h z ′ ) - - - ( 4 )
载体绕Y轴旋转的基准角速率ωy为:
ω y = Δγ T s - - - ( 5 )
(5)固定俯仰角θ=0°,倾斜角γ=0°,载体绕Z轴顺时针缓慢旋转,根据坐标旋转公式,求得Ts时间内绕Z轴的旋转角度
Figure BDA00003492146700045
为:
Figure BDA00003492146700046
载体绕Z轴旋转的基准角速率ωz为:
Figure BDA00003492146700047
(6)利用递推最小二乘法拟合(1)式,得到速率陀螺X、Y、Z轴的标度因数误差和零偏误差分别为:
N x = 1.1910 , B ωx = 0.0042 N y = 1.2450 , B ωy = 0.0027 N z = 1.2198 , B ωz = 0.0015 - - - ( 8 )
将误差系数代入误差方程并验证,静态条件下,航姿测量系统测量的航向角误差由补偿前十分钟漂移一百多度减小到补偿后漂移在2.5°以内,动态条件下,三个轴在输入角速率小于25°/s时,误差均能保持在±2°以内。

Claims (1)

1.一种速率陀螺校准方法,其特征在于包括以下步骤:
步骤一、建立速率陀螺的误差模型为:
ω xm ω ym ω zm = N x 0 0 0 N y 0 0 0 N z ω x ω y ω z + B ωx B ωy B ωz - - - ( 1 ) 式中,ωxm、ωym、ωzm分别为MEMS陀螺X、Y、Z轴输出角速度的测量值;ωx、ωy、ωz分别为陀螺X、Y、Z轴的实际输入角速度值;Nx、Ny、Nz分别为陀螺X、Y、Z轴的标度因数;Bωx、Bωy、Bωz分别为陀螺X、Y、Z轴的零偏误差;
步骤二、将安装有三轴磁传感器和三轴速率陀螺的航向姿态测量系统固定在水平的三轴无磁转台上,航姿测量系统的Z轴指天,X轴和Y轴与转台内框中轴线重合;
步骤三、固定航向角
Figure FDA00003492146600012
倾斜角γ=0°,载体绕X轴逆时针缓慢旋转;定义hx、hy、hz分别为T时刻三轴磁传感器输出的地磁场矢量在X、Y、Z轴上的分量,h'x、h'y、h'z为T+Ts时刻磁传感器的三轴输出数据;根据坐标旋转公式,求得Ts时间内载体绕X轴的旋转角度Δθ为:
Δθ = arctan ( h y ′ h z - h y h z ′ h y h y ′ + h z h z ′ ) - - - ( 2 )
选择磁传感器的数据更新频率为旋转时间差Ts,载体绕X轴旋转的基准角速率ωx为:
ω x = Δθ T s - - - ( 3 )
步骤四、固定航向角
Figure FDA00003492146600015
俯仰角θ=0°,载体绕Y轴逆时针缓慢旋转,根据坐标旋转公式,求得Ts时间内载体绕Y轴的旋转角度Δγ为:
Δγ = arctan ( - h x ′ h z - h x h z ′ h x h x ′ + h z h z ′ ) - - - ( 4 )
载体绕Y轴旋转的基准角速率ωy为:
ω y = Δγ T s
步骤五、固定俯仰角θ=0°,倾斜角γ=0°,载体绕Z轴顺时针缓慢旋转,根据坐标旋转公式,求得Ts时间内绕Z轴的旋转角度
Figure FDA00003492146600018
为:
Figure FDA00003492146600021
载体绕Z轴旋转的基准角速率ωz为:
Figure FDA00003492146600022
步骤六、根据步骤三至步骤五得到的三个基准速率与陀螺输出的三轴角速率的测量值拟合(1)式,得到速率陀螺X、Y、Z轴的标度因数和零偏误差。
CN201310289732.2A 2013-07-10 2013-07-10 速率陀螺校准方法 Expired - Fee Related CN103411623B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310289732.2A CN103411623B (zh) 2013-07-10 2013-07-10 速率陀螺校准方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310289732.2A CN103411623B (zh) 2013-07-10 2013-07-10 速率陀螺校准方法

Publications (2)

Publication Number Publication Date
CN103411623A true CN103411623A (zh) 2013-11-27
CN103411623B CN103411623B (zh) 2015-12-23

Family

ID=49604651

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310289732.2A Expired - Fee Related CN103411623B (zh) 2013-07-10 2013-07-10 速率陀螺校准方法

Country Status (1)

Country Link
CN (1) CN103411623B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104101363A (zh) * 2014-07-28 2014-10-15 中国电子科技集团公司第二十六研究所 一种用于旋转载体横向姿态测量的陀螺仪动态标定方法
CN105445764A (zh) * 2015-11-13 2016-03-30 上海华测导航技术股份有限公司 一种用于gnss接收机空间姿态精度检验方法
CN106033131A (zh) * 2015-03-20 2016-10-19 阿里巴巴集团控股有限公司 一种地磁传感器校准方法、装置及智能设备
CN106840205A (zh) * 2017-01-19 2017-06-13 北京小鸟看看科技有限公司 陀螺仪校准补偿方法及装置、虚拟现实头戴设备
CN107024674A (zh) * 2017-05-26 2017-08-08 北京航空航天大学 一种基于递推最小二乘法的磁强计现场快速标定方法
CN107727115A (zh) * 2017-09-27 2018-02-23 青岛海信移动通信技术股份有限公司 陀螺仪校正方法及装置
CN110530396A (zh) * 2019-07-31 2019-12-03 苏州浪潮智能科技有限公司 一种陀螺仪传感器误差的处理方法及设备
CN113008271A (zh) * 2019-08-15 2021-06-22 深圳市瑞立视多媒体科技有限公司 校准3d旋转差异的数学模型构建方法、校准方法及其装置
CN114166246A (zh) * 2021-11-19 2022-03-11 西安航天精密机电研究所 一种速率光纤陀螺仪标度因数标定方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101029833A (zh) * 2007-03-12 2007-09-05 北京航空航天大学 一种捷联mems陀螺动态误差标定方法
CN101246023A (zh) * 2008-03-21 2008-08-20 哈尔滨工程大学 微机械陀螺惯性测量组件的闭环标定方法
US20090089001A1 (en) * 2007-08-14 2009-04-02 American Gnc Corporation Self-calibrated azimuth and attitude accuracy enhancing method and system (SAAAEMS)
CN101568801A (zh) * 2006-12-20 2009-10-28 萨甘安全防护公司 校准轴对称振动陀螺测试仪中的比例因子的方法
CN102564451A (zh) * 2011-12-09 2012-07-11 北京理工大学 一种平台式惯性导航系统的系统级补偿方法
CN102564459A (zh) * 2012-01-17 2012-07-11 北京理工大学 一种单轴旋转调制捷联式惯性导航系统海上校准方法
CN102680004A (zh) * 2012-05-30 2012-09-19 北京航空航天大学 一种挠性陀螺位置姿态测量系统pos的标度因数误差标定与补偿方法
CN102927994A (zh) * 2012-10-23 2013-02-13 北京航空航天大学 一种斜置冗余捷联惯性导航系统的快速标定方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101568801A (zh) * 2006-12-20 2009-10-28 萨甘安全防护公司 校准轴对称振动陀螺测试仪中的比例因子的方法
CN101029833A (zh) * 2007-03-12 2007-09-05 北京航空航天大学 一种捷联mems陀螺动态误差标定方法
US20090089001A1 (en) * 2007-08-14 2009-04-02 American Gnc Corporation Self-calibrated azimuth and attitude accuracy enhancing method and system (SAAAEMS)
CN101246023A (zh) * 2008-03-21 2008-08-20 哈尔滨工程大学 微机械陀螺惯性测量组件的闭环标定方法
CN102564451A (zh) * 2011-12-09 2012-07-11 北京理工大学 一种平台式惯性导航系统的系统级补偿方法
CN102564459A (zh) * 2012-01-17 2012-07-11 北京理工大学 一种单轴旋转调制捷联式惯性导航系统海上校准方法
CN102680004A (zh) * 2012-05-30 2012-09-19 北京航空航天大学 一种挠性陀螺位置姿态测量系统pos的标度因数误差标定与补偿方法
CN102927994A (zh) * 2012-10-23 2013-02-13 北京航空航天大学 一种斜置冗余捷联惯性导航系统的快速标定方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
李建利等: "微小型捷联惯性测量单元标定及补偿方法", 《宇航学报》, vol. 29, no. 3, 31 May 2008 (2008-05-31) *
袁赣南等: "MEMS陀螺随机漂移的状态空间模型分析及应用", 《传感技术学报》, vol. 24, no. 6, 30 June 2011 (2011-06-30) *
郭鹏飞等: "斜置惯性测量单元的一体化标定技术", 《中国惯性技术学报》, vol. 15, no. 3, 30 June 2007 (2007-06-30) *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104101363A (zh) * 2014-07-28 2014-10-15 中国电子科技集团公司第二十六研究所 一种用于旋转载体横向姿态测量的陀螺仪动态标定方法
CN106033131A (zh) * 2015-03-20 2016-10-19 阿里巴巴集团控股有限公司 一种地磁传感器校准方法、装置及智能设备
CN105445764A (zh) * 2015-11-13 2016-03-30 上海华测导航技术股份有限公司 一种用于gnss接收机空间姿态精度检验方法
CN106840205A (zh) * 2017-01-19 2017-06-13 北京小鸟看看科技有限公司 陀螺仪校准补偿方法及装置、虚拟现实头戴设备
CN107024674A (zh) * 2017-05-26 2017-08-08 北京航空航天大学 一种基于递推最小二乘法的磁强计现场快速标定方法
CN107024674B (zh) * 2017-05-26 2019-04-26 北京航空航天大学 一种基于递推最小二乘法的磁强计现场快速标定方法
CN107727115A (zh) * 2017-09-27 2018-02-23 青岛海信移动通信技术股份有限公司 陀螺仪校正方法及装置
CN110530396A (zh) * 2019-07-31 2019-12-03 苏州浪潮智能科技有限公司 一种陀螺仪传感器误差的处理方法及设备
CN113008271A (zh) * 2019-08-15 2021-06-22 深圳市瑞立视多媒体科技有限公司 校准3d旋转差异的数学模型构建方法、校准方法及其装置
CN114166246A (zh) * 2021-11-19 2022-03-11 西安航天精密机电研究所 一种速率光纤陀螺仪标度因数标定方法
CN114166246B (zh) * 2021-11-19 2024-04-16 西安航天精密机电研究所 一种速率光纤陀螺仪标度因数标定方法

Also Published As

Publication number Publication date
CN103411623B (zh) 2015-12-23

Similar Documents

Publication Publication Date Title
CN103411623B (zh) 速率陀螺校准方法
CN108168574B (zh) 一种基于速度观测的8位置捷联惯导系统级标定方法
CN110160554B (zh) 一种基于寻优法的单轴旋转捷联惯导系统标定方法
CN102980577B (zh) 一种微型捷联航姿系统及其工作方法
CN107655493B (zh) 一种光纤陀螺sins六位置系统级标定方法
CN102221372B (zh) 使用离心机和转台对惯性测量单元进行误差标定的方法
CN106052595B (zh) 基于激光陀螺捷联惯导的三轴转台轴线垂直度检测方法
CN101290326B (zh) 石英挠性加速度计测量组件的参数辨识标定方法
CN102879832B (zh) 用于地磁要素测量系统的非对准误差校正方法
CN104459828B (zh) 基于绕轴法的地磁矢量系统非对准校正方法
CN103323625B (zh) 一种mems-imu中加速度计动态环境下的误差标定补偿方法
CN104807475B (zh) 动中通卫星天线倾角校准过程中零点漂移值的测量方法
CN103697918B (zh) 一种三轴正交一轴斜置构型光纤陀螺惯测装置的标定方法
CN105180968A (zh) 一种imu/磁强计安装失准角在线滤波标定方法
CN103196462A (zh) 一种mimu中mems陀螺仪的误差标定补偿方法
CN101246023A (zh) 微机械陀螺惯性测量组件的闭环标定方法
CN102636183B (zh) 基于光纤监测和双轴转台离心机的挠性陀螺二次过载项测试方法
CN106969783A (zh) 一种基于光纤陀螺惯性导航的单轴旋转快速标定技术
CN112595350B (zh) 一种惯导系统自动标定方法及终端
CN101963512A (zh) 船用旋转式光纤陀螺捷联惯导系统初始对准方法
CN104697521B (zh) 一种采用陀螺冗余斜交配置方式测量高速旋转体姿态和角速度的方法
CN106940193A (zh) 一种基于Kalman滤波的船舶自适应摇摆标定方法
CN107402007A (zh) 一种提高微型ahrs模块精度的方法和微型ahrs模块
CN105628025A (zh) 一种恒速偏频/机抖激光陀螺惯导系统导航方法
CN102636184B (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
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: 20151223

Termination date: 20160710