-
星库记录了恒星位置、星等、光谱等天文参数。由于编制星库使用仪器、观测条件和处理方法不一样,不同星库中恒星的颗数不一样,恒星的位置精度也会有差异。星库是设计制作导航星库的基础,也是决定星敏感器定姿精度的关键因素之一。表 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}}$ ,v1,v2,v3为3颗定姿星的惯性坐标系下的单位矢量,矩阵V一定程度上反映了恒星的空间几何分布[9-10]。假设惯性坐标i的原点为O,3颗星对应的单位矢量为Ov1,Ov2,Ov3。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颗星。第2步:随机生成光轴指向,并搜索当前光轴指向所在天区视场内的星,任取3颗星计算其星矢量组成的四面体的体积,分别计算每种组合中恒星矢量组成四面体体积。
第3步:基于恒星几何分布选取导航星。由2.1所述可知,用于定姿的恒星矢量空间构成体积越大,定姿误差越小,定姿精度越高。以(120°,30°)天区为例,在基础星库中该天区共有30颗星,任取3颗导航星,则有4 060种组合。图 3给出了姿态误差与四面体体积的变化趋势图,从图中可以看出,四面体体积越大,姿态角的误差越小,即姿态角度的精度越高。图 4给出了精简后的星库在(120°,30°)天区的星点分布图,图中红色的3颗星为空间几何分布体积最小的3颗星,可以看出这3星在平面坐标上基本处于一条直线上,其星矢量形成的空间四面体的体积最小,解算出的姿态数据误差也是最大的,可先删除定姿误差较大的3颗星。
第4步:在减少恒星分布对定姿误差影响后,再次遍历全天球,每个视场保留最亮的15颗,以进一步精简星库,精简之后星库共有导航星4 065颗。
第5步:为了保证导航星库在视场中的星分布尽可能均匀,本文将内切圆视场划分成4个面积相等的扇形区域,如图 4所示。遍历每个视场,统计每个扇形区域的星颗数,若出现有星颗数为0的扇形区,则将基础星库在该扇形区最亮的星添加到导航星库中。仿真结果表明,共有2 260个天区出现星颗数为0的区域,共添加126颗星,即星库星颗数共为4 191颗。图 5~7给出了扇区为0区域加添星的过程。图 8为星库的全天球分布图。
-
为了验证星库的有效性和完备性,对制作的导航星库进行全天球识别验证,输入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 从表 2和表 3可以看出,在测试的全天球轨道中,本文方法制作的星库可以较好地满足全天球捕获的要求,且颗数大于10颗的天区占到97.64%,覆盖性较好。
为了进一步验证本文方法制作的星库对提升姿态精度的有效性,在全天球随机生成了10万个光轴指向,计算得到理想姿态与实际姿态的误差,如图 9所示。
从图中可以看出,利用本文方法制作的星库,删除了定姿误差较大的恒星,提升了星敏的定姿精度,极大地减小了X、Y、Z3个方向的理想姿态与实际误差。与经典的导航星库制作方法相比,利用本文方法制作的导航星库计算出的姿态精度更高,姿态输出也更为稳定,且不会出现如图 2所示姿态误差变化较大的情况。
Establishment of a Navigation Star Database Based on Star Distribution
-
摘要:
为了进一步提升星敏感器的定姿精度,分析了恒星的空间几何分布对姿态精度的影响,提出了一种基于恒星分布的星敏感器导航星库制作方法。在分析对比了现有星库的基础上,选择了精度较高的依巴谷星库作为基础星库,并剔除了暗星和双星以精简星库。遍历全天球对每个视场的星库中的空间几何分布进行分析对比,删除定姿较差的3颗星。对星库进行遍历补偿,保证了星库中恒星分布的均匀性。实验结果表明,全天球轨道测试中,本文方法制作星库的全天识别率为100%,且星颗分布数10颗以上天区占到97.64%,具有良好的覆盖性,同时有效提升了姿态精度。
Abstract:In order to further improve attitude precision of star sensors, the influence of star spatial geometric distribution to attitude precision is analyzed and a method of establishing a navigation star database based on star distribution is proposed. On the basis of analyzing the current star database, Hipparcos is selected as the basic star database, while, dark-star and double-star are deleted in the star database, and then star spatial geometric distribution in every FOV (Field of View) is analyzed and compared while traversing all-sky, three stars which determine the worst attitude precision are deleted. Finally, Traversing compensationof star database is executed to ensure the uniformity of star distribution. The experimental results show that all-sky recognition rate is 100% by using the star database made by the proposed method in all celestial orbit tests, and the areas with more than 10 stars distribution r account for 97.64%. The porposed method is of good coverage, and can effectively improve the attitude accuracy.
-
Key words:
- navigation star table/
- star distribution/
- attitude precision/
- star sensor
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常见的几种星库
Table 1Several 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 表 2筛选后星库的全天识别情况
Table 2The 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 3The 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 -
[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.