发明名称 利用背景天文信息对于航天器进行相对导航的系统
摘要 本发明公开了一种利用背景天文信息对于航天器进行相对导航的系统,该系统包括有目标航天器及参考天体的识别单元、星图匹配识别单元、焦平面上夹角关系确定单元、相对轨道运动关系单元和卡尔曼滤波解算单元;该系统保存在导航计算机的处理器内。本发明借助追踪航天器的硬件平台,通过利用目标天文背景的信息作为参考,最后结合目标航天器相对于追踪航天器的距离,采用滤波方法计算出目标航天器在追踪航天器的惯性坐标系OCXCYCZC下的相对运动状态,最终实现对目标航天器的跟踪定位。
申请公布号 CN102607563A 申请公布日期 2012.07.25
申请号 CN201210044323.1 申请日期 2012.02.22
申请人 北京航空航天大学 发明人 陈统;王楷;徐世杰
分类号 G01C21/24(2006.01)I 主分类号 G01C21/24(2006.01)I
代理机构 北京永创新实专利事务所 11121 代理人 李有浩
主权项 1.一种利用背景天文信息对于航天器进行相对导航的系统,其特征在于:该系统包括有目标航天器及参考天体的识别单元(231)、星图匹配识别单元(232)、焦平面上夹角关系确定单元(233)、相对轨道运动关系单元(234)和卡尔曼滤波解算单元(235);目标航天器及参考天体的识别单元(231)第一方面接收可见光敏感器(201)获得的目标天文背景(300)中的图像信息PC;第二方面采用基于局部熵的星体位置提取方法对图像信息PC,再通过线性内插方法获得目标航天器(100)在光学系统焦平面坐标系O<sub>S</sub>X<sub>S</sub>Y<sub>S</sub>Z<sub>S</sub>上的目标航天器-像位置(x<sub>100</sub>,y<sub>100</sub>,0)和参考天体(400)在光学系统焦平面坐标系O<sub>S</sub>X<sub>S</sub>Y<sub>S</sub>Z<sub>S</sub>上的参考天体-像位置<img file="FSA00000675312700011.GIF" wi="259" he="58" />第三方面将目标航天器-像位置(x<sub>100</sub>,y<sub>100</sub>,0)和参考天体-像位置<img file="FSA00000675312700012.GIF" wi="232" he="59" />输出给焦平面上夹角关系确定单元(233);第四方面将参考天体-像位置<img file="FSA00000675312700013.GIF" wi="233" he="59" />输出给星图匹配识别单元(232);星图匹配识别单元(232)第一方面采用多边形角距匹配算法方法对参考天体-像位置<img file="FSA00000675312700014.GIF" wi="230" he="58" />进行处理,得到参考天体(400)在惯性坐标系O<sub>C</sub>X<sub>C</sub>Y<sub>C</sub>Z<sub>C</sub>下的当前时刻t<sub>k</sub>的参考天体-惯性单位矢量<maths num="0001"><![CDATA[<math><mrow><mover><mi>S</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>=</mo><mo>{</mo><msub><mover><mi>S</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>,</mo><msub><mover><mi>S</mi><mo>&OverBar;</mo></mover><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mover><mi>S</mi><mo>&OverBar;</mo></mover><mi>n</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>}</mo><mo>,</mo></mrow></math>]]></maths>第二方面将当前时刻t<sub>k</sub>的参考天体-惯性单位矢量<maths num="0002"><![CDATA[<math><mrow><mover><mi>S</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>=</mo><mo>{</mo><msub><mover><mi>S</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>,</mo><msub><mover><mi>S</mi><mo>&OverBar;</mo></mover><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mover><mi>S</mi><mo>&OverBar;</mo></mover><mi>n</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>}</mo></mrow></math>]]></maths>输出给卡尔曼滤波解算单元(235);焦平面上夹角关系确定单元(233)第一方面将接收到的目标航天器-像位置(x<sub>100</sub>,y<sub>100</sub>,0)和参考天体-像位置<img file="FSA00000675312700017.GIF" wi="256" he="59" />第二方面计算目标航天器(200)视线与参考天体(400)视线之间在当前时刻t<sub>k</sub>的夹角<maths num="0003"><![CDATA[<math><mrow><mrow><mover><mi>&gamma;</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>=</mo><mo>{</mo><msub><mi>&gamma;</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>,</mo><msub><mi>&gamma;</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>&gamma;</mi><mi>n</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>}</mo></mrow><mo>;</mo></mrow></math>]]></maths>第三方面将夹角<maths num="0004"><![CDATA[<math><mrow><mover><mi>&gamma;</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>=</mo><mo>{</mo><msub><mi>&gamma;</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>,</mo><msub><mi>&gamma;</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>&gamma;</mi><mi>n</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>}</mo></mrow></math>]]></maths>输出给卡尔曼滤波解算单元(235);相对轨道运动关系单元(234)依据前一时刻t<sub>k-1</sub>的位置与速度估计值<maths num="0005"><![CDATA[<math><mrow><msub><mover><mi>M</mi><mo>&OverBar;</mo></mover><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>&Delta;x</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>&Delta;y</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>&Delta;z</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd><mtd><mi>&Delta;</mi><msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd><mtd><mi>&Delta;</mi><msub><mover><mi>y</mi><mo>&CenterDot;</mo></mover><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd><mtd><mi>&Delta;</mi><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd></mtr></mtable></mfenced><mi>T</mi></msup></mrow></math>]]></maths>和位置与速度递推关系<img file="FSA000006753127000111.GIF" wi="477" he="74" />获得当前时刻t<sub>k</sub>的位置与速度递推值<maths num="0006"><![CDATA[<math><mrow><msub><mover><mi>M</mi><mo>&OverBar;</mo></mover><mrow><mi>k</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>&Delta;x</mi><mrow><mi>k</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>&Delta;y</mi><mrow><mi>k</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>&Delta;z</mi><mrow><mi>k</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd><mtd><mi>&Delta;</mi><msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mrow><mi>k</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd><mtd><mi>&Delta;</mi><msub><mover><mi>y</mi><mo>&CenterDot;</mo></mover><mrow><mi>k</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd><mtd><mi>&Delta;</mi><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mrow><mi>k</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>;</mo></mrow></math>]]></maths><img file="FSA000006753127000113.GIF" wi="134" he="74" />表示当前时刻t<sub>k</sub>的位置与速度的递推值,<img file="FSA000006753127000114.GIF" wi="106" he="69" />表示前一时刻t<sub>k-1</sub>的位置与速度估计值,<img file="FSA000006753127000115.GIF" wi="69" he="56" />表示追踪航天器200在惯性坐标系O<sub>C</sub>X<sub>C</sub>Y<sub>C</sub>Z<sub>C</sub>下的位置,<img file="FSA000006753127000116.GIF" wi="39" he="66" />表示递推关系函数;卡尔曼滤波解算单元(235)依据滤波关系<maths num="0007"><![CDATA[<math><mrow><msub><mover><mi>M</mi><mo>&OverBar;</mo></mover><mi>k</mi></msub><mo>=</mo><mover><mi>G</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><msub><mover><mi>M</mi><mo>&OverBar;</mo></mover><mrow><mi>k</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo><mi>&rho;</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>,</mo><mover><mi>S</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>,</mo><mover><mi>&gamma;</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow></mrow></math>]]></maths>对当前时刻t<sub>k</sub>的测量距离ρ(t<sub>k</sub>)、当前时刻t<sub>k</sub>的参考天体-惯性单位矢量<maths num="0008"><![CDATA[<math><mrow><mover><mi>S</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>=</mo><mo>{</mo><msub><mover><mi>S</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>,</mo><msub><mover><mi>S</mi><mo>&OverBar;</mo></mover><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mover><mi>S</mi><mo>&OverBar;</mo></mover><mi>n</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>}</mo><mo>,</mo></mrow></math>]]></maths>当前时刻t<sub>k</sub>的夹角<maths num="0009"><![CDATA[<math><mrow><mover><mi>&gamma;</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>=</mo><mo>{</mo><msub><mi>&gamma;</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>,</mo><msub><mi>&gamma;</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>&gamma;</mi><mi>n</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>}</mo></mrow></math>]]></maths>和递推值<img file="FSA00000675312700023.GIF" wi="133" he="74" />进行滤波处理,得到当前时刻t<sub>k</sub>的目标相对位置与速度的估计值<img file="FSA00000675312700024.GIF" wi="185" he="69" />表示卡尔曼滤波关系函数。
地址 100191 北京市海淀区学院路37号