发明名称 一种基于横坐标系的极区航行船舶捷联惯性导航系统重调方法
摘要 本发明的目的在于提供一种基于横坐标系的极区航行船舶捷联惯性导航系统重调方法,在t<sub>n</sub>时刻,获得横经度位置误差和横纬度位置误差,并依此对系统位置误差进行重调,获得经纬度位置误差和航向误差与t<sub>n</sub>时刻平台漂移角之间的关系;在t<sub>n+1</sub>时刻再次观测,获取第二点的横经度位置误差和横纬度位置误差,并对系统位置误差再次进行重调,得到t<sub>n+1</sub>时刻位置航向误差与陀螺漂移之间的关系;在t<sub>n+2</sub>时刻第三次观测,获得第三点的横经度位置误差和横纬度位置误差,得到t<sub>n+2</sub>时刻位置航向误差与陀螺漂移之间的关系;计算得到t<sub>n+2</sub>时刻陀螺漂移和航向误差值,并将它们补偿到系统中,完成系统重调。本发明可有效解除传统综合校正方法中对舰船速度和航行方式的限制。
申请公布号 CN103940446A 申请公布日期 2014.07.23
申请号 CN201410143315.1 申请日期 2014.04.10
申请人 哈尔滨工程大学 发明人 奔粤阳;尹冬寒;李倩;赵维珩;李明宇
分类号 G01C25/00(2006.01)I 主分类号 G01C25/00(2006.01)I
代理机构 代理人
主权项 一种基于横坐标系的极区航行船舶捷联惯性导航系统重调方法,其特征是:(1)在第一个位置点t<sub>n</sub>时刻,将外部测量设备测得的位置信息与惯性导航系统的位置信息输出值作差,得到第一点的横经度位置误差<img file="FDA0000489277520000011.GIF" wi="149" he="84" />和横纬度位置误差<img file="FDA0000489277520000012.GIF" wi="148" he="83" />后,并依此对系统位置误差进行重调,应用横坐标系下的P方程计算得到重调后第一位置点横经纬度位置航向误差与t<sub>n</sub>时刻平台漂移角之间的关系式:P(t<sub>n</sub>)=T(t<sub>n</sub>)ψ<sup>i</sup>(t<sub>n</sub>)其中i表示地心惯性坐标系,P(t<sub>n</sub>)为t<sub>n</sub>时刻的横经纬度位置误差和航向误差矩阵,ψ<sup>i</sup>(t<sub>n</sub>)为t<sub>n</sub>时刻地心惯性坐标系下的平台漂移角,T(t<sub>n</sub>)为P(t<sub>n</sub>)和ψ<sup>i</sup>(t<sub>n</sub>)二者之间的坐标转换矩阵;(2)在第二个位置点t<sub>n+1</sub>时刻再次观测,将外部测量设备测得的位置信息与惯性导航系统的位置信息输出值作差,得到第二点的横经度位置误差<img file="FDA0000489277520000013.GIF" wi="175" he="81" />和横纬度位置误差<img file="FDA0000489277520000014.GIF" wi="204" he="84" />并依此对系统位置误差进行重调,结合步骤(1),得到t<sub>n+1</sub>时刻位置航向误差与陀螺漂移之间的关系P(t<sub>n+1</sub>)=T(t<sub>n+1</sub>)([T(t<sub>n</sub>)]<sup>‑1</sup>·P(t<sub>n</sub>)+U(t<sub>n+1</sub>‑t<sub>n</sub>)·ε<sup>b</sup>)其中,<img file="FDA0000489277520000015.GIF" wi="595" he="139" />为矩阵<img file="FDA0000489277520000016.GIF" wi="58" he="74" />在时刻t<sub>n</sub>到t<sub>n+1</sub>的积分,这里惯性坐标系下的ψ方程为平台漂移角和陀螺漂移的关系方程,其表达式为<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msup><mi>&psi;</mi><mi>i</mi></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msup><mi>&psi;</mi><mi>i</mi></msup><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mi>t</mi></msubsup><msubsup><mi>C</mi><mi>b</mi><mi>i</mi></msubsup><mi>dt</mi><mo>)</mo></mrow><msup><mi>&epsiv;</mi><mi>b</mi></msup></mrow>]]></math><img file="FDA0000489277520000017.GIF" wi="562" he="107" /></maths>其中ψ<sup>i</sup>(t)表示t时刻惯性坐标系下的平台漂移角,ψ<sup>i</sup>(0)为初始平台漂移角,<img file="FDA0000489277520000018.GIF" wi="68" he="75" />为载体坐标系与惯性坐标系的转换矩阵,ε<sup>b</sup>为载体坐标系下的陀螺漂移;(3)在第三个位置点t<sub>n+2</sub>时刻进行第三次观测,由外部测量设备测得的位置信息与惯性导航系统的位置信息输出值作差,得到第三点的横经度位置误差<img file="FDA0000489277520000019.GIF" wi="181" he="85" />和横纬度位置误差<img file="FDA00004892775200000110.GIF" wi="210" he="84" />由步骤(1)和步骤(2)得到的结果,得到t<sub>n+2</sub>时刻位置航向误差与陀螺漂移之间的关系P(t<sub>n+2</sub>)=T(t<sub>n+2</sub>)([T(t<sub>n+1</sub>)]<sup>‑1</sup>·P(t<sub>n+1</sub>)+U(t<sub>n+2</sub>‑t<sub>n+1</sub>)·ε<sup>b</sup>)(4)由步骤(2)得到的t<sub>n+1</sub>时刻横经纬度位置航向误差与陀螺漂移之间的关系以及由步骤(3)得到的t<sub>n+2</sub>时刻横经纬度位置航向误差与陀螺漂移之间的关系,得到t<sub>n+2</sub>时刻系统的陀螺漂移及航向误差:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open='[' close=']'><mtable><mtr><mtd><mi>&delta;K</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>n</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msubsup><mi>&epsiv;</mi><mi>x</mi><mi>b</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>&epsiv;</mi><mi>y</mi><mi>b</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>&epsiv;</mi><mi>z</mi><mi>b</mi></msubsup></mtd></mtr></mtable></mfenced><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>l</mi><mn>11</mn></msub></mtd><mtd><msub><mi>l</mi><mn>12</mn></msub></mtd><mtd><msub><mi>l</mi><mn>13</mn></msub></mtd><mtd><msub><mi>l</mi><mn>14</mn></msub></mtd></mtr><mtr><mtd><msub><mi>l</mi><mn>21</mn></msub></mtd><mtd><msub><mi>l</mi><mn>22</mn></msub></mtd><mtd><msub><mi>l</mi><mn>23</mn></msub></mtd><mtd><msub><mi>l</mi><mn>24</mn></msub></mtd></mtr><mtr><mtd><msub><mi>m</mi><mn>11</mn></msub><msub><mi>l</mi><mn>31</mn></msub></mtd><mtd><msub><mi>m</mi><mn>11</mn></msub><msub><mi>l</mi><mn>32</mn></msub><mo>+</mo><msub><mi>m</mi><mn>12</mn></msub></mtd><mtd><msub><mi>m</mi><mn>11</mn></msub><msub><mi>l</mi><mn>33</mn></msub><mo>+</mo><msub><mi>m</mi><mn>13</mn></msub></mtd><mtd><msub><mi>m</mi><mn>11</mn></msub><msub><mi>l</mi><mn>34</mn></msub><mo>+</mo><msub><mi>m</mi><mn>14</mn></msub></mtd></mtr><mtr><mtd><msub><mi>m</mi><mn>21</mn></msub><msub><mi>l</mi><mn>31</mn></msub></mtd><mtd><msub><mi>m</mi><mn>21</mn></msub><msub><mi>l</mi><mn>32</mn></msub><mo>+</mo><msub><mi>m</mi><mn>22</mn></msub></mtd><mtd><msub><mi>m</mi><mn>21</mn></msub><msub><mi>l</mi><mn>33</mn></msub><mo>+</mo><msub><mi>m</mi><mn>23</mn></msub></mtd><mtd><msub><mi>m</mi><mn>21</mn></msub><msub><mi>l</mi><mn>34</mn></msub><mo>+</mo><msub><mi>m</mi><mn>24</mn></msub></mtd></mtr></mtable></mfenced><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>&CenterDot;</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>&delta;</mi><mover><mi>&lambda;</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><msub><mi>t</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>&delta;</mi><mover><mi>L</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><msub><mi>t</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>&delta;</mi><mover><mi>&lambda;</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><msub><mi>t</mi><mrow><mi>n</mi><mo>+</mo><mn>2</mn></mrow></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>&delta;</mi><mover><mi>L</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><msub><mi>t</mi><mrow><mi>n</mi><mo>+</mo><mn>2</mn></mrow></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000489277520000021.GIF" wi="1661" he="356" /></maths>其中δK(t<sub>n</sub>)为t<sub>n</sub>时刻的航向误差;<img file="FDA0000489277520000022.GIF" wi="245" he="84" />为陀螺漂移;l<sub>ij</sub>,(i=1,2,3,4;j=1,2,3,4)为步骤(2)得到的位置航向误差与陀螺漂移之间的转换矩阵的元素;m<sub>uv</sub>,(u=1,2,3,4;v=1,2,3,4)为步骤三得到的位置航向误差与陀螺漂移之间的转换矩阵的元素;将得到的横坐标系下的航向误差δK(t<sub>n</sub>)以及陀螺漂移<img file="FDA0000489277520000023.GIF" wi="238" he="84" />补偿到惯性导航系统中,即完成横坐标系下的捷联惯性导航系统的系统重调。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室