CN110657821A - 一种轨道检测惯性测量装置的标定方法 - Google Patents

一种轨道检测惯性测量装置的标定方法 Download PDF

Info

Publication number
CN110657821A
CN110657821A CN201810700844.5A CN201810700844A CN110657821A CN 110657821 A CN110657821 A CN 110657821A CN 201810700844 A CN201810700844 A CN 201810700844A CN 110657821 A CN110657821 A CN 110657821A
Authority
CN
China
Prior art keywords
calibration
temperature
inertia
zero position
along
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
CN201810700844.5A
Other languages
English (en)
Other versions
CN110657821B (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.)
Beijing Automation Control Equipment Institute BACEI
Original Assignee
Beijing Automation Control Equipment Institute BACEI
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 Beijing Automation Control Equipment Institute BACEI filed Critical Beijing Automation Control Equipment Institute BACEI
Priority to CN201810700844.5A priority Critical patent/CN110657821B/zh
Publication of CN110657821A publication Critical patent/CN110657821A/zh
Application granted granted Critical
Publication of CN110657821B publication Critical patent/CN110657821B/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
    • G01C25/00Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
    • G01C25/005Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass initial alignment, calibration or starting-up of inertial devices

Landscapes

  • Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Gyroscopes (AREA)

Abstract

本发明属于惯性测量检测技术领域,具体涉及一种轨道检测惯性测量装置的标定方法。本方法分两个步骤,第一步为惯性测量装置的温度标定,第二步为惯性测量装置的温度梯度标定。首先是对温度进行标定,采用常规的六位置、六速率标定试验,解算出惯性器件的零位和标度因数,用其对惯性器件数据进行修正;然后对温度梯度进行标定,将温度梯度对惯性器件的影响等效成惯性器件误差随通电时间的变化规律。本标定方法首次将惯性器件随温度梯度的变化规律等效成随通电时间的变化规律,并利用惯性器件的标定参数简化惯性测量装置的标定方法,技术简单实用、有效,具有创新性,应用该新型标定技术能够大大提高标定精度,同时又能高效完成测试任务,达到了研究设计目标。

Description

