CN109752568A - 基于主成分分析的微电子机械系统加速度计标定方法 - Google Patents
基于主成分分析的微电子机械系统加速度计标定方法 Download PDFInfo
- Publication number
- CN109752568A CN109752568A CN201910077271.XA CN201910077271A CN109752568A CN 109752568 A CN109752568 A CN 109752568A CN 201910077271 A CN201910077271 A CN 201910077271A CN 109752568 A CN109752568 A CN 109752568A
- Authority
- CN
- China
- Prior art keywords
- axis
- accelerometer
- amount
- error model
- calibrated
- 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
Landscapes
- Micromachines (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开一种基于主成分分析的微电子机械系统加速度计标定方法,包括:(10)采集包括加速度计三轴原始数据和二轴转台码盘数据;(20)构建待标定加速度计的线性误差模型;(30)将二轴转台码盘数据代入线性误差模型,优化求解线性误差模型各项系数;(40)将所述各轴理论输出与各轴对应的原始数据比较,得到各轴残差量;(50)分析待标定加速度计各轴的残差量,找到该轴影响最大非线性误差量;(60)构成完整的加速度计误差模型,得到各项系数值,用所述整体误差模型及各项系数值对微电子机械系统加速度计进行标定。本发明的微电子机械系统加速度计标定方法,成本低、效率高、适应性好。
Description
技术领域
本发明属于微电子机械系统设备标定技术领域,特别是一种成本低、效率高、适应性好的基于主成分分析的微电子机械系统加速度计标定方法。
背景技术
惯性导航(Inertial Navigation)技术是于上世纪中叶开始逐步发展起来的一项完全自主式的导航技术。这项技术自助性强、不易受外界环境干扰、输出信息量较大、隐蔽性较强、实时性高等特点,因此应用范围和领域非常广。微电子机械系统(micro-electromechanical system,MEMS)加速度计是微小型惯性导航系统中非常重要的元件之一,其性能的优劣程度直接决定了惯性导航系统的导航精度。它能利用感受加速度并将其转换为电信号的方式来测量加速力,从而感测运载体的线运动信息。据统计在由各种因素引起的总的导航误差中,有70%以上的误差是由惯性仪表引起的,因此提高加速度计精度有着非常重要的意义。
加速度计误差的误差源主要分为确定性误差和非确定性误差。确定性误差主要有零偏、刻度因子、轴系失准等系统误差;非确定性误差主要有量化噪声、白噪声、偏差不稳定性、随机游走等随机误差。
加速度计的标定方法根据观测量、标定场所、标定层次的不同可分为系统级标定及分立标定、外场标定及实验室标定、系统标定及器件标定。可以通过转台对加速度计进行分立标定,标定各个相关误差参数。
因此现有的标定方法,对转台的依赖性较高,为了达到标定精度,通常选用精密线振动台。精密线振动台是一种重要的用来模拟高加速度环境的惯导测试设备,它可以为被测惯性仪表提供足够高的输入比力,对加速度计的高阶非线性误差模型系数产生充足的激励,以提高其所属惯导系统的工作精度。但是这种设备体积大,价格昂贵,也加大了标定难度,限制了标定精度。而且针对不同原理的加速度计,影响其非线性误差的原因会存在差异,统一的非线性误差模型会限制其标定精度。
总之,现有技术存在的问题是:微电子机械系统加速度计标定成本高、效率低、适应性差。
发明内容
本发明的目的在于提供一种基于主成分分析的微电子机械系统加速度计标定方法,成本低、效率高、适应性好。
实现本发明目的的技术解决方案为:
一种基于主成分分析的微电子机械系统加速度计标定方法,包括如下步骤:
(10)实验数据采集:将待标定加速度计置于包含能够进行码盘读数的二轴转台的标定平台上进行角速度计标定实验,采集包括加速度计三轴原始数据和二轴转台码盘数据;
(20)线性误差模型建立:根据加速度计的零偏、三轴比例因子和三轴之间的非正交量关系,构建待标定加速度计的线性误差模型;
(30)线性误差模型求解:将二轴转台码盘数据代入线性误差模型,优化求解线性误差模型各项系数;
(40)各轴残差量获取:将线性误差模型各项系数代入线性误差模型,求得待标定加速度计各轴理论输出,将所述各轴理论输出与各轴对应的原始数据比较,得到各轴残差量;
(50)各轴最大非线性误差量获取:分析待标定加速度计各轴的残差量,找到该轴影响最大非线性误差量;
(60)加速度计标定:联立所找出的各轴的最大的非线性误差量,结合线性模型,构成完整的加速度计误差模型,所有系数一起参与优化,得到各项系数值,用所述整体误差模型及各项系数值对微电子机械系统加速度计进行标定。
本发明与现有技术相比,其显著优点为:
1、成本低:传统的标定方法需要借助高精度,高成本的精密转台等实验设备,而本发明方法通过计算机辅助处理实验数据,降低了对原始数据的精度要求,能将那些采集于精度不高的实验设备的实验数据进行修正,并最终将所有参数一起优化寻找到最优的系统误差模型,因此降低了对实验设备的依赖性,从而有效的控制了实验成本;
2、适应性好、效率高:本发明方法可以针对不同种的MEMS加速度计设计不同的非线性误差模型,因为不同加速度计的非线性特征不一样,也就是说本发明可以适应无论是内部结构导致的误差模型的差异性还是由于安装装配过程中导致的误差模型差异性。因此适用范围广,标定过程不繁琐,这也大大提高了标定效率。
下面结合附图和具体实施方式对本发明作进一步的详细描述。
附图说明
图1是本发明基于主成分分析的微电子机械系统加速度计标定方法的主流程图。
图2是图1中各轴非线性误差影响量获取步骤的流程图。
具体实施方式
如图1所示,本发明基于主成分分析的微电子机械系统加速度计标定方法,包括如下步骤:
(10)实验数据采集:将待标定加速度计置于包含能够进行码盘读数的二轴转台的标定平台上进行角速度计标定实验,采集包括加速度计三轴原始数据和二轴转台码盘数据;
例如,将MEMS加速度计放置在能够进行码盘读数的二轴转台的标定平台上进行角速度计标定实验,采集包括加速度计三轴原始数据和二轴转台码盘数据;
(20)线性误差模型建立:根据加速度计的零偏、三轴比例因子和三轴之间的非正交量关系,构建待标定加速度计的线性误差模型;
例如,针对上述MEMS加速度计,可以建立线性误差模型如下式:
式中为理论值,为实验值,(Sx,Sy,Sz)为三轴比例因子,(k12,k12,k21,k23,k31,k32)为三轴之间的非正交量,(bx,by,bz)为三轴零偏;
(30)线性误差模型求解:将二轴转台码盘数据代入线性误差模型,优化求解线性误差模型各项系数;
(40)各轴残差量获取:将线性误差模型各项系数代入线性误差模型,求得待标定加速度计各轴理论输出,将所述各轴理论输出与各轴对应的原始数据比较,得到各轴残差量;
(50)各轴最大非线性误差量获取:分析待标定加速度计各轴的残差量,找到该轴影响最大非线性误差量;
每个轴分开讨论,即使出现由于操作不当而引起的实验误差也只需要将那某一步重新实验,采集数据即可,不需要将整个实验重新做,提高了标定效率;且此方法适用于不同原理的加速度计,具有普适性。
如图2所示,所述(50)各轴非线性误差影响量获取步骤包括:
(51)特征样例获取:根据待标定加速度计各轴残差量,列出各轴非线性特征模型,利用加速度计三轴原始数据粗拟合出各轴对应的特征系数,构成该轴特征样例如下,
AA=[ai];
式中,ai表示该轴的加速度计的第i个非线性误差特征;
以x轴为例,构成该轴特征样例如下,
其中,
p01=0.0005809;p02=3.327e-06;p03=-0.0003614;
p04=-0.0004825;p05=4.906e-05;p06=0.001881;p07=-1.842e-05;
(52)均值分析:按下式得到均值归一化值,
其中,
(53)协方差矩阵计算:根据下式得到用以衡量两个变量总体误差的协方差矩阵,
当cov(ai,aj)>0时,表明ai与aj正相关,
当cov(ai,aj)<0时,表明ai与aj负相关,
当cov(ai,aj)=0时,表明ai与aj不相关;
例中,对7组变量,构建出7*7的协方差矩阵
(54)特征值特征向量计算:求解矩阵A的齐次方程组:
|A-λE|=0
式中,E是n维单位阵,方程组的解[λi]为矩阵A的i个特征值,对应的基础解系[xi]为对应的特征向量,对向量做归一化,得到对应的归一化后的特征向量;
(55)特征向量矩阵计算:
将特征值按数值从大到小排序,并将对应的特征向量按列向量的方式组合成特征向量矩阵,如下式,
M=[xi];
(56)最大影响非线性误差量选取:
对特征向量矩阵做投影,得到投影数据,
MM=ΔAA'×M,
选取投影数据MM中最小的一个作为最大影响非线性误差量。
例中,选择的最大影响非线性误差量为
(60)加速度计标定:联立所找出的各轴的最大的非线性误差量,结合线性模型,构成完整的加速度计误差模型,所有系数一起参与优化,得到各项系数值,用所述整体误差模型及各项系数值对微电子机械系统加速度计进行标定。此步的整体优化降低了对实验原始数据精度的依赖性,因此本方法不需要依赖高精度的实验设备。
例中,结合线性模型,构成完整的加速度计误差模型,如下
所有系数一起参与优化,得到各项系数值,用所述整体误差模型及各项系数值对微电子机械系统加速度计进行标定。
本发明方法通过计算机辅助处理实验数据,降低了对原始数据的精度要求,能将那些采集于精度不高的实验设备的实验数据进行修正,并最终将所有参数一起优化寻找到最优的系统误差模型,因此降低了对实验设备的依赖性,从而有效的控制了实验成本。同时,本发明方法可以针对不同种的MEMS加速度计设计不同的非线性误差模型,因为不同加速度计的非线性特征不一样,也就是说本发明可以适应无论是内部结构导致的误差模型的差异性还是由于安装装配过程中导致的误差模型差异性。因此适用范围广,标定过程不繁琐,这也大大提高了标定效率。
Claims (2)
1.一种基于主成分分析的微电子机械系统加速度计标定方法,其特征在于,包括如下步骤:
(10)实验数据采集:将待标定加速度计置于包含能够进行码盘读数的二轴转台的标定平台上进行角速度计标定实验,采集包括加速度计三轴原始数据和二轴转台码盘数据;
(20)线性误差模型建立:根据加速度计的零偏、三轴比例因子和三轴之间的非正交量关系,构建待标定加速度计的线性误差模型;
(30)线性误差模型求解:将二轴转台码盘数据代入线性误差模型,优化求解线性误差模型各项系数;
(40)各轴残差量获取:将线性误差模型各项系数代入线性误差模型,求得待标定加速度计各轴理论输出,将所述各轴理论输出与各轴对应的原始数据比较,得到各轴残差量;
(50)各轴最大非线性误差量获取:分析待标定加速度计各轴的残差量,找到该轴影响最大非线性误差量;
(60)加速度计标定:联立所找出的各轴的最大的非线性误差量,结合线性模型,构成完整的加速度计误差模型,所有系数一起参与优化,得到各项系数值,用所述整体误差模型及各项系数值对微电子机械系统加速度计进行标定。
2.根据权利要求1所述的加速度计标定方法,其特征在于,所述(50)各轴非线性误差影响量获取步骤包括:
(51)特征样例获取:根据待标定加速度计各轴残差量,列出各轴非线性特征模型,利用加速度计三轴原始数据粗拟合出各轴对应的特征系数,构成该轴特征样例如下,
AA=[ai];
式中,ai表示该轴的加速度计的第i个非线性误差特征;
(52)均值分析:按下式得到均值归一化值,
其中,
(53)协方差矩阵计算:根据下式得到用以衡量两个变量总体误差的协方差矩阵,
当cov(ai,aj)>0时,表明ai与aj正相关,
当cov(ai,aj)<0时,表明ai与aj负相关,
当cov(ai,aj)=0时,表明ai与aj不相关;
(54)特征值特征向量计算:求解矩阵A的齐次方程组:
|A-λE|=0
式中,E是n维单位阵,方程组的解[λi]为矩阵A的i个特征值,对应的基础解系[xi]为对应的特征向量,对向量做归一化,得到对应的归一化后的特征向量;
(55)特征向量矩阵计算:
将特征值按数值从大到小排序,并将对应的特征向量按列向量的方式组合成特征向量矩阵,如下式,
M=[xi];
(56)最大影响非线性误差量选取:
对特征向量矩阵做投影,得到投影数据,
MM=ΔAA'×M,
选取投影数据MM中最小的一个作为最大影响非线性误差量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910077271.XA CN109752568B (zh) | 2019-01-28 | 2019-01-28 | 基于主成分分析的微电子机械系统加速度计标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910077271.XA CN109752568B (zh) | 2019-01-28 | 2019-01-28 | 基于主成分分析的微电子机械系统加速度计标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109752568A true CN109752568A (zh) | 2019-05-14 |
CN109752568B CN109752568B (zh) | 2020-12-04 |
Family
ID=66406275
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910077271.XA Active CN109752568B (zh) | 2019-01-28 | 2019-01-28 | 基于主成分分析的微电子机械系统加速度计标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109752568B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110927822A (zh) * | 2019-12-03 | 2020-03-27 | 吉林大学 | 一种评估希尔伯特变换算法获得磁梯度张量准确性的方法 |
DE102022126970A1 (de) | 2022-10-14 | 2024-04-25 | Deutsches Zentrum für Luft- und Raumfahrt e.V. | Kalibrierung eines Beschleunigungssensors |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6360602B1 (en) * | 1999-07-29 | 2002-03-26 | Litton Systems, Inc. | Method and apparatus reducing output noise in a digitally rebalanced accelerometer |
US7720661B2 (en) * | 2005-02-28 | 2010-05-18 | Honeywell International Inc. | Low vibration rectification in a closed-loop, in-plane MEMS device |
US20110184645A1 (en) * | 2010-01-28 | 2011-07-28 | Sirf Technology Holdings, Inc. | Use of accelerometer only data to improve gnss performance |
WO2014001470A1 (fr) * | 2012-06-29 | 2014-01-03 | Movea | Procede de calibration continue d'un capteur |
CN103591965A (zh) * | 2013-09-12 | 2014-02-19 | 哈尔滨工程大学 | 一种舰载旋转式捷联惯导系统在线标定的方法 |
CN103759742A (zh) * | 2014-01-22 | 2014-04-30 | 东南大学 | 基于模糊自适应控制技术的捷联惯导非线性对准方法 |
CN103941042A (zh) * | 2014-04-28 | 2014-07-23 | 北京航天控制仪器研究所 | 一种陀螺加速度计多位置误差系数标定方法 |
CN105938149A (zh) * | 2016-06-24 | 2016-09-14 | 南京理工大学 | 一种基于加速度传感器的加速度记录仪校准装置及方法 |
CN105956609A (zh) * | 2016-04-23 | 2016-09-21 | 北京化工大学 | 一种基于svm的加速度计混合动态模型建立的方法 |
CN106959121A (zh) * | 2017-02-28 | 2017-07-18 | 北京联合大学 | 一种自主反向调优的超限学习算法在磁罗盘误差补偿中的应用方法 |
CN108982918A (zh) * | 2018-07-27 | 2018-12-11 | 北京航天控制仪器研究所 | 基准不确定情况下加速度计组合误差系数分离与标定方法 |
-
2019
- 2019-01-28 CN CN201910077271.XA patent/CN109752568B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6360602B1 (en) * | 1999-07-29 | 2002-03-26 | Litton Systems, Inc. | Method and apparatus reducing output noise in a digitally rebalanced accelerometer |
US7720661B2 (en) * | 2005-02-28 | 2010-05-18 | Honeywell International Inc. | Low vibration rectification in a closed-loop, in-plane MEMS device |
US20110184645A1 (en) * | 2010-01-28 | 2011-07-28 | Sirf Technology Holdings, Inc. | Use of accelerometer only data to improve gnss performance |
WO2014001470A1 (fr) * | 2012-06-29 | 2014-01-03 | Movea | Procede de calibration continue d'un capteur |
CN103591965A (zh) * | 2013-09-12 | 2014-02-19 | 哈尔滨工程大学 | 一种舰载旋转式捷联惯导系统在线标定的方法 |
CN103759742A (zh) * | 2014-01-22 | 2014-04-30 | 东南大学 | 基于模糊自适应控制技术的捷联惯导非线性对准方法 |
CN103941042A (zh) * | 2014-04-28 | 2014-07-23 | 北京航天控制仪器研究所 | 一种陀螺加速度计多位置误差系数标定方法 |
CN105956609A (zh) * | 2016-04-23 | 2016-09-21 | 北京化工大学 | 一种基于svm的加速度计混合动态模型建立的方法 |
CN105938149A (zh) * | 2016-06-24 | 2016-09-14 | 南京理工大学 | 一种基于加速度传感器的加速度记录仪校准装置及方法 |
CN106959121A (zh) * | 2017-02-28 | 2017-07-18 | 北京联合大学 | 一种自主反向调优的超限学习算法在磁罗盘误差补偿中的应用方法 |
CN108982918A (zh) * | 2018-07-27 | 2018-12-11 | 北京航天控制仪器研究所 | 基准不确定情况下加速度计组合误差系数分离与标定方法 |
Non-Patent Citations (5)
Title |
---|
HOUZHOU 等: ""The application of ICP accelerometer sensor signal in a Data Record"", 《ICCASM2010》 * |
张纯: ""加速度传感器自动标定系统的研究与实现"", 《中国优秀硕士学位论文全文数据库》 * |
李永光: ""高精度加速度计标定与补偿技术研究"", 《中国优秀硕士学位论文全文数据库》 * |
王海龙: ""惯性平台系统自标定实验设计与辨识"", 《中国优秀硕士学位论文全文数据库》 * |
陈林: ""低成本MEMS惯性测量组件复合标定与误差补偿方法研究"", 《中国优秀硕士学位论文全文数据库》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110927822A (zh) * | 2019-12-03 | 2020-03-27 | 吉林大学 | 一种评估希尔伯特变换算法获得磁梯度张量准确性的方法 |
CN110927822B (zh) * | 2019-12-03 | 2021-05-18 | 吉林大学 | 一种评估希尔伯特变换算法获得磁梯度张量准确性的方法 |
DE102022126970A1 (de) | 2022-10-14 | 2024-04-25 | Deutsches Zentrum für Luft- und Raumfahrt e.V. | Kalibrierung eines Beschleunigungssensors |
Also Published As
Publication number | Publication date |
---|---|
CN109752568B (zh) | 2020-12-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100368774C (zh) | 光纤陀螺惯测装置快速启动和精度保证的工程实现方法 | |
Zhang et al. | Allan variance analysis on error characters of MEMS inertial sensors for an FPGA-based GPS/INS system | |
US9250300B2 (en) | Dynamic magnetometer calibration | |
CN109752568A (zh) | 基于主成分分析的微电子机械系统加速度计标定方法 | |
CN101576489A (zh) | 光干涉检测甲烷或二氧化碳装置及精度自动补偿检测方法 | |
CN103808349B (zh) | 矢量传感器的误差校正方法和装置 | |
CN112762962B (zh) | 基于温度滞后模型的微机电系统加速度计零偏补偿方法 | |
CN110553641A (zh) | 一种基于相关系数的提高惯性制导精度的方法 | |
CN108225374A (zh) | 一种融合遗传算法的Allan方差分析法 | |
CN103033198A (zh) | 一种设置光纤陀螺仿真信号中随机误差参数的方法 | |
CN110375772B (zh) | 自适应卡尔曼滤波的环形激光器随机误差建模与补偿方法 | |
CN101661047B (zh) | 一种捷联惯导加速度计输出转换装置 | |
CN112082547A (zh) | 组合导航系统优化方法、装置、电子设备及存储介质 | |
KR101331956B1 (ko) | 아날로그 mems 센서를 이용하는 고정밀 ins 모듈 및 그 구동 방법 | |
CN109902357B (zh) | 一种复变差分的光滑非线性结构动响应灵敏度分析方法 | |
CN110726852A (zh) | 一种mems加速度计温度补偿方法 | |
CN112067844B (zh) | Mems传感器阵列式高精度输出控制方法 | |
CN116026328A (zh) | 微惯导的零偏滞回效应补偿模型的构建方法和补偿方法 | |
CN108398576B (zh) | 一种静态误差标定系统及方法 | |
CN110888142A (zh) | 基于mems激光雷达测量技术的航天器隐藏目标点测量方法 | |
CN111707288B (zh) | 一种无水平台的微型惯性传感器零位测量方法 | |
CN113865621B (zh) | 任意六位置MEMS陀螺仪及其g值敏感系数标定方法 | |
CN113639769B (zh) | 一种基于散斑相关性提高传感器灵敏度的方法 | |
US20040179208A1 (en) | Structure for sophisticated surveying instrument with coordinate board for position identification | |
CN116151172B (zh) | Mems器件模型构建方法、装置及设计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |