半参数估计的自然样条函数法.pdfVIP

第29卷第5期 武汉大学学报·信息科学版 VoI．29No．5 2004年5月 GeomaticsandInformationScienceofWuhan 2004 University May 文章编号：167l一8860(2004)05。0398—04 文献标识码：A 半参数估计的自然样条函数法 吴 云1 孙海燕1 马学忠1 (1武汉大学测绘学院，武汉市珞喻路129号，430079) 摘 要：用补偿最小二乘原理，得到了参数和非参数分量的惟一解，并通过模拟计算，对半参数回归模型和参 数模型的计算结果进行了比较。结果表明．半参数回归方法能较好地将观测值中具有连续光滑特性的系统误 差分离出来。 关键词：半参数回归；系统误差；自然样条函数；平滑参数；补偿最小二乘原理 中图法分类号：P207．2 在数据处理中，观测值中的系统误差或参数 必要观测数；4为行维偶然误差向量；B为列满秩 模型与实际情况的偏差是影响参数估计精度的主 设计矩阵；非参数向量S一[s，，S。，…，S。]7，是描 要因素。例如，在GPS长基线测量中，双差观测 述系统误差的挖维向量，其中S，是某些特定量t， 值中仍然存在电离层和对流层折射误差的影 的函数。式(1)的误差方程为： 响…，如果在参数模型中不顾及这些系统误差，将 v一黜+§一L (2) 会严重地影响到平差精度，甚至会导致错误的结 式中，§一(j。，j。，…，；。)1，即系统误差s的估值。 果『2】。观测值中的系统误差往往随着某些因素的 设s(f)为区间[f。，t。]上的自然样条插值函 连续变化，有时其函数特性十分复杂，难以用少量 数，t，(i一1，2，…，门)为节点，且tl￡，f。。s(￡)满 参数准确地表达。自然样条函数[31属于无穷维函 足插值条件： 数空问的光滑曲线．用来描述连续变化的模型误 s(t，)一s， (3) 差有其独特的灵活性和简洁性。本文用自然样条 可以找到惟一的满足上述条件的自然样条插值函 函数来表示随时问连续变化的系统误差，即用参 数[3]，因此，观测方程为： 数在表示观测值中可线性化部分的基础上，增加 用自然样条函数表示随时间连续变化的系统误差 1．2补偿最Ih--乘原理及其解 部分，也就是非参数分量部分，从而构成半参数模 补偿最小二乘原理是在最--乘法的目标函 型。为了求得参数与半参数分量的惟一解，引入 数上增加一个包含非参数部分的补偿项，即 least 了补偿最小二乘原理(penalisedsquares， PI，s)14】。本文的算例验证了自然样条函数的半 sum=y(L，一Bx—s(￡，))2+ i=l 参数模型和补偿最小二乘原理在数据处理中的可 r， 行性。 aI”(s”(f))2dt—mi