欢迎来到天天文库
浏览记录
ID:43842639
大小:211.37 KB
页数:3页
时间:2019-10-15
《三种粗差检测方法的比较及分析》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、第20卷·专辑·淮海工学院学报(自然科学版)Vol.20S.I.2011年12月JournalofHuaihaiInstituteofTechnology(NaturalScienceEdition)Dec.2011DOI:10.3969/j.issn.1672-6685.2011.S0.050*三种粗差检测方法的比较及分析李海军,何丽媛,李岩(连云港市勘察测绘院有限公司,江苏连云港222001)摘要:阐述了粗差检测的原理,结合具体算例,运用三种不同的数学模型,对同一算例进行了粗差检测,比较了数据探测法、QUAD法和LEGE法等三种粗差检测方法的优缺点。关键词:最小二乘法;粗差;数据
2、探测中图分类号:文献标识码:A文章编号:1672-6685(2011)S0-0129-030引言一种较为直观且简单可行的方法,它不仅能确定k个粗差的位置(k≤r-1,r为多余观测数),而且可经典的平差方法是按最小二乘法来进行的,观以同时求得各个粗差的数值大小,故称为“粗差定测数据难免会存在粗差,如果粗差不排除,所得结果位定值法”,简称LEGE法。则不是最优的。近年来,随着各种电子观测仪、GPS等现代测量技术的应用,对外业观测大量数据中的2算例及分析粗差进行探测显得尤为重要。虽然在数据处理方面的各种理论已经相当成熟,但当观测数据由于含有图1为一水准网,高差观测数n=19,必要观测粗差而
3、偏离正态分布时,估计结果与真实值偏离很数t=10,多余观测数r=n-t=9。分别模拟了一组远,精度也很差,这样,即使是经典的最小二乘法对含有1、2、3、4个粗差的观测真误差,见表1中第4、5、粗差也无能为力。因此在进行最小二乘平差前,对6、7列。以L珟表示高差真值向量(单位为mm),S表示观测数据进行认真检测,定位并修正粗差观测,以提各水准路线距离(单位为km),ε表示真误差向量,见高参数估计的准确度和精度,一直是测绘工程领域表1。(本例原始数据取自文献[1]中的例2)的热点问题。本文运用三种不同的数学模型,对同一算例进行了粗差检测,比较了三种方法的优缺点。1粗差检测原理巴尔达教授从
4、已知单位权方差出发导出了以服从正态分布的标准残差为统计量的数据探测法,这一类方法的主要问题是不能自动准确地给出粗差的定位。欧吉坤等人从真误差入手,利用真误差与观测值之间的解析关系,提出拟准观测的概念。借鉴图1水准网拟稳平差思想,附加“拟准观测的真误差范数极小”的条件,解决了关于真误差的秩亏方程组求确定解2.1数据探测法的问题。这种方法称为“拟准检定法”,简称QUAD由模拟四个粗差数据计算而得,结果见表2,3。法。於宗俦等人提出的多维粗差同时定位定值法是该方法每次只能测出一个粗差,若要测另一粗差,*收稿日期:2011-10-11;修订日期:2011-11-10作者简介:李海军(1982
5、-),男,江苏连云港人,连云港市勘察测绘院有限公司助理工程师,主要从事城市测绘生产和研究。130淮海工学院学报(自然科学版)2011年12月须先剔除所发现的粗差,重新计算,逐次重复。此例第检测出的粗差与实际值有较大差别,一个粗差可能需一次计算得出L2含有粗差,对其进行改正后再计算,要若干次的检验改正才能剔除干净。此例中探测到的依次发现L7,L11,再次探测时重新发现L7仍有粗差。所有粗差值:ε2=11.662857,ε7=-(4.865570+2.907684)=-7.773254,ε11=4.902005;与实际表1数据计算结果模拟粗差(分别为15.9,-11.6,12.4)相差较
6、大。No.L珟/mmS/kmε/mmε/mmε/mmε/mm由于每次只考虑一个粗差,未顾及各改正数相关1171302.9680.50.50.50.5性,方法并不是绝对可靠,对某些粗差无法发现;同27513.3410.90.915.9338637.635-1.4-1.4-1.4-1.4时,也可能对某些没有粗差的观测值造成误判,当作4117336.547-0.5-0.5-0.5-0.5粗差剔除。在闭合环11—10—12—17中,L11和L175100117.320-0.8-0.8-0.8-0.8同时存在符号相同的粗差,造成当改正L11后,再次6309268.9221.11.114.11.
7、17409386.218-1.6-1.6-1.6-11.6检验得出不再有粗差的结论,以致无法检验出L17中881197.7310.2-8.20.20.2的粗差。同时,L3的u值却越来越大,这在较严格的时9156885.623-9.3-0.3-0.3-0.3候还会认为其存在粗差。这个过程见表3。10149372.3290.80.812.80.811110734.7260.40.40.412.42.2QUAD法12134363.915-0.7-0.7-0.7-
此文档下载收益归作者所有