中文核心期刊

中国科学引文数据库(CSCD)来源期刊

中国高校优秀科技期刊

中国宇航学会深空探测技术专业委员会会刊

高级检索

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

一种基于恒星分布的星敏感器导航星库制作方法

叶志龙 孙朔冬 陈纾 吴迪 武斌

叶志龙, 孙朔冬, 陈纾, 吴迪, 武斌. 一种基于恒星分布的星敏感器导航星库制作方法[J]. 深空探测学报(中英文), 2018, 5(1): 90-96. doi: 10.15982/j.issn.2095-7777.2018.01.013
引用本文: 叶志龙, 孙朔冬, 陈纾, 吴迪, 武斌. 一种基于恒星分布的星敏感器导航星库制作方法[J]. 深空探测学报(中英文), 2018, 5(1): 90-96. doi: 10.15982/j.issn.2095-7777.2018.01.013
YE Zhilong, SUN Shuodong, CHEN Shu, WU Di, WU Bin. Establishment of a Navigation Star Database Based on Star Distribution[J]. Journal of Deep Space Exploration, 2018, 5(1): 90-96. doi: 10.15982/j.issn.2095-7777.2018.01.013
Citation: YE Zhilong, SUN Shuodong, CHEN Shu, WU Di, WU Bin. Establishment of a Navigation Star Database Based on Star Distribution[J]. Journal of Deep Space Exploration, 2018, 5(1): 90-96. doi: 10.15982/j.issn.2095-7777.2018.01.013

一种基于恒星分布的星敏感器导航星库制作方法

doi: 10.15982/j.issn.2095-7777.2018.01.013
详细信息
    作者简介:

    叶志龙(1988– ),男,工程师,硕士,主要研究方向:天文导航及光电敏感器。通信地址:上海市闵行区中春路1555号(201109)电话:(021)24183216 E-mail:zhilongye@163.com

  • ● Star distribution is considered to evaluate navigation star database. ● 100% all-sky recognition rate and good star coverage can be achieved. ● Better attitude precision can be obtained by proposed star database.

