测量中病态问题的一种迭代算法 an iteration method of ill-conditioned problems in surveying.pdfVIP

测量中病态问题的一种迭代算法 张 丽1’2，张书毕1’2，张秋昭1’2，牛作鹏1’2 (1．中国矿业大学环境与测绘学院．江苏徐州221116； 2．中国矿业大学江苏省资源环境信息工程重点实验室．江苏徐州221116) 摘 要：系统病态性问题在测量数据处理领域是常见的，具有许多种解算方法。迭代法 是其中的一种。通过分析病态矩阵的性质，基于谱修正迭代法，构造了一种迭代算法，实例表 明：该算法是有效的，提高了迭代速度，具有实际意义。 关键词：病态；最小二乘估计；迭代法 中图分类号：P207文献标志码：A 文章编号：1008—9268(2010)06—0027—03 N宕=W (1) 0 引 言 式中：N=B7PB；W=B丁PL；P—D21． 在参数的最小二乘估计中，经常会遇到病态问 将式(1)两边同时加上又，得 题。这是由于系统结构不好，观测信息量少导致估 (N+J)文一w+支 (2) 计不稳定，从而使结果严重偏离真值。例如，在 式中，J为t×1的单位阵。 GPS快速定位中，由于要求测量时间短，GPS卫星 因为N是对称正定的满秩方阵，所以可以构 与测站之间的几何关系变化不明显，导致观测结构 造如下迭代格式 很差，要同时解算坐标分量和相位模糊度参数，就 支(K十1)=(N+J)一1(w+宕‘K’)， 会遇到严重的病态问题【1]。不采取有效的措施，就 (K=1，2，3，…) (3) 得不到精确的定位结果。 当法矩阵N为病态矩阵，对式(1)直接求逆时 当前，解决病态问题的主要方法是以岭估计为 所得结果与真值相距甚远，采用此谱修正迭代法能 主的有偏估计法和基于奇异值分解的直接解算方 得到较为精确的解。 法。文献E2-1提出了谱修正迭代法。谱修正迭代法 2 迭代算法 是一种直接解算方法，在一定程度上改善了最小二 乘估计。但是算例表明，该方法的收敛速度并不十 2．1 谱修正迭代法的分析 分理想。分析产生问题的原因，在此基础上构造了 为说明问题，首先将迭代式(3)式改写如下 一种迭代算法，大大提高了收敛速度。 (4) 为简化，令(N+J)_1=G，所以式(4)变为 1谱修正迭代法 文眯+I)=G宕Ⅲ’+Gw 设观测方程为 由数值分析可知，此迭代格式收敛的充要条件 L—BX+厶 式中：B为以×t的设计矩阵；L为，z×l的观测向敛越快，p(G)越接近1，收敛越慢，因而所需迭代次 量；△为姐×1的观测误差向量，且有E(△)=0，D 数越多。引入谱半径的定义：设A∈R积“，At，Az， (厶)=D。；x为￡×l的未知参数向量。 令x的估值为文，由间接平差得误差方程为 半径‘3|。根据定义，谱修正迭代法中，迭代矩阵G y=．救一L 的谱半径为 根据最／b乘原理有法方程 1 』D(G)5—1