发明名称 基于INS修正的机载长航时天文导航系统的方法
摘要 本发明涉及基于INS修正的机载长航时天文导航系统的方法,用于长航时航空飞行器导航技术领域。通过纯惯导输出的位置信息,引导具有二维转动自由度物理平台上的天文观测子系统,对天空中的恒星目标进行观测,并进行恒星点目标图像在地理坐标系的位置信息解算,利用解算出的位置信息修正惯导输出的含有随时间积累的数据。本方法实现了基于INS的含误差项修正的机载长航时高精度天文定位系统,适用于机载物理平台上的天文定位导航。
申请公布号 CN104567868A 申请公布日期 2015.04.29
申请号 CN201410850164.3 申请日期 2014.12.30
申请人 中国科学院西安光学精密机械研究所 发明人 邓小国;杨晓许;谢梅林;魏宇;江波;唐慧君
分类号 G01C21/02(2006.01)I 主分类号 G01C21/02(2006.01)I
代理机构 西安智邦专利商标代理有限公司 61211 代理人 杨亚婷
主权项 一种基于INS修正以提高机载长航时天文导航系统的方法,其特征在于:当天文导航系统只能观测到一颗导航恒星时,该方法包括以下步骤:1】载机的惯导系统提供载机的的位置信息<img file="FDA0000648582530000011.GIF" wi="44" he="85" />和<img file="FDA0000648582530000012.GIF" wi="75" he="85" />2】天文导航系统根据接收到的惯导系统信息确定导航恒星;当确定的导航恒星仅有一颗时,按照步骤2.1‑2.4计算惯导系统的输出误差;2.1】跟踪观测到导航恒星;得到导航恒星在地理坐标系下的观测方位角A和高度角h;2.2】通过惯导系统提供的经纬度,以及从星历表获得的最佳导航恒星的赤纬δ<sub>A</sub>和时角t<sub>G</sub>,通过式(1)可以计算得到导航恒星的计算高度角<img file="FDA0000648582530000013.GIF" wi="46" he="76" />和方位角<img file="FDA0000648582530000014.GIF" wi="73" he="74" />即<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>sin</mi><mover><mi>h</mi><mo>^</mo></mover><mo>=</mo><mi>sin</mi><msub><mi>&delta;</mi><mi>A</mi></msub><mi>sin</mi><mover><mi>L</mi><mo>^</mo></mover><mo>+</mo><mi>cos</mi><msub><mi>&delta;</mi><mi>A</mi></msub><mi>cos</mi><mover><mi>L</mi><mo>^</mo></mover><mi>cos</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>G</mi></msub><mo>+</mo><mover><mi>&lambda;</mi><mo>^</mo></mover><mo>)</mo></mrow></mrow>]]></math><img file="FDA00006485825300000117.GIF" wi="879" he="78" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>cos</mi><mover><mi>A</mi><mo>^</mo></mover><mo>=</mo><mfrac><mrow><mi>cos</mi><mover><mi>L</mi><mo>^</mo></mover><mi>sin</mi><msub><mi>&delta;</mi><mi>A</mi></msub><mo>-</mo><mi>sin</mi><mover><mi>L</mi><mo>^</mo></mover><mi>cos</mi><msub><mi>&delta;</mi><mi>A</mi></msub><mi>cos</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>G</mi></msub><mo>+</mo><mover><mi>&lambda;</mi><mo>^</mo></mover><mo>)</mo></mrow></mrow><mrow><mi>cos</mi><mover><mi>h</mi><mo>^</mo></mover></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000648582530000015.GIF" wi="1290" he="162" /></maths>其中,<img file="FDA0000648582530000016.GIF" wi="50" he="80" />和<img file="FDA0000648582530000017.GIF" wi="47" he="74" />表示INS输入的经度和纬度;2.3】计算导航恒星的观测方位角A与高度角h和其计算方位角<img file="FDA0000648582530000018.GIF" wi="45" he="71" />与高度角<img file="FDA0000648582530000019.GIF" wi="46" he="72" />之差,如下:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>&Delta;h</mi><mo>=</mo><mover><mi>h</mi><mo>^</mo></mover><mo>-</mo><mi>h</mi></mrow>]]></math><img file="FDA00006485825300000110.GIF" wi="224" he="85" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>&Delta;A</mi><mo>=</mo><mover><mi>A</mi><mo>^</mo></mover><mo>-</mo><mi>A</mi></mrow>]]></math><img file="FDA00006485825300000111.GIF" wi="242" he="86" /></maths>2.4】设△λ和△L为惯导系统的输出的经纬度误差,根据下式(2)计算△λ和△L:<maths num="0005" id="cmaths0005"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><mi>&Delta;&lambda;</mi><mo>=</mo><mfrac><mrow><msub><mi>&alpha;</mi><mn>1</mn></msub><mrow><mo>(</mo><mo>-</mo><mi>cos</mi><mi>A</mi><mi>sinh</mi><mo>&CenterDot;</mo><mi>&Delta;h</mi><mo>-</mo><mi>sin</mi><mi>A</mi><mi>cosh</mi><mo>&CenterDot;</mo><mi>&Delta;A</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&beta;</mi><mn>2</mn></msub><mo>&CenterDot;</mo><mi>cosh</mi><mo>&CenterDot;</mo><mi>&Delta;h</mi></mrow><mrow><msub><mi>&alpha;</mi><mn>1</mn></msub><msub><mi>&beta;</mi><mn>2</mn></msub><mo>-</mo><msub><mi>&alpha;</mi><mn>2</mn></msub><msub><mi>&beta;</mi><mn>1</mn></msub></mrow></mfrac></mtd></mtr><mtr><mtd><mi>&Delta;L</mi><mo>=</mo><mfrac><mrow><msub><mi>&beta;</mi><mn>1</mn></msub><mrow><mo>(</mo><mo>-</mo><mi>cos</mi><mi>A</mi><mi>sinh</mi><mo>&CenterDot;</mo><mi>&Delta;h</mi><mo>-</mo><mi>sin</mi><mi>A</mi><mi>cosh</mi><mo>&CenterDot;</mo><mi>&Delta;A</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&beta;</mi><mn>2</mn></msub><mo>&CenterDot;</mo><mi>cosh</mi><mo>&CenterDot;</mo><mi>&Delta;h</mi></mrow><mrow><msub><mi>&alpha;</mi><mn>2</mn></msub><msub><mi>&beta;</mi><mn>1</mn></msub><mo>-</mo><msub><mi>&alpha;</mi><mn>1</mn></msub><msub><mi>&beta;</mi><mn>2</mn></msub></mrow></mfrac></mtd></mtr></mtable></mfenced>]]></math><img file="FDA00006485825300000112.GIF" wi="1174" he="300" /></maths>其中:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>&alpha;</mi><mn>1</mn></msub><mo>=</mo><mi>cos</mi><mover><mi>L</mi><mo>^</mo></mover><mi>sin</mi><msub><mi>&delta;</mi><mi>A</mi></msub><mo>-</mo><mi>sin</mi><mover><mi>L</mi><mo>^</mo></mover><mi>cos</mi><msub><mi>&delta;</mi><mi>A</mi></msub><mi>cos</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>G</mi></msub><mo>+</mo><mover><mi>&lambda;</mi><mo>^</mo></mover><mo>)</mo></mrow></mrow>]]></math><img file="FDA00006485825300000113.GIF" wi="841" he="93" /></maths><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>&beta;</mi><mn>1</mn></msub><mo>=</mo><mo>-</mo><mi>cos</mi><mover><mi>L</mi><mo>^</mo></mover><mi>cos</mi><msub><mi>&delta;</mi><mi>A</mi></msub><mi>sin</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>G</mi></msub><mo>+</mo><mover><mi>&lambda;</mi><mo>^</mo></mover><mo>)</mo></mrow></mrow>]]></math><img file="FDA00006485825300000114.GIF" wi="601" he="92" /></maths><maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>&alpha;</mi><mn>2</mn></msub><mo>=</mo><mo>-</mo><mi>sin</mi><mover><mi>L</mi><mo>^</mo></mover><mi>sin</mi><msub><mi>&delta;</mi><mi>A</mi></msub><mo>-</mo><mi>cos</mi><mover><mi>L</mi><mo>^</mo></mover><mi>cos</mi><msub><mi>&delta;</mi><mi>A</mi></msub><mi>cos</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>G</mi></msub><mo>+</mo><mover><mi>&lambda;</mi><mo>^</mo></mover><mo>)</mo></mrow></mrow>]]></math><img file="FDA00006485825300000115.GIF" wi="869" he="92" /></maths><maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msub><mi>&beta;</mi><mn>2</mn></msub><mo>=</mo><mi>sin</mi><mover><mi>L</mi><mo>^</mo></mover><mi>cos</mi><msub><mi>&delta;</mi><mi>A</mi></msub><mi>sin</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>G</mi></msub><mo>+</mo><mover><mi>&lambda;</mi><mo>^</mo></mover><mo>)</mo></mrow></mrow>]]></math><img file="FDA00006485825300000116.GIF" wi="557" he="88" /></maths>3】通过下式计算经天文导航系统修正后的载机的地理经纬度信息如下:<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><mi>L</mi><mo>=</mo><mover><mi>L</mi><mo>^</mo></mover><mo>-</mo><mi>&Delta;L</mi></mrow>]]></math><img file="FDA0000648582530000021.GIF" wi="236" he="74" /></maths><maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><mi>&lambda;</mi><mo>=</mo><mover><mi>&lambda;</mi><mo>^</mo></mover><mo>-</mo><mi>&Delta;&lambda;</mi><mo>.</mo></mrow>]]></math><img file="FDA0000648582530000022.GIF" wi="300" he="95" /></maths>
地址 710119 陕西省西安市高新区新型工业园信息大道17号