初轨计算的参考矢量法.pdfVIP

第38卷 第4期 天 文 学 报 V0l 38．No4 1997年 11月 ACTA ASTRONOM1CA SINICA Nov．1997 L ． 彤 一3 初轨计算的参考矢量法 尸／弓，初轨计算的参考矢量法 ， j I’ 重 泣 壅： 吴 连 大 (中国科学院系坑科学研究所 北京 00O80) (中国科学院紫金山天文台 南京 210008) 摘 要 本文从最优估计角度研究短弧资料定初轨问题，提出一种最小二乘估计算法，同时提出 一 种符合实用的参考矢量法．仿真结果显示，参考矢量法克服了现有初轨算法所承受病态影 响，把定轨精度提高到接近最佳水平． 糕词盟1监， 皿 绔 迤 l／引 吉 在文[1】中，我们通过病态分析与误差分析，看到由短弧观测引起的法方程的病态性， 对现有初轨计算方法的定轨精度有不同程度的影响；在某些情况下，影响还比较严重．本 文将从最优估计的角度研究初轨计算的最优算法，提出了一种最小二乘估计算法，并给出 了一种实用的参考矢量法，该算法克服了观测误差对法方程病态的影响，将定轨提高到接 近最优估计的水平． 按照统计理论，对一般的参数估计问题，其估计误差的方差，存在Cramer-RaoT界． 对线性模型，采用最优加权的最小二乘估计，达到了这个下界．但对初轨这样的非线性系 统，无一算法可达到这个下界，只能接近这个下界．本文的算法，其估计误差的方差就十 分接近这个下界．有关Cramer-RaoT界及本文算法的精度的比较，我们将另文给出． 2 初轨的最小二乘估计算法 与文【1]一样，本文仅讨论拉普拉斯类型的方法，即采用下述的二体动力学模型 = r0+g ／o， - 其中I0 为历元时刻to天体的坐标和速度矢量． l r =p I +Rl 譬 1996年7月I2日收到原稿，1996年I2月i3日收到修改藉 ． 国家自然科学基金项目 ， 中科院紫盒山天文台客座研究员 354 天 文 学 报 38卷 式中rl和 ．分别为历元时刻 天体和观测站的坐标矢量，P．为站星之间的斜距，t为站到星 的方向矢量． 假定只有赤经和赤纬的( ，69观测量，且假设观测误差(d珥)，(d6 )相互独立，均值 为零，标准差分别为 与 ，在下述条件下： (3) 由观测量( ，d 求r0， 的最小二乘估计，应当是满足下述残差平方和达到极小 O=∑cos (晦～ +∑((1— ) ， (4) 式中( ， 为计算值，由估计 、 通过(1)、(2)H-g，其中 =