CN103323625A - 一种mems-imu中加速度计动态环境下的误差标定补偿方法 - Google Patents
一种mems-imu中加速度计动态环境下的误差标定补偿方法 Download PDFInfo
- Publication number
- CN103323625A CN103323625A CN2013102332878A CN201310233287A CN103323625A CN 103323625 A CN103323625 A CN 103323625A CN 2013102332878 A CN2013102332878 A CN 2013102332878A CN 201310233287 A CN201310233287 A CN 201310233287A CN 103323625 A CN103323625 A CN 103323625A
- Authority
- CN
- China
- Prior art keywords
- accelerometer
- turntable
- inside casing
- housing
- output
- 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
Links
Images
Abstract
本发明给出了一种MEMS-IMU中三轴加速度计动态环境下的误差标定补偿方法,其包括:使转台外框按速率模式运行;其次让Z轴回到零位,转台内框按照正弦模式运行;改变IMU的安装,分别让X轴、Y轴加速度计的敏感轴方向与天向平行,和Z轴相同的运行模式让加速度计敏感动态加速度和动态角速度;记录下三个轴向加速度计的输出数据,进行分析处理,可以得到MEMS-IMU中三轴加速度计的安装误差,标定系数,零点漂移以及动态角速度对加速度计的影响因子。本方法简单易实现,对MEMS-IMU系统精度提高明显,特别适用于应用低成本、低精度MEMS加速度计的系统中。
Description
技术领域
本发明给出了一种MEMS-IMU中加速度计动态环境下的误差标定补偿方法,可用于MEMS-IMU中三轴加速度计的误差标定和补偿,特别适用于低成本微机电加速度计构成的惯性测量单元。
背景技术
MEMS(Micro-Electro-MechanicalSystems,MEMS)微惯性测量单元一般由陀螺仪和加速度计构成,由于其体积小、质量轻、成本低等优点越来越多的被应用到惯性导航的各个领域,给惯性测量带来全新的变革和发展机遇。
微惯性测量单元是整个惯性测量系统的核心传感器,决定着整个系统的性能和测量精度。加速度计作为两个核心元件之一,由于其精度不是很高,所以对其进行主要误差的分析和标定就显得尤为重要。
通常情况下,标定加速度计的安装误差和比例因子方法都是在静态或者准静态条件下进行,这样得到的误差补偿模型对于加速度计测量动态加速度时的误差补偿结果不理想。
发明内容
本发明提出的目的是通过双轴转台模拟的动态环境下,标定MEMS-IMU中三轴加速度计主要误差,并给出误差补偿模型。
实现本发明目的的解决方案及步骤如下:
将MEMS-IMU安装在双轴转台的基座上,Z加速度计、X加速度计敏感轴分别指向转台的外框和内框,转台外框指天向,内框指北向,则Y加速度计敏感轴自动指西向。
待IMU(Inertial Measurement Unit,IMU)预热一段时间之后,固定转台内框,将转台外框设定为速率模式,速率为ω,此时Z加速度计敏感动态加速度a=g·cos(ωt),其中g为当地地球重力加速度、t为转台运行时间。同步采集Y加速度计输出aayz和Z加速度计输出aazz。
转台内框转动90度后保持静止,将转台外框设定为速率模式,速率同为ω,此时Z加速度计敏感动态加速度a=g·cos(ωt),其中g为当地地球重力加速度、t为转台运行时间,同步采集X加速度计输出aaxz。
转台外框指天向且保持静止,将转台内框设定同为幅值为A,频率为f的正弦速率模式,同步采集转台内框角速度ωz、X加速度计输出aωxz、Y加速度计输出aωyz和Z加速度计输出aωzz。
调整IMU三个加速计的指向,使Y加速度计指向外框,X加速度计指向内框,Z加速度计指向东向,固定转台内框,将转台外框设定为速率模式,速率同为ω,此时Y加速度计敏感动态加速度a=g·cos(ωt),其中g为当地地球重力加速度、t为转台运行时间,同步采集Y加速度计输出aayy和Z加速度计输出aazy。
转台内框转动90度后保持静止,将转台外框设定为速率模式,速率同为ω,此时Y加速度计敏感动态加速度a=g·cos(ωt),其中g为当地地球重力加速度、t为转台运行时间,同步采集X加速度计输出aaxy。
转台外框指天向且保持静止,将转台内框设定同为幅值为A,频率为f的正弦速率模式,同步采集转台内框角速度ωy、X加速度计输出aωxy、Y加速度计输出aωyy和Z加速度计输出aωzy。
调整IMU三个加速计的指向,使X加速度计指向外框,Y加速度计指向内框,Z加速度计指向东向.固定转台内框,将转台外框设定为速率模式,速率同为ω,此时X加速度计敏感动态加速度a=g·cos(ωt),其中g为当地地球重力加速度、t为转台运行时间,同步采集X加速度计输出aaxx和Z加速度计输出aazx。
转台内框转动90度后保持静止,将转台外框设定为速率模式,速率同为ω,此时X加速度计敏感动态加速度a=g·cos(ωt),其中g为当地地球重力加速度、t为转台运行时间,同步采集Y加速度计输出aayx。
转台外框指天向且保持静止,将转台内框设定同为幅值为A,频率为f的正弦速率模式,同步采集转台内框角速度ωx、X加速度计输出aωxx、Y加速度计输出aωyx和Z加速度计输出aωzx。
本发明的创新点在于:MEMS-IMU中三轴加速度计敏感双轴转台模拟的动态环境条件下加速度和角速度,通过加速度计实际输出,分析加速度计主要误差受到的动态加速度和角速度的影响,从而求出加速度计动态误差补偿模型,实现误差的补偿。
附图说明
图1为MEMS-IMU中三轴加速度计安装示意图;
图2为与本发明具体实施方案一致的方法流程图。
具体实施方案
结合附图介绍本发明详细实施方案和步骤:
将如图1所示三轴加速度计安装在双轴转台的基座上,Z加速度计、X加速度计敏感轴分别指向转台的外框和内框,转台外框指天向,内框指北向,则Y加速度计敏感轴自动指西向,如图2-(a)所示。
待IMU预热一段时间之后,固定转台内框,将转台外框设定为速率模式运行,速率为ω,此时Z加速度计敏感动态加速度a=g·cos(ωt),其中g为当地地球重力加速度、t为转台运行时间,同步采集Y加速度计输出aayz和Z加速度计输出aazz,由采集到的Z加速度计输出数据拟合得到曲线aazz=hz·cos(ωt)+azz0,其中hz、azz0分别为由Z加速度计采集数据计算得到正弦曲线幅值和零点漂移。通过分析采集到的数据aayz、aazz和a之间的关系,拟合得到Czz=g/hz,aayz=Cyz·a+ayz0,其中Czz、Cyz分别为Z轴方向有动态加速度输入时Z加速度计、Y加速度计的输出误差系数,ayz0为此时计算得到的Y加速度计零点漂移。
转台内框顺时针转动90度后保持静止,如图2-(b)所示。将定转台外框设定为速率模式,速率同为ω,此时Z加速度计敏感动态加速度a=g·cos(ωt),其中g为当地地球重力加速度、t为转台运行时间,同步采集X加速度计输出aaxz。通过分析采集到的数据aaxz和a之间的关系,拟合得到aaxz=Cxz·a+axz0,其中Cxz为Z轴方向有动态加速度输入时X加速度计的输出误差系数,axz0为此时计算得到的X加速度计零点漂移。
转台外框指天向且保持静止,将转台内框设定为幅值为A,频率为f的正弦速率模式,如图2-(c)所示,同步采集转台内框角速度ωz、X加速度计输出aωxz、Y加速度计输出aωyz和Z加速度计输出aωzz;并分析ωz别和aωxz、aωyz、aωzz之间的关系,拟合出关系式aωxz=Txz·ωz+aωxz0、aωyz=Tyz·ωz+aωyz0、aωzz=Tzz·ωz+aωzz0,其中Txz、Tyz、Tzz分别是Z轴加速度计方向上有动态角速度输入时X、Y、Z三个加速度计输出与动态角速度的比例因子,aωxz0、aωyz0、aωzz0分别表示此时三个加速度计的零点漂移。
调整IMU三个加速计的指向,使Y加速度计指向外框,X加速度计指向内框,Z加速度计指向东向,如图2-(d)所示。固定转台内框,将转台外框设定为速率模式,速率同为ω,此时Y加速度计敏感动态加速度a=g·cos(ωt),其中g为当地地球重力加速度、t为转台运行时间,同步采集Y加速度计输出aayy和Z加速度计输出aazy,由采集到的Y加速度计输出数据拟合得到曲线aayy=hy·cos(ωt)+ayy0,其中hy、ayy0分别为由Y加速度计采集数据计算得到正弦曲线幅值和零点漂移。通过分析采集到的数据aayy、aazy和a之间的关系,拟合得到Cyy=g/hy,aazy=Czy·a+azy0,其中Cyy、Czy分别为Y轴方向有动态加速度输入时Y加速度计、Z加速度计的输出误差系数,azy0为此时计算得到的Z加速度计零点漂移。
转台内框顺时针转动90度后保持静止,如图2-(e)所示。将转台外框设定为速率模式运行,速率同为ω,此时Y加速度计敏感动态加速度a=g·cos(ωt),其中g为当地地球重力加速度、t为转台运行时间,同步采集X加速度计输出aaxy。通过分析采集到的数据aaxy和a之间的关系,拟合得到aaxy=Cxy·a+axy0,其中Cxy为Y轴方向有动态加速度输入时X加速度计的输出误差系数,axy0为此时计算得到的X加速度计零点漂移。
转台外框指天向且保持静止,将转台内框设定同为幅值为A,频率为f的正弦速率模式,如图2-(f)所示。同步采集转台内框角速度ωy、X加速度计输出aωxy、Y加速度计输出aωyy和Z加速度计输出aωzy;并分析ωy分别和aωxy、aωyy、aωzy之间的关系,拟合出关系式aωxy=Txy·ωy+aωxy0、aωyy=Tyy·ωy+aωyy0、aωzy=Tzy·ωy+aωzy0,其中Txy、Tyy、Tzy分别是Y轴加速度计方向上有动态角速度输入时X、Y、Z三个加速度计输出与动态角速度的比例因子,aωxy0、aωyy0、aωzy0分别表示此时三个加速度计的零点漂移。
调整IMU三个加速计的指向,使X加速度计指向外框,Y加速度计指向内框,Z加速度计指向东向,如图2-(g)所示。固定转台内框,将转台外框设定为速率模式,速率同为ω,此时X加速度计敏感动态加速度a=g·cos(ωt),其中g为当地地球重力加速度、t为转台运行时间,同步采集X加速度计输出aaxx和Z加速度计输出aazx,由采集到的X加速度计输出数据拟合得到曲线aaxx=hx·cos(ωt)+axx0,其中hx、axx0分别为由X加速度计采集数据计算得到正弦曲线幅值和零点漂移。通过分析采集到的数据aaxx、aazx和a之间的关系,拟合得到Cxx=g/hx,aazx=Czx·a+azx0,其中Cxx、Czx分别为X轴方向有动态加速度输入时X加速度计、Z加速度计的输出误差系数,azx0为此时计算得到的Z加速度计零点漂移。
转台内框顺时针转动90度后保持静止,如图2-(h)所示。将转台外框设定为速率模式运行,速率同为ω,此时X加速度计敏感动态加速度a=g·cos(ωt),其中g为当地地球重力加速度、t为转台运行时间,同步采集Y加速度计输出aayx。通过分析采集到的数据aayx和a之间的关系,拟合得到aayx=Cyx·a+ayx0,其中Cyx为X轴方向有动态加速度输入时Y加速度计的输出误差系数,ayx0为此时计算得到的Y加速度计零点漂移。
转台外框指天向且保持静止,将转台内框设定同为幅值为A,频率为f的正弦速率模式,如图2-(i)所示。同步采集转台内框角速度ωx、X加速度计输出aωxx、Y加速度计输出aωyx和Z加速度计输出aωzx;并分析ωx分别和aωxx、aωyx、aωzx之间的关系,拟合出关系式aωxx=Txx·ωx+awxx0、aωyx=Tyx·ωx+awyx0、aωzx=Tzz·ωx+awzx0,其中Txx、Tyx、Tzx分别是X轴加速度计方向上有动态角速度输入时X、Y、Z三个加速度计输出与动态角速度的比例因子,aωxy0、aωyy0、aωzy0分别表示此时三个加速度计的零点漂移。
结合以上具体实施方案确定出来的各个误差参数和比例因子,确定出三轴加速度计的误差补偿模型如式(1)所示:
其中
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (1)
1.一种基于微电子机械系统MEMS惯性测量单元IMU中加速度计动态环境下的误差标定补偿方法,其特征在于包括以下步骤:
(1)将三轴加速度计安装在双轴转台的基座上,Z加速度计和X加速度计的敏感轴分别指向转台的外框和内框,转台外框指天向,内框指北向,Y加速度计的敏感轴自动指西向;
(2)待IMU预热一段时间之后,固定转台内框,将转台外框设定为速率模式,速率为ω,此时Z加速度计敏感动态加速度a=g·cos(ωt),其中g为当地地球重力加速度、t为转台运行时间,同步采集Y加速度计输出aayz和Z加速度计输出aazz;
(3)转台内框转动90度后保持静止,将转台外框设定为速率模式,角速率同为ω,此时Z加速度计敏感动态加速度a=g·cos(ωt),其中g为当地地球重力加速度、t为转台运行时间,同步采集X加速度计输出aaxz;
(4)转台外框指天向且保持静止,将转台内框设定为幅值为A,频率为f的正弦速率模式,同步采集转台内框角速率ωz、X加速度计输出aωxz、Y加速度计输出aωyz和Z加速度计输出aωzz;
(5)调整IMU三个加速计的指向,使Y加速度计指向外框,X加速度计指向内框,Z加速度计指向东向,固定转台内框,将转台外框设定为速率模式,角速率同为ω,此时Y加速度计敏感动态加速度a=g·cos(ωt),其中g为当地地球重力加速度、t为转台运行时间;同步采集Y加速度计输出aayy和Z加速度计输出aazy;
(6)转台内框转动90度后保持静止,将转台外框设定为速率模式,角速率同为ω,此时Y加速度计敏感动态加速度a=g·cos(ωt),其中g为当地地球重力加速度、t为转台运行时间,同步采集X加速度计输出aaxy;
(7)转台外框指天向且保持静止,将转台内框设定同为幅值为A,频率为f的正弦速率模式,同步采集转台内框角速率ωy、X加速度计输出aωxy、Y加速度计输出aωyy和Z加速度计输出aωzy;
(8)调整IMU三个加速计的指向,使X加速度计指向外框,Y加速度计指向内框,Z加速度计指向东向,固定转台内框,将转台外框设定为速率模式,角速率同为ω,此时X加速度计敏感动态加速度a=g·cos(ωt),其中g为当地地球重力加速度、t为转台运行时间,同步采集X加速度计输出aaxx和Z加速度计输出aazx;
(9)转台内框转动90度后保持静止,将转台外框设定为速率模式,角速率同为ω,此时X加速度计敏感动态加速度a=g·cos(ωt),其中g为当地地球重力加速度、t为转台运行时间,同步采集Y加速度计输出aayx;
(10)转台外框指天向且保持静止,将转台内框设定同为幅值为A,频率为f的正弦速率模式,同步采集转台内框角速率ωx、X加速度计输出aωxx、Y加速度计输出aωyx和Z加速度计输出aωzx。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310233287.8A CN103323625B (zh) | 2013-06-13 | 2013-06-13 | 一种mems-imu中加速度计动态环境下的误差标定补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310233287.8A CN103323625B (zh) | 2013-06-13 | 2013-06-13 | 一种mems-imu中加速度计动态环境下的误差标定补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103323625A true CN103323625A (zh) | 2013-09-25 |
CN103323625B CN103323625B (zh) | 2014-10-15 |
Family
ID=49192494
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310233287.8A Active CN103323625B (zh) | 2013-06-13 | 2013-06-13 | 一种mems-imu中加速度计动态环境下的误差标定补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103323625B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103675352A (zh) * | 2013-12-19 | 2014-03-26 | 中北大学 | 一种弹用捷联三轴加速度计组合静动态参数综合标定方法 |
CN103884870A (zh) * | 2014-03-13 | 2014-06-25 | 工业和信息化部电子第五研究所 | 提高加速度计校准精度的方法和装置 |
CN103983274A (zh) * | 2014-04-11 | 2014-08-13 | 湖北航天技术研究院总体设计所 | 一种适用于低精度无方位基准双轴转位设备的惯性测量单元标定方法 |
CN105842481A (zh) * | 2016-03-18 | 2016-08-10 | 重庆邮电大学 | 一种惯性导航系统中基于十字叉的三轴加速度计标定方法 |
CN103994775B (zh) * | 2014-04-11 | 2017-01-04 | 湖北航天技术研究院总体设计所 | 一种适用于低精度有方位基准双轴转位设备的惯性测量单元标定方法 |
CN106840100A (zh) * | 2017-03-13 | 2017-06-13 | 蒋海涛 | 一种数字式倾角传感器及测量方法 |
CN107505662A (zh) * | 2017-09-25 | 2017-12-22 | 华中科技大学 | 一种三轴mems重力仪 |
CN109613183A (zh) * | 2019-01-03 | 2019-04-12 | 北京英视睿达科技有限公司 | 低精度监测设备数据的修定方法及系统 |
CN110530395A (zh) * | 2019-05-31 | 2019-12-03 | 北京航天时代光电科技有限公司 | 小量程加速度计标定方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101059384A (zh) * | 2007-05-18 | 2007-10-24 | 南京航空航天大学 | 一种捷联mems惯性测量单元及安装误差标定方法 |
CN101290326A (zh) * | 2008-06-13 | 2008-10-22 | 哈尔滨工程大学 | 石英挠性加速度计测量组件的参数辨识标定方法 |
CN101598545A (zh) * | 2009-07-01 | 2009-12-09 | 哈尔滨工程大学 | 捷联系统加速度计安装误差角直接测量方法 |
-
2013
- 2013-06-13 CN CN201310233287.8A patent/CN103323625B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101059384A (zh) * | 2007-05-18 | 2007-10-24 | 南京航空航天大学 | 一种捷联mems惯性测量单元及安装误差标定方法 |
CN101290326A (zh) * | 2008-06-13 | 2008-10-22 | 哈尔滨工程大学 | 石英挠性加速度计测量组件的参数辨识标定方法 |
CN101598545A (zh) * | 2009-07-01 | 2009-12-09 | 哈尔滨工程大学 | 捷联系统加速度计安装误差角直接测量方法 |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103675352A (zh) * | 2013-12-19 | 2014-03-26 | 中北大学 | 一种弹用捷联三轴加速度计组合静动态参数综合标定方法 |
CN103884870A (zh) * | 2014-03-13 | 2014-06-25 | 工业和信息化部电子第五研究所 | 提高加速度计校准精度的方法和装置 |
CN103884870B (zh) * | 2014-03-13 | 2016-08-24 | 工业和信息化部电子第五研究所 | 提高加速度计校准精度的方法和装置 |
CN103983274A (zh) * | 2014-04-11 | 2014-08-13 | 湖北航天技术研究院总体设计所 | 一种适用于低精度无方位基准双轴转位设备的惯性测量单元标定方法 |
CN103994775B (zh) * | 2014-04-11 | 2017-01-04 | 湖北航天技术研究院总体设计所 | 一种适用于低精度有方位基准双轴转位设备的惯性测量单元标定方法 |
CN105842481A (zh) * | 2016-03-18 | 2016-08-10 | 重庆邮电大学 | 一种惯性导航系统中基于十字叉的三轴加速度计标定方法 |
CN106840100A (zh) * | 2017-03-13 | 2017-06-13 | 蒋海涛 | 一种数字式倾角传感器及测量方法 |
CN107505662A (zh) * | 2017-09-25 | 2017-12-22 | 华中科技大学 | 一种三轴mems重力仪 |
CN109613183A (zh) * | 2019-01-03 | 2019-04-12 | 北京英视睿达科技有限公司 | 低精度监测设备数据的修定方法及系统 |
CN110530395A (zh) * | 2019-05-31 | 2019-12-03 | 北京航天时代光电科技有限公司 | 小量程加速度计标定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103323625B (zh) | 2014-10-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103323625B (zh) | 一种mems-imu中加速度计动态环境下的误差标定补偿方法 | |
CN101290326B (zh) | 石英挠性加速度计测量组件的参数辨识标定方法 | |
CN101949710B (zh) | Gnss辅助mems惯性传感器零偏的快速在线动态标定方法 | |
CN101424534B (zh) | 惯性/重力组合导航半实物模拟装置 | |
CN102692239B (zh) | 一种基于旋转机构的光纤陀螺八位置标定方法 | |
CN103808331A (zh) | 一种mems三轴陀螺仪误差标定方法 | |
CN103196462A (zh) | 一种mimu中mems陀螺仪的误差标定补偿方法 | |
CN101178313A (zh) | 适合于光纤陀螺捷联惯性导航系统的地速检测方法 | |
CN102589573A (zh) | 微型组合导航系统中的传感器野外标定方法 | |
KR101658473B1 (ko) | Mems자이로스코프의 가속도 민감도 보정 방법 | |
CN104931046A (zh) | 一种微型惯性测量系统 | |
CN103411623A (zh) | 速率陀螺校准方法 | |
CN102636184B (zh) | 无角运动环境下基于离心机的挠性陀螺比力敏感项标定方法 | |
CN110383002A (zh) | 空间取向的确定 | |
CN102455183A (zh) | 三轴姿态传感器 | |
CN104931047A (zh) | 一种基于稳压电路的微型惯性测量系统 | |
CN108917788B (zh) | 一种惯性平台系统加速度计动态精度的测试方法和系统 | |
Yingbo et al. | Calibration method of quartz accelerometer on dynamic centrifuge | |
Li et al. | Status quo and developing trend of MEMS-gyroscope technology | |
CN104931053A (zh) | 一种具有供电优化的微型惯性测量系统 | |
CN103868527A (zh) | 一种标定捷联惯性组合加速度计组合的方法 | |
CN103940443B (zh) | 一种陀螺仪误差标定的方法 | |
CN102589568B (zh) | 车辆捷联惯性导航系统的三轴陀螺常值漂移快速测量方法 | |
CN102607590B (zh) | 一种基于光纤陀螺监测的挠性陀螺过载项抗干扰测试装置 | |
CN105758422A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |