CN113295887A - 一种微机电三轴加速度计交叉耦合系数快速标定方法 - Google Patents

一种微机电三轴加速度计交叉耦合系数快速标定方法 Download PDF

Info

Publication number
CN113295887A
CN113295887A CN202110528736.6A CN202110528736A CN113295887A CN 113295887 A CN113295887 A CN 113295887A CN 202110528736 A CN202110528736 A CN 202110528736A CN 113295887 A CN113295887 A CN 113295887A
Authority
CN
China
Prior art keywords
accelerometer
axis
axis accelerometer
follows
electromechanical
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
Application number
CN202110528736.6A
Other languages
English (en)
Other versions
CN113295887B (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.)
Kaifeng Navigation Control Technology Co ltd
Original Assignee
Beijing Kaikai Hangyu Navigation Control 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 Beijing Kaikai Hangyu Navigation Control Technology Co ltd filed Critical Beijing Kaikai Hangyu Navigation Control Technology Co ltd
Priority to CN202110528736.6A priority Critical patent/CN113295887B/zh
Publication of CN113295887A publication Critical patent/CN113295887A/zh
Application granted granted Critical
Publication of CN113295887B publication Critical patent/CN113295887B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P21/00Testing or calibrating of apparatus or devices covered by the preceding groups

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Force Measurement Appropriate To Specific Purposes (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)

Abstract

本发明公开了一种微机电三轴加速度计交叉耦合系数快速标定测试方法,主要实现步骤为:将微机电三轴加速度计通过高精度标定测试工装,安装到三轴转台上,设计加速度计标定位置编排方案为六位置标定测试方法;建立准确的加速度计误差模型,通过误差模型计算加速度计的标度因数、零位和交叉耦合系数。针对现有技术标定时间长、操作繁琐、计算复杂的缺点,该方法简单、快速,极大提高生产效率。

Description

一种微机电三轴加速度计交叉耦合系数快速标定方法
技术领域
本发明涉及一种微机电三轴加速度计交叉耦合系数快速标定测试方法,特别适用于小体积MEMS三轴加速度计交叉耦合系数的标定,属于惯性测量单元技术领域,可适用于快速标定三轴加速度计的场合。
背景技术
MEMS加速度计是采用微机电(Micro Electro-Mechanical System,MEMS)技术设计的惯性器件,由它们构成的惯性测量单元测量载体相对惯性空间的角速度和加速度,能够自主、连续、实时提供位置、速度和姿态信息,具有隐蔽性好,不受气候条件限制,无信号丢失,不受干扰等优点,在国民经济和国防安全等众多领域有重要的用途。
微机电三轴加速度计是一个MEMS芯片,里面包含了三个正交设置的加速度计,微机电三轴加速度计是惯性测量单元的核心部件,它的精度决定了MIMU的测量精度,MIMU误差源中确定性的系统误差占总误差的绝大部分,因此,对确定性的系统误差进行分离、标定、补偿对提高惯性测量单元精度有决定性的作用。
目前,三轴加速度计的交叉耦合系数的标定基本采用多位置静态标定,比如8位置、12位置、24位置等多种标定方法,可以比较精确的标定加速度计交叉耦合系数,但是标定时间长,操作繁琐,计算复杂。
发明内容
本发明的技术解决问题是:针对现有技术标定时间长、操作繁琐、计算复杂的缺点,提出一种微机电三轴加速度计交叉耦合系数快速标定测试方法,该方法简单、快速,极大提高生产效率。
本发明的技术解决方案是:利用三轴转台,采用六位置标定下特有的误差模型来计算三轴加速度计交叉耦合系数。具体步骤如下:
【1】将微机电三轴加速度计通过高精度标定测试工装,安装到三轴转台上,设计加速度计标定位置编排方案为六位置标定测试方法,具体位置如下:
(1)位置1是X轴加速度计向上,Y轴加速度计和Z轴加速度计处于水平位置。
(2)位置2是X轴加速度计向下,Y轴加速度计和Z轴加速度计处于水平位置。
(3)位置3是Y轴加速度计向上,X轴加速度计和Z轴加速度计处于水平位置。
(4)位置4是Y轴加速度计向下,X轴加速度计和Z轴加速度计处于水平位置。
(5)位置5是Z轴加速度计向上,X轴加速度计和Y轴加速度计处于水平位置。
(6)位置6是Z轴加速度计向下,X轴加速度计和Y轴加速度计处于水平位置。
【2】建立准确的加速度计误差模型,通过误差模型计算加速度计的标度因数、零位和交叉耦合系数,具体公式如下:
每个位置下的加速度计方程如下:
Figure BDA0003066242530000021
式中:
l表示加速度计输出数据的平均值,k表示加速度计标度因数,m表示加速度计交叉耦合系数,b表示加速度计零位,a表示三轴转台输入的重力加速度,下标x,y,z表示三个轴向的加速度计。
六个位置下的输入重力加速度表示如下:
a1=[g 0 0]T
a2=[-g 0 0]T
a3=[0 g 0]T
a4=[0 -g 0]T
a5=[0 0 g]T
a6=[0 0 -g]T
六位置下的最小二乘矩阵设计为如下形式
Figure BDA0003066242530000031
相应的,六个位置下三轴加速度计输出平均值表示如下:
u1=[lx1 ly1 lz1]T X轴向上时
u2=[lx2 ly2 lz2]T X轴向下时
u3=[lx3 ly3 lz3]T Y轴向上时
u4=[lx4 ly4 lz4]T Y轴向下时
u5=[lx5 ly5 lz5]T Z轴向上时
u6=[lx6 ly6 lz6]T Z轴向下时
将这六个位置的加速度计输出数据表示为一个矩阵形式,如下:
U=[u1 u2 u3 u4 u5 u6]
使用最小二乘方法,计算得到加速度计的标度因数、零位、交叉耦合系数,方程如下:
M=U·AT·(AAT)-1
本发明的优点在于:
本发明根据简单的六位置标定方法得到的加速度计测试数据,利用一种特定的最小二乘矩阵计算获得准确的微机电三轴加速度计交叉耦合系数,该方法简单、快速,极大提高生产效率。
附图说明
图1为本发明的流程图。
图2为本发明的六个位置示意图。
具体实施方式
下面将结合本发明中的附图,对本发明中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请参阅图1-2,本发明提供了一种微机电三轴加速度计交叉耦合系数快速标定测试方法,下面以某种三轴微机电加速度计为例进行说明。
三轴微机电加速度计标定方法的具体流程按图1进行。
【1】将微机电三轴加速度计通过高精度标定测试工装,安装到三轴转台上,设计加速度计标定位置编排方案为六位置标定测试方法,具体位置如图2所示,这六个翻转位置的具体含义是:
(1)位置1是X轴加速度计向上,Y轴加速度计和Z轴加速度计处于水平位置。
(2)位置2是X轴加速度计向下,Y轴加速度计和Z轴加速度计处于水平位置。
(3)位置3是Y轴加速度计向上,X轴加速度计和Z轴加速度计处于水平位置。
(4)位置4是Y轴加速度计向下,X轴加速度计和Z轴加速度计处于水平位置。
(5)位置5是Z轴加速度计向上,X轴加速度计和Y轴加速度计处于水平位置。
(6)位置6是Z轴加速度计向下,X轴加速度计和Y轴加速度计处于水平位置。
【2】建立准确的加速度计误差模型,通过误差模型计算加速度计的标度因数、零位和交叉耦合系数,具体公式如下:
每个位置下的加速度计方程如下:
Figure BDA0003066242530000051
六位置下转台输入的最小二乘矩阵设计为如下形式
Figure BDA0003066242530000052
相应的,六个位置下三轴加速度计输出计算得到平均值,将这六个位置的加速度计输出数据表示为一个矩阵形式,如下:
Figure BDA0003066242530000061
使用最小二乘方法,计算得到加速度计的标度因数、零位、交叉耦合系数,方程如下:
M=U·AT·(AAT)-1
以上对本发明所提供的一种微机电加速度计交叉耦合系数快速标定测试方法进行了详细介绍,对于本领域的一般技术人员,依据本发明的具体实施过程及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。本发明未详细描述内容为本领域技术人员公知技术。

