粗差的定位与估计.docx

精品文档 精品文档 芁 螅芅 芆肈 蒁芀 莈蚃 蒇薆 肅芆 薀衿 蝿薃 腿蒃 袄袈 薀蒈 膀蒃 蚇螄 薃膆 蚀莁 薁莄 荿羅 蚆蚈 螀羀 螈芄 螆膆 莅羆 袀蕿 腿芄 薈膃 膃薈 芄腿 蕿袀 羆莅 膆螆 芄螈 羀螀 蚈蚆 羅荿 莄薁 莁蚀 膆薃 螄蚇 蒃膀 蒈薀 袈袄 蒃腿 薃蝿 衿薀 芆肅 薆蒇 蚃莈 芀蒁 肈芆 芅螅 螃芁 蚁莁 蒅芃 肄芇 袃袀 螈芁 膇蒄 袂薅 袃蒀 膈膁 蚅蚆 袅肈 羃聿 蕿肁 莇蚃 蚄莆 肃蚈 羀羁 袅袄 蒃羈 膂蒁 莁羁 薇膅 蒆芆 节肀 薈袁 艿莆 芅螈 莂虿 罿螂 蚇莃 羄肆 莂芈 莀蚂 葿薄 螃薈 蒂袇 螁薂 袆螅PAGE PAGE7 精品文档 PAGE 粗差的定位与估计 实验目的 今有一模拟的立体像对（ f为100.5mm，摄影比例尺为1:12000，像点坐标量 测中误差为2.8μm），在y坐标中人为的加入4бo~100бo的粗差，请用本课程的 知识将这些粗差一一找出。 粗差的定位方法 2.1概述 可靠性理论给出了平差系统发现粗差的能力和不可发现的粗差对平差结果 的影响，同时也给出了检测和发现粗差的统计检验量。 但是可靠性研究的一个最 终目标是如何在平差过程中自动的发现粗差的存在，并正确的指出粗差的位置， 从而将它从平差中剔除，这就是所谓的粗差定位问题。它不仅仅是一个理论问题， 而更主要的是算法上的问题。 2.2选权迭代法的基本思路 选权迭代法的基本思想是：由于粗差未知，平差仍从惯常的最小二乘法开 始，但在每次平差后，根据其残差和有关其他参数，按所选择的权函数，计算每 一个观测值在下一步迭代平差中的权。 如果权函数选择得当，且粗差可定位，则 含粗差观测值的权将愈来愈小， 直到趋近于零。迭代终止时，相应的残差将直接 指出粗差的值，而平差的结果将不受粗差的影响。 这样便实现了粗差的自动定位 和改正。 该方法从下列最小条件出发： pivi2 min (1) 式中权函数 (v1) (v) （ ） pi f(vi ,...) v 1,2,... (2) 对于目前已提出的各种权函数，可以按其内容区分为： 残差vi的函数 pi f1(vi) (3) 标准化残差wi的函数 pif1(wi) (4) 方差估计?li的函数 pif1(?li) (5) 若按其形式，则可以区分为幂函数和指数函数。 已见到的一些权函数可列举如下： 1）Lq迭代法（最小范数迭代法） 该法由最小范数法vi q min演变而来，当q=1和q=0时，相应的权函数为 pi(v1) 1 c （6） vi(v) pi(v1) 1 （7） 2(v) c （vi ） 式中，c为一适当小的正数。 引入它是为了在改正数为零时不至引起迭代求解的困难。 由于此解法中的第一步仍为惯常的最小二乘平差， 所以它与严格的最小范数 解法有明显的区别。 2）丹麦法 丹麦法有多种形式的权函数，由Krarup等人提出的一种适合于摄影测量平差的权函数为： 第一次迭代 pi 1 0.05 vi 4.4 第二三次迭代 pi exp (8) 3.0 第二三次迭代 piexp vi Kubik教授提出的权函数则为：  0.05 1, 若 vp0 c f(v) v p0 m0 , exp cm0 若 vp0 c m0 其中，p0 权系数，m0 为观测中误差而常数c=3。 3）带权数据探测法 此法由加拿大EL-HaKing提出并试验过，其权函数为 (v1) 1 ,当wi 4.13 pi wi 2 ，当 wi 4.13 1 vi 式中，wi ?0 qvii 2.2粗差的估计 可以根据残差向量估计残差的大小，若设 l?i为第i 个观测值粗差的估值， 则有： （9） l?i vi li 而该估值的精度 l?i可按误差传播定律导出，其结果为 l?i l i （10） i ri 3算法步骤 (1) 列立误差方程，令各权因子初值均为 1，即令w1w2 wn1， WL，则P(0) P，p为观测权阵； (2) T ? T PL ? 和残差V的第一次估值： 解算法方程B PBX B 0，得出参数X （） T -1T ?1 PL X （B PB）B （1） （1） ? -L VBX (3) （1） （Vi） wi确定各观测值新的权因子，按 pipiwi构造新的 由V 按 Vi 等价权 p (1) ，再解算法方程 T ? T PL ? 和残差的第二次 B PBX B 0，得出参数X 估值； ?(2) (B T P (1) B) 1 B T（1） X PL V (2) ?(2) L BX 由V(2)构造新的等价权P(1)，再解算法方程，类似迭代计算，直至前后两次解的差值符合限差要求为止； 最后结果为 V (k) ?(k1) L BX 由于pi piwi，而wi (vi)，（vi） ，故随着 函数的选取不同，