主权项 |
一种基于近距离地标测距的组合导航方法,其特征在于,所述方法包括步骤: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>δθh</mi><mi>δ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>δθh</mi><mi>δ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>Δ</mi><msubsup><mover><mi>ρ</mi><mo>^</mo></mover><mi>INS</mi><mn>2</mn></msubsup><mo>-</mo><mi>Δ</mi><msubsup><mover><mi>ρ</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、利用所述状态方程和量测方程,进行卡尔曼滤波,以实时修正惯性导航系统参数误差和器件参数误差,完成组合导航。 |