发明名称 一种适用于极区的惯导系统综合校正方法
摘要 本发明公开了能够实现极区航行运载器在无阻尼条件下进行综合校正的一种适用于极区的惯导系统综合校正方法。本方法通过以横地理坐标系作为导航坐标系,在极区范围内通过横坐标系转换,将导航参数转换到横地理坐标系下。以横地理坐标系下的惯导输出误差量作为状态量,速度作为外部观测量,建立卡尔曼滤波模型,估测出水平误差角并引入到传统两点校正算法中,在无阻尼状态条件下完成极区惯导系统综合校正过程。本发明在无阻尼条件下实现极区航行运载器综合校正,解除了传统综合校正方案对运载器航行方式的限制,提高了极区惯导系统工作的灵活性。
申请公布号 CN103983277A 申请公布日期 2014.08.13
申请号 CN201410206180.9 申请日期 2014.05.16
申请人 哈尔滨工程大学 发明人 孙枫;尹冬寒;赵立莎;奔粤阳;李敬春;杨晓龙;赵维珩;杨祥龙;李明宇;阮双双
分类号 G01C25/00(2006.01)I 主分类号 G01C25/00(2006.01)I
代理机构 代理人
主权项 一种适用于极区的惯导系统综合校正方法,其特征在于,包括以下几个步骤:步骤一:将船用惯导系统输出的导航参数量转换到横地理坐标系下,导航参数包括:经纬度值、航向、速度和姿态信息;地理坐标系与横地理坐标系的转换矩阵为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>C</mi><mi>n</mi><mover><mi>n</mi><mo>&OverBar;</mo></mover></msubsup><mo>=</mo><msubsup><mi>C</mi><mi>i</mi><mover><mi>n</mi><mo>&OverBar;</mo></mover></msubsup><mo>&CenterDot;</mo><msubsup><mi>C</mi><mi>n</mi><mi>i</mi></msubsup></mrow>]]></math><img file="FDA0000506079740000011.GIF" wi="248" he="77" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>C</mi><mi>i</mi><mover><mi>n</mi><mo>&OverBar;</mo></mover></msubsup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>cos</mi><mover><mi>&lambda;</mi><mo>&OverBar;</mo></mover><mi>cos</mi><mi>&Omega;t</mi></mtd><mtd><mi>cos</mi><mover><mi>&lambda;</mi><mo>&OverBar;</mo></mover><mi>sin</mi><mi>&Omega;t</mi></mtd><mtd><mo>-</mo><mi>sin</mi><mover><mi>&lambda;</mi><mo>-</mo></mover></mtd></mtr><mtr><mtd><mo>-</mo><mi>sin</mi><mover><mi>L</mi><mo>&OverBar;</mo></mover><mi>sin</mi><mover><mi>&lambda;</mi><mo>&OverBar;</mo></mover><mi>cos</mi><mi>&Omega;t</mi><mo>-</mo><mi>cos</mi><mover><mi>L</mi><mo>&OverBar;</mo></mover><mi>sin</mi><mi>&Omega;t</mi></mtd><mtd><mi>cos</mi><mover><mi>L</mi><mo>&OverBar;</mo></mover><mi>cos</mi><mi>&Omega;t</mi><mo>-</mo><mi>sin</mi><mover><mi>L</mi><mo>&OverBar;</mo></mover><mi>sin</mi><mover><mi>&lambda;</mi><mo>&OverBar;</mo></mover><mi>sin</mi><mi>&Omega;t</mi></mtd><mtd><mo>-</mo><mi>sin</mi><mover><mi>L</mi><mo>&OverBar;</mo></mover><mi>cos</mi><mover><mi>&lambda;</mi><mo>&OverBar;</mo></mover></mtd></mtr><mtr><mtd><mi>cos</mi><mover><mi>L</mi><mo>&OverBar;</mo></mover><mi>sin</mi><mover><mi>&lambda;</mi><mo>&OverBar;</mo></mover><mi>cos</mi><mi>&Omega;t</mi><mo>-</mo><mi>sin</mi><mover><mi>L</mi><mo>&OverBar;</mo></mover><mi>sin</mi><mi>&Omega;t</mi></mtd><mtd><mi>cos</mi><mover><mi>L</mi><mo>&OverBar;</mo></mover><mi>sin</mi><mover><mi>&lambda;</mi><mo>&OverBar;</mo></mover><mi>sin</mi><mi>&Omega;t</mi><mo>+</mo><mi>sin</mi><mover><mi>L</mi><mo>&OverBar;</mo></mover><mi>cos</mi><mi>&Omega;t</mi></mtd><mtd><mi>cos</mi><mover><mi>L</mi><mo>-</mo></mover><mi>cos</mi><mover><mi>&lambda;</mi><mo>&OverBar;</mo></mover></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000506079740000012.GIF" wi="1782" he="243" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>C</mi><mi>n</mi><mi>i</mi></msubsup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mo>-</mo><mi>sin</mi><mrow><mo>(</mo><mi>&lambda;</mi><mo>+</mo><mi>&Omega;t</mi><mo>)</mo></mrow></mtd><mtd><mo>-</mo><mi>sin</mi><mi>L</mi><mi>cos</mi><mrow><mo>(</mo><mi>&lambda;</mi><mo>+</mo><mi>&Omega;t</mi><mo>)</mo></mrow></mtd><mtd><mi>cos</mi><mi>L</mi><mi>cos</mi><mrow><mo>(</mo><mi>&lambda;</mi><mo>+</mo><mi>&Omega;t</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>cos</mi><mrow><mo>(</mo><mi>&lambda;</mi><mo>+</mo><mi>&Omega;t</mi><mo>)</mo></mrow></mtd><mtd><mo>-</mo><mi>sin</mi><mi>L</mi><mi>sin</mi><mrow><mo>(</mo><mi>&lambda;</mi><mo>+</mo><mi>&Omega;t</mi><mo>)</mo></mrow></mtd><mtd><mi>cos</mi><mi>L</mi><mi>sin</mi><mrow><mo>(</mo><mi>&lambda;</mi><mo>+</mo><mi>&Omega;t</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>cos</mi><mi>L</mi></mtd><mtd><mi>sin</mi><mi>L</mi></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000506079740000013.GIF" wi="1223" he="232" /></maths><img file="FDA0000506079740000014.GIF" wi="71" he="76" />表示惯性坐标系到横地理坐标系的转换矩阵,<img file="FDA0000506079740000015.GIF" wi="64" he="78" />表示地理坐标系到惯性坐标系的转换矩阵,Ω为地球自转角速度,t为时间,λ为地理坐标系下的经度值,横地理坐标系下的经纬度值为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mover><mi>L</mi><mo>&OverBar;</mo></mover><mo>=</mo><mi>arcsin</mi><mrow><mo>(</mo><mi>cos</mi><mi>L</mi><mi>sin</mi><mi>&lambda;</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000506079740000016.GIF" wi="473" he="82" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mover><mi>&lambda;</mi><mo>&OverBar;</mo></mover><mo>=</mo><mi>arctan</mi><mrow><mo>(</mo><mfrac><mrow><mi>cos</mi><mi>L</mi><mi>cos</mi><mi>&lambda;</mi></mrow><mrow><mi>sin</mi><mi>L</mi></mrow></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000506079740000017.GIF" wi="502" he="140" /></maths>L为地理坐标系下的纬度值,<img file="FDA0000506079740000018.GIF" wi="48" he="67" />为横地理坐标系下的纬度值,<img file="FDA0000506079740000019.GIF" wi="50" he="66" />为横地理坐标系下的经度值,横地理坐标系下航向与速度为:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mover><mi>K</mi><mo>&OverBar;</mo></mover><mo>=</mo><mi>K</mi><mo>+</mo><mi>arctan</mi><mfrac><mrow><mi>sin</mi><mi>&lambda;</mi></mrow><mrow><mi>sin</mi><mi>L</mi><mi>cos</mi><mi>&lambda;</mi></mrow></mfrac><mo>-</mo><mi>&pi;</mi></mrow>]]></math><img file="FDA00005060797400000110.GIF" wi="615" he="128" /></maths><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mover><mi>V</mi><mo>&OverBar;</mo></mover><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mover><mi>V</mi><mo>&OverBar;</mo></mover><mi>x</mi></msub></mtd></mtr><mtr><mtd><msub><mover><mi>V</mi><mo>&OverBar;</mo></mover><mi>y</mi></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>V</mi><mi>sin</mi><mover><mi>K</mi><mo>&OverBar;</mo></mover></mtd></mtr><mtr><mtd><mi>V</mi><mi>cos</mi><mover><mi>K</mi><mo>&OverBar;</mo></mover></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA00005060797400000111.GIF" wi="444" he="169" /></maths>K表示地理坐标系下的航向,<img file="FDA00005060797400000112.GIF" wi="57" he="67" />为横地理坐标系下的航向,π为圆周率,V为地理坐标系下的速度,<img file="FDA00005060797400000113.GIF" wi="49" he="67" />为横地理坐标系下的速度,<img file="FDA00005060797400000114.GIF" wi="54" he="79" />为横坐标系下的东向速度,<img file="FDA00005060797400000115.GIF" wi="56" he="86" />为横坐标系下的北向速度;步骤二:利用得到的横坐标系下的导航参数进行横坐标系下的导航机械编排,确定系统的状态变量和量测量;步骤三:根据系统的状态变量和量测量建立卡尔曼滤波器,由卡尔曼滤波方法得到横坐标系下的导航系统的x轴姿态误差角<img file="FDA00005060797400000116.GIF" wi="54" he="81" />和y轴姿态误差角<img file="FDA00005060797400000117.GIF" wi="79" he="83" />步骤四:将得到的x轴姿态误差角<img file="FDA0000506079740000021.GIF" wi="50" he="91" />和y轴姿态误差角<img file="FDA0000506079740000022.GIF" wi="52" he="84" />引入到导航系统的两点校正中,完成极区两点校正。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室