稳健估计.ppt

是用来判断估计统计量对异常值敏感程度的指标，反映了在不同位置上异常数据对估值所造成的相对影响的大小 一定要掌握等价权的思想！ 该方法的关键是确定等价权。选择不同的 Rou函数，就构成不同的权函数，通常权函数是一个在平差过程中随改正数变化的量，经过多次迭代，从而使含有粗差的异常观测的权函数为零（或接近于零）。 这样一种通过在平差过程中的变权实现参数估计的稳健性的方法，称之为选权迭代法。 数学中的稳健估计，假定观测值是等权的，而在测量中引入权阵的稳健估计理论由周江文教授提出，称之为等价权抗差估计。 当所有改正数均在-c和c之间时，Huber估计就是经典的最小二乘估计。而当改正数大于c时，改正数越大，权越小，从而对参数估计的影响越小 将改正数分了四段 参考文献（部分） 杨元喜. 抗差估计理论及其应用. 八一出版社 1993 黄维彬. 近代平差理论及其应用.解放军出版社.1992 李德仁，袁修孝.误差处理与可靠性理论. 武汉大学出版社. 2002 周江文，欧吉坤，杨元喜等. 测量误差理论新探. 地震出版社.1999 周江文. 经典误差理论与抗差估计. 测绘学报，18（2），1989：115-120 刘大杰，陶本藻. 实用测量数据处理方法. 测绘出版社，2000 周秋生. 测量控制网优化设计.测绘出版社1992 陶本藻. 测量数据统计分析. 测绘出版社1992 陶本藻. 测量数据处理的统计理论和方法. 测绘出版社2007 武汉大学测绘学院测量平差学科组. 误差理论与测量平差基础. 武 汉大学出版社. 2003 编程与练习 * * 第四章　　粗差估计（数据探测和稳健估计） 、　模型误差，粗差 函数模型 随机模型 模型：描述客观实际的数学式 所建模型与实际模型之差——模型误差 最小二乘具有良好的均衡误差特性，不具备抗粗差干扰能力， 对粗差相当敏感。估计失实；例如：x真值为5，进行8次观测： L1=5.001,L2=5.002,L3=4.998, L4=4.993,L5=5.001,L6=5.008, L7=5.500,L8=4.997,采用最小二乘（平均值） 2　．切尾法：去掉一个最小值和一个最大值，期于平均。 其他估计法： 中位数法 观测值按大小排列，取中，偶数取中间两数平均： 4.993，4.997，4.998，5.001，5.002，5.008，5.500， 随机误差——最小二乘 系统误差——附有参数的 粗差 粗差是一种异常大误差 不能用最小二乘 其他方法见“实用测量数据处理”P52例1 粗差估计：1）识别法——数据探测 探测粗差数据的位置，剔除，保干净数据 2）调节法——稳健估计 采用抗干扰性强的估计法，削弱粗差对平差结果的影响 识别——把粗差归入函数模型 调节——把粗差归入随机模型 平均漂移模型 方差膨胀模型 无模型误差 可靠性：多余观测、误差、可靠性的关系。 如果粗差、系统误差不能在平差正确发现、消除，损害结果，不可靠 可靠性任务：1、研究平差系统发现、区分模型误差的能力，以及不可发现、 区分的模型误差对平差的影响 2、平差过程中发现区分模型误差的方法 不可靠 有一定可靠行，不可定位 可靠性高，可定位 可靠性：衡量成果可靠成度的指标 粗差对残差的影响 残差（改正数） ：发现粗差的重要信号 观测值残差大 粗差（可能） 设： 最小二乘： 二、残差理论（可靠性理论） 平差的几何条件反映了观测误差对改正数的作用程度。 的性质 2)幂等阵的秩等于迹： 3)R为降秩阵 4)r为多余观测值 首先讨论 的性质。 1)R为幂等阵 1、 某一观测值改正数 ，粗差 对观测值的影响。 其中， 受所有观测值的影响，大小取决于 ， 。 = - 某一观测值的改正数 受所有观测值的影响，大小取决于 ⒈ = - 2、 某一观测粗差 对所有改正数的影响，观测值含有的影响为： ， 。影响大小由 决定。 3、 对自身的影响 多余观测分量 ， ，多余观测分量为零，可靠性为零。 ， ，纯多余观测，对求参数无作用。 一般 ： ，粗差不能完全反映在本观测值 所以难以正确定位，只能用假设检验法 不能简单地直接由残差来检测粗差。 1、基本思想（Baarda数据探测法） 通过假设检验的方法 改正数 三、数据探测法 在平差中检测