一种惯性测量单元的无定向多位置高精度标定方法技术

一种惯性测量单元的无定向多位置高精度标定方法，本发明专利技术涉及一种利用精密转台精确标定惯性测量单元（ＩＭＵ）误差系数的方法。该方法可在精密转台不指北的情况下，通过ＩＭＵ在六个位置上的１２次旋转，精确标定出陀螺仪标度因数、陀螺仪常值漂移、陀螺仪安装误差、与ｇ有关误差项和加速度计标度因数、加速度计常值偏置、加速度计安装误差共７类３３个误差系数。本发明专利技术具有精度高且操作简单的特点，不仅提高了ＩＭＵ的使用精度，同时也大大提高了标定的效率。

【技术实现步骤摘要】
一种惯性测量单元的无定向多位置高精度标定方法
本专利技术涉及一种利用精密转台精确标定惯性测量单元(InertialMeasurement Unit，IMU)误差系数的方法，可用于标定挠性陀螺IMU、液浮陀螺IMU或MEMS陀螺IMU等。
技术介绍
捷联惯性导航系统(Strapdown Inertial Navigation System，SINS)是一种测量运动载体位置、速度和姿态的精密仪器，它将陀螺仪和加速度计固联在载体上，用计算机平台代替物理平台，具有结构简单、成本低、可靠性高等优点，已经成为了惯性导航技术的主要发展方向。惯性测量单元(InertialMeasurement Unit，IMU)由陀螺仪和加速度计及相关电路组成，是捷联惯导系统的核心，其误差包含确定性误差和随机误差两部分，其中确定性误差即系统误差约占总误差的90％左右，是捷联惯导系统最主要的误差源。因此，捷联惯导在使用前必须通过标定试验确定出IMU的各项误差系数，并在系统中进行补偿。传统的标定方法包括静态多位置试验方法和角速率试验方法两种。角速度率试验方法精度较高，但是只能标定出IMU角速度通道的标度因数和安装误差共9个误差系数。静态多位置试验方法的基本原理是利用转台提供的方位基准(即需要转台精确指北)和水平基准，将地球自转角速度和重力加速度作为输入IMU的标称量，并与IMU中陀螺仪和加速度计的输出进行比较，建立IMU的误差方程，然后将精密旋转基准(转台)转动到多个不同位置，当位置数与输入输出方程中的未知数即误差系数的个数相等时，即可通过联立的方程组求解出各项误差系数。静态多位置试验方法可标定出角速-->度通道和加速度通道的全部33个误差系数，但是其精度不如角速率试验的精度高，因此工程实践中经常利用角速率试验标定角速度通道的标度因数和安装误差，然后通过多位置试验标定其余各项误差系数，整个标定过程非常复杂繁琐，效率较低。尤其是静态多位置试验方法，需要对惯性测量单元进行精确定向，即确定其与北向的夹角。通常利用经纬仪等高精度测量仪器来实现惯性测量单元的精确定向，这个过程十分复杂，进一步降低了标定的效率。
技术实现思路
本专利技术的技术解决问题是：克服现有技术的不足，提供一种惯性测量单元的无定向多位置精确标定方法，该方法可在精密转台不指北的情况下，通过IMU在六个位置上的12次旋转，精确标定出陀螺仪标度因数、陀螺仪常值漂移、陀螺仪安装误差、与g有关误差项和加速度计标度因数、加速度计常值偏置、加速度计安装误差共7类33个误差系数。本专利技术具有精度高且操作简单的特点，不仅提高了IMU的使用精度，同时也大大提高了标定的效率。本专利技术的技术解决方案为：一种惯性测量单元的无定向多位置高精度标定方法，其特征在于包括以下步骤：(1)将IMU安装在三轴转台上，使IMU的Z轴向上，与转台的Z轴重合，IMU的X轴和Y轴处于水平面内；(2)使转台绕Z轴逆时针旋转360，然后保持IMU静止1-3分钟，记录旋转及静止过程中IMU输出的角速度和加速度；(3)使转台绕Z轴顺时针旋转360，然后保持IMU静止1-3分钟，记录旋转及静止过程中IMU输出的角速度和加速度；(4)旋转转台使IMU的Z轴向下，X轴和Y轴处于水平面内，然后重复试验步骤(2)和(3)；(5)旋转转台使IMU的X轴向上，Z轴和Y轴处于水平面内，重复试-->验步骤(2)和(3)；(6)旋转转台使IMU的X轴向下，Z轴和Y轴处于水平面内，重复试验步骤(2)和(3)；(7)旋转转台使IMU的X轴向上，Z轴和X轴处于水平面内，重复试验步骤(2)和(3)；(8)旋转转台使IMU的X轴向下，Z轴和X轴处于水平面内，重复试验步骤(2)和(3)；(9)利用六个位置上12次旋转的12组数据，根据IMU的误差模型，采用IMU误差系数计算公式标定出IMU的全部33个误差系数。本专利技术的原理是：当IMU的Z轴向上进行逆时针和顺时针旋转时，IMU的误差方程如下：SzPz1+t=(ωz1++ωez1)+Dz+Mzyωey1+Mzxωex1+Dzzg---(1)]]>SxPx1+t=ωex1+Dx+Mxyωey1+Mxz(ωz1++ωez1)+Dxzg---(2)]]>SyPy1+t=ωey1+Dy+Myxωex1+Myz(ωz1++ωez1)+Dyzg---(3)]]>SzPz1-t=(ωz1-+ωez1)+Dz+Mzyωey1+Mzxωex1+Dzzg---(4)]]>SxPx1-t=ωex1+Dx+Mxyωey1+Mxz(ωz1-+ωez1)+Dxzg---(5)]]>SyPy1-t=ωey1+Dy+Myxωex1+Myz(ωz1-+ωez1)+Dyzg---(6)]]>kzNz1t=g+Bz---(7)]]>kxNx1t=Ixzg+Bx---(8)]]>kyNy1t=Iyzg+By---(9)]]>其中，下标中的“1”表示第1位置，即Z轴向上的位置，上标“+”和“-”-->分别表示逆时针和顺时针旋转，ωilk(i＝x，y，z；k＝+，-)表示IMU相对地球的旋转角速度在i轴的投影，ωeilk(i＝x，y，z；k＝+，-)表示IMU系统在第1位置逆(顺)时针旋转时，地球自转角速度ωe在i轴的投影。g表示重力加速度。S为陀螺仪的标度因数，P为陀螺仪的输出，t为数据采集时间，M为陀螺仪的安装误差，ω为陀螺仪的输入角速度，D为陀螺仪的常值漂移，G为陀螺仪的与g有关误差系数，k为加速度计的标度因数，N为加速度计的输出，C为加速度计的安装误差，f为加速度计的输入加速度，B为加速度计的常值偏置。由于IMU顺时针和逆时针都旋转了360度，将方程(10)～(15)积分，则可得：SxPx1+=∫0t[Dx+Mxz(ωz1++ωez1)+Dxzg]dτ---(11)]]>SyPy1+=∫0t[Dy+Myz(ωz1++ωez1)+Dyzg]dτ---(12)]]>SxPx1-=∫0t[Dx+Mxz(ωz1-+ωez1)+Dxzg]dτ---(14)]]>SyPy1-=∫0t[Dy+Myz(ωz1-+ωez1)+Dyzg]dτ---(15)]]>对于第2～6位置(即Z轴向下、X轴向上、X轴向下、Y轴向上及Y轴向下5个位置)，根据上述方法，可建立每个位置上的7个方程，将第1～6位置建立方程组联立求解，可以得到IMU的误差系数的计算公式如式(16)～式(23)所示：Si＝1440°/(Pij+-Pij--Pi(j+1)+-Pi(j+1)-)    (16)Di＝Si(Pij++Pij-+Pi(j+1)++Pi(j+1)-)dτ/4t   (17)Mni＝Sn(Pnj+-Pnj本文档来自技高网...

【技术保护点】
一种惯性测量单元的无定向多位置高精度标定方法，其特征在于包括以下步骤：（１）将ＩＭＵ安装在三轴转台上，使ＩＭＵ的Ｚ轴向上，与转台的Ｚ轴重合，ＩＭＵ的Ｘ轴和Ｙ轴处于水平面内；（２）使转台绕Ｚ轴逆时针旋转３６０°，然后保持ＩＭＵ 静止１－３分钟，记录旋转及静止过程中ＩＭＵ输出的角速度和加速度；（３）使转台绕Ｚ轴顺时针旋转３６０°，然后保持ＩＭＵ静止１－３分钟，记录旋转及静止过程中ＩＭＵ输出的角速度和加速度；（４）旋转转台使ＩＭＵ的Ｚ轴向下，Ｘ轴和Ｙ轴 处于水平面内，然后重复试验步骤（２）和（３）；（５）旋转转台使ＩＭＵ的Ｘ轴向上，Ｚ轴和Ｙ轴处于水平面内，重复试验步骤（２）和（３）；（６）旋转转台使ＩＭＵ的Ｘ轴向下，Ｚ轴和Ｙ轴处于水平面内，重复试验步骤（２）和（３）；　 　（７）旋转转台使ＩＭＵ的Ｘ轴向上，Ｚ轴和Ｘ轴处于水平面内，重复试验步骤（２）和（３）；（８）旋转转台使ＩＭＵ的Ｘ轴向下，Ｚ轴和Ｘ轴处于水平面内，重复试验步骤（２）和（３）；（９）利用六个位置上１２次旋转的１２组数据，根据Ｉ ＭＵ的误差模型，采用ＩＭＵ误差系数计算公式标定出ＩＭＵ的全部３３个误差系数。...

【技术特征摘要】
1、一种惯性测量单元的无定向多位置高精度标定方法，其特征在于包括以下步骤：(1)将IMU安装在三轴转台上，使IMU的Z轴向上，与转台的Z轴重合，IMU的X轴和Y轴处于水平面内；(2)使转台绕Z轴逆时针旋转360°，然后保持IMU静止1-3分钟，记录旋转及静止过程中IMU输出的角速度和加速度；(3)使转台绕Z轴顺时针旋转360°，然后保持IMU静止1-3分钟，记录旋转及静止过程中IMU输出的角速度和加速度；(4)旋转转台使IMU的Z轴向下，X轴和Y轴处于水平面内，然后重复试验步骤(2)和(3)；(5)旋转转台使IMU的X轴向上，Z轴和Y轴处于水平面内，重复试验步骤(2)和(3)；(6)旋转转台使IMU的X轴向下，Z轴和Y轴处于水平面内，重复试验步骤(2)和(3)；(7)旋转转台使IMU的X轴向上，Z轴和X轴处于水平面内，重复试验步骤(2)和(3)；(8)旋转转台使IMU的X轴向下，Z轴和X轴处于水平面内，重复试验步骤(2)和(3)；(9)利用六个位置上12次旋转的12组数据，根据IMU的误差模型，采用IMU误差系数计算公式标定出IMU的全部33个误差系数。2、根据权利要求1所述的惯性测量单元的无定向多位置高精度标定方法，其特征在于：所述的步骤(9)中IMU的误差模型包括角速度通道误差模型和加速度通道误差模型，分别如下：SPt=Mω+D+Gf]]>kN＝Cf+B其中，S为陀螺仪的标度因数，P为陀螺仪的输出，t为数据采集时间，M为陀螺仪的安装误差，ω为陀螺仪的输入角速度，D为陀螺仪的常值漂移，G为陀螺仪的与g有关误差系数，g为重力加速度，k为加速度计的标度因数，N为加速度计的输出，C为加速度计的安装误差，f为加速度计的输入加速度，B为加...

【专利技术属性】
技术研发人员：房建成，刘百奇，张钰，张海鹏，杨胜，高翌春，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：11[中国|北京]