Claims (1)

1.一种微机电三轴加速度计交叉耦合系数快速标定测试方法,其特征在于步骤如下:
【1】将微机电三轴加速度计通过高精度标定测试工装,安装到三轴转台上,设计加速度计标定位置编排方案为六位置标定测试方法,具体位置如下:
(1)位置1是X轴加速度计向上,Y轴加速度计和Z轴加速度计处于水平位置。
(2)位置2是X轴加速度计向下,Y轴加速度计和Z轴加速度计处于水平位置。
(3)位置3是Y轴加速度计向上,X轴加速度计和Z轴加速度计处于水平位置。
(4)位置4是Y轴加速度计向下,X轴加速度计和Z轴加速度计处于水平位置。
(5)位置5是Z轴加速度计向上,X轴加速度计和Y轴加速度计处于水平位置。
(6)位置6是Z轴加速度计向下,X轴加速度计和Y轴加速度计处于水平位置。
【2】建立准确的加速度计误差模型,通过误差模型计算加速度计的标度因数、零位和交叉耦合系数,具体公式如下:
每个位置下的加速度计方程如下:
Figure FDA0003066242520000011
式中:
l表示加速度计输出数据的平均值,k表示加速度计标度因数,m表示加速度计交叉耦合系数,b表示加速度计零位,a表示三轴转台输入的重力加速度,下标x,y,z表示三个轴向的加速度计。
六个位置下的输入重力加速度表示如下:
a1=[g 0 0]T
a2=[-g 0 0]T
a3=[0 g 0]T
a4=[0 -g 0]T
a5=[0 0 g]T
a6=[0 0 -g]T
六位置下的最小二乘矩阵设计为如下形式
Figure FDA0003066242520000021
相应的,六个位置下三轴加速度计输出平均值表示如下:
u1=[lx1 ly1 lz1]T X轴向上时
u2=[lx2 ly2 lz2]T X轴向下时
u3=[lx3 ly3 lz3]T Y轴向上时
u4=[lx4 ly4 lz4]T Y轴向下时
u5=[lx5 ly5 lz5]T Z轴向上时
u6=[lx6 ly6 lz6]T Z轴向下时
将这六个位置的加速度计输出数据表示为一个矩阵形式,如下:
U=[u1 u2 u3 u4 u5 u6]
使用最小二乘方法,计算得到加速度计的标度因数、零位、交叉耦合系数,方程如下:
M=U·AT·(AAT)-1
CN202110528736.6A 2021-05-14 2021-05-14 一种微机电三轴加速度计交叉耦合系数快速标定方法 Active CN113295887B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110528736.6A CN113295887B (zh) 2021-05-14 2021-05-14 一种微机电三轴加速度计交叉耦合系数快速标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110528736.6A CN113295887B (zh) 2021-05-14 2021-05-14 一种微机电三轴加速度计交叉耦合系数快速标定方法

