05线性参数的最小二乘法处理.pptVIP

非线性误差方程的近似值应当尽可能准确，只有当参数近似值与平差值很接近时才成立；若无法得到可靠的近似值时，可迭代计算。 五、非线性最小二乘法 （1）法方程不可逆 （2）有偏估计 （3）观测值稳定性不同 （4）权比例不固定 （5）观测值含有噪声 （6）观测值不独立 （7）观测值含有粗差 （8）非最小二乘准则 ……… 问题思考： 加权秩亏自由网平差 有偏估计 拟稳平差 验后方差分量估计 Kalman滤波 相关平差 粗差探测 贝叶斯估计 ……… 现代 平差 总结 * 最小二乘法（Least Square Method）的历史 最小二乘法的基本原理 线性参数的最小二乘法 等精度、不等精度 非线性参数的最小二乘法 等精度、不等精度 精度估计 直接测量量的精度估计 最小二乘估计量的精度估计 组合测量的最小二乘原理 ，非线性最小二乘 * 6.626 * * * 权阵P的含义 观测值不独立时 二、最小二乘法解算 2.2 不等精度最小二乘法解算 这是一个典型的多元函数求极值的问题，其解法主要取决于目标函数自变量的选取。 如果自变量之间存在约束条件，则需采用条件极值法； 如果自变量之间函数独立，则可采用自由极值法； 如何求函数V TPV 的最小值？ 最小二乘本质上就是求函数的极值问题。 二、最小二乘法解算 2.2 不等精度最小二乘法解算 极值方程 由误差方程 极值方程转化为： 误差方程代入上式： 得 二、最小二乘法解算 2.2 不等精度最小二乘法解算 法方程 N——系数阵 U——自由项 参数的解向量： 二、最小二乘法解算 2.2 不等精度最小二乘法解算 说明: 1. 法方程系数阵N 的维数t×t ; 2. 法方程的秩r(N)=r(ATPA)=r(A)=t 时，即N 满秩，法方程才有唯一解; 3. 当r(A)=t 时，称矩阵A列满秩，为保证 A列满秩，未知参数选取是关键; 二、最小二乘法解算 2.2 不等精度最小二乘法解算 命题即得证。 思路： 设另有解 ，改正数向量V1，只要 如何证明： 是VTPV 唯一的最小值点。 二、最小二乘法解算 2.2 不等精度最小二乘法解算 因为V1和V均是误差方程的解，所以 两式相减，得 所以， 是VTPV 唯一的最小值点。 二、最小二乘法解算 2.2 不等精度最小二乘法解算 参数平差的解向量为： 两边取数学期望，得 如何证明： 就是未知参数的无偏估计。 二、最小二乘法解算 2.2 不等精度最小二乘法解算 * 二、最小二乘法解算 算例3：已知某测量过程的误差方程和标准差，试求解未知参数的最小二乘解。 ?=0.24 * 二、最小二乘法解算 算例3 如图所示沿三条水准路线测定结点P的高程，观测结果如表所示。试求P点的高程的最或然值。（注：箭头指向高程增加的方向） 已知点 已知点高程 H/m 水准路线 观测高差 h/m 路线距离 A 12.431 1 7.421 2 B 10.252 2 5.237 5 Ｃ 1.567 3 3.453 10 二、最小二乘法解算 算例4 解法1： ① 定权，选足数且独立的未知参数 取10km观测高差的权为单位权，即S0=10，由水准路线定权公式 ②′列出观测方程： 设P点的高程的最或然值为 得各段水准路线高差的权： 二、最小二乘法解算 算例4 ③ 组成法方程 ④ 解法方程，求得未知向量估值： ⑤ 回代误差方程，求出改正数和观测值平差值； ② 列出误差方程式： 二、最小二乘法解算 算例4 解法2： ① 定权，设P点高程的平差值为 得各段水准路线高差的权： ② 沿不同路径计算P点高程 权： 5 2 1 ③ 利用权中数计算P点高程 二、最小二乘法解算 算例4 最小二乘法数据处理的目的 （1）给出待求量的最可信赖估计值 （2）确定可信赖程度，即求解估计量的精度 三、最小二乘法精度估算 精度估计思路 （1）为求出最小二乘平差值的精度，需要估算测量数据精度； （2）测量数据的精度由标准差?表示， ?无真值，只能依据测量数据计算其估计值； * 可以证明 是自由度为 n-t 的 变量。根据 变量的性质： 将不是 的无偏估计量 三、最小二乘法精度估算 3.1 测量数据精度估算 将等精度测量中残余误差平方和变为加权的残余误差平方和，就可以借用等精度测量数据的精度估计方法，对不等精度测量数据的精度进行估计。 * 三、最小二乘法精度估算 3.1 测量数据精度估算 单位权标准差写成矩阵形式即为： * 三、最小二乘法精度估算 3.1 测量数据精度估算 * 3.2 平差量精度估算 三、最小二乘法精度估算 X的协方差阵