M—split估计与LS关系及在粗差探测中的应用.pdfVIP

第2期 矿 山 测 量 No．2 2013年 4月 MINESURVEYING Apr．2013 doi：10．3969／j．issn．1001—358X．2013．02．20 M—split估计与 LS关系及在粗差探测中的应用木 李 娟 ，王爱生 ，王鑫渊 (1．中国矿业大学(北京)煤炭资源与安全开采国家重点实验室，北京 100083； 2．江苏师范大学测绘学院，徐州 221116；3．徐州矿业有限公司，徐州 221000) 摘要：文中介绍 了一种新的平差方法——M—split估计，针对观测值 中无粗差的情况，分析 了该方法 与最小二乘法(Ls)的关系，并给出了两组结果的融合算法。如果观测值 中有粗差，则通过算例说明 该方法能够进行抗差估计 ，同时探测粗差 。 关键词 ：M—split估计 ；抗差 ；最小二乘平差；粗差 中图分类号 ：P207．2 文献标识码 ：B 文章编号 ：1001—358X(2013)02—0057—03 自从 Huber提 出抗差估计 的概念…以来 ，抗差 于是 ，估值 和 就是下面最优 问题的解 估计的理论和应用得到迅速 的发展 J。抗差 M一估 maxK (y；O，0B)=Ks，B(y；O，0B) (2) 计也早已被应用于大地测量平差计算 中 。杨元喜 ， B 等人以M一估计原理为基础发展 了抗差卡尔曼滤 这就是 M—split估计的基本原理。 波 、抗差 配置 、方差 一协方差矩 阵 的抗差估 使用 M—split估计进行数据处理时发现 ，该算 计 、抗 差 贝叶斯 估 计 。在建 立 中国最新 的 法存在的一个关键 问题是，它使用牛顿迭代法将观 CGCS2000坐标系统时对天文大地 网的平差也采用 测值分别按不同的准则进行计算，得到两个不同的 了抗差处理的方法 。 结果 ，而不管观测值 中有没有粗差。本文首先推导 2009年 ，Z．Wisniewski提 出了一种 同时估计正 了M—split估计与最小二乘估计的关系 ，然后提出 平差结果取舍 的判别准则 ，目的是对观测值 中有无 常参数 和系统误差或粗差 的方法 [91。这种方 粗差进行 区分 。如果没有粗差，则要对两组结果进 法基于下面的假设 ：观测 Y；值属于不 同的随机变量 行融合。如果有粗差，则只取其 中单位权 中误差小 中的一种 (以二元分布为例)，ya—P P ∈或 ye～ 的一组结果。 P。。∈P ，观测值的分布是二者的复合 。新理论 的基 本思路是给每一个观测值赋予一个属性 ，称为初级 1 平方型 M—split估计与最dx--乘的关 系 分解能(elementarysplitpotentia1)，该属性描述 了观 令 ：V=AX一 (3) 测值是 —P。P ∈P 或 YB～P0BPPB的可能性。 按照[9]，平方型M—split估计的最优化 问题可 假设初级分解能与概率 P =P (Yi；0)或 PB；PB