Publications (2)

Publication Number Publication Date
CN113295887A true CN113295887A (zh) 2021-08-24
CN113295887B CN113295887B (zh) 2024-04-05

Family

ID=77322335

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110528736.6A Active CN113295887B (zh) 2021-05-14 2021-05-14 一种微机电三轴加速度计交叉耦合系数快速标定方法

Country Status (1)

Country Link
CN (1) CN113295887B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101832782A (zh) * 2010-03-26 2010-09-15 中北大学 一种微惯性测量组合现场快速标定方法
CN104677381A (zh) * 2015-01-29 2015-06-03 中国空空导弹研究院 一种微惯性测量单元的测试系统
CN108592952A (zh) * 2018-06-01 2018-09-28 北京航空航天大学 基于杆臂补偿与正反倍速率同时标定多mimu误差的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101832782A (zh) * 2010-03-26 2010-09-15 中北大学 一种微惯性测量组合现场快速标定方法
CN104677381A (zh) * 2015-01-29 2015-06-03 中国空空导弹研究院 一种微惯性测量单元的测试系统
CN108592952A (zh) * 2018-06-01 2018-09-28 北京航空航天大学 基于杆臂补偿与正反倍速率同时标定多mimu误差的方法

Also Published As

Publication number Publication date
CN113295887B (zh) 2024-04-05

Similar Documents

Publication Publication Date Title
CN111678538B (zh) 一种基于速度匹配的动态水平仪误差补偿方法
CN102980577B (zh) 一种微型捷联航姿系统及其工作方法
JP5474938B2 (ja) 加速度計および誤差補償方法
CN106767917B (zh) 一种斜置冗余惯导系统标定误差模型建模方法
CN201561759U (zh) 惯性姿态方位测量装置
CN109709628B (zh) 一种旋转加速度计重力梯度仪标定方法
CN107202578B (zh) 一种基于mems技术的捷联式垂直陀螺仪解算方法
CN106969783A (zh) 一种基于光纤陀螺惯性导航的单轴旋转快速标定技术
CN103868648A (zh) 三轴气浮仿真实验平台的质心测量方法
CN111121820B (zh) 基于卡尔曼滤波的mems惯性传感器阵列融合方法
CN112683274A (zh) 一种基于无迹卡尔曼滤波的无人机组合导航方法和系统
CN113916257B (zh) 三轴mems加计组合惯性测量单元标定方法
CN112906185B (zh) 一种基于人工智能的mems惯性传感器异构阵列及其设计方法
CN105136142A (zh) 一种基于微惯性传感器的室内定位方法
CN113295887B (zh) 一种微机电三轴加速度计交叉耦合系数快速标定方法
CN116972875A (zh) 一种基于陀螺仪的物体运动轨迹监测方法
CN105842481A (zh) 一种惯性导航系统中基于十字叉的三轴加速度计标定方法
CN112729280B (zh) 一种基于多面体式阵列结构的微惯性姿态测量装置及方法
CN108398576A (zh) 一种静态误差标定系统及方法
CN110954131A (zh) 一种光纤陀螺输入轴失准角标定工装
CN113945227B (zh) 一种基于离心台的惯性测量单元全量程加速度标校方法
CN112336341A (zh) 基于惯性、位置传感器的人体跌倒检测系统
CN113865621B (zh) 任意六位置MEMS陀螺仪及其g值敏感系数标定方法
CN110824583A (zh) 一种航空重力仪用在线自主标定方法
CN101509934B (zh) 微加速度与微角速率单片集成传感器

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
TA01 Transfer of patent application right

Effective date of registration: 20240305

Address after: 230088 1056, North building, Hefei original animation Park Management Co., Ltd., No. 19, Tianzhi Road, high tech Zone, Hefei, Anhui Province

Applicant after: Kaifeng Navigation Control Technology Co.,Ltd.

Country or region after: Zhong Guo

Address before: 100176 room 716, building 4, yard 5, Liangshuihe 2nd Street, Daxing District, Beijing

Applicant before: Beijing Kaikai Hangyu Navigation Control Technology Co.,Ltd.

Country or region before: Zhong Guo

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant