发明名称 一种利用车体非完整性约束的惯性自主导航方法
摘要 本发明公开了一种利用车体非完整性约束的惯性自主导航方法。本方法利用在运动过程中车体非完整性约束,通过平滑惯性自主导航系统输出的侧向和垂向速度累积,构建卡尔曼滤波观测方程,并利用卡尔曼滤波方程周期性计算出载车的导航参数误差并修正导航参数,从而实现了利用车体非完整性约束惯性导航系统的自主导航。本发明有效的解决了惯性导航系统自主导航精度发散的问题,有效的保证了惯性定位定向设备的精度。该方法使用方便,无需增加额外辅助设备。
申请公布号 CN103712620B 申请公布日期 2016.03.30
申请号 CN201310609301.X 申请日期 2013.11.27
申请人 北京机械设备研究所 发明人 高文劭
分类号 G01C21/16(2006.01)I;G01C21/20(2006.01)I 主分类号 G01C21/16(2006.01)I
代理机构 中国航天科工集团公司专利中心 11024 代理人 岳洁菱
主权项 一种利用车体非完整性约束的惯性自主导航方法,其特征在于具体步骤为:第一步 构建利用车体非完整性约束的惯性自主导航系统利用车体非完整性约束的惯性自主导航系统,包括:三轴硅微陀螺仪、三轴硅微加速度计、导航计算机、惯性导航解算模块、导航修正模块和通讯模块;三轴硅微陀螺仪和三轴硅微加速度计分别与导航计算机实现SPI总线连接;惯性导航解算模块、导航修正模块和通讯模块分别置于导航计算机中;惯性导航解算模块用于利用采集的角速率和加速度通过数学运算得到载体在地理坐标系内的姿态、速度和位置;导航修正模块用于利用车体非完整性约束构建观测方程,采用卡尔曼滤波方程周期性修正利用车体非完整性约束的惯性自主导航系统的导航参数,即载车的姿态、速度和位置;通讯模块用于利用车体非完整性约束的惯性自主导航系统内外的数据采集和通讯;第二步 惯性导航解算模块对导航参数初始化载车静止条件下,惯性导航解算模块接收来自三轴硅微陀螺仪和三轴硅微加速度计的数据,对三轴硅微陀螺仪和三轴硅微加速度计的零偏进行补偿,同时接收预存信息完成对惯性自主导航系统姿态、速度和位置信息的初始化;第三步 惯性导航解算模块进行导航解算载车启动后,惯性导航解算模块接收三轴硅微陀螺仪和三轴硅微加速度计的数据,计算出载车在地理坐标系内的姿态、速度和位置;第四步 导航修正模块进行修正准备导航修正模块接收导航解算模块输出的速度,并利用姿态转换矩阵将接收到的速度转换为车体的侧向速度和垂向速度;转换公式为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mover><mi>v</mi><mo>^</mo></mover><mi>b</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><msubsup><mover><mi>C</mi><mo>^</mo></mover><mi>b</mi><mi>n</mi></msubsup><mo>)</mo></mrow><mi>T</mi></msup><msub><mover><mi>v</mi><mo>^</mo></mover><mi>n</mi></msub><mo>;</mo></mrow>]]></math><img file="FDA0000880013500000011.GIF" wi="302" he="103" /></maths>( )<sup>T</sup>表示向量转置,<img file="FDA0000880013500000012.GIF" wi="66" he="70" />为三维向量,其第一项和第三项为侧向<img file="FDA0000880013500000013.GIF" wi="58" he="75" />和垂向速度<img file="FDA0000880013500000014.GIF" wi="78" he="78" /><img file="FDA0000880013500000015.GIF" wi="86" he="86" />为惯性自主导航系统的姿态转换矩阵;<img file="FDA0000880013500000016.GIF" wi="62" he="70" />为惯性导航解算模块解算的速度,为三维向量,表示为<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mover><mi>v</mi><mo>^</mo></mover><mi>n</mi></msub><mo>=</mo><msup><mfenced open = "[" close = "]"><mtable><mtr><mtd><msubsup><mi>v</mi><mi>E</mi><mi>n</mi></msubsup></mtd><mtd><msubsup><mi>v</mi><mi>N</mi><mi>n</mi></msubsup></mtd><mtd><msubsup><mi>v</mi><mi>U</mi><mi>n</mi></msubsup></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>;</mo></mrow>]]></math><img file="FDA0000880013500000017.GIF" wi="421" he="94" /></maths>导航修正模块同时计算载车运动情况下的系统状态转移矩阵;第五步 导航修正模块利用车体非完整性约束计算观测向量和观测矩阵车体在运动情况下,受非完整性约束的影响,其侧向和垂向速度为零,<img file="FDA0000880013500000018.GIF" wi="63" he="77" />和<img file="FDA0000880013500000019.GIF" wi="60" he="77" />就是车体的速度误差;将<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>Z</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msup><mfenced open = "[" close = "]"><mtable><mtr><mtd><msubsup><mi>v</mi><mi>x</mi><mi>b</mi></msubsup></mtd><mtd><msubsup><mi>v</mi><mi>z</mi><mi>b</mi></msubsup></mtd></mtr></mtable></mfenced><mi>T</mi></msup></mrow>]]></math><img file="FDA0000880013500000021.GIF" wi="343" he="91" /></maths>作为系统的观测向量;导航修正模块在修正时刻,计算从上一次修正时刻到本次修正时刻累积接收的车体的侧向速度和垂向速度之和,并计算这两项的平均值,将这两项作为<img file="FDA0000880013500000022.GIF" wi="55" he="79" />和<img file="FDA0000880013500000023.GIF" wi="78" he="77" />计算滤波时刻的观测矩阵H(k),<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>H</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><msub><mi>C</mi><mn>13</mn></msub><msubsup><mi>v</mi><mi>N</mi><mi>n</mi></msubsup><mo>-</mo><msub><mi>C</mi><mrow><mn>1</mn><mi>2</mi></mrow></msub><msubsup><mi>v</mi><mi>U</mi><mi>n</mi></msubsup></mrow></mtd><mtd><mrow><msub><mi>C</mi><mn>11</mn></msub><msubsup><mi>v</mi><mi>N</mi><mi>n</mi></msubsup><mo>-</mo><msub><mi>C</mi><mn>13</mn></msub><msubsup><mi>v</mi><mi>E</mi><mi>n</mi></msubsup></mrow></mtd><mtd><mrow><msub><mi>C</mi><mn>12</mn></msub><msubsup><mi>v</mi><mi>E</mi><mi>n</mi></msubsup><mo>-</mo><msub><mi>C</mi><mn>11</mn></msub><msubsup><mi>v</mi><mi>N</mi><mi>n</mi></msubsup></mrow></mtd><mtd><msub><mi>C</mi><mn>11</mn></msub></mtd><mtd><msub><mi>C</mi><mn>12</mn></msub></mtd><mtd><msub><mi>C</mi><mn>13</mn></msub></mtd><mtd><msub><mn>0</mn><mrow><mn>1</mn><mo>&times;</mo><mn>9</mn></mrow></msub></mtd></mtr><mtr><mtd><mrow><msub><mi>C</mi><mn>33</mn></msub><msubsup><mi>v</mi><mi>N</mi><mi>n</mi></msubsup><mo>-</mo><msub><mi>C</mi><mi>32</mi></msub><msubsup><mi>v</mi><mi>U</mi><mi>n</mi></msubsup></mrow></mtd><mtd><mrow><msub><mi>C</mi><mn>31</mn></msub><msubsup><mi>v</mi><mi>N</mi><mi>n</mi></msubsup><mo>-</mo><msub><mi>C</mi><mn>33</mn></msub><msubsup><mi>v</mi><mi>E</mi><mi>n</mi></msubsup></mrow></mtd><mtd><mrow><msub><mi>C</mi><mn>32</mn></msub><msubsup><mi>v</mi><mi>E</mi><mi>n</mi></msubsup><mo>-</mo><msub><mi>C</mi><mn>31</mn></msub><msubsup><mi>v</mi><mi>N</mi><mi>n</mi></msubsup></mrow></mtd><mtd><msub><mi>C</mi><mn>31</mn></msub></mtd><mtd><msub><mi>C</mi><mn>32</mn></msub></mtd><mtd><msub><mi>C</mi><mn>33</mn></msub></mtd><mtd><msub><mn>0</mn><mrow><mn>1</mn><mo>&times;</mo><mn>9</mn></mrow></msub></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000880013500000024.GIF" wi="1556" he="158" /></maths>其中C<sub>ij</sub>为计算姿态矩阵<img file="FDA0000880013500000025.GIF" wi="69" he="77" />的(i,j)项;<img file="FDA0000880013500000026.GIF" wi="67" he="71" />和<img file="FDA0000880013500000027.GIF" wi="56" he="70" />为惯性导航解算模块计算的速度分量的第一项和第三项;第六步 导航修正模块计算导航参数误差导航修正模块利用计算得到的运动误差方程、观测向量、观测矩阵,采用卡尔曼滤波方程计算得到利用车体非完整性约束的惯性自主导航系统的导航参数误差;第七步 惯性导航解算模块修正导航参数惯性导航解算模块接收导航修正模块传递的导航参数误差,修正利用车体非完整性约束的惯性自主导航系统的姿态、速度和位置误差,将导航参数误差和观测向量Z(k)值归零;经过上述各步骤实现了利用车体非完整性约束的惯性自主导航。
地址 100854 北京市海淀区北京142信箱208分箱