发明名称 一种基于近距离地标测距的组合导航方法
摘要 本发明公开了一种基于近距离地标测距的组合导航方法,方法包括以下步骤:S1、将捷联惯导系统和两个测距装置安装在载体上,采集捷联惯导系统惯性测量单元正常工作状态下输出的载体运动信息以进行惯导解算;S2、载体行驶至预设地标区域内,采集两个测距装置到地标点的距离以及地标点的位置信息;S3、建立状态方程;S4、建立量测方程,将两个测距装置测量距离的平方差作为量测信息;S5、利用状态方程和量测方程,进行卡尔曼滤波,以实时修正惯性导航系统参数误差和器件参数误差,实现组合导航。实施本发明方法可有效避免信号干扰,具备较强的隐蔽性和自主性,且消除了误差累积的问题。
申请公布号 CN104864868A 申请公布日期 2015.08.26
申请号 CN201510288879.9 申请日期 2015.05.29
申请人 湖北航天技术研究院总体设计所 发明人 刘明;胡华峰;周海;李旦;罗伟;杨元侠
分类号 G01C21/16(2006.01)I;G01C21/20(2006.01)I 主分类号 G01C21/16(2006.01)I
代理机构 武汉东喻专利代理事务所(普通合伙) 42224 代理人 宋业斌
主权项 一种基于近距离地标测距的组合导航方法,其特征在于,所述方法包括步骤:S1、将捷联惯导系统和两个测距装置安装在载体上,采集捷联惯导系统惯性测量单元正常工作状态下输出的载体运动信息并进行惯导解算;S2、载体行驶至预设地标区域内,采集两个测距装置到地标点的距离以及地标点的位置信息;S3、建立状态方程<img file="FDA0000728026490000011.GIF" wi="263" he="62" />其中,F为状态转移矩阵,w为系统噪声;状态向量<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>x</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mi>INSSys</mi></msub></mtd></mtr><mtr><mtd><msub><mi>x</mi><mi>INSSens</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000728026490000012.GIF" wi="297" he="160" /></maths>所述x<sub>INSSys</sub>为惯性导航系统参数误差状态向量,所述x<sub>INSSens</sub>为加速度计和陀螺参数误差状态向量;S4、建立量测方程z=Hx+v,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>z</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>z</mi><mrow><mi>Obs</mi><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>z</mi><mrow><mi>Obs</mi><mn>2</mn></mrow></msub></mtd></mtr></mtable></mfenced><mo>,</mo><mi>H</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>H</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>H</mi><mn>2</mn></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000728026490000013.GIF" wi="537" he="163" /></maths>z<sub>Obs1</sub>对应的量测矩阵<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>H</mi><mn>1</mn></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>2</mn><msup><mrow><mo>(</mo><msubsup><mi>C</mi><mi>B</mi><mi>N</mi></msubsup><mrow><mo>(</mo><msubsup><mi>l</mi><mn>1</mn><mi>B</mi></msubsup><mo>-</mo><msubsup><mi>l</mi><mn>2</mn><mi>B</mi></msubsup><mo>)</mo></mrow><mo>)</mo></mrow><mi>T</mi></msup><msubsup><mi>D</mi><mi>&delta;&theta;h</mi><mi>&delta;R</mi></msubsup></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000728026490000014.GIF" wi="689" he="132" /></maths>z<sub>Obs2</sub>对应的量测矩阵H<sub>2</sub>=[0 1 0];所述<img file="FDA0000728026490000015.GIF" wi="68" he="83" />为惯导解算的姿态矩阵,所述<img file="FDA0000728026490000016.GIF" wi="140" he="84" />分别为两个测距装置在体坐标系下的投影<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>D</mi><mi>&delta;&theta;h</mi><mi>&delta;R</mi></msubsup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mi>R</mi></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mo>-</mo><mi>R</mi></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000728026490000017.GIF" wi="446" he="235" /></maths>R为地球半径;观测量<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>z</mi><mrow><mi>Obs</mi><mn>1</mn></mrow></msub><mo>=</mo><mi>&Delta;</mi><msubsup><mover><mi>&rho;</mi><mo>^</mo></mover><mi>INS</mi><mn>2</mn></msubsup><mo>-</mo><mi>&Delta;</mi><msubsup><mover><mi>&rho;</mi><mo>~</mo></mover><mi>Distmt</mi><mn>2</mn></msubsup><mo>,</mo></mrow>]]></math><img file="FDA0000728026490000018.GIF" wi="476" he="78" /></maths><img file="FDA0000728026490000019.GIF" wi="131" he="78" />为惯导计算的两个测距装置到地标点距离的平方差,所述<img file="FDA00007280264900000110.GIF" wi="154" he="78" />为所述两个测距装置到地标点测量距离的平方差;观测量<img file="FDA00007280264900000111.GIF" wi="380" he="87" />所述<img file="FDA00007280264900000112.GIF" wi="94" he="98" />为惯导计算的载体高度,所述<img file="FDA00007280264900000113.GIF" wi="96" he="87" />为对应地标点的实际高度;v为量测噪声;S5、利用所述状态方程和量测方程,进行卡尔曼滤波,以实时修正惯性导航系统参数误差和器件参数误差,完成组合导航。
地址 430040 湖北省武汉市金山大道9号