(第三讲秩亏自由平差.pptVIP

上节广义最小二乘准则： 1、基本模型为： 2、平差准则： 秩亏自由网平差 一、问题的提出 自由网： 当控制网中没有必要的起算数据时，通常称为自由网。 附合网、独立网： 当控制网中除必要起算数据时外，还有多余的起算数据的网，称为附合网；等于必要起算数据，称独立网。 自由网平差方法分为： 经典自由网平差和秩亏自由网平差两种。 一些特殊用途的控制网，如变形观测网、沉降监测网等，一般为自由网。 1、经典自由网平差 选定x3的高程为已知，则可列出误差方程为： 法方程： 系数阵的行列式不为零，即R（N）=2，非奇异，方程有唯一解： 经典平差法的条件： 是在控制网中必需设定足够的坐标起算数据； 也可设定各点的高程近似值时，取x3的已知高程为近似值，但 。（即设一点的高程为已知）其函数模型为： 这就成为附有条件的间接平差了。 2、秩亏自由网平差 组法方程： 法方程系数阵： 可见，系数阵的行列式等于零,即是一个奇异阵，方程有无穷多组解。 产生秩亏的原因：就是平差网形中缺少的必要起算数据个数。 秩亏数d：就是秩亏自由网中的基准亏损数， d=R（B）-R（B） （ R‘（B）是B的列满秩数，R（B）是实际秩数。） 秩亏自由网平差： 如果网中不设起始数据或没有必要的起算数据，而且又设所有网点坐标为参数，这样的平差问题称为秩亏自由网平差。 思考： 在没有起算数据的网中，秩亏数和什么个数 相等？ 水准网、测角网、测边网、边角网以及GPS网的秩亏数各是多少？ 二、秩亏自由网平差原理 秩亏自由网平差的函数模型为 相应的误差方程为 随机模型为 法方程为 问题的提出:在秩亏自由网平差中,如果像经典平差平差那样,只要求遵循最小二乘原则求未知参数的解,将不可能取得唯一确定的估计量; 解决方法:为了得唯一确定的估计量,需要在遵循最小二乘原则基础上附加另外条件; 附加条件的前提:该条件的确定应保证所求得的未知数的估计量是最优的. 这样的最优解是唯一存在的,它就是法方程的最小范数解! 1、秩亏法方程的最小范数解（直接解法） 设满足法方程的一个解为X,取其平方和的开方为 称为向量X的范数,几何意义是向量的长度。 最小范数满足条件,称为最小范数条件,其表达式为 法方程若有一解X满足其范数最小,这个解就称为最小范数解。 求最小范数的法方程解过程： 即求下列数学解： 得： 广义逆矩阵的概念 1）广义逆A- 1、定义：设A的秩R（A）=r≤min(n,m)，满足下列矩阵方程的A-定义为A的广义逆 2、广义逆A-的计算 A是非奇异方阵，凯利逆A-1就是A的广义逆。 A是列满秩时 A是行满秩时 A是降秩矩阵时：秩分解法、降阶法。 降阶法： 在秩亏的方阵A中，删去d个某一行和相应的某一列降阶求逆，删去位置均以“0”代之，即得奇异方阵的广义逆A-。 可见A-不唯一。 例如： 可以验证： 2）广义逆A+（Moore-Penrose广义逆、伪逆） 1、定义：满足下列四个条件，即 2、 A+的计算 当A为对称方阵时： 值得说明的是： 1）因广义逆不唯一，但可以证明，用不同的广义逆（NN）-代入上式后，求得的X向量却是相同的，故X有唯一解！ 2）以上解法又称为“直接解法”。 精度估计 参数估值的协因数阵： 单位权方差估值仍为： 2、附加条件法（是一种实用算法） 自由网误差方程为 为消除秩亏，附加条件 按最小二乘原则，作函数 得法方程 解法方程，得X解 或者，整理得： 3. 伪观测值法 以上介绍可知： 直接算法需要求广义逆； 附加条件法或伪观测法需要确定两项：Px和S。 按平差基准不同可将自由网平差分为三类: 1）以全部网点重心为基准(简称重心基准) 的秩亏自由网平差；(PX=E） 2）以网中部份相对稳定点重心为基准(简称拟稳基准) 的拟稳自由网平差(简称拟稳平差，PX1=0,PX2=E,) ； 3）即网中存在d 个起始数据, 这就是固定基准下的经典自由网平差。 秩亏问题解决：经典平差（附加固定的基准条件）和伪逆平差（直接利用广义逆求解 ）； 优缺点：解法简捷 ，但没考虑到解法物理意义，不能反映真实情况。 提出：拟稳平差理论。 “拟稳平差”的基本思想：考虑到监测网中的点，处于不同的地质构造和地球物理环境，随着时间的延伸，都可能发生变动，但是总存在相对变化小的，即相对稳定的点。 4、自由网拟稳平差 我国周江文研究员针对变形监测网提出了一种所谓拟稳平差的自由网平差方法。 1、自由网拟稳平差原理 设不稳定未知数X1，稳定未知数为X