非线性最小二乘平差

《非线性最小二乘平差》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、o6-2非线性模型平...o6-3非线性模型平...o6-4非线性模型...o6-1问题的提出非线性最小二乘平差6-1问题的提出经典平差是基于线性模型的平差方法。然而在现实世界中，严格的线性模型并不多见。测量上大量的数学模型也是非线性模型。传统的线性模型平差中的很多理论在非线性模型平差中就不一定适用；线性模型平差中的很多结论在非线性模型平差中就不一定成立；线性模型平差中的很多优良统计性质在非线性模型平差中就不一定存在。例如，在线性模型平差中，当随机误差服从正态分布时，未知参数X的最小二乘估计具有一致无偏性和方差最小性。但在非线性模型平差中，

2、即使随机误差严格服从正态分布，未知参数X的非线性最小二乘估计也是有偏的。其方差一般都不能达到最小值。 对于测量中大量的非线性模型，在经典平差中总是进行线性近似（经典的测量平差中称之为线性化），即将其展开为台劳级数，并取至一次项，略去二次以上各项。如此线性近似，必然会引起模型误差。过去由于测量精度不高，线性近似所引起的模型误差往往小于观测误差，故可忽略不计。随着科学技术的不断发展，现在的观测精度已大大提高，致使因线性近似所产生的模型误差与观测误差相当，有些甚至还会大于观测误差。例如，GPS载波相位观测值的精度很高，往往小于因线性近似所产生的模

3、型误差。因此，用近似的理论、模型、方法去处理具有很高精度的观测结果，从而导致精度的损失，这显然是不合理的。现代科学技术要求估计结果的精度尽可能高。这样，传统线性近似的方法就不一定能满足当今科学技术的要求。另外，有些非线性模型对参数的近似值十分敏感，若近似值精度较差，则线性化会产生较大的模型误差。由于线性近似后，没有顾及因线性近似所引起的模型误差，而用线性模型的精度评定理论去评定估计结果的精度，从而得到一些虚假的优良统计性质，人为地拔高了估计结果的精度。 鉴于上述各种原因，对非线性模型平差进行深入的研究是很有必要的。非线性模型的平差和精度估计

4、以及相应的误差理论研究也是当前国内外测绘界研究的前沿课题之一。电子教材 > 第六章非线性模型平差 > 6-2非线性模型平差原理一、非线性误差方程测量中大量的观测方程是非线性方程。比如导线测量中，以待定点坐标为未知参数的角度观测方程和边长观测方程分别为：（6-2-1）式中：为待定点坐标的真值，分别为角度观测值和边长观测值的真误差。角度观测值和边长观测值的观测方程（6-2-1）式是待定点坐标真值（）的非线性函数。又如在GPS伪距测量中，第j颗卫星至测站k的几何距离的观测方程为：也是测站点k的待定坐标真值（）的非线性函数。一般地，用L表示的观测向

5、量，用表示的未知参数向量的真值，用△表示的真误差向量，则非线性观测方程可写为：（6-2-2）式中:，是由n个的非线性函数组成的的向量;。（6-2-2）式就是我们所要讨论的一般的非线性模型。在一般的非线性模型（6-2-2）式中，用未知参数向量和真误差向量的估计值代替其真值，得非线性误差方程如下：（6-2-3）式中：V为观测值的改正数向量（残差向量）；为参数向量的估值。  二、非线性模型平差由非线性误差方程（6-2-3）式知，非线性误差方程（6-2-3）式中仅有n个方程，而有n+t个未知数（n个观测值的改正数和t个参数）。因此非线性误差方程（6

6、-2-3）式是非线性不定方程组，有无穷组解。在这无穷组解中，必然有一组解能使（6-2-4）我们将满足（6-2-4）式的一组解作为最优解，并称（6-2-4）式所确定的为的一个非线性最小二乘[23]估计。本书中将求解非线性最小二乘估计的过程称为非线性模型平差。可见，非线性模型平差与线性模型平差的是完全一致的。（6-2-4）式的几何意义就是观测空间至解空间的距离最短，或者说是解轨迹π上离观测值L最近的点（见图6-1）。L到π的距离就是‖V‖。图6-1在非线性模型（6-2-3）式中，若存在一阶连续偏导数，且的非线性最小二乘估计量存在，则残差向量V在

7、处垂直于切空间T（见图6-1）[24]。一、非线性最小二乘估计的近似解当非线性模型（6-2-3）式的非线性强度[24]较弱时，可以将非线性模型在处线性近似，并用线性模型的求解理论和方法来近似地求解非线性模型（6-2-3）式。这也就是我们大家所熟悉的传统方法——线性化方法，即将非线性模型（6-2-3）式在处用台劳级数展开，取至一次项，得：（6-3-1）令（6-3-2）（6-3-3）则（6-3-1）式可写为：（6-3-4）（6-3-4）式就是我们熟悉的间接平差的误差方程。由间接平差知，根据最小二乘原理可解得：（6-3-5）于是参数X的非线性平差

8、结果为：（6-3-6）例6-1（本例取自参考文献[24]）已知非线性模型为。其中参数和的真值为。的5个真值（用参数的真值X算得）和相应的5个同精度独立观测值列于表6-1。表6-1