验拟合后剩余残差的统计检和}一,r待估参数的选择.pdfVIP

第 40卷 第 4期 天 文 学 报 Vol40，No4 1999年 11月 ACTA ASTRONOMICA SINICA N0v ．1999 拟合后剩余残差的统计检验与 } 一，r 待估参数 的选择 IJ f 胡小工 黄 城 廖新浩 (中国科学院上{孽天文台 上{孽200030) 摘要 在空闻技术的实际资料处理中，经常遇到不同持估参数组音解算的RMS水平大 翦相当而估值结果相差甚远 的情况 车文提 出采用残差统计捡验 的方法 ．结音 RMS的大小 ． 来判断精密定轨和参数解算 的好坏 将拟 台后 的剩余残差序列看作一十 随机变量样本 ，并假 设其满足正态分布，采用高阶矩和皮尔逊]c2统计量来度量残差序刊偏离正态分布的程度 根 据统计检验 的结果 ，提 出选择台适 的持怙参数 的标准 ， 美羹词 统计检验 。残差序列，侍估参数 薤 迄 f 一 — — — — 一 — — — — 一 — — — 引言 残差平方和的平方根 即RMS，是评价精密定轨和参数解算 的重要统计量 ．在实际资 料处理中，它几乎 是说 明解 算好 坏 的唯一指标 ．倒如 ，在 美 国 Texas大学 空 间研 究 中心 CSR发展 的激光测距资料处理软件 UTOPIA中，对解算是否收敛提 出的判据是 J：(1) 本次迭代与上次迭代解算 的初始历元卫星位置差小于 lena；(2)lRMS—LRMSl0．01 RMS．其中RMS=√∑ ，(P一 )，Pi0， 分别是第f个观测值与计算值，，是第i次 厂——I_ —— 观测的权重．LRMS= ／∑ 。f △ 1是所谓线性预报残差平方和的平方根，△z是待估 Y 、d ， 参数的改正量．LRMS由GivensGenderman正交变换算法直接给出…．第一个判据是为了 防止 由法方程严重病态引起 的假收敛而 引进的，而第二个判据 的意义是要求待估参数误 差的高阶项 不再对 RMS有显著 的贡献 ． 在通常的残差分析理论中，残差被看成是一个时间序列 ．通过对其进行滤波 、谱分析 或更复杂 的协方差分析 (这时需残差偏导文件)，人们期望能发现资料处理 中采用 的力学 模型和测量模型所存在 的(系统性)误差 ．这种误差反 映在残差序列 中就是一个较 明显 的 变化趋势分量 ．例如 ，地球 引力场模型中低阶带谐项系数 的时变就是通过这种残差分析方 法发现的_2j．不 同于上述作法 ，本文假定 。在解算较好 的情况下 ，拟台后剩余残差 中不应 再包台明显的信号(因为力学模型和测量模型已足够准确)．因此可 以将残差看成随机变 1998 08—19收到原稿、1998 10 26收到謦改藕 ·国家 自然科学基盘 (批准号 和 中科院 “九五 重大壹助(批准号 K1951—304) 天 文 学 报 40卷 量的一个样本而不再认 为不 同历元 的残差 间存 在显著 的动 力学联 系．本 文更进 一步假定 残差应大致服从零均值 的正态分布 ，即高斯 白噪声 有理 由相信 ，相反的情况 ，即残差严重 偏离高期 白噪声 的情况 ，即使不算解算失败 ，至少也给理解解算结果及其精度带来不便 ． 注 意到常用 的统计量 RMS就是 随机样本 的标 准差 ，而本文试 图说 明，残 差 的高阶矩 和其 他统计量对评价解算结果也有重要的作用 ，甚至在某些情况下是关键 的作用． 促使 我们提 出残差统计检验的原 因是在利用 UTO