2009-现代测量平差与半参数估计[1].ppt

现代测量平差与半参数估计 陶本藻教授武汉大学测绘学院地球空间环境与大地测量教育部重点实验室 函数模型是描述观测量与待求参数间的数学函数关系的模型 是确定客观实际的本质或特征的模型 随机模型是描述平差问题中的随机量 如观测量 及其相互间统计相关性质的模型 1 测量平差数学模型 经典平差模型 R A UR Q nX为非随机参数 经典平差公式 R A t ud u tR Q nX非随机 秩亏自由网平差 R Q g nR A uX非随机 具有奇异协方差的平差模型 配置 似合推估 模型 广义测量平差原理 最小二乘统一理论 2 广义高斯 马尔柯夫模型 最小二乘统一理论 Rao在文中提出的最小二小乘准则是 R A tR D gn g tu tX非随机 陶本藻 刘大杰 5 5 1990 从奇异正态分布的密度函数 3 平差系统的模型误差 模型误差分为函数模型误差和随机模型误差两类最小二乘平差参数X的估值具有最优无偏性 单位权方差的估值具有无偏性和渐进最优性 这些良好的统计性质都是基于模型误差不显著的情况 但在实际平差系统中 由于种种原因的建模近似 例如非线性观测方程的线性化 未顾及或近似考虑某种系统误差影响 观测值的先验协方差阵不尽合理等原因都会造成函数模型和随机模型产生误差 模型近似在回归拟合模型中则更为突出 4 测量中的半参数回归模型及估计方法4 1 测量中的半参数回归模型和平差估计准则 补偿最小二乘原理 半参数回归模型 平差估计准则 是一个适当给定的矩阵 称为正规化矩阵 是一个给定的纯量因子 在极小化过程中对和起平衡作用 称为平滑因子 4 2 正规化矩阵正定时半参数模型的估计 由拉格朗日乘数法 构造函数 参数估计 非参数估计 4 3 最小二乘核估计与附加系统参数平差 附加系统参数估计 核估计 偏核光滑估计 4 4 最小二乘配置与补偿最小二乘 半参数 配置 4 5 正规化矩阵的一个选择 如果观测值是在时刻得到的一个时间序列 若假设相邻时刻的模型误差与的差别不应太大 因此可令 4 6 确定的周江文方法时间序列系统误差模型设为 例如的协方差阵为 5 半参数估计的自然样条函数法 5 1半参数回归模型中的自然样条插值函数 设为区间上的自然样条插值函数 为节点 且 满足插值条件 满足上述条件的插值函数中 自然样条函数是最光滑的 Green Silverman1987 5 2补偿最小二乘原理及其解 Fessler 1991 G矩阵为严格的对角占优矩阵 即正定矩阵 设 是半正定矩阵 补偿最小二乘准则项表示为 按照求条件极值的拉格朗日乘数法 构造函数 得到法方程为 法方程系数阵满足 时可逆 5 3平滑参数的确定 交叉核实法 去除时刻的观测值 也就是以作为观测值 对于给定的平滑参数得到观测值估计和系统误差估计 通过三次样条函数 可以内插出时刻的系统误差估计 记为 从而得到时刻的观测值估计如果平滑参数选取的比较合适的话 那么 若选取的值 使 Green Silverman1994 等价模型 6 半参数模型假设检验 1 假设检验 2 1 2 联合平差拒绝域 7 时间序列分析与残差分析 7 1 AR P 模型 求解 平差求解可得 F检验定阶 给定显著性水平 例如 0 05 0 01 查F分布表可得临界值若 则是合适的 检验不显著 若 则是不合适的 检验显著 7 2 半参数模型 模型 Y X S A式中 A S X 7 3 两模型的比较 分别用AR p 模型和半参数模型对香河台东西向长约700米的短水准观测数据进行计算比较 AR p 模型求得的平差值图 半参数模型求得的平差值图 分别用AR p 模型和半参数模型求得的平差值比较图 上面的图中虚线表示用AR p 模型求得的平差值 实线表示半参数模型求得的平差值从图中可以看出用AR p 模型求得的平差值的变化较大 而且快 曲线也不光滑 而用半参数模型求得的平差值的变化小 而且慢 曲线也非常光滑 由此可以