发明名称 一种基于速度加比力匹配的传递对准时间延迟估计方法
摘要 本发明公开一种基于速度加比力匹配的传递对准时间延迟估计算法。其实现方法是:利用已对准的主惯导对子惯导进行粗对准,再利用子惯导相对主惯导的速度差、比力差作为滤波观测量,并将主惯导信息的时间延迟量扩充为滤波系统状态变量,结合惯导系统的误差模型,利用卡尔曼滤波算法估计出时间延迟以及子惯导相对于主惯导的平台误差角。本发明适用于主、子惯导均为平台式惯导系统时,延迟时间较短的速度加比力匹配传递对准。
申请公布号 CN103344251B 申请公布日期 2015.09.30
申请号 CN201310227454.8 申请日期 2013.06.08
申请人 哈尔滨工程大学 发明人 徐博;郭宇;池姗姗;陈春;肖永平;王文佳;田学林;金辰
分类号 G01C25/00(2006.01)I;G01C21/20(2006.01)I 主分类号 G01C25/00(2006.01)I
代理机构 代理人
主权项 一种基于速度加比力匹配的传递对准时间延迟估计算法,其特征是:(1)子惯导预热,并利用已对准的平台式主惯导系统对子惯导系统进行一步传递对准,完成子惯导系统的粗对准;(2)主、子惯导进行导航解算,并采集主惯导速度、比力信息传输至子惯导解算计算机;(3)在子惯导计算机中,构造卡尔曼滤波基本方程,并利用主、子惯导速度差、比力差进行标准离散卡尔曼滤波迭代解算;其中,主、子惯导均为平台式惯导系统,所选取的滤波系统状态变量为:X=[δv<sub>x</sub> δv<sub>y</sub> φ<sub>x</sub> φ<sub>y</sub> φ<sub>z</sub> ▽<sub>x</sub> ▽<sub>y</sub> ε<sub>x</sub> ε<sub>y</sub> ε<sub>z</sub> Δt]<sup>T</sup>式中,δv<sub>x</sub>、δv<sub>y</sub>分别为子惯导东向、北向速度误差;φ<sub>x</sub>、φ<sub>y</sub>、φ<sub>z</sub>分别为子惯导相对于主惯导的平台误差角;▽<sub>x</sub>、▽<sub>y</sub>分别为子惯导加速度计零偏;ε<sub>x</sub>、ε<sub>y</sub>、ε<sub>z</sub>分别为子惯导陀螺漂移;Δt为主惯导信息的时间延迟;结合子惯导系统误差微分方程,建立滤波状态方程;速度加比力匹配法的系统观测量为:Z=[δv<sub>x</sub> δv<sub>y</sub> δf<sub>x</sub> δf<sub>y</sub>]<sup>T</sup>式中,δf<sub>x</sub>、δf<sub>y</sub>分别为主、子惯导比力之差;滤波观测方程为:Z=HX+V<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>H</mi><mo>=</mo></mrow><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><mi>Dv</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><mi>Dv</mi><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><msub><mi>f</mi><mi>z</mi></msub></mtd><mtd><msub><mi>f</mi><mi>y</mi></msub></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><mi>Da</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><msub><mi>f</mi><mi>z</mi></msub></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><msub><mi>f</mi><mi>x</mi></msub></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><mi>Da</mi><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000697057580000011.GIF" wi="1149" he="309" /></maths>式中,<maths num="0002" id="cmaths0002"><math><![CDATA[<mi>Dv</mi><mrow><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mover><mi>V</mi><mo>&CenterDot;</mo></mover><mi>m</mi><mi>n</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><mi>Da</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mover><mi>f</mi><mo>&CenterDot;</mo></mover><mi>m</mi><mi>n</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000697057580000012.GIF" wi="649" he="83" /></maths>其中,<img file="FDA0000697057580000013.GIF" wi="298" he="85" />分别为主惯导速度与比力的微分;(4)卡尔曼滤波结束后,得到延迟时间估计值以及时间延迟补偿后的平台误差角。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室