Establishment of a Navigation Star Database Based on Star Distribution

  • 摘要:

    为了进一步提升星敏感器的定姿精度,分析了恒星的空间几何分布对姿态精度的影响,提出了一种基于恒星分布的星敏感器导航星库制作方法。在分析对比了现有星库的基础上,选择了精度较高的依巴谷星库作为基础星库,并剔除了暗星和双星以精简星库。遍历全天球对每个视场的星库中的空间几何分布进行分析对比,删除定姿较差的3颗星。对星库进行遍历补偿,保证了星库中恒星分布的均匀性。实验结果表明,全天球轨道测试中,本文方法制作星库的全天识别率为100%,且星颗分布数10颗以上天区占到97.64%,具有良好的覆盖性,同时有效提升了姿态精度。

    Highlights
    ● Star distribution is considered to evaluate navigation star database. ● 100% all-sky recognition rate and good star coverage can be achieved. ● Better attitude precision can be obtained by proposed star database.
  • 图  1  精简前后的星库三维图

    Fig.  1  3-D image of star datebase before and after simplification

    图  2  姿态误差与体积变化关系

    Fig.  2  Relationship between attitude error and volume

    图  3  精简星库恒星分布图

    Fig.  3  The distribution map of simplified star datebase

    图  4  视场划分图

    Fig.  4  Split graph of FOV

    图  5  删除误差较大的定姿星

    Fig.  5  Delete large error star for attitude determination

    图  6  选取视场中最亮的15颗星

    Fig.  6  Select the most 15 brightest star in FOV

    图  7  扇形区中添加星后

    Fig.  7  Star is added in section

    图  8  导航星库的二维分布图

    Fig.  8  2-D distribution map of navigation star table

    图  9  理想姿态与实际姿态的误差

    Fig.  9  Error between ideal attitude and actual attitude

    表  1  常见的几种星库

    Table  1  Several common star datebase

    星库名称 星库历元 总星数 星库精度/(″) 分布密度(star/sep.deg)
    FK5 J2000.0 1 535 0.02 0.037
    FK4 Exl J2000.0 3 117 0.04 0.075
    FK6 J2000.0 4 150 0.000 4 0.100
    SAO J2000.0 258 997 1.0 6.267
    PPM(North) J2000.0 181 731 0.27 4.398
    PPM(South) J2000.0 197 179 0.11 4.771
    GSC 2.2 J2000.0 455 851 237 0.3 11 031.060
    Hipparcos J2000.0 118 218 0.000 8 2.861
    Tychσ1 J2000.0 1 058 332 0.007 25.610
    TRC J2000.0 990 182 0.040 23.961
    Tychσ2 J2000.0 2 539 913 0.007 61.463
    下载: 导出CSV

    表  2  筛选后星库的全天识别情况

    Table  2  The Recognition of all-sky by navigation star datebase

    轨道参数
    /(°)
    最少识
    别星数
    是否
    识别
    覆盖
    率/%
    轨道参数
    /(°)
    最少识
    别星数
    是否
    识别
    覆盖
    率/%
    赤纬0 7 100 赤纬–10 8 100
    赤纬10 7 100 赤纬–20 8 100
    赤纬20 6 100 赤纬–30 9 100
    赤纬30 7 100 赤纬–40 9 100
    赤纬40 4 100 赤纬–50 6 100
    赤纬50 8 100 赤纬–60 9 100
    赤纬60 8 100 赤纬–70 8 100
    赤纬70 9 100 赤纬–80 11 100
    赤纬80 14 100 100
    下载: 导出CSV

    表  3  星库中星数的分布情况

    Table  3  The distribution of star number in star datebase

    星数情况/颗 天区数 百分比/%
    小于5 2 0.02
    5~10 234 2.34
    10~15 2 324 23.24
    15~20 4 764 47.64
    大于20 2 676 26.76
    下载: 导出CSV
  • [1] 朱长征,居永忠,杜晓辉. 导航星库制定方法研究[J]. 宇航学报,2010,31(5):1327-1330.

    ZHU C Z,JU Y Z,DU X H. Approach to founding a guide star catalogue[J]. Journal of Astronautics,2010,31(5):1327-1330.
    [2] 郑万波,刘智,郝志航. 二维精简索引分层导航星库的构造[J]. 吉林大学学报(信息科学版),2003,21(2):123-127.

    ZHENG W B,LIU Z,HAO Z H. Construction of delaminating guide star database with 2-D condensed index[J]. Journal of Changchun Post and Telecommunication Institute,2003,21(2):123-127.
    [3] 许士文,李葆华,张迎春,等. 用字符匹配进行星图识别的导航星库存储方法[J]. 哈尔滨工业大学学报[J]. 2005,37(6):819-821.

    XU S W,LI B H,ZHANG Y C,et al. A store navigation star database method for star map recognition using character match[J]. Journal of Harbin Institute of Technology,2005,37(6):819-821.
    [4] 蔡志武,韩春好,陈金平. 一种高精度星载导航星库的构建方法[J]. 测绘科学技术学报,2006,23(1):29-32.

    CAI Z W,HAN C H,CHEN J P. A approach of construction high precision guide star catalogue on board[J]. Journal of Zhengzhou Institute of Surveying and Mapping,2006,23(1):29-32.
    [5] ZHANG H, SANG H. A double-FOV star sensor for high dynamic spacecraft navigation[C]//Proceedings of the SPIE, MIPPR 2009. Yichang, China: [s.n.], 2009.
    [6] LIEBE C C. Star tracker for attitude determination[J]. IEEE Aerospace and Electronic Systems Magazine,1995,10(6):10-16.
    [7] SUN Y, XIAO Y, GENG Y. On-obrit calibration of star sensor based on a new lens distortion model[C]// The 32 nd Chinese Control Conference. Xi'an, China: IEEE, 2013: 4989-4994.
    [8] 张磊,魏仲慧,何昕. 快速全天自主星图识别[J]. 光学精密工程,2009,17(4):909-915.

    ZHANG L. WEI Z H,HE X. Fast all-sky autonmous star identification[J]. Optics and Precision Engineering,2009,17(4):909-915
    [9] 王丽娜,赵慧,熊智,等. 基于导航恒星几何分布的天文导航定姿系统误差建模及误差特性分析[J]. 兵工学报,2015,36(10):1933-1942.

    WANG L N,ZHAO H,XIONG Z,et al. Modeling and error characteristics analysis of celestial attitude determination error based on geometric configuration of guide stars[J]. Acta Armamentarii,2015,36(10):1933-1942.
    [10] 赵慧,熊智,王丽娜,等. 基于恒星几何构型分布的天文定位误差建模及误差特性分析[J]. 兵工学报,2015,36(5):813-822.

    ZHAO H,XIONG Z,WANG L N,et al. Modeling of celestial positioning error and analysis of error characteristics based on distribution of guide stars[J]. Acta Armamentarii,2015,36(5):813-822.
  • [1] 王卓, 徐瑞.  基于多目标优化的深空探测器姿态组合规划方法 . 深空探测学报(中英文), 2021, 8(2): 147-153. doi: 10.15982/j.issn.2096-9287.2021.20200069
    [2] 裴福俊, 严鸿, 朱明君.  太阳敏感器辅助的分布式EKF-SLAM火星车自主导航方法 . 深空探测学报(中英文), 2020, 7(2): 191-196. doi: 10.15982/j.issn.2095-7777.2020.20171117001
    [3] 徐浩, 裴福俊, 蒋宁.  一种基于李群描述的深空探测器姿态估计方法 . 深空探测学报(中英文), 2020, 7(1): 102-108. doi: 10.15982/j.issn.2095-7777.2020.20171117002
    [4] 段建锋, 刘勇, 李勰, 陈明, 王兆魁.  “嫦娥4号”中继星任务轨道确定问题初探 . 深空探测学报(中英文), 2018, 5(6): 531-538. doi: 10.15982/j.issn.2095-7777.2018.06.005
    [5] 杨成伟, 郑建华.  脉冲星X射线数据处理与分析 . 深空探测学报(中英文), 2018, 5(3): 219-225. doi: 10.15982/j.issn.2095-7777.2018.3.003
    [6] 王元超, 郑建华, 潘之辰, 李明涛.  脉冲星候选样本分类方法综述 . 深空探测学报(中英文), 2018, 5(3): 203-211,218. doi: 10.15982/j.issn.2095-7777.2018.3.001
    [7] 李艳丽, 郑建华, 高东.  基于FPGA+DSP的X射线脉冲星导航原理样机的设计与实现 . 深空探测学报(中英文), 2018, 5(3): 226-234. doi: 10.15982/j.issn.2095-7777.2018.3.004
    [8] 易韦韦, 偶晓娟, 许静文, 李晶, 李冰.  脉冲星导航试验卫星观测数据处理与分析 . 深空探测学报(中英文), 2018, 5(3): 241-245,261. doi: 10.15982/j.issn.2095-7777.2018.3.006
    [9] 喻子原, 刘劲, 宁晓琳, 马辛, 桂明臻, 康志伟.  面向编队飞行的天文多普勒差分/脉冲星组合导航 . 深空探测学报(中英文), 2018, 5(3): 212-218. doi: 10.15982/j.issn.2095-7777.2018.3.002
    [10] 李亮, 王广利, 郭丽, 胡小工.  脉冲星导航的天体测量考虑 . 深空探测学报(中英文), 2018, 5(3): 235-240. doi: 10.15982/j.issn.2095-7777.2018.3.005
    [11] 康珅, 单家元.  航天器分布式有限时间编队方法 . 深空探测学报(中英文), 2017, 4(4): 390-394. doi: 10.15982/j.issn.2095-777.2017.04.0013
    [12] 邵巍, 顾天昊.  一种基于星表特征直线匹配的着陆器位姿估计算法 . 深空探测学报(中英文), 2017, 4(3): 281-286. doi: 10.15982/j.issn.2095-7717.2017.03.013
    [13] 郑伟, 张璐, 王奕迪.  基于星联网的深空自主导航方案设计 . 深空探测学报(中英文), 2017, 4(1): 31-37. doi: 10.15982/j.issn.2095-7777.2017.01.005
    [14] 程宇峰, 润一, 王密.  深空探测光学导航敏感器在轨几何定标方法 . 深空探测学报(中英文), 2016, 3(3): 228-236. doi: 10.15982/j.issn.2095-7777.2016.03.006
    [15] 刘瑞霞, 张剑桥.  基于测速测角敏感器的火星探测器自主导航方法研究 . 深空探测学报(中英文), 2016, 3(3): 219-224. doi: 10.15982/j.issn.2095-7777.2016.03.004
    [16] 平劲松, 简念川, 张添翼, 张建辉.  47亿千米距离上的追踪——研究团队成功测量到即将飞掠冥王星的"新地平线"探测器多普勒速度 . 深空探测学报(中英文), 2015, 2(2): 192-192. doi: 10.15982/j.issn.2095-7777.2015.02.015
    [17] 武长青, 徐瑞, 朱圣英.  基于对数势函数的深空探测器姿态规划与控制方法 . 深空探测学报(中英文), 2015, 2(4): 365-370. doi: 10.15982/j.issn.2095-7777.2015.04.011
    [18] 刘也, 曹建峰.  一种基于星载GPS的月球探测器导航算法 . 深空探测学报(中英文), 2015, 2(4): 313-317. doi: 10.15982/j.issn.2095-7777.2015.04.003
    [19] 王大铁, 李骥, 黄翔宇, 张洪华.  月球软着陆过程高精度自主导航避障方法 . 深空探测学报(中英文), 2014, 1(1): 44-51.
    [20] 魏二虎, 杨洪洲, 张帅, 刘经南, 易慧.  脉冲星非实时平差的火星探测自主导航模型 . 深空探测学报(中英文), 2014, 1(4): 298-302. doi: 10.15982/j.issn.2095-7777.2014.04.009
  • 加载中
