发明名称 一种双轴旋转惯导系统多位置自主标定方法
摘要 本发明提供一种双轴旋转惯导系统多位置自主标定方法,适用于双轴旋转惯导系统,避免系统组件定期拆卸,能够提高长航时导航自主性的在线标定方法。包括以下步骤:步骤1:建立双轴旋转惯导系统的器件误差模型和导航误差方程;步骤2、预热陀螺仪和加速度计组件,进行基于卡尔曼滤波的单位置精对准;步骤3、根据精对准结果转动载体,调整载体位置到近似与导航坐标系重合;步骤4、按照十位置转位方法转动环架,在每个位置采集导航解算得到的速度误差并计算姿态误差变化量,得到观测量;步骤5、根据每个位置的导航结果利用最小二乘计算出需要标定的误差。
申请公布号 CN104165638B 申请公布日期 2017.02.08
申请号 CN201410386915.0 申请日期 2014.08.07
申请人 北京理工大学 发明人 付梦印;邓志红;王博;任乾
分类号 G01C25/00(2006.01)I 主分类号 G01C25/00(2006.01)I
代理机构 北京理工大学专利中心 11120 代理人 仇蕾安
主权项 一种双轴旋转惯导系统多位置自主标定方法,其特征在于,包括以下步骤:步骤1:建立双轴旋转惯导系统的器件误差模型和导航误差方程;步骤2、预热陀螺仪和加速度计组件,进行基于卡尔曼滤波的单位置精对准;步骤3、根据精对准结果转动载体,调整载体位置到近似与导航坐标系重合;步骤4、按照十位置转位方法转动环架,在每个位置采集导航解算得到的速度误差并计算姿态误差变化量,得到观测量;步骤5、根据每个位置的导航结果利用最小二乘计算出需要标定的误差;所述步骤一中所述的建立器件误差模型包括以下步骤:第1步:建立双轴旋转惯导系统的器件误差模型;其中:加速度计的误差模型为:<img file="FDA0001070182790000011.GIF" wi="998" he="72" />其中<img file="FDA0001070182790000012.GIF" wi="59" he="62" />为加速度计输出误差,<img file="FDA0001070182790000013.GIF" wi="590" he="231" />为从IMU坐标系即s系到加速度计坐标系的转换矩阵,S<sub>aij</sub>为i和j方向上加速度计的安装误差角,S<sub>aij</sub>中的i=x、y、z,j=x、y、z且i≠j,K<sub>ax</sub>、K<sub>ay</sub>和K<sub>az</sub>分别为x、y和z方向上的加速度计刻度因数误差,f<sup>s</sup>为速度计的输入比力,<img file="FDA0001070182790000014.GIF" wi="398" he="78" /><img file="FDA0001070182790000015.GIF" wi="166" he="63" />和<img file="FDA0001070182790000016.GIF" wi="50" he="46" />为x、y和z方向上的加速度计零偏;陀螺仪的误差模型为:<img file="FDA0001070182790000017.GIF" wi="1005" he="70" />其中ε<sub>n</sub>为陀螺仪输出误差,<img file="FDA0001070182790000018.GIF" wi="612" he="230" />为从IMU坐标系到陀螺仪坐标系的转换矩阵,S<sub>gij</sub>为为i和j方向上陀螺仪的安装误差角,S<sub>gij</sub>中的<img file="FDA0001070182790000019.GIF" wi="470" he="47" />且i≠j,K<sub>gx</sub>、K<sub>gy</sub>.和K<sub>gz</sub>分别为x、y和z方向上的陀螺仪刻度因数误差,ω<sub>s</sub>为输入角速度,ε=[ε<sub>x</sub> ε<sub>y</sub> ε<sub>z</sub>]<sup>T</sup>为陀螺仪零偏;第2步:建立双轴旋转惯导系统的导航误差方程:<img file="FDA0001070182790000021.GIF" wi="1414" he="498" />其中δV<sub>E</sub>,δV<sub>N</sub>,δV<sub>U</sub>分别为东向、北向和天向速度误差,Δφ<sub>E</sub>、Δφ<sub>N</sub>和Δφ<sub>U</sub>分别为系统东向、北向和天向的失准角变化量,<img file="FDA0001070182790000022.GIF" wi="197" he="55" />和<img file="FDA0001070182790000023.GIF" wi="62" he="55" />分别为<img file="FDA0001070182790000024.GIF" wi="66" he="61" />东向、北向和天向的分量,ε<sub>E</sub>、ε<sub>N</sub>和ε<sub>U</sub>分别为ε<sub>n</sub>东向、北向和天向的分量,g为重力加速度,L为当地纬度,ω<sub>ie</sub>为地球自转角速度,φ<sub>E0</sub>、φ<sub>N0</sub>、φ<sub>U0</sub>为系统东向、北向和天向的初始失准角;第3步:由于自标定过程需要单位置精对准,单位置精对准后,采用陀螺仪和加速度计的误差模型参数表示系统的初始失准角φ<sub>E0</sub>、φ<sub>N0</sub>、φ<sub>U0</sub>,则有:<img file="FDA0001070182790000025.GIF" wi="1430" he="444" />其中<img file="FDA0001070182790000026.GIF" wi="427" he="62" />和ε<sub>0</sub>=[ε<sub>E0</sub> ε<sub>N0</sub> ε<sub>U0</sub>]<sup>T</sup>分别为初始对准位置处等效的加速度计和陀螺仪误差,即:<img file="FDA0001070182790000027.GIF" wi="1446" he="183" /><img file="FDA0001070182790000028.GIF" wi="1462" he="206" />第4步:提取观测量;在双轴旋转惯导系统中,有如下关系成立:<img file="FDA0001070182790000029.GIF" wi="1220" he="79" />其中,<img file="FDA00010701827900000210.GIF" wi="67" he="63" />为从外环架坐标系到载体坐标系的转换矩阵,<img file="FDA00010701827900000211.GIF" wi="76" he="70" />为从IMU坐标系到内环架坐标系的转换矩阵,<img file="FDA00010701827900000212.GIF" wi="62" he="55" />为实时导航解算出的姿态矩阵,<img file="FDA00010701827900000213.GIF" wi="61" he="71" />为单位置精对准 结束时刻的姿态矩阵,Δφ×为姿态误差变化量Δφ=[Δφ<sub>E</sub> Δφ<sub>N</sub> Δφ<sub>U</sub>]<sup>T</sup>构成的反对称矩阵,因此有<img file="FDA0001070182790000031.GIF" wi="1118" he="403" />对式(3)中的姿态误差方程进行拉普拉斯变换并移项,可得<img file="FDA0001070182790000032.GIF" wi="1398" he="134" />其中<img file="FDA0001070182790000033.GIF" wi="796" he="222" />对式(9)进行拉普拉斯反变换可得Δφ(t)=B(t)*φ<sub>0</sub>‑B(t)*ε<sup>n</sup>(t)      (10)其中*代表卷积,<img file="FDA0001070182790000034.GIF" wi="865" he="221" />对速度误差方程进行积分可得<img file="FDA0001070182790000035.GIF" wi="1278" he="471" />
地址 100081 北京市海淀区中关村南大街5号