一种轨道检测惯性测量装置的标定方法
技术领域
本发明属于惯性测量检测技术领域,具体涉及一种轨道检测惯性测量装置的标定方法。
背景技术
惯性测量装置中的惯性器件(陀螺仪和加速度计)在在装入系统后,由于残余应力释放、温度环境影响、安装结构变形等原因,其参数会发生一定的改变,因此都需要对惯性测量装置进行标定。目前国内对惯性测量装置常用的标定方法是在全温度范围内选取若干个(一般为八个或十个)温度点,每个温度点通过六位置、六速率等常规方法进行标定,然后将每个温度点标定出的惯性器件(陀螺仪和加速度计)参数采用分段线性拟合或者二阶非线性拟合等方式进行处理。用于轨道检测的惯性测量装置与一般惯性测量装置不同,具有工作时间长、无需正交等特点,采用常规标定方法不仅费时费力,而且标定精度不高,不足以满足轨道检测对惯性器件的需求。
本成果提供的标定方法就是在这种背景下产生的,该方法不仅节约了标定时间,而且标定精度有很大提高,完全满足要求。
惯性测量装置中的惯性器件(陀螺仪和加速度计)在生产过程中,都要对零位和标度因数的温度系数进行详细标定,但是在装入惯性测量装置后,由于温度梯度等方面的影响,惯性器件的零位和标度因数会发生变化,因此需要在惯性测量装置中重新进行标定。
发明内容
本成果发明了一种利用惯性器件的温度标定系数,在惯性测量装置中对温度梯度进行标定同时对惯性器件误差进行修正的方法,不仅提高惯性测量装置的精度,而且提高了标定效率。
为了实现这一目的,本发明采取的技术方案是:
一种轨道检测惯性测量装置的标定方法,本方法分两个步骤,第一步为惯性测量装置的温度标定,第二步为惯性测量装置的温度梯度标定。首先是对温度进行标定,采用常规的六位置、六速率标定试验,解算出惯性器件的零位和标度因数,用其对惯性器件数据进行修正;然后对温度梯度进行标定,将温度梯度对惯性器件的影响等效成惯性器件误差随通电时间的变化规律。
一种轨道检测惯性测量装置的标定方法,所述第一步惯性测量装置的温度标定,即在全温度范围内选取三个温度点,分别为低温、常温、高温点,在每个温度点首先进行常规的六位置、六速率标定试验,通过常规的标定解算方法解算出惯性器件的零位和标度因数,将惯性器件的零位和标度因数与单器件标定时的数值进行对比,得出差值,三个温度点共得三个差值,然后取平均做为修正值对惯性器件的标定参数进行修正,从而实现惯性测量装置的标定;
X轴加速度计的零位的修正公式如下:
式中:ΔKax0为X轴加速度计零位修正量,ΔKax0g、ΔKax0c、ΔKax0d分别为高温、常温、低温时X轴加速度计零位修正量。用ΔKax0作为修正量对X轴加速度计的器件级标定参数进行修正,作为系统级标定参数写入惯性测量装置的flash中,实现误差补偿;
利用上述的方法可以得出Y、Z加速度计零位,X、Y、Z加速度计标度因数的修正量;
对于陀螺仪,也用上述同样的方法修正,X轴陀螺仪的零位的公式如下:
式中:ΔKgx0为X轴陀螺仪零位修正量,ΔKgx0g、ΔKgx0c、ΔKgx0d分别为高温、常温、低温时X轴陀螺仪零位修正量。用ΔKgx0作为修正量对X轴陀螺仪的器件级标定参数进行修正,作为系统级标定参数写入惯性测量装置的flash中,实现温度误差补偿;
利用上述的方法可以得出Y、Z陀螺仪零位,X、Y、Z陀螺仪标度因数的修正量。
一种轨道检测惯性测量装置的标定方法,所述第二步为惯性测量装置的温度梯度标定温度误差补偿完成后,接下来进行温度梯度误差补偿;
将惯性器件随温度梯度的变化规律等效成随通电时间的变化规律,选取三个位置,每个位置采集惯性器件的输出数据,根据三个位置分离出陀螺仪和加速度计的零位、标度因数随通电时间的变化规律,然后采用如下公式进行拟合;
X轴加速度计的零位的拟合公式如下:
Kax=Kax0+Kax1×t+Kax2×t2 (3)
式中:Kax为X轴加速度计零位,Kax0、Kax1、Kax2为X轴加速度计零位的随时间变化的系数,t为通电时间。采用最小二乘法便可计算出Kax0、Kax1、Kax2,将Kax0、Kax1、Kax2代入到公式(3)中即可得到X轴加速度计零位随通电时间的变化关系,也就是X轴加速度计零位随温度梯度的变化关系;
用上述同样的方法可以得到X加速度计的标度因数、Y和Z加速度计的零位和标度因数随时间变化的关系,也就是随温度梯度的变化关系;
对于陀螺仪,采用用同样的方法,X轴陀螺仪的零位计算公式如下:
Kgx=Kgx0+Kgx1×t+Kgx2×t2 (4)
式中:Kgx为X轴陀螺仪零位,Kgx0、Kgx1、Kgx2为X轴陀螺仪零位的随时间变化的系数,t为通电时间。采用最小二乘法便可计算出Kgx0、Kgx1、Kgx2,将Kgx0、Kgx1、Kgx2代入到公式(4)中即可得到X轴陀螺仪零位随通电时间的变化关系,也就是X轴陀螺仪零位随温度梯度的变化关系;
同样的方法计算得到X轴陀螺仪的标度因数、Y和Z陀螺仪的零位和标度因数随时间变化的关系,随温度梯度的变化关系;
将以上补偿公式和补偿参数写入惯性测量装置的flash中,从而实现惯性器件误差随温度梯度的变化关系。
一种轨道检测惯性测量装置的标定方法,所述每个位置采集惯性器件的输出数据的采集时间为40min。
本发明的有益效果为:本标定方法首次将惯性器件随温度梯度的变化规律等效成随通电时间的变化规律,并利用惯性器件的标定参数简化惯性测量装置的标定方法,技术简单实用、有效,具有创新性,应用该新型标定技术能够大大提高标定精度,同时又能高效完成测试任务,达到了研究设计目标。
附图说明
图1为轨道检测惯性测量装置的标定方法。
具体实施方式
下面结合附图及具体实施例对本发明详细说明。
下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整的叙述,显然,所描述的实施例仅是本发明的一部分实施例,而不是全部的实施例。基于本发明的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明的保护范围。
本方法分两个步骤,分别实现惯性测量装置的温度标定和温度梯度标定。首先是对温度进行标定,采用常规的六位置、六速率标定试验,解算出惯性器件的零位和标度因数,用其对惯性器件数据进行修正;然后对温度梯度进行标定,将温度梯度对惯性器件的影响等效成惯性器件误差随通电时间的变化规律。
首先为温度标定。全温度范围内选取三个温度点,分别为低温、常温、高温点,在每个温度点首先进行常规的六位置、六速率标定试验,通过常规的标定解算方法解算出惯性器件(三个陀螺仪和三个加速度计)的零位和标度因数,将惯性器件的零位和标度因数与单器件标定时的数值进行对比,得出差值,三个温度点共得三个差值,然后取平均做为修正值对惯性器件的标定参数进行修正,从而实现惯性测量装置的标定。
以X轴加速度计的零位为例,修正公式如下:
Figure BDA0001714331320000051
式中:ΔKax0为X轴加速度计零位修正量,ΔKax0g、ΔKax0c、ΔKax0d分别为高温、常温、低温时X轴加速度计零位修正量。用ΔKax0作为修正量对X轴加速度计的器件级标定参数进行修正,作为系统级标定参数写入惯性测量装置的flash中,实现误差补偿。
同样的方法可以得出Y、Z加速度计零位,X、Y、Z加速度计标度因数的修正量。
对于陀螺仪来说,也用同样的方法修正,以X轴陀螺仪的零位为例,公式如下:
Figure BDA0001714331320000052
式中:ΔKgx0为X轴陀螺仪零位修正量,ΔKgx0g、ΔKgx0c、ΔKgx0d分别为高温、常温、低温时X轴陀螺仪零位修正量。用ΔKgx0作为修正量对X轴陀螺仪的器件级标定参数进行修正,作为系统级标定参数写入惯性测量装置的flash中,实现温度误差补偿。利用上述的方法可以得出Y、Z陀螺仪零位,X、Y、Z陀螺仪标度因数的修正量。
温度误差补偿完成后,接下来进行温度梯度误差补偿。
实施方式是将惯性器件随温度梯度的变化规律等效成随通电时间的变化规律,具体做法为:选取三个位置,每个位置采集惯性器件的输出数据,采集时间为40min,根据三个位置很容易分离出陀螺仪和加速度计的零位、标度因数随通电时间的变化规律,然后采用如下公式进行拟合。
同样以X轴加速度计的零位为例,拟合公式如下:
Kax=Kax0+Kax1×t+Kax2×t2 (3)
式中:Kax为X轴加速度计零位,Kax0、Kax1、Kax2为X轴加速度计零位的随时间变化的系数,t为通电时间。采用最小二乘法便可计算出Kax0、Kax1、Kax2,将Kax0、Kax1、Kax2代入到公式(3)中即可得到X轴加速度计零位随通电时间的变化关系,也就是X轴加速度计零位随温度梯度的变化关系。
同样的方法可以得到X加速度计的标度因数、Y和Z加速度计的零位和标度因数随时间变化的关系,也就是随温度梯度的变化关系。
对于陀螺仪来说,也用同样的方法,以X轴陀螺仪的零位为例,公式如下:
Kgx=Kgx0+Kgx1×t+Kgx2×t2 (4)
式中:Kgx为X轴陀螺仪零位,Kgx0、Kgx1、Kgx2为X轴陀螺仪零位的随时间变化的系数,t为通电时间。采用最小二乘法便可计算出Kgx0、Kgx1、Kgx2,将Kgx0、Kgx1、Kgx2代入到公式(4)中即可得到X轴陀螺仪零位随通电时间的变化关系,也就是X轴陀螺仪零位随温度梯度的变化关系。
同样的方法可以得到X轴陀螺仪的标度因数、Y和Z陀螺仪的零位和标度因数随时间变化的关系,随温度梯度的变化关系。
将以上补偿公式和补偿参数写入惯性测量装置的flash中,从而实现惯性器件误差随温度梯度的变化关系。

Claims (4)

1.一种轨道检测惯性测量装置的标定方法,其特征在于:本方法分两个步骤,第一步为惯性测量装置的温度标定,第二步为惯性测量装置的温度梯度标定。首先是对温度进行标定,采用常规的六位置、六速率标定试验,解算出惯性器件的零位和标度因数,用其对惯性器件数据进行修正;然后对温度梯度进行标定,将温度梯度对惯性器件的影响等效成惯性器件误差随通电时间的变化规律。
2.如权利要求1所述的一种轨道检测惯性测量装置的标定方法,其特征在于:所述第一步惯性测量装置的温度标定,即在全温度范围内选取三个温度点,分别为低温、常温、高温点,在每个温度点首先进行常规的六位置、六速率标定试验,通过常规的标定解算方法解算出惯性器件的零位和标度因数,将惯性器件的零位和标度因数与单器件标定时的数值进行对比,得出差值,三个温度点共得三个差值,然后取平均做为修正值对惯性器件的标定参数进行修正,从而实现惯性测量装置的标定;
X轴加速度计的零位的修正公式如下:
Figure FDA0001714331310000011
式中:ΔKax0为X轴加速度计零位修正量,ΔKax0g、ΔKax0c、ΔKax0d分别为高温、常温、低温时X轴加速度计零位修正量。用ΔKax0作为修正量对X轴加速度计的器件级标定参数进行修正,作为系统级标定参数写入惯性测量装置的flash中,实现误差补偿;
利用上述的方法可以得出Y、Z加速度计零位,X、Y、Z加速度计标度因数的修正量;
对于陀螺仪,也用上述同样的方法修正,X轴陀螺仪的零位的公式如下:
Figure FDA0001714331310000021
式中:ΔKgx0为X轴陀螺仪零位修正量,ΔKgx0g、ΔKgx0c、ΔKgx0d分别为高温、常温、低温时X轴陀螺仪零位修正量。用ΔKgx0作为修正量对X轴陀螺仪的器件级标定参数进行修正,作为系统级标定参数写入惯性测量装置的flash中,实现温度误差补偿;
利用上述的方法可以得出Y、Z陀螺仪零位,X、Y、Z陀螺仪标度因数的修正量。
3.如权利要求1所述的一种轨道检测惯性测量装置的标定方法,其特征在于:所述第二步为惯性测量装置的温度梯度标定温度误差补偿完成后,接下来进行温度梯度误差补偿;
将惯性器件随温度梯度的变化规律等效成随通电时间的变化规律,选取三个位置,每个位置采集惯性器件的输出数据,根据三个位置分离出陀螺仪和加速度计的零位、标度因数随通电时间的变化规律,然后采用如下公式进行拟合;
X轴加速度计的零位的拟合公式如下:
Kax=Kax0+Kax1×t+Kax2×t2 (3)
式中:Kax为X轴加速度计零位,Kax0、Kax1、Kax2为X轴加速度计零位的随时间变化的系数,t为通电时间。采用最小二乘法便可计算出Kax0、Kax1、Kax2,将Kax0、Kax1、Kax2代入到公式(3)中即可得到X轴加速度计零位随通电时间的变化关系,也就是X轴加速度计零位随温度梯度的变化关系;
用上述同样的方法可以得到X加速度计的标度因数、Y和Z加速度计的零位和标度因数随时间变化的关系,也就是随温度梯度的变化关系;
对于陀螺仪,采用用同样的方法,X轴陀螺仪的零位计算公式如下:
Kgx=Kgx0+Kgx1×t+Kgx2×t2 (4)
式中:Kgx为X轴陀螺仪零位,Kgx0、Kgx1、Kgx2为X轴陀螺仪零位的随时间变化的系数,t为通电时间。采用最小二乘法便可计算出Kgx0、Kgx1、Kgx2,将Kgx0、Kgx1、Kgx2代入到公式(4)中即可得到X轴陀螺仪零位随通电时间的变化关系,也就是X轴陀螺仪零位随温度梯度的变化关系;
同样的方法计算得到X轴陀螺仪的标度因数、Y和Z陀螺仪的零位和标度因数随时间变化的关系,随温度梯度的变化关系;
将以上补偿公式和补偿参数写入惯性测量装置的flash中,从而实现惯性器件误差随温度梯度的变化关系。
4.如权利要求3所述的一种轨道检测惯性测量装置的标定方法,其特征在于:所述每个位置采集惯性器件的输出数据的采集时间为40min。
CN201810700844.5A 2018-06-29 2018-06-29 一种轨道检测惯性测量装置的标定方法 Active CN110657821B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810700844.5A CN110657821B (zh) 2018-06-29 2018-06-29 一种轨道检测惯性测量装置的标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810700844.5A CN110657821B (zh) 2018-06-29 2018-06-29 一种轨道检测惯性测量装置的标定方法

Publications (2)

Publication Number Publication Date
CN110657821A true CN110657821A (zh) 2020-01-07
CN110657821B CN110657821B (zh) 2021-08-10

Family

ID=69026847

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810700844.5A Active CN110657821B (zh) 2018-06-29 2018-06-29 一种轨道检测惯性测量装置的标定方法

Country Status (1)

Country Link
CN (1) CN110657821B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112082572A (zh) * 2020-08-24 2020-12-15 中国石油天然气集团有限公司 一种标定钻具磁干扰的装置及方法
CN114964223A (zh) * 2022-04-19 2022-08-30 北京自动化控制设备研究所 一种微惯导系统测量误差补偿方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6175807B1 (en) * 1999-06-09 2001-01-16 Litton Systems, Inc. Temperature compensation method for strapdown inertial navigation systems
CN1314946C (zh) * 2005-11-04 2007-05-09 北京航空航天大学 一种消除陀螺常值漂移影响的惯性测量单元混合标定方法
CN102853850B (zh) * 2012-09-11 2015-07-01 中国兵器工业集团第二一四研究所苏州研发中心 基于单轴转台的三轴mems陀螺旋转积分标定方法
CN107621269B (zh) * 2016-07-15 2020-08-04 北京计算机技术及应用研究所 光纤陀螺温度漂移误差补偿方法
CN106500728B (zh) * 2016-10-26 2019-08-02 北京小鸟看看科技有限公司 一种陀螺仪的温度漂移补偿方法和装置

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112082572A (zh) * 2020-08-24 2020-12-15 中国石油天然气集团有限公司 一种标定钻具磁干扰的装置及方法
CN112082572B (zh) * 2020-08-24 2023-04-25 中国石油天然气集团有限公司 一种标定钻具磁干扰的装置及方法
CN114964223A (zh) * 2022-04-19 2022-08-30 北京自动化控制设备研究所 一种微惯导系统测量误差补偿方法
CN114964223B (zh) * 2022-04-19 2023-09-12 北京自动化控制设备研究所 一种微惯导系统测量误差补偿方法