图(9) / 表 (3)
计量
  • 文章访问数:  1618
  • HTML全文浏览量:  64
  • PDF下载量:  920
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-05-11
  • 修回日期:  2017-06-20
  • 网络出版日期:  2020-09-10
  • 刊出日期:  2018-02-01

一种基于恒星分布的星敏感器导航星库制作方法

doi: 10.15982/j.issn.2095-7777.2018.01.013
    作者简介:

    叶志龙(1988– ),男,工程师,硕士,主要研究方向:天文导航及光电敏感器。通信地址:上海市闵行区中春路1555号(201109)电话:(021)24183216 E-mail:zhilongye@163.com

  • ● Star distribution is considered to evaluate navigation star database. ● 100% all-sky recognition rate and good star coverage can be achieved. ● Better attitude precision can be obtained by proposed star database.

摘要: 

为了进一步提升星敏感器的定姿精度,分析了恒星的空间几何分布对姿态精度的影响,提出了一种基于恒星分布的星敏感器导航星库制作方法。在分析对比了现有星库的基础上,选择了精度较高的依巴谷星库作为基础星库,并剔除了暗星和双星以精简星库。遍历全天球对每个视场的星库中的空间几何分布进行分析对比,删除定姿较差的3颗星。对星库进行遍历补偿,保证了星库中恒星分布的均匀性。实验结果表明,全天球轨道测试中,本文方法制作星库的全天识别率为100%,且星颗分布数10颗以上天区占到97.64%,具有良好的覆盖性,同时有效提升了姿态精度。

