CN104236533B - 一种陀螺仪数据融合方法 - Google Patents

一种陀螺仪数据融合方法 Download PDF

Info

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
Application number
CN201410476399.0A
Other languages
English (en)
Other versions
CN104236533A (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.)
Guangzhou Boguan Intelligent Technology Co., Ltd.
Original Assignee
GUANGZHOU BOSMA PHOTOELECTRIC TECHNOLOGY Co Ltd
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 GUANGZHOU BOSMA PHOTOELECTRIC TECHNOLOGY Co Ltd filed Critical GUANGZHOU BOSMA PHOTOELECTRIC TECHNOLOGY Co Ltd
Priority to CN201410476399.0A priority Critical patent/CN104236533B/zh
Publication of CN104236533A publication Critical patent/CN104236533A/zh
Application granted granted Critical
Publication of CN104236533B publication Critical patent/CN104236533B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C19/00Gyroscopes; Turn-sensitive devices using vibrating masses; Turn-sensitive devices without moving masses; Measuring angular rate using gyroscopic effects
    • G01C19/02Rotary 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>&amp;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。
CN201410476399.0A 2014-09-17 2014-09-17 一种陀螺仪数据融合方法 Active CN104236533B (zh)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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
CN109163721B (zh) 姿态测量方法及终端设备
CN105300379B (zh) 一种基于加速度的卡尔曼滤波姿态估计方法及系统
CN104898681B (zh) 一种采用三阶近似毕卡四元数的四旋翼飞行器姿态获取方法
CN104236533B (zh) 一种陀螺仪数据融合方法
CN103940425B (zh) 一种磁‑惯性组合捷联测量方法
CN105157705B (zh) 一种半捷联雷达导引头视线角速度提取方法
CN108592943B (zh) 一种基于opreq方法的惯性系粗对准计算方法
CN106979780A (zh) 一种无人车实时姿态测量方法
CN109724602A (zh) 一种基于硬件fpu的姿态解算系统及其解算方法
CN110044385B (zh) 一种大失准角情况下的快速传递对准方法
CN108225370A (zh) 一种运动姿态传感器的数据融合与解算方法
CN106370178B (zh) 移动终端设备的姿态测量方法及装置
CN106767798B (zh) 一种无人机导航用的位置及速度的实时估计方法及系统
CN110986928B (zh) 光电吊舱三轴陀螺仪漂移实时修正方法
CN116817896B (zh) 一种基于扩展卡尔曼滤波的姿态解算方法
CN103712598A (zh) 一种小型无人机姿态确定系统与确定方法
CN107063254A (zh) 一种陀螺地磁组合的姿态解算方法
CN108344413A (zh) 一种水下滑翔器导航系统及其低精度与高精度转换方法
CN109556606A (zh) 基于模糊补偿和卡尔曼滤波的传感器角度测量优化方法
CN112798021A (zh) 基于激光多普勒测速仪的惯导系统行进间初始对准方法
CN111189474A (zh) 基于mems的marg传感器的自主校准方法
CN111307114B (zh) 基于运动参考单元的水面舰船水平姿态测量方法
CN104121930B (zh) 一种基于加表耦合的mems陀螺漂移误差的补偿方法
CN105277210B (zh) 一种任意安装多轴陀螺组合标定方法
CN109674480B (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