CN105241474B - 一种斜置构型惯导***标定方法 - Google Patents

一种斜置构型惯导***标定方法 Download PDF

Info

Publication number
CN105241474B
CN105241474B CN201410326080.XA CN201410326080A CN105241474B CN 105241474 B CN105241474 B CN 105241474B CN 201410326080 A CN201410326080 A CN 201410326080A CN 105241474 B CN105241474 B CN 105241474B
Authority
CN
China
Prior art keywords
accelerometer
axis
gyro
gyroscope
coordinate system
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
CN201410326080.XA
Other languages
English (en)
Other versions
CN105241474A (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 CN201410326080.XA priority Critical patent/CN105241474B/zh
Publication of CN105241474A publication Critical patent/CN105241474A/zh
Application granted granted Critical
Publication of CN105241474B publication Critical patent/CN105241474B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Navigation (AREA)
  • Gyroscopes (AREA)

Abstract

本发明属于导航技术领域,涉及一种敏感轴斜置的惯导***的标定方法。本发明技术方案通过采用归一正交化的方法,把斜置的陀螺、加速度计敏感轴通过方向余弦矩阵转换到弹体轴上,使得虚拟的陀螺坐标系、加速度计坐标系与真实的弹体系重合,从而法解决了敏感轴斜置的惯导***的参数辨识问题。

Description

一种斜置构型惯导***标定方法
技术领域
本发明属于导航技术领域,涉及一种敏感轴斜置的惯导***的标定方法。
背景技术
捷联惯导***通常采用余度技术提高***的可靠性和精度,惯性敏感器件多采用倾斜配置。而对于非冗余的惯性导航***,其敏感器件也可以采用斜置,其优点在于可以扩大载体三个轴方向的角速度和加速度测量范围,节省惯性器件的安装空间。但是这种安装方式给惯测装置的标定带来了困难。常规的标定方法要求惯性器件沿正交的载体坐标系安装,通过转台上的位置试验、速率试验,辨识惯性器件的各项参数。
因此,亟需研制一种斜置构型惯导***标定方法,既可以在捷联惯导***或非冗余的惯性导航***中倾斜配置惯性敏感器件,又可以准确、快速辨识敏感轴斜置的惯导***的各项参数,方便对惯导***进行标定。
发明内容
本发明的目的在于提供一种斜置构型惯导***标定方法,从而准确、快速辨识敏感轴斜置的惯导***的各项参数。
为了实现这一目的,本发明采取的技术方案如下:
一种斜置构型惯导***标定方法,具体包括以下步骤:
一、对惯导***陀螺仪、加速度计进行配置
设置惯导***三个非正交的陀螺仪分别为G1、G2、G3,三个正交的弹体轴分别为OXb、OYb、OZb;陀螺仪的三个输入轴OG1、OG2、OG3均匀分布在以+OYb为中心,与+OYb轴的夹角为α的圆锥面上,相邻的陀螺仪的两个输入轴在OXbZb平面内的投影夹角为120°,其中,OG1在OXbZb平面内的投影与+OXb重合;
三个非正交的加速度计分别为A1、A2、A3,加速度计的三个输入轴OA1、OA2、OA3均匀分布在以+OYb为中心,与+OYb轴夹角为β的圆锥面上,相邻的加速度计的两个输入轴在OXbZb平面内的投影夹角为120°,其中,OA1在OXbZb平面内的投影与+OXb重合;
其中,α+β=90°;在本具体实施例中:α为54.73°,β为35.27°。
二、分别对陀螺仪、加速度计输出归一化
设定陀螺仪G1、G2、G3,加速度计A1、A2、A3输出的量纲均为LSB/s;
在标定之前,将陀螺仪的输出量纲转化为°/s,加速度计的输出量纲转化g;
陀螺仪、加速度计的输出归一化公式如下:
式中:
——陀螺仪G1、G2、G3原始输出脉冲数,单位:LSB/s;
——陀螺仪G1、G2、G3零位,单位:LSB/s;
——陀螺仪G1、G2、G3标度因数,单位:(LSB/s)/(°/s);
Ngx1、Ng2、Ng3——陀螺仪G1、G2、G3归一化后输出,单位:°/s;
——加速度计A1、A2、A3原始输出脉冲数,单位:LSB/s;
——加速度计A1、A2、A3零位,单位:LSB/s;
——加速度计A1、A2、A3标度因数,单位:(LSB/s)/g;
Na1、Na2、Na3——加速度计A1、A2、A3归一化后输出,单位:g;
归一化处理之后,陀螺仪G1、G2、G3加速度计A1、A2、A3输出仍然是非正交的;
三、陀螺仪、加速度计输出正交化
(3.1)对陀螺仪输出正交化
确定O-XbYbZb与O-G1G2G3之间的转换矩阵:
式中:
O-Ng1Ng2Ng3——步骤二中得到的陀螺仪G1、G2、G3归一化后输出;
O-NgxNgyNgz——与弹体系O-XbYbZb重合的、虚拟陀螺坐标系;
经过如上公式,非正交的陀螺坐标系O-G1G2G3通过转换矩阵Tg转化到成正交的虚拟陀螺坐标系O-NgxNgyNgz
(3.2)对加速度计输出正交化
确定O-XbYbZb与O-A1A2A3之间的转换矩阵:
式中:
O-Na1Na2Na3——步骤二中得到的加速度计A1、A2、A3归一化后输出;
O-NaxNayNaz——与弹体系O-XbYbZb重合的、虚拟的加速度计坐标系;
经过如上公式,非正交的加速度计坐标系O-A1A2A3通过转换矩阵Ta转化到成正交的虚拟加速度计坐标系O-NaxNayNaz
四、确定惯导***数学模型
归一正交化后,惯导***数学模型如下:
式中:
Ngx,Ngy,Ngz——陀螺通道在虚拟陀螺坐标系中各坐标轴上的输出;
Dfx,Dfy,Dfz——陀螺通道在虚拟陀螺坐标系中各坐标轴上的常值漂移;
Sgx,Sgy,Sgz——陀螺通道在虚拟陀螺坐标系中各坐标轴上的标度因数;
Kgij——i轴方向对j陀螺仪的安装误差系数;
Dix——X轴向线运动对X轴上陀螺输出的影响;
Diy——X轴向线运动对Y轴上陀螺输出的影响;
Diz——X轴向线运动对Z轴上陀螺输出的影响;
Dox——Y轴向线运动对X轴上陀螺输出的影响;
Doy——Y轴向线运动对Y轴上陀螺输出的影响;
Doz——Y轴向线运动对Z轴上陀螺输出的影响;
Dsx——Z轴向线运动对X轴上陀螺输出的影响;
Dsy——Z轴向线运动对Y轴上陀螺输出的影响;
Dsz——Z轴向线运动对Z轴上陀螺输出的影响;
式中:
Nax,Nay,Naz——加速度计通道在虚拟的加速度计坐标系各坐标轴上的脉冲输出;
Kax0,Kay0,Kaz0——加速度计通道在虚拟的加速度计坐标系各坐标轴上的偏值;
Kaij——i轴对j加速度计通道的安装误差系数;
Ka1x,Ka1y,Ka1z——加速度计通道在虚拟的加速度计坐标系各坐标轴上的标度因数。
进一步的,如上所述的一种斜置构型惯导***标定方法,其中:α为54.73°,β为35.27°。
本发明技术方案的有益效果在于,通过采用归一正交化的方法,把斜置的陀螺、加速度计敏感轴通过方向余弦矩阵转换到弹体轴上,使得虚拟的陀螺坐标系、加速度计坐标系与真实的弹体系重合,从而法解决了敏感轴斜置的惯导***的参数辨识问题。
具体实施方式
下面结合具体实施例对本发明技术方案进行详细说明。
一种斜置构型惯导***标定方法,具体包括以下步骤:
一、对惯导***陀螺仪、加速度计进行配置
设置惯导***三个非正交的陀螺仪分别为G1、G2、G3,三个正交的弹体轴分别为OXb、OYb、OZb;陀螺仪的三个输入轴OG1、OG2、OG3均匀分布在以+OYb为中心,与+OYb轴的夹角为α的圆锥面上,相邻的陀螺仪的两个输入轴在OXbZb平面内的投影夹角为120°,其中,OG1在OXbZb平面内的投影与+OXb重合;
三个非正交的加速度计分别为A1、A2、A3,加速度计的三个输入轴OA1、OA2、OA3均匀分布在以+OYb为中心,与+OYb轴夹角为β的圆锥面上,相邻的加速度计的两个输入轴在OXbZb平面内的投影夹角为120°,其中,OA1在OXbZb平面内的投影与+OXb重合;
其中,α+β=90°;
二、分别对陀螺仪、加速度计输出归一化
设定陀螺仪G1、G2、G3,加速度计A1、A2、A3输出的量纲均为LSB/s;
在标定之前,将陀螺仪的输出量纲转化为°/s,加速度计的输出量纲转化g;
陀螺仪、加速度计的输出归一化公式如下:
式中:
——陀螺仪G1、G2、G3原始输出脉冲数,单位:LSB/s;
——陀螺仪G1、G2、G3零位,单位:LSB/s;
——陀螺仪G1、G2、G3标度因数,单位:(LSB/s)/(°/s);
Ngx1、Ng2、Ng3——陀螺仪G1、G2、G3归一化后输出,单位:°/s;
——加速度计A1、A2、A3原始输出脉冲数,单位:LSB/s;
——加速度计A1、A2、A3零位,单位:LSB/s;
——加速度计A1、A2、A3标度因数,单位:(LSB/s)/g;
Na1、Na2、Na3——加速度计A1、A2、A3归一化后输出,单位:g;
归一化处理之后,陀螺仪G1、G2、G3加速度计A1、A2、A3输出仍然是非正交的;
三、陀螺仪、加速度计输出正交化
(3.1)对陀螺仪输出正交化
确定O-XbYbZb与O-G1G2G3之间的转换矩阵:
式中:
O-Ng1Ng2Ng3——步骤二中得到的陀螺仪G1、G2、G3归一化后输出;
O-NgxNgyNgz——与弹体系O-XbYbZb重合的、虚拟陀螺坐标系;
经过如上公式,非正交的陀螺坐标系O-G1G2G3通过转换矩阵Tg转化到成正交的虚拟陀螺坐标系O-NgxNgyNgz
(3.2)对加速度计输出正交化
确定O-XbYbZb与O-A1A2A3之间的转换矩阵:
式中:
O-Na1Na2Na3——步骤二中得到的加速度计A1、A2、A3归一化后输出;
O-NaxNayNaz——与弹体系O-XbYbZb重合的、虚拟的加速度计坐标系;
经过如上公式,非正交的加速度计坐标系O-A1A2A3通过转换矩阵Ta转化到成正交的虚拟加速度计坐标系O-NaxNayNaz
四、确定惯导***数学模型
归一正交化后,惯导***数学模型如下:
式中:
Ngx,Ngy,Ngz——陀螺通道在虚拟陀螺坐标系中各坐标轴上的输出;
Dfx,Dfy,Dfz——陀螺通道在虚拟陀螺坐标系中各坐标轴上的常值漂移;
Sgx,Sgy,Sgz——陀螺通道在虚拟陀螺坐标系中各坐标轴上的标度因数;
Kgij——i轴方向对j陀螺仪的安装误差系数;
Dix——X轴向线运动对X轴上陀螺输出的影响;
Diy——X轴向线运动对Y轴上陀螺输出的影响;
Diz——X轴向线运动对Z轴上陀螺输出的影响;
Dox——Y轴向线运动对X轴上陀螺输出的影响;
Doy——Y轴向线运动对Y轴上陀螺输出的影响;
Doz——Y轴向线运动对Z轴上陀螺输出的影响;
Dsx——Z轴向线运动对X轴上陀螺输出的影响;
Dsy——Z轴向线运动对Y轴上陀螺输出的影响;
Dsz——Z轴向线运动对Z轴上陀螺输出的影响;
式中:
Nax,Nay,Naz——加速度计通道在虚拟的加速度计坐标系各坐标轴上的脉冲输出;
Kax0,Kay0,Kaz0——加速度计通道在虚拟的加速度计坐标系各坐标轴上的偏值;
Kaij——i轴对j加速度计通道的安装误差系数;
Ka1x,Ka1y,Ka1z——加速度计通道在虚拟的加速度计坐标系各坐标轴上的标度因数。

Claims (2)

1.一种斜置构型惯导***标定方法,其特征在于,具体包括以下步骤:
(1)对惯导***陀螺仪、加速度计进行配置
设置惯导***三个非正交的陀螺仪分别为G1、G2、G3,三个正交的弹体轴分别为OXb、OYb、OZb;陀螺仪的三个输入轴OG1、OG2、OG3均匀分布在以+OYb为中心,与+OYb轴的夹角为α的圆锥面上,相邻的陀螺仪的两个输入轴在OXbZb平面内的投影夹角为120°,其中,OG1在OXbZb平面内的投影与+OXb重合;
三个非正交的加速度计分别为A1、A2、A3,加速度计的三个输入轴OA1、OA2、OA3均匀分布在以+OYb为中心,与+OYb轴夹角为β的圆锥面上,相邻的加速度计的两个输入轴在OXbZb平面内的投影夹角为120°,其中,OA1在OXbZb平面内的投影与+OXb重合;
其中,α+β=90°;
(2)分别对陀螺仪、加速度计输出归一化
设定陀螺仪G1、G2、G3,加速度计A1、A2、A3输出的量纲均为LSB/s;
在标定之前,将陀螺仪的输出量纲转化为°/s,加速度计的输出量纲转化为g;
陀螺仪、加速度计的输出归一化公式如下:
式中:
——陀螺仪G1、G2、G3原始输出脉冲数,单位:LSB/s;
——陀螺仪G1、G2、G3零位,单位:LSB/s;
——陀螺仪G1、G2、G3标度因数,单位:(LSB/s)/(°/s);
Ngx1、Ng2、Ng3——陀螺仪G1、G2、G3归一化后输出,单位:°/s;
——加速度计A1、A2、A3原始输出脉冲数,单位:LSB/s;
——加速度计A1、A2、A3零位,单位:LSB/s;
——加速度计A1、A2、A3标度因数,单位:(LSB/s)/g;
Na1、Na2、Na3——加速度计A1、A2、A3归一化后输出,单位:g;
归一化处理之后,陀螺仪G1、G2、G3加速度计A1、A2、A3输出仍然是非正交的;
(3)陀螺仪、加速度计输出正交化
(3.1)对陀螺仪输出正交化
确定O-XbYbZb与O-G1G2G3之间的转换矩阵:
式中:
Ng1Ng2Ng3——步骤(2)中得到的陀螺仪G1、G2、G3归一化后输出;
O-NgxNgyNgz——与弹体系O-XbYbZb重合的、虚拟陀螺坐标系;
经过如上公式,非正交的陀螺坐标系O-G1G2G3通过转换矩阵Tg转化到成正交的虚拟陀螺坐标系O-NgxNgyNgz
(3.2)对加速度计输出正交化
确定O-XbYbZb与O-A1A2A3之间的转换矩阵:
式中:
Na1Na2Na3——步骤(2)中得到的加速度计A1、A2、A3归一化后输出;
O-NaxNayNaz——与弹体系O-XbYbZb重合的、虚拟的加速度计坐标系;
经过如上公式,非正交的加速度计坐标系O-A1A2A3通过转换矩阵Ta转化到成正交的虚拟加速度计坐标系O-NaxNayNaz
(4)确定惯导***数学模型
归一正交化后,惯导***数学模型如下:
式中:
Ngx,Ngy,Ngz——陀螺通道在虚拟陀螺坐标系中各坐标轴上的输出;
Dfx,Dfy,Dfz——陀螺通道在虚拟陀螺坐标系中各坐标轴上的常值漂移;
Sgx,Sgy,Sgz——陀螺通道在虚拟陀螺坐标系中各坐标轴上的标度因数;
Kgij——i轴方向对j陀螺仪的安装误差系数;
Dix——X轴向线运动对X轴上陀螺输出的影响;
Diy——X轴向线运动对Y轴上陀螺输出的影响;
Diz——X轴向线运动对Z轴上陀螺输出的影响;
Dox——Y轴向线运动对X轴上陀螺输出的影响;
Doy——Y轴向线运动对Y轴上陀螺输出的影响;
Doz——Y轴向线运动对Z轴上陀螺输出的影响;
Dsx——Z轴向线运动对X轴上陀螺输出的影响;
Dsy——Z轴向线运动对Y轴上陀螺输出的影响;
Dsz——Z轴向线运动对Z轴上陀螺输出的影响;
式中:
Nax,Nay,Naz——加速度计通道在虚拟的加速度计坐标系各坐标轴上的脉冲输出;
Kaox,Kaoy,Kaoz——加速度计通道在虚拟的加速度计坐标系各坐标轴上的偏值;
Kaij——i轴对j加速度计通道的安装误差系数;
Ka1x,Ka1y,Ka1z——加速度计通道在虚拟的加速度计坐标系各坐标轴上的标度因数。
2.如权利要求1所述的一种斜置构型惯导***标定方法,其特征在于,步骤(1)中:α为54.73°,β为35.27°。
CN201410326080.XA 2014-07-10 2014-07-10 一种斜置构型惯导***标定方法 Active CN105241474B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410326080.XA CN105241474B (zh) 2014-07-10 2014-07-10 一种斜置构型惯导***标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410326080.XA CN105241474B (zh) 2014-07-10 2014-07-10 一种斜置构型惯导***标定方法

Publications (2)

Publication Number Publication Date
CN105241474A CN105241474A (zh) 2016-01-13
CN105241474B true CN105241474B (zh) 2018-06-26

Family

ID=55039203

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410326080.XA Active CN105241474B (zh) 2014-07-10 2014-07-10 一种斜置构型惯导***标定方法

Country Status (1)

Country Link
CN (1) CN105241474B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106705995A (zh) * 2016-11-23 2017-05-24 极翼机器人(上海)有限公司 一种MEMS陀螺仪g值敏感系数的标定方法
CN106767917B (zh) * 2016-12-08 2019-09-27 中国人民解放军国防科学技术大学 一种斜置冗余惯导***标定误差模型建模方法
CN107607116A (zh) * 2017-10-30 2018-01-19 北京理工大学 一种高动态惯性测量单元
CN111561948B (zh) * 2019-12-05 2023-07-28 北京计算机技术及应用研究所 四轴冗余捷联惯导的***级标定方法
CN114061620B (zh) * 2021-11-09 2024-05-10 武汉华中天易星惯科技有限公司 一种四冗余度惯导分立式标定方法及标定***

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102997932A (zh) * 2011-09-15 2013-03-27 北京自动化控制设备研究所 一种消除高精度惯导***标定中转台抖动影响的方法
CN103697918A (zh) * 2013-12-26 2014-04-02 北京航天时代光电科技有限公司 一种三轴正交一轴斜置构型光纤陀螺惯测装置的标定方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030033850A1 (en) * 2001-08-09 2003-02-20 Challoner A. Dorian Cloverleaf microgyroscope with electrostatic alignment and tuning
US8265826B2 (en) * 2003-03-20 2012-09-11 Hemisphere GPS, LLC Combined GNSS gyroscope control system and method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102997932A (zh) * 2011-09-15 2013-03-27 北京自动化控制设备研究所 一种消除高精度惯导***标定中转台抖动影响的方法
CN103697918A (zh) * 2013-12-26 2014-04-02 北京航天时代光电科技有限公司 一种三轴正交一轴斜置构型光纤陀螺惯测装置的标定方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
光学捷联惯组系数基准一致性转换方法的研究;戚红向等;《航天控制》;20110430;第29卷(第2期);第15-22页 *
捷联惯性传感器多余度配置的误差标定技术研究;华冰等;《传感器技术》;20051231;第24卷(第5期);第31-33页 *

Also Published As

Publication number Publication date
CN105241474A (zh) 2016-01-13

Similar Documents

Publication Publication Date Title
CN105241474B (zh) 一种斜置构型惯导***标定方法
CN105203129B (zh) 一种惯导装置初始对准方法
JP5043358B2 (ja) 傾斜角演算方法及び傾斜角演算装置
JP2008002992A (ja) 姿勢角検出装置と姿勢角検出方法
JP2007040763A (ja) 加速度センサの補正装置
CN103323625B (zh) 一种mems-imu中加速度计动态环境下的误差标定补偿方法
CN105910606A (zh) 一种基于角速度差值的方向修正方法
JP6739662B2 (ja) 空間的な向きの決定
KR20140025319A (ko) 근 자기장의 동적 추적 및 보상을 위한 장치 및 방법
JP7025215B2 (ja) 測位システム及び測位方法
CN103776450A (zh) 适用于高速旋转飞行体的半捷联式惯性测量与导航算法
TWI407102B (zh) 三維空間運動感測方法
Xu et al. A robust incremental-quaternion-based angle and axis estimation algorithm of a single-axis rotation using MARG sensors
CN106708088B (zh) 坐标计算方法及装置、飞行控制方法及***、无人机
JP2017102030A (ja) ピッチ角速度補正値算出装置、姿勢角算出装置およびピッチ角速度補正値算出方法
CN107576334B (zh) 惯性测量单元的标定方法及装置
RU2013136361A (ru) Интегрированная бесплатформенная система навигации средней точности для беспилотного летательного аппарата
CN103411625B (zh) 动力调谐陀螺测斜仪的标定补偿方法及标定补偿模型
JP2008241676A (ja) 電子コンパスシステム
CN116125789A (zh) 一种基于四元数的姿态算法参数自动匹配***及方法
CN105928544A (zh) 微惯性测量组合单元的快速自标定方法及装置
JP2009053039A (ja) 車両姿勢推定装置、車両姿勢推定方法
EP3669759A1 (en) System and method for dynamic posturography
Lima et al. Performance evaluation of attitude estimation algorithms in the design of an ahrs for fixed wing uavs
JP6454857B2 (ja) 姿勢検出装置及び姿勢検出方法

Legal Events

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