本发明涉及加速度误差标定技术领域，特别是涉及一种静态误差标定系统及方法。
背景技术：
加速度传感器是一种非常重要的惯性传感器，在惯性导航、姿态测量、振动监测等领域发挥着核心作用。随着微机电系统(microelectro-mechanicalsystems)技术的飞速发展，mems三轴加速度传感器以其体积小、重量轻、功耗低、成本低和可靠性高等优点被大量应用于手机、平衡车、无人机、运动捕捉、虚拟现实等众多领域。其中应用最广泛的是mems三轴加速度传感器组合，能够实现载体的姿态感知和辅助导航等功能。
由于受工作原理和制造工艺的影响，mems三轴加速度传感器的测量精度较低。特别是出厂后的mems三轴加速度传感器存在一系列静态误差，其静态误差包括三轴零偏、三轴刻度因子和三轴轴间交叉误差共9个误差参数。对mems三轴加速度传感器的标定是指采用一定的技术手段确定这9个误差参数或其中几个参数的值，如果不通过标定对这些误差的值进行确定进而补偿，mems三轴加速度传感器的测量精度将大打折扣，甚至影响正常的使用。
目前mems三轴加速度传感器标定绝大多数采用多位置的方法。具体包括无依托的六位置法和依托转台的六位置法。
无依托的六位置法的具体步骤为：首先设计需要采集三轴加速度数据的六个位置；然后依次将mems三轴加速度传感器按六个位置放置，并采集其输出的加速度数据；最后对采集到的六个位置的数据进行处理，计算获得加速度传感器的零偏和刻度因子误差。目前大疆等无人机厂家的产品都采用无依托六位置法进行现场标定。虽然无依托六位置法对环境要求不高，仅需要一个平整的台面即可。但由于六次测量只能建立起6个包含误差参数的方程式，无法求解出全部9个误差参数，因此无依托六位置法仅能够标定三轴零偏和三轴刻度因子6个静态误差参数，无法标定轴间交叉误差，而轴间交叉误差是三轴mems三轴加速度传感器用于感知载体姿态时的主要误差来源，特别是对很大一部分物理正交的mems三轴加速度传感器而言，轴间正交性很难保证，其轴间交叉误差会比较严重。因此无依托六位置法无法标定mems三轴加速度传感器的全部静态误差参数，因此具有标定效果较差的问题。
依托转台的六位置法也比较常见，其步骤与无依托六位置法类似，不同在于需要将mems三轴加速度传感器固定于转台上，并将转台移动至设计的六个位置，采集相应数据并处理可获得加速度传感器的零偏、刻度因子和轴间交叉误差。依托转台的六位置法能够将mems三轴加速度传感器的每个轴与当地重力矢量对齐，能够从建立起的6个方程直接得到9个静态误差参数的解析表达式。因此依托转台的六位置法能够标定包括三轴零偏、刻度因子和轴间交叉误差在内的全部9个静态误差参数，但转台的成本非常高，且体积大、不便移动，无法用于现场标定。
技术实现要素：
本发明的目的是提供一种静态误差标定系统及方法，以实现全面标定静态误差，提高标定静态误差的精度，降低标定系统的体积和成本。
为实现上述目的，本发明提供一种静态误差标定系统，所述标定系统包括：三轴加速度传感器、辅助装置、水平台面；
三轴加速度传感器放置在水平台面上，获得六个位置；
三轴加速度传感器放置在辅助装置上，辅助装置放置在水平台面上，获得三个位置。
可选的，所述辅助装置包括凹槽，所述凹槽是由长方体向下凹陷形成的。
可选的，所述凹槽为“v”字形凹槽，所述“v”字形凹槽的角度为30°至60°之间的任一角度。
可选的，所述三轴加速度传感器放置在水平台面上，获得六个位置，具体包括：三轴加速度传感器的底面与水平台面接触放置；三轴加速度传感器的左面与水平台面接触放置；三轴加速度传感器的右面与水平台面接触放置；三轴加速度传感器的后面与水平台面接触放置；三轴加速度传感器的前面与水平台面接触放置；三轴加速度传感器的正面与水平台面接触放置；
所述三轴加速度传感器放置在辅助装置上，辅助装置放置在水平台面上，获得三个位置，具体包括：三轴加速度传感器的底面和右面分别与辅助装置的凹槽接触放置；三轴加速度传感器的底面和后面与辅助装置的凹槽接触放置；三轴加速度传感器的后面和左面分别与辅助装置的凹槽接触放置。
本发明还提供一种静态误差标定方法，所述方法包括：
分别获取三轴加速度传感器在各位置的m组测量数据；
根据三轴加速度传感器在各位置的m组测量数据确定各位置的均值向量；
获取初始静态误差参数向量；
根据初始静态误差参数向量和各位置的均值向量确定各位置的测量误差函数；
根据各位置的测量误差函数确定当前静态误差参数向量；
根据当前静态误差参数向量和各位置的测量误差函数确定代价函数值；
判断代价函数值是否小于设定阈值；如果代价函数值小于设定阈值，则将当前静态误差参数向量作为最优静态误差参数向量，输出最优静态误差参数向量；如果代价函数值大于等于设定阈值，则将当前静态误差参数向量作为初始静态误差参数向量，重新确定当前静态误差参数向量。
可选的，所述根据三轴加速度传感器各位置的m组测量数据确定各位置的均值向量，具体包括：
对所述在各位置的m组测量数据进行剔除野值处理，获得各位置的n组预处理数据；其中，n为小于等于m的整数；
根据各位置的n组预处理数据确定各位置的均值向量。
可选的，所述根据初始静态误差参数向量和各位置的均值向量确定各位置的测量误差函数，具体包括：
根据初始静态误差参数向量和各位置的均值向量确定各位置的静态误差参数模型；
根据各位置的静态误差参数模型确定各位置的测量误差函数。
可选的，所述根据各位置的测量误差函数确定当前静态误差参数向量，具体包括：
根据各位置的测量误差函数确定测量误差向量；
根据初始静态误差参数向量和测量误差向量确定当前静态误差参数向量。
可选的，所述各个位置具体包括：
三轴加速度传感器的底面与水平台面接触放置；三轴加速度传感器的左面与水平台面接触放置；三轴加速度传感器的右面与水平台面接触放置；三轴加速度传感器的后面与水平台面接触放置；三轴加速度传感器的前面与水平台面接触放置；三轴加速度传感器的正面与水平台面接触放置；三轴加速度传感器的底面和右面分别与辅助装置的凹槽接触放置；三轴加速度传感器的底面和后面与辅助装置的凹槽接触放置；三轴加速度传感器的后面和左面分别与辅助装置的凹槽接触放置。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
本发明设计辅助装置，先将三轴加速度传感器放置在水平台面上，获得六个位置；然后将三轴加速度传感器放置在辅助装置上，辅助装置放置在水平台面上，获得三个位置；获取三轴加速度传感器在9个位置处于不同姿态的多组测量数据，根据9组不同测量数据实现全方面标定三轴加速度传感器的静态误差，不仅提高了标定静态误差的精度，还摆脱了依托转台和精确水平面进行静态误差标定，所需设备仅为一水平台面和本发明设计的标定辅助装置即可，降低标定系统的体积和成本。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例静态误差标定系统结构图；
图2为本发明实施例静态误差标定系统坐标分析图；
图3为本发明实施例辅助装置结构示意图；
图4为本发明实施例静态误差标定方法流程图。
其中，1、三轴加速度传感器，2、水平台面，3、辅助装置。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明的目的是提供一种静态误差标定系统及方法，以实现全面标定静态误差，提高标定静态误差的精度，降低标定系统的体积和成本。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
受工作原理和制造工艺的影响，mems惯性器件的测量精度受到制约，在测量输出中总是包含着各种误差项，具体包括静态误差和随机误差两类。其中随机误差包括mems器件受环境(如温度、湿度)、电气噪声等因素影响产生的不确定误差，需要在实际使用过程中利用滤波等方法进行抑制。静态误差包括三轴的零偏、刻度因子和轴间交叉误差，零偏误差是由于mems加速度传感器在零输入的情况下输出不为零而引起的误差，刻度因子误差是实际输出与说明书上不同而引起的误差，轴间交叉误差是由于工艺原因测量的三轴之间不成交而引起的误差，通过标定可确定这些误差参数的值，进而进行补偿，使加速度传感器的输出更接近真实值。
图1为本发明实施例静态误差标定系统结构图，其中，图1中(a)-(i)为三轴加速度传感器放置的9个不同方式；图2为本发明实施例静态误差标定系统坐标分析图，其中，图2中(a)-(i)为三轴加速度传感器不同位置的坐标分析图；如图1-图2所示，本发明提供一种静态误差标定系统，所述标定系统包括：三轴加速度传感器1、辅助装置3、水平台面2。
图3为本发明实施例辅助装置结构示意图；(a)为辅助装置立体结构示意图；(b)为辅助装置主视图；(c)为辅助装置俯视图；(d)为辅助装置右视图；如图3所示，所述辅助装置3包括凹槽，所述凹槽是由长方体向下凹陷形成的；所述凹槽为“v”字形凹槽，所述“v”字形凹槽的角度为30°至60°之间的任一角度，本发明以45°为例构件的凹槽。
三轴加速度传感器1出厂时，厂家会标出三轴加速度传感器1的敏感轴，将三轴加速度传感器1的三个敏感轴记分别为x轴，y轴，z轴，其关系满足图2，在此基础上记x轴向为三轴加速度传感器1的前面，y轴向为右面，z轴向为底面，与x轴反向为后面，与y轴反向为左面，与z轴反向为正面。
如图1所示，三轴加速度传感器1放置在水平台面2上，获得六个位置；具体包括：三轴加速度传感器1的底面与水平台面2接触放置；三轴加速度传感器1的左面与水平台面2接触放置；三轴加速度传感器1的右面与水平台面2接触放置；三轴加速度传感器1的后面与水平台面2接触放置；三轴加速度传感器1的前面与水平台面2接触放置；三轴加速度传感器1的正面与水平台面2接触放置；三轴加速度传感器1放置在辅助装置3上，辅助装置3放置在水平台面2上，获得三个位置。所述三轴加速度传感器1放置在辅助装置3上，辅助装置3放置在水平台面2上，获得三个位置，具体包括：三轴加速度传感器1的底面和右面分别与辅助装置3的凹槽接触放置；三轴加速度传感器1的底面和后面与辅助装置3的凹槽接触放置；三轴加速度传感器1的后面和左面分别与辅助装置3的凹槽接触放置。
图4为本发明实施例静态误差标定方法流程图；如图4所示，本发明还提供一种静态误差标定方法，所述方法具体包括以下步骤：
步骤10：分别获取三轴加速度传感器1在各位置的m组测量数据；其中，m为大于等于10的整数。
步骤20：根据三轴加速度传感器1在各位置的m组测量数据确定各位置的均值向量。
步骤30：获取初始静态误差参数向量。
步骤40：根据初始静态误差参数向量和各位置的均值向量确定各位置的测量误差函数。
步骤50：根据各位置的测量误差函数确定当前静态误差参数向量。
步骤60：根据当前静态误差参数向量和各位置的测量误差函数确定代价函数值。
步骤70：判断代价函数值是否小于设定阈值；如果代价函数值小于设定阈值，则将当前静态误差参数向量作为最优静态误差参数向量，输出最优静态误差参数向量；如果代价函数值大于等于设定阈值，则将当前静态误差参数向量作为初始静态误差参数向量，重新确定各位置的测量误差函数。
下面对给个步骤进行详细论述：
步骤10：分别获取三轴加速度传感器在各位置的m组测量数据；其中，m为大于等于10的整数。本发明获取测量数据的个数以及获取测量数据的频率根据实际需求确定。
步骤20：所述根据三轴加速度传感器1在各位置的m组测量数据确定各位置的均值向量
步骤201：对所述在各位置的m组测量数据进行剔除野值处理，获得各位置的n组预处理数据；其中，n为小于等于m的整数。
步骤202：根据各位置的n组预处理数据确定各位置的均值向量
步骤30：获取初始静态误差参数向量β(s)，其中s为大于等于0的整数。
当s＝0时，则β(0)＝(111000000)t，β(0)为理想无误差情况；当s＝1时，则β(1)＝(sx(1)sy(1)sz(1)kxy(1)kxz(1)kyz(1)bx(1)by(1)bz(1))t；当s＝2时，则β(2)＝(sx(2)sy(2)sz(2)kxy(2)kxz(2)kyz(2)bx(2)by(2)bz(2))t；依次类推，初始静态误差参数向量β(s)的具体公式为：β(s)＝(sx(s)sy(s)sz(s)kxy(s)kxz(s)kyz(s)bx(s)by(s)bz(s))t，其中，sx(s)、sy(s)、sz(s)分为第s步x轴、y轴、z轴向的零偏；kxy(s)、kxz(s)、kyz(s)分为第s步xy轴、xz轴、yz轴的轴间交叉误差；bx(s)、by(s)、bz(s)分别为第s步x轴、y轴、z轴向的刻度因子。
步骤40：所述根据初始静态误差参数向量β(s)和各位置的均值向量
步骤401：根据初始静态误差参数向量β(s)和各位置的均值向量
步骤402：根据各位置的静态误差参数模型
其中，
步骤50：所述根据各位置的测量误差函数
步骤501：根据各位置的测量误差函数
步骤502：根据初始静态误差参数向量β(s)和测量误差向量r(s)确定当前静态误差参数向量β(s+1)；具体公式为：
β(s+1)＝β(s)-(j)-1r(s)；
其中，β(s)为初始静态误差参数向量，r(s)为测量误差向量，(j)-1为雅可比矩阵的逆矩阵，
步骤60：根据当前静态误差参数向量β(s+1)和各位置的测量误差函数
步骤70：判断代价函数值s(β(s+1))是否小于设定阈值；如果代价函数值s(β(s+1))小于设定阈值，则将当前静态误差参数向量β(s+1)作为最优静态误差参数向量，输出最优静态误差参数向量；如果代价函数值s(β(s+1))大于等于设定阈值，则将当前静态误差参数向量β(s+1)作为初始静态误差参数向量β(s)，并执行步骤40。
本发明设计辅助装置3，先将三轴加速度传感器1放置在水平台面2上，获得六个位置；然后将三轴加速度传感器1放置在辅助装置3上，辅助装置3放置在水平台面2上，获得三个位置；获取三轴加速度传感器1在9个位置处于不同姿态的多组测量数据，根据9组不同测量数据实现全方面标定三轴加速度传感器1的静态误差，不仅提高了标定静态误差的精度，还摆脱了依托转台和精确水平面进行静态误差标定，所需设备仅为一水平台面2和本发明设计的标定辅助装置即可，降低标定系统的体积和成本。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。