Also Published As

Publication number Publication date
CN110657821B (zh) 2021-08-10

Similar Documents

Publication Publication Date Title
CN110006450B (zh) 一种激光捷联惯导系统在卧式三轴转台上的标定方法
KR102008597B1 (ko) 자세 측정 시스템의 정밀 보정 방법
CN109631940B (zh) 一种四环惯性稳定平台框架零位标校方法
CN103389114B (zh) 基于g敏感性的三轴MEMS陀螺非正交误差标定方法
CN102564455B (zh) 星敏感器安装误差四位置标定与补偿方法
CN103823084A (zh) 一种三轴加速度传感器的校准方法
CN110657821B (zh) 一种轨道检测惯性测量装置的标定方法
CN101246023A (zh) 微机械陀螺惯性测量组件的闭环标定方法
CN104121928B (zh) 一种适用于低精度有方位基准单轴转位设备的惯性测量单元标定方法
CN109827545B (zh) 一种基于双mems加速度计的在线倾角测量方法
CN102679999A (zh) 星敏感器安装误差四位置标定与补偿方法
CN115855046A (zh) 一种捷联惯导系统加速度计六位置温补标定方法
CN102654515B (zh) 三轴加速度传感器的z敏感轴的校准方法
CN111220817B (zh) 一种捷联惯组三轴加速度计的标定方法
CN103994775B (zh) 一种适用于低精度有方位基准双轴转位设备的惯性测量单元标定方法
CN110631605B (zh) 一种陀螺阵列标定方法及系统
CN106897480B (zh) 一种抑制光纤陀螺标度因数模型迟滞效应的建模方法
CN104635193B (zh) 一种电表用的soc芯片实时时钟高精度补偿方法
CN106767620B (zh) 一种高精度位移测量系统的传感器安装检测方法
CN113959464B (zh) 一种陀螺仪辅助的加速度计现场校准方法和系统
CN110186479B (zh) 一种惯性器件误差系数确定方法
CN114413933B (zh) 一种加速度计动态较准方法、系统及存储介质
CN110823255A (zh) 一种基于比力观测无需调平引北的系统级自标定方法
CN113984090B (zh) 一种轮式机器人imu误差在线标定与补偿方法及装置
CN111982151B (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