主权项 |
一种双惯导联合旋转调制导航与在线相对性能评估方法,其特征在于包括以下步骤:步骤一:坐标系定义,定义导航坐标系(n系)为当地水平地理坐标系,坐标轴分别指向北向‑东向‑地向(N‑E‑D),载体坐标系(b系)坐标轴分别沿载体的横滚轴‑俯仰轴‑偏航轴(前‑右‑下),惯导系统1、2的惯性测量单元(IMU)坐标系分别为s<sub>1</sub>、s<sub>2</sub>,坐标轴指向同载体坐标系定义;步骤二:联合旋转调制策略编排,惯导系统1、2分别绕方位轴进行4位置8次序旋转调制转位,转位规律不同;步骤三:联合系统状态方程确定,取惯导系统1与惯导系统2的姿态误差、速度误差、位置误差的状态差为联合系统状态,16个误差状态为:x=[(φ<sub>N1</sub>‑φ<sub>N2</sub>) (φ<sub>E1</sub>‑φ<sub>E2</sub>) (φ<sub>D1</sub>‑φ<sub>D2</sub>) (δv<sub>N1</sub>‑δv<sub>N2</sub>) (δv<sub>E1</sub>‑δv<sub>E2</sub>)(δL<sub>1</sub>‑δL<sub>2</sub>) (δλ<sub>1</sub>‑δλ<sub>2</sub>) ε<sub>x1</sub> ε<sub>y1</sub> ε<sub>x2</sub> ε<sub>y2</sub> (ε<sub>z1</sub>‑ε<sub>z2</sub>) ▽<sub>x1</sub> ▽<sub>y1</sub> ▽<sub>x2</sub> ▽<sub>y2</sub>]<sup>T</sup> (1)即x=[φ<sub>N12</sub> φ<sub>E12</sub> φ<sub>D12</sub> δv<sub>N12</sub> δv<sub>E12</sub> δL<sub>12</sub> δλ<sub>12</sub>ε<sub>x1</sub> ε<sub>y1</sub> ε<sub>x2</sub> ε<sub>y2</sub> ε<sub>z12</sub> ▽<sub>x1</sub> ▽<sub>y1</sub> ▽<sub>x2</sub> ▽<sub>y2</sub>]<sup>T</sup>其中,上标T表示向量或矩阵的转置,φ<sub>N12</sub>=(φ<sub>N1</sub>‑φ<sub>N2</sub>)、φ<sub>E12</sub>=(φ<sub>E1</sub>‑φ<sub>E2</sub>)、φ<sub>D12</sub>=(φ<sub>D1</sub>‑φ<sub>D2</sub>)分别为惯导系统1与惯导系统2的姿态误差向量的差值,δv<sub>N12</sub>=(δv<sub>N1</sub>‑δv<sub>N2</sub>)、δv<sub>E12</sub>=(δv<sub>E1</sub>‑δv<sub>E2</sub>)分别为惯导系统1与惯导系统2的北向、东向速度误差的差值,δL<sub>12</sub>=(δL<sub>1</sub>‑δL<sub>2</sub>)、δλ<sub>12</sub>=(δλ<sub>1</sub>‑δλ<sub>2</sub>)分别为惯导系统1与惯导系统2的纬度、经度误差的差值,ε<sub>x1</sub>、ε<sub>x2</sub>、ε<sub>y1</sub>、ε<sub>y2</sub>分别为惯导系统1与惯导系统2的IMU对应水平坐标轴的陀螺常值漂移,ε<sub>z12</sub>=(ε<sub>z1</sub>‑ε<sub>z2</sub>)为惯导系统1与惯导系统2的IMU对应天向坐标轴的陀螺常值漂移之差(绕方位轴单轴旋转时,水平加表零偏、陀螺常值漂移可分离,但天向的陀螺常值漂移不可分离),▽<sub>x1</sub>、▽<sub>x2</sub>、▽<sub>y1</sub>、▽<sub>y2</sub>分别为惯导系统1与惯导系统2的IMU对应水平坐标轴的加表常值零偏,根据联合系统状态确定联合系统状态方程;步骤四:量测方程确定,扣除安装位置误差以及杆臂造成的速度误差,在缺少外界参考信息条件下,取1s更新一次的惯导系统1、2的速度和位置之差为相应的量测量,量测量为:z(t)=[δv<sub>N1</sub>‑δv<sub>N2</sub> δv<sub>E1</sub>‑δv<sub>E2</sub> δL<sub>1</sub>‑δL<sub>2</sub> δλ<sub>1</sub>‑δλ<sub>2</sub>]<sup>T</sup> (2)其中,<img file="FDA0000753676810000011.GIF" wi="395" he="70" />为惯导系统1、2的北向速度差值,<img file="FDA0000753676810000012.GIF" wi="387" he="70" />为惯导系统1、2的东向速度差值,<img file="FDA0000753676810000013.GIF" wi="330" he="78" />为惯导系统1、2的纬度差值,<img file="FDA0000753676810000014.GIF" wi="317" he="73" />为惯导系统1、2的经度差值,根据量测状态确定量测方程;步骤五:卡尔曼滤波,根据联合系统状态方程、量测方程构建卡尔曼滤波器,1s进行一次量测更新,对两套惯导系统各自的陀螺、加表零偏进行估计;步骤六:惯导性能在线评估,根据陀螺、加表零偏的估计值对惯导系统1、2的相对性能进行评估,零偏小的系统作为优选系统。 |