注释:
1)  ● Star distribution is considered to evaluate navigation star database. ● 100% all-sky recognition rate and good star coverage can be achieved. ● Better attitude precision can be obtained by proposed star database.

English Abstract

叶志龙, 孙朔冬, 陈纾, 吴迪, 武斌. 一种基于恒星分布的星敏感器导航星库制作方法[J]. 深空探测学报(中英文), 2018, 5(1): 90-96. doi: 10.15982/j.issn.2095-7777.2018.01.013
引用本文: 叶志龙, 孙朔冬, 陈纾, 吴迪, 武斌. 一种基于恒星分布的星敏感器导航星库制作方法[J]. 深空探测学报(中英文), 2018, 5(1): 90-96. doi: 10.15982/j.issn.2095-7777.2018.01.013
YE Zhilong, SUN Shuodong, CHEN Shu, WU Di, WU Bin. Establishment of a Navigation Star Database Based on Star Distribution[J]. Journal of Deep Space Exploration, 2018, 5(1): 90-96. doi: 10.15982/j.issn.2095-7777.2018.01.013
Citation: YE Zhilong, SUN Shuodong, CHEN Shu, WU Di, WU Bin. Establishment of a Navigation Star Database Based on Star Distribution[J]. Journal of Deep Space Exploration, 2018, 5(1): 90-96. doi: 10.15982/j.issn.2095-7777.2018.01.013
    • 星敏感器作为一种高精度的姿态测量器件,已广泛应用于卫星,深空探测等航天飞行器中,其通过匹配探测恒星信息与导航星库中的恒星信息,并采用QUEST(Quaternion Estimator)算法解算姿态。导航星库是星敏感器进行星图识别和跟踪的关键信息,其合理性和完备性对星敏感器工作的鲁棒性和精度都有着重要的影响。现有的导航星库制作方法主要有星等阈值法、星等加权法、正交网格法、回归选取法、二维精简索引分层构造法等[1-4]。星等阈值法是通过选取星等小于一定阈值的星构成导航星库,该方法虽然简单,也具有良好的匹配性,但是恒星分布的不均匀性会导致星图识别准确性下降。若采用提高星等阈值来改善此问题,则会导致星库容量增大,不利于星敏感器存储使用,也增加星库搜索的时间。星等加权法通过对星等较低的星赋予较大权值来提高导航星库恒星较少天区的星颗数,但依旧存在与星等阈值法类似的问题。正交网格法则是通过映射法将候选导航星映射从球面到平面上,并利用互不重叠的等面积正交网格将其划分,选取每个网格中最亮的星作为导航星。此方法可通过调节网格大小筛选导航,有较强的灵活性。但是未考虑到星敏感器的实际探测能力,导致一些网格中的导航星不可用。回归选取法则是将星等阈值法中的固定阈值变为动态阈值,可有效地筛选出导航星,但动态阈值的自适应性仍不够完善。二维精简索引分层构造法把导航星分为局部亮星和暗星。以局部亮星构造导航三角形,节约了导航星库空间,大大缩短了搜索时间,然而该方法对星敏感器的视场有一定的要求,仅适用于探测星等较低的第二代星敏感器。

      由此可见,现有的导航星库多以星等和分布均匀性为依据选取导航星,此类方法虽然可以较为完备地筛选出导航星库,但未考虑星库中恒星的空间分布对于定姿的影响,往往会出现星敏感器视场划过某一天区的时候,其输出姿态的误差有一定幅度增大[5-8]。因此,本文考虑利用恒星的空间分布与定姿精度的关系,并结合星等和恒星分布的均匀性来选取制作导航星库。

    • 星库记录了恒星位置、星等、光谱等天文参数。由于编制星库使用仪器、观测条件和处理方法不一样,不同星库中恒星的颗数不一样,恒星的位置精度也会有差异。星库是设计制作导航星库的基础,也是决定星敏感器定姿精度的关键因素之一。表 1给出了天文领域实际使用的几种星库。

      表 1  常见的几种星库

      Table 1.  Several common star datebase

      星库名称 星库历元 总星数 星库精度/(″) 分布密度(star/sep.deg)
      FK5 J2000.0 1 535 0.02 0.037
      FK4 Exl J2000.0 3 117 0.04 0.075
      FK6 J2000.0 4 150 0.000 4 0.100
      SAO J2000.0 258 997 1.0 6.267
      PPM(North) J2000.0 181 731 0.27 4.398
      PPM(South) J2000.0 197 179 0.11 4.771
      GSC 2.2 J2000.0 455 851 237 0.3 11 031.060
      Hipparcos J2000.0 118 218 0.000 8 2.861
      Tychσ1 J2000.0 1 058 332 0.007 25.610
      TRC J2000.0 990 182 0.040 23.961
      Tychσ2 J2000.0 2 539 913 0.007 61.463

      表 1中SAO星库和Hipparcos星库较容易获得,其中Hipparcos星库的精度较高,可达到0.000 8″,因此,本文选择Hipparcos星库为基础星库。从表 1中可以看出,Hipparcos星库中恒星的数目近12万颗,数目过于庞大,且星敏感器光学系统的敏感能力有限。因此,可根据星敏感器的敏感极限来删减星库。另一方面,星库中存在的双星和变星,对姿态解算有一定影响,可删除双星和变星进一步精简星库。

    • 星敏感器的定姿误差来源主要有星点质心定位误差,光学系统误差以及定姿算法的转换误差等。为便于分析,将星点定位误差和光学系统误差认为是天文观测误差。星敏感器通过获取星点在星敏感器测量坐标系的坐标,以及星点对应于惯性坐标的坐标,通过Quest算法解算得到姿态信息。

    • 假设恒星在测量坐标系b下的三维坐标为${ {w}_k} =$ ${\left[ {\begin{array}{*{20}{c}} {{x_{bk}}} & {{y_{bk}}} & {{z_{bk}}} \end{array}} \right]^{\rm T}}$$k = 1,2,3 \ldots n$,在惯性坐标系i下的三维坐标为${ {v}_k} = {\left[ {\begin{array}{*{20}{c}} {{x_{ik}}} & {{y_{ik}}} & {{z_{ik}}} \end{array}} \right]^{\rm T}}$$k = 1,2,3 \ldots n$,则

      $$ { {v}_k} = {A}_b^i{w_k} $$ (1)

      其中:${A}_b^i$为观测坐标系b到惯性坐标系i的姿态转移矩阵。若记${W} = {\left[ {\begin{array}{*{20}{c}} \!\!\!\!\!{w_1^{\rm T}}\!\!\!\!\!\!\!\! & {w_2^{\rm T}}\! & \!{\begin{array}{*{20}{c}} \!\!\!\!\!\!\!\!\!\! \cdots\!\!\!\!\!\!\!\! & {w_n^{\rm T}} \!\!\!\!\!\!\! \end{array}} \end{array}} \right]^{\rm{T}}}$${V} = {\left[ {\begin{array}{*{20}{c}} {v_1^{\rm T}} & {v_2^{\rm T}} & {\begin{array}{*{20}{c}} \!\!\!\!\!\!\!\!\!\! \cdots\!\!\!\!\!\!\!\! & {v_n^{\rm T}}\!\!\!\!\!\!\! \end{array}} \end{array}} \right]^{\rm{T}}}$,依据上式可知

      $$ {W} = {VA}_b^i $$ (2)

      星敏感器在观测星数不小于3时,可用最小二乘法求解姿态转移矩阵,即

      $$ {A}_b^i = {\left( {{{V}^{{\rm T}}}{V}} \right)^{ - 1}}{{V}^{{\rm T}}}{W} = {QW} $$ (3)

      其中:${Q} = {\left( {{{V}^{{\rm T}}}{V}} \right)^{ - 1}}{{V}^{{\rm T}}}$

      由于存在天文观测误差,恒星的实际观测矢量为$\tilde{ W} = {W} + \Delta {W}$,其中∆W为天文观测误差。则

      $$ \widetilde {{A}_b^i} = {Q} \tilde{ W} = {QW} + {Q}\Delta {W} = {A}_b^i + \Delta {A}_b^i $$ (4)

      若以观测坐标系b为参考,记定姿误差矢量为$ {\delta} \varphi = {\left[ {\begin{array}{*{20}{c}} {\delta {\varphi _1}} & {\delta {\varphi _2}} & {\delta {\varphi _3}} \end{array}} \right]^{\rm{T}}}$,当定姿误差角为一小量时,$\widetilde {{A}_b^i}$可表示为

      $$ \widetilde {{A}_b^i} = {A}_{b'}^i = {A}_b^i{A}_{b'}^b $$ (5)

      式中${A}_{b'}^b = \left[ {\begin{array}{*{20}{c}} 1 & {\delta {\varphi _3}} & { - \delta {\varphi _2}}\\ { - \delta {\varphi _3}} & 1 & {\delta {\varphi _1}}\\ {\delta {\varphi _2}} & { - \delta {\varphi _1}} & 1 \end{array}} \right]$b’为包含误差计算的测量坐标系。依据式(4)和式(5)可得

      $$ \Delta {A}_b^i = \widetilde {{A}_b^i} - {A}_b^i = {A}_b^i\left[ {\begin{array}{*{20}{c}} 0 & {\delta {\varphi _3}} & { - \delta {\varphi _2}}\\ { - \delta {\varphi _3}} & 0 & {\delta {\varphi _1}}\\ {\delta {\varphi _2}} & { - \delta {\varphi _1}} & 0 \end{array}} \right] = {A}_b^i{M} $$ (6)

      式中:${M} = \left[ {\begin{array}{*{20}{c}} 0 & {\delta {\varphi _3}} & { - \delta {\varphi _2}}\\ { - \delta {\varphi _3}} & 0 & {\delta {\varphi _1}}\\ {\delta {\varphi _2}} & { - \delta {\varphi _1}} & 0 \end{array}} \right] = {A}_i^b\Delta {A}_b^i$,记M的协方差矩阵为PM,误差矩阵$\Delta {A}_b^i$的协方差矩阵为P,定姿误差矢量δφ的协方差矩阵为Pφ

      $$ \begin{array}{l} {{P}_{{M}}} = {M}{{M}^{\rm T}} = \left( {{A}_i^b\Delta {A}_b^i} \right){\left( {{A}_i^b\Delta {A}_b^i} \right)^{\rm T}} = {A}_i^b{{P}_\Delta }{A}_b^i\\ = \left[ {\begin{array}{*{20}{c}} {{{\left( {\delta {\varphi _2}} \right)}^2} + {{\left( {\delta {\varphi _3}} \right)}^2}} & { - \delta {\varphi _1}\delta {\varphi _2}} & { - \delta {\varphi _1}\delta {\varphi _3}}\\ { - \delta {\varphi _1}\delta {\varphi _2}} & {{{\left( {\delta {\varphi _1}} \right)}^2} + {{\left( {\delta {\varphi _3}} \right)}^2}} & { - \delta {\varphi _2}\delta {\varphi _3}}\\ { - \delta {\varphi _1}\delta {\varphi _3}} & { - \delta {\varphi _2}\delta {\varphi _3}} & {{{\left( {\delta {\varphi _1}} \right)}^2} + {{\left( {\delta {\varphi _2}} \right)}^2}} \end{array}} \right] \end{array} $$ (7)
      $$ {{P}_\varphi } = \delta \varphi {\left( {\delta \varphi } \right)^{\rm{T}}} = \left[ {\begin{array}{*{20}{c}} {{{\left( {\delta {\varphi _1}} \right)}^2}} & {\delta {\varphi _1}\delta {\varphi _2}} & {\delta {\varphi _1}\delta {\varphi _3}}\\ {\delta {\varphi _1}\delta {\varphi _2}} & {{{\left( {\delta {\varphi _2}} \right)}^2}} & {\delta {\varphi _2}\delta {\varphi _3}}\\ {\delta {\varphi _1}\delta {\varphi _3}} & {\delta {\varphi _2}\delta {\varphi _3}} & {{{\left( {\delta {\varphi _3}} \right)}^2}} \end{array}} \right] $$ (8)

      由式(7)和式(8)可知

      $$ \begin{array}{l} {\rm{tr}}{{P}_\varphi } = {\left( {\delta {\varphi _1}} \right)^2} + {\left( {\delta {\varphi _2}} \right)^2} + {\left( {\delta {\varphi _3}} \right)^2} = {\rm{tr}}{{P}_{{M}}}\\ = \frac{1}{2}{\rm{tr}}\left( {{A}_i^b{{P}_\Delta }{A}_b^i} \right) = \frac{1}{2}{\rm{tr}}\left( {{{P}_\Delta }{A}_i^b{A}_b^i} \right) = \frac{1}{2}{\rm{tr}}\left( {{{P}_\Delta }} \right) \end{array} $$ (9)

      天文观测误差∆W主要取决于敏感器件及星点提取等误差,与选星无关。因此,假设∆W一定,其方差为$\sigma _W^2$。由式(3)、式(4)和式(9)可得

      $$ \begin{array}{l} {{P}_\Delta } = {\rm{E}}\left( {\Delta {A}_b^i{{\left( {\Delta {A}_b^i} \right)}^{\rm T}}} \right) = {\rm{E}}\left( {{Q}\Delta {W}{{\left( {{Q}\Delta {W}} \right)}^{\rm T}}} \right)\\ = {Q}{\rm{E}}\left( {\Delta {W}{{\left( {\Delta {W}} \right)}^{\rm T}}} \right){{Q}^{\rm T}} = {Q}{{Q}^{\rm T}}\sigma _W^2 = {\left( {{{V}^{\rm T}}{V}} \right)^{ - 1}}\sigma _W^2 \end{array} $$ (10)

      由式(9)和式(10)可知

      $$ {\rm{tr}}{{P}_\varphi } = \frac{1}{2}{\rm{tr}}\left( {{{P}_\Delta }} \right) = \frac{{{\rm{tr}}\left( {{{\left( {{{V}^{\rm T}}{V}} \right)}^{ - 1}}} \right)}}{2}\sigma _W^2 = \frac{{{\rm{tr}}\left( {{{\left( {{{\rm{V}}^{\rm{T}}}{\rm{V}}} \right)}^{\rm{*}}}} \right)}}{{2{\rm{det}}\left( {{{\rm{V}}^{\rm{T}}}{\rm{V}}} \right)}}\sigma _W^2 $$ (11)

      以Quest算法的三星定姿为例,其中${V} = \left[ v\!\!{v_1^{\rm T}} \quad {v_2^{\rm T}} \cdots \right.$$\left. {v_n^{\rm T}} \right]^{\rm{T}}$v1v2v3为3颗定姿星的惯性坐标系下的单位矢量,矩阵V一定程度上反映了恒星的空间几何分布[9-10]。假设惯性坐标i的原点为O,3颗星对应的单位矢量为Ov1Ov2Ov3。3颗恒星与原点构成的四面体体积Tvol

      $$ {T_{vol}} = \frac{1}{6}\left( {{v_1} \times {v_2}} \right) \cdot {v_3} = \frac{1}{6}{\rm{det}}\left( {V} \right) $$ (12)

      将上式代入式(10),式(9),令${B} = {\left( {{{V}^{\rm T}}{V}} \right)^*}$

      $$ \begin{array}{l} {\left( {\delta {\varphi _1}} \right)^2} + {\left( {\delta {\varphi _2}} \right)^2} + {\left( {\delta {\varphi _3}} \right)^2} = \displaystyle\frac{{{\rm{tr}}\left( {B} \right)}}{{72{T_{vol}}^2}}\sigma _W^2{\rm{tr}}{{P}_\varphi }\\ {\rm{ = }}\frac{{{\rm{tr}}\left( {{{\left( {{{V}^{\rm T}}{V}} \right)}^*}} \right)}}{{2{\rm{det}}\left( {{{V}^{\rm T}}{V}} \right)}}\sigma _W^2{\rm{ = }}\displaystyle\frac{{{\rm{tr}}\left( {B} \right)}}{{72{T_{vol}}^2}}\sigma _W^2 \end{array} $$ (13)

      由式(13)可知,定姿误差的方差与星矢量构成的体积以及矩阵B的迹相关。为了说明星矢量构成体积和B的迹对定姿精度的影响量级,以(120°,30°)天区为例,在天区探测到的恒星共有30颗,共有4 060种定姿组合,其中$\displaystyle\frac{{\max \left( {{T_{vol}}^2} \right)}}{{\min \left( {{T_{vol}}^2} \right)}} = 1.95 \times {10^7}$$\displaystyle\frac{{\max \left( {{\rm{tr}}\left( {B} \right)} \right)}}{{\min \left( {{\rm{tr}}\left( {B} \right)} \right)}} = 34.57$,由此可见,恒星的空间几何分布是决定定姿误差的主要因素。因此,可考虑恒星的空间几何分布因素作为导航星库选取策略,以提升星敏感器的定姿精度。

    • 综上所述,本文选取Hipparcos星库作为基础星库来筛选导航星库。

      第1步:星库精简,根据星敏感器的敏感极限星等精简星库,删除原始星库中星等大于6.2的恒星,同时删除双星,光学头部的视场为14.5° × 14.5°,探测器尺寸为2 048 × 2 048,若星点成像为3 × 3的高斯分布,则能分辨两颗恒星的最小间距为6个像素,且星跟踪模式下的星跟踪半径最大设为30,则剔除星对角矩小于$\displaystyle\frac{{14.5^\circ }}{{2048}} \times 30 \approx 0.212^\circ $,即删星对角矩小于0.212°的双星。图 1为精简前后的星库三维图示。精简后的星库共有4 975颗星。

      图  1  精简前后的星库三维图

      Figure 1.  3-D image of star datebase before and after simplification

      第2步:随机生成光轴指向,并搜索当前光轴指向所在天区视场内的星,任取3颗星计算其星矢量组成的四面体的体积,分别计算每种组合中恒星矢量组成四面体体积。

      第3步:基于恒星几何分布选取导航星。由2.1所述可知,用于定姿的恒星矢量空间构成体积越大,定姿误差越小,定姿精度越高。以(120°,30°)天区为例,在基础星库中该天区共有30颗星,任取3颗导航星,则有4 060种组合。图 3给出了姿态误差与四面体体积的变化趋势图,从图中可以看出,四面体体积越大,姿态角的误差越小,即姿态角度的精度越高。图 4给出了精简后的星库在(120°,30°)天区的星点分布图,图中红色的3颗星为空间几何分布体积最小的3颗星,可以看出这3星在平面坐标上基本处于一条直线上,其星矢量形成的空间四面体的体积最小,解算出的姿态数据误差也是最大的,可先删除定姿误差较大的3颗星。

      第4步:在减少恒星分布对定姿误差影响后,再次遍历全天球,每个视场保留最亮的15颗,以进一步精简星库,精简之后星库共有导航星4 065颗。

      图  2  姿态误差与体积变化关系

      Figure 2.  Relationship between attitude error and volume

      图  3  精简星库恒星分布图

      Figure 3.  The distribution map of simplified star datebase

      图  4  视场划分图

      Figure 4.  Split graph of FOV

      第5步:为了保证导航星库在视场中的星分布尽可能均匀,本文将内切圆视场划分成4个面积相等的扇形区域,如图 4所示。遍历每个视场,统计每个扇形区域的星颗数,若出现有星颗数为0的扇形区,则将基础星库在该扇形区最亮的星添加到导航星库中。仿真结果表明,共有2 260个天区出现星颗数为0的区域,共添加126颗星,即星库星颗数共为4 191颗。图 57给出了扇区为0区域加添星的过程。图 8为星库的全天球分布图。

      图  5  删除误差较大的定姿星

      Figure 5.  Delete large error star for attitude determination

      图  6  选取视场中最亮的15颗星

      Figure 6.  Select the most 15 brightest star in FOV

      图  7  扇形区中添加星后

      Figure 7.  Star is added in section

    • 为了验证星库的有效性和完备性,对制作的导航星库进行全天球识别验证,输入15组动力学参数进行全天识别测试全天球的覆盖性,即赤经0°到360°,赤纬分别为–90°到90°的轨道。识别情况如表 2所示。为了进一步评价星库的完备性,随机生成1万个光轴指向,统计每个视场内的星颗数情况,如表 3所示。

      表 2  筛选后星库的全天识别情况

      Table 2.  The Recognition of all-sky by navigation star datebase

      轨道参数
      /(°)
      最少识
      别星数
      是否
      识别
      覆盖
      率/%
      轨道参数
      /(°)
      最少识
      别星数
      是否
      识别
      覆盖
      率/%
      赤纬0 7 100 赤纬–10 8 100
      赤纬10 7 100 赤纬–20 8 100
      赤纬20 6 100 赤纬–30 9 100
      赤纬30 7 100 赤纬–40 9 100
      赤纬40 4 100 赤纬–50 6 100
      赤纬50 8 100 赤纬–60 9 100
      赤纬60 8 100 赤纬–70 8 100
      赤纬70 9 100 赤纬–80 11 100
      赤纬80 14 100 100

      表 3  星库中星数的分布情况

      Table 3.  The distribution of star number in star datebase

      星数情况/颗 天区数 百分比/%
      小于5 2 0.02
      5~10 234 2.34
      10~15 2 324 23.24
      15~20 4 764 47.64
      大于20 2 676 26.76

      图  8  导航星库的二维分布图

      Figure 8.  2-D distribution map of navigation star table

      表 2表 3可以看出,在测试的全天球轨道中,本文方法制作的星库可以较好地满足全天球捕获的要求,且颗数大于10颗的天区占到97.64%,覆盖性较好。

      为了进一步验证本文方法制作的星库对提升姿态精度的有效性,在全天球随机生成了10万个光轴指向,计算得到理想姿态与实际姿态的误差,如图 9所示。

      图  9  理想姿态与实际姿态的误差

      Figure 9.  Error between ideal attitude and actual attitude

      从图中可以看出,利用本文方法制作的星库,删除了定姿误差较大的恒星,提升了星敏的定姿精度,极大地减小了XYZ 3个方向的理想姿态与实际误差。与经典的导航星库制作方法相比,利用本文方法制作的导航星库计算出的姿态精度更高,姿态输出也更为稳定,且不会出现如图 2所示姿态误差变化较大的情况。

    • 本文分析了恒星分布对定姿精度的影响,进而提出了一种基于恒星几何分布的星敏感器导航星库的制作方法。首先根据星敏感器的敏感极限以及跟踪半径的大小精简星库,计算每个视场内不同恒星分布组合的定姿误差,删除误差较大的星以进一步精简星库。同时对于视场内恒星分布的均匀性进行了优化,添加了星颗数为0的扇区,使得导航星库进一步完善,经全天球识别测试和天区星数统计结果表明,本文方法制作的星库完备性和均匀性较好,同时能够有效提升星敏感器的定姿精度。

参考文献 (10)

目录

    /

    返回文章
    返回