CN104236533B - 一种陀螺仪数据融合方法 - Google Patents
一种陀螺仪数据融合方法 Download PDFInfo
- Publication number
- CN104236533B CN104236533B CN201410476399.0A CN201410476399A CN104236533B CN 104236533 B CN104236533 B CN 104236533B CN 201410476399 A CN201410476399 A CN 201410476399A CN 104236533 B CN104236533 B CN 104236533B
- Authority
- CN
- China
- Prior art keywords
- angle
- data fusion
- accelerometer
- modulus value
- attitude angle
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C19/00—Gyroscopes; Turn-sensitive devices using vibrating masses; Turn-sensitive devices without moving masses; Measuring angular rate using gyroscopic effects
- G01C19/02—Rotary gyroscopes
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Gyroscopes (AREA)
Abstract
本发明公开了一种陀螺仪数据融合方法,包括以下步骤:A、在陀螺仪设备开启时,利用静止状态下加速度计输出的三个轴矢量,计算三个轴矢量和的模值;B、在运动状态下加速度计输出的三个轴矢量,计算三个轴矢量和的模值;C、利用上一时刻的姿态角、根据加速度计的输出得到的姿态角、根据角速度计得到的特定轴的角速度、上一时刻与当前时刻的间隔时间以及数据融合参考系数,计算当前的姿态角。本发明通过在数据融合过程中引入数据融合参考系数,在不增加计算复杂度的情况下实现高精度数据融合,同时解决了现有技术在设备存在加速度的情况下不能进行姿态定位的问题。本发明作为一种陀螺仪数据融合方法可广泛应用于陀螺仪技术领域。
Description
技术领域
本发明涉及陀螺仪技术领域,尤其是一种陀螺仪数据融合方法。
背景技术
现有设备中采用陀螺仪来做智能控制,大多设备采用双轴或三轴陀螺仪来控制设备的姿态角,由于加速度计在静止时数据是精确的,但是当设备运动起来后,数据就不再精确了,角加速度计计算物体的旋转运动,至此只用一轴的数据,无论是单独利用加速度计还是单独利用角加速度计都无法得到设备的角度,而且由于陀螺仪本身的精度和制造工艺等问题,使得由硬件直接输出的数据具有漂移和累积误差,还有数据耦合过程中也会存在误差,因此输出的角度数据无法长时间保持高精度,因此利用陀螺仪控制设备的技术难点在于陀螺仪数据的融合运算。
当前主要有以下几种陀螺仪数据融合的方法,基于四元数融合法,基于卡尔曼算法的数据融合,基于互补滤波的方式进行数据融合。虽然基于卡尔曼算法的数据融合的精度较高,但是计算量巨大,过程复杂,基于四元数法计算量比卡尔曼算法小一些,互补滤波计算量最小但是精度要低一些,而且目前上面的这些融合算法都是在设备没有做加速运动的情况,在加速运动的情况下,加速度计的数据量融合了重力还有设备本身的加速度,无法再利用它进行姿态定位。
发明内容
为了解决上述技术问题,本发明的目的是:提供一种适用于设备存在加速度的陀螺仪数据融合方法。
本发明所采用的技术方案是:一种陀螺仪数据融合方法,包括有以下步骤:
A、在陀螺仪设备开启时,利用静止状态下加速度计输出的三个轴矢量,计算三个轴矢量和的模值;
B、在运动状态下加速度计输出的三个轴矢量,计算三个轴矢量和的模值;
C、利用上一时刻的姿态角、根据加速度计的输出得到的姿态角、根据角速度计得到的特定轴的角速度、上一时刻与当前时刻的间隔时间以及数据融合参考系数,计算当前的姿态角。
进一步,所述步骤A中对三个轴矢量进行零值偏移误差补偿后再计算三个轴矢量和的模值。
进一步,所述数据融合参考系数根据上述步骤A、步骤B中的模值计算所得。
进一步,所述步骤C中当前姿态角的计算公式为:
其中anglet-1为上一时刻的姿态角,angle_accelerate为根据加速度计的输出得到的姿态角,angle_gyro_rate为根据角速度计得到的特定轴的角速度,T为上一时刻与当前时刻的间隔时间,K为数据融合参考系数,anglet为当前的姿态角。
进一步,所述数据融合参考系数K的计算公式为:
其中G为所述步骤A中的模值, 为所述步骤B中的模值,C和D为常数因子。
进一步,对于不同取值范围的,设置不同值的常数因子C和常数因子D。
本发明的有益效果是:本发明通过在数据融合过程中引入数据融合参考系数,在不增加计算复杂度的情况下实现高精度数据融合,同时解决了现有技术在设备存在加速度的情况下不能进行姿态定位的问题。
附图说明
图1为本发明方法的步骤流程图;
图2为某时刻陀螺仪状态在空间直角坐标系中的示意图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步说明:
参照图1,一种陀螺仪数据融合方法,包括有以下步骤:
A、在陀螺仪设备开启时,利用静止状态下加速度计输出的三个轴矢量,计算三个轴矢量和的模值;
加速度计输出值(A_x,A_y,A_z)为在设备坐标系下重力加速度在三个轴的矢量,加速度计三个轴的矢量和的模值计算公式如下:
G即重力加速度值,不同的纬度的值会有一定的误差,一般在不精确的情况下取9.8m/s2,用来做参照值。
B、在运动状态下加速度计输出的三个轴矢量,计算三个轴矢量和的模值;
当陀螺仪处在稳态时加速度值为0(静止,或者只有旋转,或者匀速运动时),三轴的和矢量的模值应该是G ,当设备做加速或减速运动时,那么瞬时加速度值G’应该大于或小于G的值,物体的加速度越大,那么G与G’的相差也就越大。
C、利用上一时刻的姿态角、根据加速度计的输出得到的姿态角、根据角速度计得到的特定轴的角速度、上一时刻与当前时刻的间隔时间以及数据融合参考系数,计算当前的姿态角。
进一步作为优选的实施方式,所述步骤A中对三个轴矢量进行零值偏移误差补偿后再计算三个轴矢量和的模值。
进一步作为优选的实施方式,所述数据融合参考系数根据上述步骤A、步骤B中的模值计算所得。
进一步作为优选的实施方式,所述步骤C中当前姿态角的计算公式为:
其中anglet-1为上一时刻的姿态角,angle_accelerate为根据加速度计的输出得到的姿态角,参照图2,具体得到姿态角的公式在下面的公式中:
tan(Axz) = Rx/Rz => Angle_y = atan2(Rx,Rz);
tan(Ayz)=Ry/Rz => Angle_x=atan2(Ry,Rz);
tan(Axy)=Ry/Rx => Angle_z=atan2(Ry,Rx);
angle_gyro_rate为根据角速度计得到的特定轴的角速度,T为上一时刻与当前时刻的间隔时间,K为数据融合参考系数,anglet为当前的姿态角,anglet的值为特定轴所转过的角度,如图2中所示的Axz角、Ayz角。
进一步,所述数据融合参考系数K的计算公式为:
其中G为所述步骤A中的模值, 为所述步骤B中的模值,C和D为常数因子。
进一步,对于不同取值范围的,设置不同值的常数因子C和常数因子D;即根据实践中的具体情况设置几个动态范围的阈值,比如当在[0,0.2]、[0.2,0.4]、[0.4,0.5]设置不同的常数因子C和常数因子D从而得到不同的K值,这样既可以提高运算的速度又可以尽量减小运算量。
以上是对本发明的较佳实施进行了具体说明,但本发明创造并不限于所述实施例,熟悉本领域的技术人员在不违背本发明精神的前提下还可以作出种种的等同变换或替换,这些等同的变形或替换均包含在本申请权利要求所限定的范围内。
Claims (3)
1.一种陀螺仪数据融合方法,其特征在于:包括有以下步骤:
A、在陀螺仪设备开启时,利用静止状态下加速度计输出的三个轴矢量,计算三个轴矢量和的模值;
B、在运动状态下加速度计输出的三个轴矢量,计算三个轴矢量和的模值;
C、利用上一时刻的姿态角、根据加速度计的输出得到的姿态角、根据角速度计得到的特定轴的角速度、上一时刻与当前时刻的间隔时间以及数据融合参考系数,计算当前的姿态角,当前姿态角的计算公式为:
anglet=anglet-1+K(angle_accelerate-anglet-1)+(1-K)(angle_gyro_rate)T其中anglet-1为上一时刻的姿态角,angle_accelerate为根据加速度计的输出得到的姿态角,angle_gyro_rate为根据角速度计得到的特定轴的角速度,T为上一时刻与当前时刻的间隔时间,K为数据融合参考系数,anglet为当前的姿态角;
所述数据融合参考系数根据上述步骤A、步骤B中的模值计算所得,数据融合参考系数K的计算公式为:
<mrow>
<mi>K</mi>
<mo>=</mo>
<mfrac>
<mi>C</mi>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mi>D</mi>
<mo>|</mo>
<msup>
<mi>G</mi>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<mi>G</mi>
<mo>|</mo>
</mrow>
</mfrac>
</mrow>
其中G为所述步骤A中的模值,G'为所述步骤B中的模值,C和D为常数因子。
2.根据权利要求1所述的一种陀螺仪数据融合方法,其特征在于:所述步骤A中对三个轴矢量进行零值偏移误差补偿后再计算三个轴矢量和的模值。
3.根据权利要求1所述的一种陀螺仪数据融合方法,其特征在于:对于不同取值范围的|G'-G|,设置不同值的常数因子C和常数因子D。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410476399.0A CN104236533B (zh) | 2014-09-17 | 2014-09-17 | 一种陀螺仪数据融合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410476399.0A CN104236533B (zh) | 2014-09-17 | 2014-09-17 | 一种陀螺仪数据融合方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104236533A CN104236533A (zh) | 2014-12-24 |
CN104236533B true CN104236533B (zh) | 2017-11-14 |
Family
ID=52225080
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410476399.0A Active CN104236533B (zh) | 2014-09-17 | 2014-09-17 | 一种陀螺仪数据融合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104236533B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107203257A (zh) * | 2016-03-17 | 2017-09-26 | 深圳多哚新技术有限责任公司 | 一种头部姿态补偿方法及相关设备 |
CN106028107B (zh) * | 2016-06-20 | 2018-12-14 | 昆明理工大学 | 一种基于加速度计和陀螺仪传感器的无线英文文本输入系统 |
CN108663044B (zh) * | 2017-03-31 | 2021-03-16 | 阿里巴巴(中国)有限公司 | 一种定位方法及装置 |
CN114440926A (zh) * | 2022-01-18 | 2022-05-06 | 武汉元生创新科技有限公司 | 一种陀螺仪检测量程扩展方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101057196A (zh) * | 2004-10-05 | 2007-10-17 | 迈克罗茵费尼蒂公司 | 控制动态装置的方法和装置 |
CN101319902A (zh) * | 2008-07-18 | 2008-12-10 | 哈尔滨工程大学 | 一种低成本组合式定位定向装置及组合定位方法 |
CN101782391A (zh) * | 2009-06-22 | 2010-07-21 | 北京航空航天大学 | 机动加速度辅助的扩展卡尔曼滤波航姿系统姿态估计方法 |
CN102289306A (zh) * | 2011-08-30 | 2011-12-21 | 江苏惠通集团有限责任公司 | 姿态感知设备及其定位、鼠标指针的控制方法和装置 |
CN202974288U (zh) * | 2012-12-05 | 2013-06-05 | 南京理工大学 | 一种微型捷联航姿系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130226505A1 (en) * | 2012-02-02 | 2013-08-29 | mCube, Incorporated | Dual Accelerometer Plus Magnetometer Body Rotation Rate Sensor-Gyrometer |
-
2014
- 2014-09-17 CN CN201410476399.0A patent/CN104236533B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101057196A (zh) * | 2004-10-05 | 2007-10-17 | 迈克罗茵费尼蒂公司 | 控制动态装置的方法和装置 |
CN101319902A (zh) * | 2008-07-18 | 2008-12-10 | 哈尔滨工程大学 | 一种低成本组合式定位定向装置及组合定位方法 |
CN101782391A (zh) * | 2009-06-22 | 2010-07-21 | 北京航空航天大学 | 机动加速度辅助的扩展卡尔曼滤波航姿系统姿态估计方法 |
CN102289306A (zh) * | 2011-08-30 | 2011-12-21 | 江苏惠通集团有限责任公司 | 姿态感知设备及其定位、鼠标指针的控制方法和装置 |
CN202974288U (zh) * | 2012-12-05 | 2013-06-05 | 南京理工大学 | 一种微型捷联航姿系统 |
Non-Patent Citations (6)
Title |
---|
一种MAV航姿估计算法及其半实物仿真;赵世峰等;《系统仿真学报》;20070930;第19卷(第18期);正文引言部分、第1-3节及图1、2 * |
动态倾角传感器及其传递特性的研究;付勇杰等;《仪表技术与传感器》;20121231(第9期);全文 * |
基于MEMS加速度计和陀螺仪的姿态检测系统;赵翔等;《铁路计算机应用》;20120331;第21卷(第3期);全文 * |
基于陀螺仪及加速度计信号融合的姿态角度测量;冯智勇等;《西南师范大学学报(自然科学版)》;20110831;第36卷(第4期);全文 * |
基于陀螺仪和加速度计的帆船运动姿态测量;苏菲等;《太赫兹科学与电子信息学报》;20140430;第12卷(第2期);全文 * |
多传感器融合在倾角测量抗振动中的应用;宋振等;《传感器世界》;20140531;第20卷(第5期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN104236533A (zh) | 2014-12-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105203129B (zh) | 一种惯导装置初始对准方法 | |
WO2020087846A1 (zh) | 基于迭代扩展卡尔曼滤波融合惯性与单目视觉的导航方法 | |
CN105953796A (zh) | 智能手机单目和imu融合的稳定运动跟踪方法和装置 | |
CN104236533B (zh) | 一种陀螺仪数据融合方法 | |
WO2021082790A1 (zh) | 一种基于IMU的uwb定位异常值处理方法 | |
CN103940425B (zh) | 一种磁‑惯性组合捷联测量方法 | |
CN106482734A (zh) | 一种用于imu多传感器数据融合的滤波方法 | |
CN110174899A (zh) | 一种基于敏捷卫星的高精度成像姿态指向控制方法 | |
CN104737205B (zh) | 自适应尺度及/或重力估计 | |
CN109724602A (zh) | 一种基于硬件fpu的姿态解算系统及其解算方法 | |
CN107063262A (zh) | 一种用于无人机姿态解算的互补滤波方法 | |
CN109682377A (zh) | 一种基于动态步长梯度下降的姿态估计方法 | |
CN110044385B (zh) | 一种大失准角情况下的快速传递对准方法 | |
CN106767798B (zh) | 一种无人机导航用的位置及速度的实时估计方法及系统 | |
CN110986928B (zh) | 光电吊舱三轴陀螺仪漂移实时修正方法 | |
CN107063254A (zh) | 一种陀螺地磁组合的姿态解算方法 | |
CN113670314B (zh) | 基于pi自适应两级卡尔曼滤波的无人机姿态估计方法 | |
CN116817896B (zh) | 一种基于扩展卡尔曼滤波的姿态解算方法 | |
CN103712598A (zh) | 一种小型无人机姿态确定系统与确定方法 | |
CN112798021A (zh) | 基于激光多普勒测速仪的惯导系统行进间初始对准方法 | |
CN111189474A (zh) | 基于mems的marg传感器的自主校准方法 | |
CN111307114B (zh) | 基于运动参考单元的水面舰船水平姿态测量方法 | |
CN108827288A (zh) | 一种基于对偶四元数的降维捷联惯性导航系统初始对准方法及系统 | |
CN105277210B (zh) | 一种任意安装多轴陀螺组合标定方法 | |
CN116125789A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CP03 | Change of name, title or address |
Address after: 510530 Second Floor, Building A5, 11 Kaiyuan Avenue, Science City, Guangzhou High-tech Industrial Development Zone, Guangdong Province Patentee after: Guangzhou Boguan Intelligent Technology Co., Ltd. Address before: 3 building, A5 building, No. 11, Kaiyuan Avenue, Guangzhou Science City, Guangdong, 510530 Patentee before: Guangzhou Bosma Photoelectric Technology Co., Ltd. |
|
CP03 | Change of name, title or address |