基于遗传算法的极短弧定轨(续)-紫金山天文台.PDFVIP

第57卷第2期 天 文 学 报 Vol.57 No.2 2016年3月 ACTA ASTRONOMICA SINICA Mar., 2016 doi: 10.15940/j.cnki.0001-5245.2016.02.005 基于遗传算法的极短弧定轨续 李鑫冉1,2,3 王 歆1,2 (1 中国科学院紫金山天文台南京210008) (2 中国科学院空间目标与碎片观测重点实验室南京210008) (3 中国科学院大学北京100049) 摘要 采用遗传算法解决极短弧定轨问题时, 由于遗传算法不同于经典方法的计算过 程, 野值剔除方法不再适用. 在遗传算法中通过在适值函数中采用不同损失函数实现了 稳健估计, 解决了极短弧定轨中的野值处理问题. 在遗传算法中不同损失函数的引入较经 典方法大大简化. 通过对多种损失函数的计算比较, 表明采用最小中值二乘(LMS, Least Median Square)和截尾最小二乘(LTS, Least Trimmed Square)估计可大幅度提高极短弧 定轨的稳健性, 具有极高的崩溃点. 关键词 航天器, 天体力学, 方法: 数值, 统计 中图分类号: P135; 文献标识码: A 引言 近年来空间目标的巡天观测采集了大量数据, 但不同于传统的跟踪观测, 巡天方式 采集到的弧段都十分短, 往往只有数十秒. 为了有效利用这些数据, 极短弧(TSA, Too [1] Short Arc)定轨越来越受到关注 . 由于短弧定轨问题本质上是病态的, 极短弧的病态性更为突出, 传统方法几乎无法 获得合理的结果. 文献[2]将遗传算法用于解决极短弧定轨问题, 采用遗传算法将定轨问 题转换为优选问题, 避免了求解线性代数方程, 从而克服了经典方法的一些不足. 遗传 算法也不同于经典优选算法, 对初值没有特殊要求, 不容易陷入局部极值, 因此是解决极 短弧定轨的一种有效手段. 文献[3]采用完全不同的变量和遗传算子, 也给出了一种遗传 算法求解极短弧定轨的具体方法. 上述研究给出了轨道计算方法, 但都没有涉及野值剔除的问题. 在实测资料的初轨 计算中, 野值问题是不可回避的, 观测数据中的少量野值会对计算结果产生很大影响, 甚 [4] 至是破坏性的 . 最常见的野值处理的方法是采用3σ准则, 在迭代过程中对于大偏差的 观测资料进行剔除, 使其不参与轨道计算, 逐步收敛后得到可靠的结果. 遗传算法不同于 2015-07-16收到原稿, 2015-08-21收到修改稿 ∗ 国家自然科学基金项目资助 †lixr@pmo.ac.cn ‡wangxin@pmo.ac.cn 182 天 文 学 报 57 卷 经典计算过程, 无法采用这样的野值剔除过程. 同时由于野值的判断与剔除需要在轨道 计算过程中逐步确定, 而剔除本身影响轨道计算结果, 因此野值剔除并不容易. 野值处理的另一种手段是采用稳健估计方法, 文献[5]在轨道计算中采用了M估计, 但其求解方法比较复杂, 特别是需要可靠初值, 可靠初值获取本身就是一个难题, 受此 局限方法的崩溃点不高. 此外方法中引入了额外的参数, 参数确定尚无可靠手段. 文 献[5]中根据历史数据的统计特性给出参数, 然而对于极短弧数据而言, 由于数据量少, 这种做法并不合适. 文献[6-7]采用了最小一乘(LAD, Least Absolute Deviation)估计, 由于最小一乘估 计无需额外参数, 避免了参数确定问题, 但其求解采用了线性规划方法, 该方法和一般计