CN105241474B - 一种斜置构型惯导***标定方法 - Google Patents
一种斜置构型惯导***标定方法 Download PDFInfo
- 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
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°。
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)
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)
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)
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 |
-
2014
- 2014-07-10 CN CN201410326080.XA patent/CN105241474B/zh active Active
Patent Citations (2)
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)
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 |