自适应拟合推估及其应用

1 地理空间信息学院2013 01 自适应拟合推估及其应用 2 常用于重力场的拟合与推估 Tscherning1978 应用非常广泛 实质是对公共点误差进行拟合 而对非公共点误差进行推估 关键是信号向量协方差函数的确定 拟合推估 一旦选定协方差函数 则通过已测点信号即可拟合出信号间的先验协方差 而且在拟合推估中视为已知常量 1背景 3 协方差函数选择得不合理用于拟合的观测信息存在异常信号向量与观测向量的协方差应协调一致 存在问题 后果 造成协方差函数的不可靠 引起拟合推估解的偏差 从随机模型误差的角度看 可部分地通过方差因子进行调整 如此 可以采用方差分量估计的方法重新调整信号向量与观测向量的先验权比 1背景 在大地测量数据处理中 方差分量估计获得了广泛研究 4 最小范数二次无偏估计 MINQUE Rao1971 最优不变二次无偏估计 BIQUE Caspary1987 Sjoeberg1984 极大似然估计 Koch1986 具有约束的极大似然估计 Searleetal 1992 最小二乘方差分量估计 TeunissenandAmiri Simkooei2007 1背景 函数模型 目标函数 5 拟合推估解 2拟合推估模型 由随机信号协方差函数计算得到得到较困难 黄维彬 1992 间通常不协调 取决于 拟合推估分布解 2拟合推估模型 目标函数 7 自适应拟合推估解 3自适应拟合推估模型 关键是构造自适应因子 3自适应拟合推估模型 基于极大似然方差分量估计的拟合推估 基于Helmert方差分量估计的拟合推估 应用方差分量估计来构造自适应因子 8 4自适应拟合推估模型高程转换算例 以某地区362个GPS 水准重合点的大地高到正常高的转换为例 在转换过程中 用332个点参加计算 另均匀选择30个点作为外部检核 高程转换算例 方案一 二次多项式模型最小二乘解方案二 标准拟合推估解方案三 基于极大似然方差分量估计的自适应拟合推估解方案四 基于Helmert方差分量估计的自适应拟合推估解 9 公共点残差 m 外部检查点检验结果 m 自适应拟合推估结果 结论 方差因子取不同值时 拟合推估解差异较大 自适应拟合推估解非常稳定 4自适应拟合推估模型高程转换算例 10 模型 准则 参数解 5自适应拟合推估水平速度场模型 1070GPS点 其中29点位CORS站 956点用来计算欧拉矢量和局部形变信号 29个CORS站和56复测站作为检核点 利用Hirvonen函数拟合信号向量协方差 水平速度场建立算例 5自适应拟合推估水平速度场模型 方案一 LS估计 方案二 拟合推估 方案三 自适应拟合推估 方案四 Hardy多面函数拟合 中国大陆欧拉矢量推估结果 观测点残差统计 mm y 2 5顾及地球物理模型的形变信息融合 85个检核点误差统计 mm y 不同方案误差分布 14 5自适应拟合推估水平速度场模型 CTRF2000形变场 欧拉矢量拟合后的剩余误差 拟合推估后的剩余误差 自适应拟合推估后的剩余误差 15 CTRF2000框架存在南东运动不同地区地壳运动规律不同 西部变形大 东部变形小 且呈右旋趋势欧拉矢量仅反映中国大陆的 整体运动部反映局部变形不同方法解得的欧拉矢量存在较大差异 考虑局部变形特征的欧拉矢量解与NNR NUVEL1A模型也有较大差异拟合推估可以充分利用地球物理信息 观测信息以及