自适应拟合推估及其应用.pptVIP

* * * 地理空间信息学院 2013.01 自适应拟合推估及其应用 * 常用于重力场的拟合与推估(Tscherning 1978)，应用非常广泛。 实质是对公共点误差进行拟合，而对非公共点误差进行推估。 关键是信号向量协方差函数的确定。 拟合推估 一旦选定协方差函数，则通过已测点信号即可拟合出信号间的先验协方差，而且在拟合推估中视为已知常量。 1 背景 * 协方差函数选择得不合理 用于拟合的观测信息存在异常 信号向量与观测向量的协方差应协调一致 存在问题 后果：造成协方差函数的不可靠，引起拟合推估解的偏差。 从随机模型误差的角度看，可部分地通过方差因子进行调整，如此，可以采用方差分量估计的方法重新调整信号向量与观测向量的先验权比 1 背景 在大地测量数据处理中，方差分量估计获得了广泛研究 * 最小范数二次无偏估计（MINQUE）(Rao 1971) 最优不变二次无偏估计（BIQUE）(Caspary 1987； Sjoeberg 1984 )、 极大似然估计（Koch 1986） 具有约束的极大似然估计（Searle et al. 1992） 最小二乘方差分量估计（Teunissen and Amiri-Simkooei 2007） 1 背景 函数模型 目标函数 * 拟合推估解 2 拟合推估模型 由随机信号协方差函数计算得到 得到较困难（黄维彬，1992 ） 、 间通常不协调 取决于 拟合推估分布解 2 拟合推估模型 目标函数 * 自适应拟合推估解 3 自适应拟合推估模型 关键是构造自适应因子 3 自适应拟合推估模型 基于极大似然方差分量估计的拟合推估 基于Helmert方差分量估计的拟合推估 应用方差分量估计来构造自适应因子 * 4 自适应拟合推估模型高程转换算例 以某地区362个GPS/水准重合点的大地高到正常高的转换为例。在转换过程中，用332个点参加计算，另均匀选择30个点作为外部检核。 高程转换算例 方案一：二次多项式模型最小二乘解 方案二：标准拟合推估解 方案三：基于极大似然方差分量估计的自适应拟合推估解 方案四：基于Helmert方差分量估计的自适应拟合推估解 * -0.329 -0.329 -0.339 -0.329 -0.329 -0.375 -0.385 最小差值 四 三 二 四 三 二 一 方案 0.085 0.296 0.085 0.296 0.088 0.296 0.1 0.085 0.296 0.085 0.296 0.100 0.307 0.105 0.304 1 RMS 最大差值 方差 因子 四 三 二 四 三 二 一 方案 -0.190 -0.194 -0.195 -0.190 -0.194 -0.197 -0.209 最小差值 0.068 0.102 0.070 0.111 0.076 0.128 0.1 0.068 0.102 0.070 0.111 0.096 0.171 0.102 0.197 1 RMS 最大差值 方差因子 公共点残差（m） 外部检查点检验结果（m） 自适应拟合推估结果 结论： 方差因子取不同值时，拟合推估解差异较大，自适应拟合推估解非常稳定。 4 自适应拟合推估模型高程转换算例 * 模型 准则 参数解 5 自适应拟合推估水平速度场模型 1070 GPS 点，其中29点位CORS站。956 点用来计算欧拉矢量和局部形变信号，29个CORS站和56复测站作为检核点。利用Hirvonen函数拟合信号向量协方差。 水平速度场建立算例 5 自适应拟合推估水平速度场模型 方案一： LS估计； 方案二：拟合推估； 方案三：自适应拟合推估； 方案四：Hardy多面函数拟合 Longitude Latitude rotation 0.234 50.6 -112.3 0.0031558 -0.0023984 -0.0009836 NNR-NUVEL1A 0.300 37.7 -88.1 0.0032090 -0.0041396 0.0001367 Scheme 3 0.302 37.2 -88.4 0.0031925 -0.0042049 0.0001159 Scheme 2 0.300 43.2 -96.6 0.0035848 -0.0037950 -0.0004415 Scheme 1 ω(°/myr ) Φ(°) λ(°) (rad/myr) (rad/myr) (rad/myr) 中国大陆欧拉矢量推估结果 观测点残差统计（mm/y） 2.25 0.00 -9.12 24.54 2.79 0.00