粗差的定位与估计

1、粗差的定位与估计1实验目的今有一模拟的立体像对(f为100.5mm,摄影比例尺为1:12000,像点坐标量测中误差为2.8口m)，在y坐标中人为的加入46。1006。的粗差，请用本课程的知识将这些粗差一一找出。2粗差的定位方法2.1概述可靠性理论给出了平差系统发现粗差的能力和不可发现的粗差对平差结果的影响,同时也给出了检测和发现粗差的统计检验量。但是可靠性研究的一个最终目标是如何在平差过程中自动的发现粗差的存在,并正确的指出粗差的位置,从而将它从平差中剔除,这就是所谓的粗差定位问题。它不仅仅是一个理论问题,而更主要的是算法上的问题。2.2选权迭代法的基本思路选权迭代法的基本思想是：由于粗差未知

2、,平差仍从惯常的最小二乘法开始,但在每次平差后,根据其残差和有关其他参数,按所选择的权函数,计算每一个观测值在下一步迭代平差中的权。如果权函数选择得当,且粗差可定位,则含粗差观测值的权将愈来愈小,直到趋近于零。迭代终止时,相应的残差将直接指出粗差的值,而平差的结果将不受粗差的影响。这样便实现了粗差的自动定位和改正。该方法从下列最小条件出发：工pv2-minii式中权函数p(v+i)=f(v(v),.)(v=1,2,.)(2)ii对于目前已提出的各种权函数,可以按其内容区分为：(1)残差Vj的函数TOC o 1-5 h zp=f(v)(3)i1i标准化残差wi的函数p=f(w)(4)i1i方差估

3、计%的函数liP=f(o)(5)i1li若按其形式，则可以区分为幂函数和指数函数。已见到的一些权函数可列举如下：1)Lq迭代法(最小范数迭代法)该法由最小范数法vfmin演变而来，当q=1和q=0时，相应的权函数为1(6)(7)p(v+1)=(6)(7)iv(V)+ci1p(V+1)=i(v2)(v)+ci式中，c为一适当小的正数。引入它是为了在改正数为零时不至引起迭代求解的困难。由于此解法中的第一步仍为惯常的最小二乘平差，所以它与严格的最小范数解法有明显的区别。2)丹麦法丹麦法有多种形式的权函数，由Krarup等人提出的一种适合于摄影测量平差的权函数为：第一次迭代第二三次迭代p第一次迭代第二

4、三次迭代p=1ifr(vA4.4305p=exp-iB(8)第二三次迭代f第二三次迭代f(v、3.0TI0.05pH，xp-J1Kubik教授提出的权函数则为:f(v)=f(v)=exp1,MVpQcmf0vp二-0-cm其中，p0权系数，m为观测中误差而常数c=3。003)带权数据探测法此法由加拿大EL-HaKing提出并试验过，其权函数为-1,当w4.13p(V+1)=jw2i11，i当w4.13i式中，wiV式中，wiioqj0Yvii2.2粗差的估计可以根据残差向量估计残差的大小，若设V/i为第i个观测值粗差的估值，则有：(9)V/=-二iIi而该估值的精度o可按误差传播定律导出，其结

5、果为Vlioo=-3(10)V心(1)3算法步骤(1)列立误差方程，令各权因子初值均为1，即令w=w=w=1，12nW=L，则P(0)=P，p为观测权阵；(2)解算法方程BtPBX-BtPL=0，得出参数X和残差V的第一次估值:X=(BtPB)-1BtPL丫=BX-L（3）由V按V2=攻确定各观测值新的权因子，按p=pw构造新的Viiii等价权P（1）,再解算法方程BtPBX-BtPL=0，得出参数X和残差，的第二次估值；X（2）二（BtP（1）B）-1BtPLV（2）=BX（2）-L（4）由口2）构造新的等价权P（1），再解算法方程，类似迭代计算，直至前后两次解的差值符合限差要求为止；（5）

6、最后结果为V（k）=BX（k-1）-LTOC o 1-5 h z由于p=pw，而w=9，（v）=史，故随着P函数的选取不同，构iiiividvii成了权函数的多种不同形式，但权函数总是一个在平差过程中随改正数变化的量，其中w与v的大小成反比，v愈大，w、p就愈小，因此经过多次迭代，iiiii从而使含有粗差的观测值的权函数为零（或接近为零），使其在平差中不起作用，而相应的观测值残差在很大程度上反映了其粗差值。这样一种通过在平差过程中变权实现参数估计的稳健性的方法，称之为选权迭代法。4计算结果相对定向中的粗差检测结果点号粗差估值（um）4306478.9443061110111017.57-143

7、.45125.22-143.45125.2248733940613321162632235666350932338445166582975324277914534575066208531076412652762215161401768946162259543464636695-345.64-433.23145.8855.32-76.7697.6365.34311.37-94.34-54.22165.43-132.54255.32-154.43-143.43-276.53134.59-156.1276.6972.23214.44-264.97-254.12167.65-132.22-252.10

8、-138.29194.3470191913701919134361340947497336835211484866124553369513396282563134047305474738476513808307621818372183715121973621933168.43-70.33144.34-152.41-168.48217.34123.13-154.25233.75-140.57194.44187.45187.33-76.45-51.35143.46-94.41-224.74-157.23207.52-159.52-243.74-475.32-173.23165.55-47.73-29.3537.727065187.357065632656.4264725.435课程学习体会与建议知识的学习是枯燥艰辛而又富有挑战的，粗差探测问题一直是理论界一直探讨的热门话题。老师的谆谆诱导、同学的出谋划策是对我的帮助很大。在此我首先感谢任课老师袁修孝导师，没有袁老师的辛勤栽培、孜孜教诲就没有我课程实验的顺利完成。然后要感谢我的同学们以及宿舍舍友，在学习过程中他们也对我有很大的帮助。在课程的学习过程中，我对于测量平差中的误差处理有了更加深入的了解，也明白了如何进行粗差的定位与估计，了解到了可靠性理论以及外部可靠性与内部可靠性的区别，这些知识对于我以后的发展进步的帮助很大，