第三章 秩亏网平差(研究生).pptVIP

在自由监测网中，假定有一部分点对于另一部分是相对稳定的。（就有稳定未知数和不稳定未知数两类） 第三章 监测网平差及参考点稳定性检验 本章主要内容 1、预备知识；2、监测网经典平差；3、秩亏自由网平差；4、监测网拟稳平差；5、参考点稳定性检验方法。 经典自由网平差算例，如图所示：选定X3的高程为已知，X1、X2的高程为未知，则可列出误差方程为： 法方程 系数阵的行列式不为零，即R（N）=2，非奇异，方程有唯一解： 经典平差法的条件： 是在控制网中必需设定足够的坐标起算数据； 3.1预备知识 3.1预备知识 如果不假设起始高程，设网中全部待定点为参数，则误差方程为： 法方程系数阵： 可见，系数阵的行列式等于零,即是一个奇异阵，方程有无穷多组解。 产生秩亏的原因：一是平差网形中缺少的必要起算数据个数，而不假设它，二是所设的未知数必须独立。 秩亏数d：就是秩亏自由网中的基准亏损数， d=R（B）-R（B）（ R‘（B）是B的列满秩数，R（B）是实际秩数。） 如果网中不设起始数据或没有必要的起算数据，而且又设所有网点坐标为参数，这样的平差问题称为秩亏自由网平差。 怎么解算秩亏自由网平差？我们首先了解广义逆矩阵。 3.1预备知识 一、广义逆矩阵 1）广义逆A- 1、定义：设A的秩R（A）=r≤min(n,m)，满足下列矩阵方程的A-定义为A的广义逆 2、广义逆A-的计算 A-广义逆的计算有多种方法，下面介绍测量中常用的简便法 当A的秩R（A）=r≤min可将(n,m)，可将A分块成 3.1预备知识 其中 ，则 3.1预备知识 算例1：设矩阵 取 3.1预备知识 2）最小范数逆 1、定义：设A的秩R（A）=r≤min(n,m)，满足下列矩阵方程的 定义为A的最小范数逆。 2、最小范数逆 的计算：在自由网平差中，最小范数逆 常取 3.1预备知识 算例2：设矩阵 ，计算最小范数逆 3）广义逆A+（Moore-Penrose广义逆、伪逆） 1、定义：满足下列四个条件，即 2、 A+的计算 当A为对称方阵时： 3.1预备知识 3.1预备知识 算例2：设矩阵 计算伪逆A+ 已知 则： 3.1预备知识 二、相容线性方程组的广义逆解 设有相容线性方程组 其解存在，设为Y，则有 按广义逆A- 有 此式称为方程组（1）的相容条件，由此可知 上式为方程组（1）的一个特解。 齐次方程组 的一般解为： 因此非齐次方程组的一般解为： 3.3秩亏自由网平差 秩亏自由网平差的函数模型为 相应的误差方程为 随机模型为 法方程为 法方程是相容方程组，其解有无穷多组 三、秩亏自由网参数的最小二乘解 3.3秩亏自由网平差 问题的提出:在秩亏自由网平差中,如果像经典平差平差那样,只要求遵循最小二乘原则求未知参数的解,将不可能取得唯一确定的估计量; 解决方法:为了得唯一确定的估计量,需要在遵循最小二乘原则基础上附加另外条件; 附加条件的前提:该条件的确定应保证所求得的未知数的估计量是最优的. 这样的最优解是唯一存在的,它就是法方程的最小范数解! 3.3秩亏自由网平差 （五）秩亏法方程的最小范数解 设满足法方程的一个解为X,取其平方和的开方为 称为向量X的范数,几何意义是向量的长度。 最小范数满足条件,称为最小范数条件,其表达式为 法方程若有一解X满足其范数最小,这个解就称为最小范数解。 3.3秩亏自由网平差 求最小范数的法方程解过程： 即求下列数学解： 得： 则 3.3秩亏自由网平差 值得说明的是：在秩亏自由网平差中，最小范数逆可取 ，与伪逆的条件相比，最小范数逆条件仅是 其中的两个条件，由于只有满足4各条件的伪逆才是唯一存在，因此最小范数逆不是唯一的，但是最小范数解却是唯一存在的。了解这一点对于研究自由网平差理论十分重要，下面作出证明。 已知最小范数逆满足两个条件 与满足下列一个条件是一致的 二、秩亏自由网平差广义逆解法 3.3秩亏自由网平差 设有两个最小范数 和 ，相应的最小范数解为： 按最小范数 的条件，则 即： 两边右乘 得 或 上式成立，必须满足 右乘任意向量Y，得： 因为 ，故有 即 可见最小范数逆不是唯一的，但是最小范数解却是唯一的 数学模型: 按照附有条件的间接平