6 第五章 条件平差.ppt

线号 高差（m） 路线长度（km） 点号 高程（m） 1 1.100 4 A 5.000 2 2.398 2 B 3.953 3 0.200 4 C 7.650 4 1.000 2 ? ? 5 3.404 2 ? ? 6 3.452 4 ? ? A ? ? ? o o o B C 1 2 3 4 5 6 P1 P2 P3 如图所示的水准网，A、B、C已知水准点，P1、P3、P3为待定点，已知水准点的高程、各水准路线的长度及观测高差列入下表： 试用条件平差法求P1、P3、P3点高程的平差值 。 作 业 建筑与城乡规划学院 * 在数学最优化问题中，拉格朗日乘数法（以数学家约瑟夫·路易斯·拉格朗日命名）是一种寻找变量受一个或多个条件所限制的多元函数的极值的方法。这种方法将一个有n 个变量与k 个约束条件的最优化问题转换为一个有n + k个变量的方程组的极值问题，其变量不受任何约束。 建筑与城乡规划学院 * §2 条件方程 右图边角网： §2 条件方程 §2 条件方程 下图独立三角网，n=? t=? r=? §2 条件方程 各类条件方程数 §2 条件方程 例 例： 1）列条件式方程 2）平差后AD的平差函数试为方程： 3）组、解法方程： §3 精度评定 测量平差 求最或然解 评定精度 观测量的平差值精度 平差值函数的精度 观测值的方差可表达为： 观测量的平差值的方差： 推广到求平差值函数的方差： 平差值函数的方差:一是求单位权方差估值；二是求该量的协因数 §3 精度评定 一、验后单位权方差的估值 即： 不管采用何种平差方法，单位权方差的公式都是不变的。 VTPV 的计算 VTPV =（QATK)TPQATK VTPV =VTPQATK =KTAQATK =KTNaaK =VTATK=(AV)TK -WTK §3 精度评定 二、协因数阵的计算 1)条件方程常数项 W（-f）的协因数 则由协因数传播率得： 列函数式： L的协因数为Q §3 精度评定 2)观测值改正数V（残差）的协因数 则由协因数传播率得： 列函数式： 二、协因数阵的计算（续） 方法一 方法二 §3 精度评定 3)联系数K的协因数 则由协因数传播率得： 列函数式： 二、协因数阵的计算（续） 则由协因数传播率得： 列函数式： §3 精度评定 二、协因数阵的计算（续） 4)观测量平差值 的协因数 ^ L 列函数式： §3 精度评定 则由协因数传播率得： §3 精度评定 即观测量平差值： 其协因数为： 在条件平差的精度评定中，这是一个非常重要的公式。 §3 精度评定 二、协因数阵的计算（续） 5)互协因数阵 函数式： 按互协因数传播率得： 可见：观测值的平差值和改正数是不相关的 §3 精度评定 二、协因数阵的计算（续） 5) 协因数阵（其他的大家自己推算） QAT Q AQ -NaaAQ -1 Q-QATNaaAQ -1 -QVV Naa -I -QAT 0 -QATNaa -1 Q-QVV 0 0 0 Q-QATNaaAQ -1 -Qvv -AQ NaaAQ -1 QATNaaAQ -1 0 -I Naa -1 QATNaa -1 0 条件平差基本向量的协因数阵 §3 精度评定 （Naa=AQAT） §3 精度评定 三、平差值函数的中误差 平差值函数： ——指根据观测值的平差值所算出来的某些量 如:平面控制网中的待定点坐标、边长、方位角等以及高程制网中的待定点平差的高程等 实际平差问题中： 观测值L 观测值的平差值 平差值的函数 平差值函数的中误差 观测值中误差σ 观测量平差值的中误差 §3 精度评定 三、平差值函数的中误差 平差值函数： ——指根据观测值的平差值所算出来的某些量 如:平面控制网中的待定点坐标、边长、方位角等以及高程制网中的待定点平差的高程等 求平差后CD 边方位角 边长 边长 §3 精度评定 设有平差值函数为 取全微分 若令 则线性化后的函数形式为 令 则 由 则 而 故 §3 精度评定 三、平差值函数的中误差 实际平差问题中： 观测值L 观测值的平差值 平差值的函数 通过平差 计算 协方差传播定律 列函数式 函数式线性化 2 §3 精度评定 解： 1）列条件方程： 其中： 2）列出用平差值的函数形式表示的CD边长： 求得： 对CD边长函数两边微分得： 令 两边除以SCD 则 ，且 3）由协方差传播定律求： 代入f 和 得： 4）计算相对中误差 由 §4 条件平差公式汇编与实例 一、条件平差公式汇编 函数模型 随机模型 条件方程: AV+W=0 法方程： 改正数方程： 观测量平差值： §4 条件平差公式汇编与实例 一