非线性最小二乘平差.docx

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

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

3、这显然是不合理的。现代科学技术要求估计结果的精度尽可能高。这样，传统线性近似的方法就不一定能满足当今科学技术的要求。另外，有些非线性模型对参数的近似值十分敏感，若近似值精度较差，则线性化会产生较大的模型误差。由于线性近似后，没有顾及因线性近似所引起的模型误差，而用线性模型的精度评定理论去评定估计结果的精度，从而得到一些虚假的优良统计性质，人为地拔高了估计结果的精度。鉴于上述各种原因，对非线性模型平差进行深入的研究是很有必要的。非线性模型的平差和精度估计以及相应的误差理论研究也是当前国内外测绘界研究的前沿课题之一。 电子教材第六章 非线性模型平差 6-2 非线性模型平差原理 一、非线性误差方程测

4、量中大量的观测方程是非线性方程。比如导线测量中，以待定点坐标为未知参数的角度观测方程和边长观测方程分别为： （6-2-1）式中：为待定点坐标的真值，分别为角度观测值和边长观测值的真误差。角度观测值和边长观测值的观测方程（6-2-1）式是待定点坐标真值（）的非线性函数。又如在GPS伪距测量中，第j颗卫星至测站k的几何距离的观测方程为： 也是测站点k的待定坐标真值（）的非线性函数。一般地，用L表示的观测向量，用表示的未知参数向量的真值，用表示的真误差向量，则非线性观测方程可写为： （6-2-2）式中: ，是由n个的非线性函数组成的的向量;。（6-2-2）式就是我们所要讨论的一般的非线性模型。在一般

5、的非线性模型（6-2-2）式中，用未知参数向量和真误差向量的估计值代替其真值，得非线性误差方程如下： （6-2-3）式中：V为观测值的改正数向量（残差向量）；为参数向量的估值。 二、非线性模型平差由非线性误差方程（6-2-3）式知，非线性误差方程（6-2-3）式中仅有n个方程，而有n + t 个未知数（n个观测值的改正数和t 个参数）。因此非线性误差方程（6-2-3）式是非线性不定方程组，有无穷组解。在这无穷组解中，必然有一组解能使 （6-2-4）我们将满足（6-2-4）式的一组解作为最优解，并称（6-2-4）式所确定的为的一个非线性最小二乘23估计。本书中将求解非线性最小二乘估计的过程称为非

6、线性模型平差。 可见，非线性模型平差与线性模型平差的是完全一致的。（6-2-4）式的几何意义就是观测空间至解空间的距离最短，或者说是解轨迹上离观测值L最近的点（见图6-1）。L到的距离就是V。 图6-1在非线性模型（6-2-3）式中，若存在一阶连续偏导数，且的非线性最小二乘估计量存在，则残差向量V在处垂直于切空间T（见图6-1）24。一、非线性最小二乘估计的近似解当非线性模型（6-2-3）式的非线性强度24较弱时，可以将非线性模型在处线性近似，并用线性模型的求解理论和方法来近似地求解非线性模型（6-2-3）式。这也就是我们大家所熟悉的传统方法线性化方法，即将非线性模型（6-2-3）式在处用台劳

7、级数展开，取至一次项，得： （6-3-1）令 （6-3-2） （6-3-3）则（6-3-1）式可写为： （6-3-4）（6-3-4）式就是我们熟悉的间接平差的误差方程。由间接平差知，根据最小二乘原理可解得： （6-3-5）于是参数X的非线性平差结果为： （6-3-6）例6-1（本例取自参考文献24）已知非线性模型为。其中参数和的真值为。的5个真值（用参数的真值X算得）和相应的5个同精度独立观测值列于表6-1。表6-1 的真值和相应的观测值i12345真值4.2028343.2589242.5270061.9594691.519394观测值4.203.252.521.951.51观测方程为： 取

8、参数X的近似值为。将观测方程在处线性近似，得误差方程：由（6-3-5）式得： 于是，由（6-3-6）式得参数X的平差值为： 参数估值的真误差为： 其范数为：二、非线性最小二乘平差的迭代解当非线性模型的非线性强度很强时，线性近似可能产生大于观测误差的模型误差，所以对于非线性模型，一般采用迭代的方法求解。求解非线性误差方程（6-2-3）式的最小二乘平差值，就是求参数X的估值，使 （6-3-7）由于是一常量，所以（6-3-7）式等价于目标函数为 （6-3-8）的非线性无约束最优化问题。因为是的非线性函数，所以对（6-3-8）式求一阶偏导数，并令其为零，得不到的显表达式。故求不出的解析解。因此，我们只

9、能设法寻找某一近似解，使 （6-3-9）成立。寻找使（6-3-9）式成立的近似解，一般只有采用迭代的方法。为此，下面介绍几种常用的迭代方法。1牛顿法设的极小值的一个近似值为，在附近将展为台劳级数，取至二次项得： （6-3-10）式中： （6-3-11） （6-3-12）称为处的Hessian矩阵。 （6-3-13）是在处的梯度方向。由于是的一个已知的近似值，故（6-3-10）式只是的函数，为了求得使（6-3-10）式成立的，将（6-3-10）式对求偏导，并令其为零，得： 移项后两边转置，顾及（6-3-12）式，得 （6-3-14）当Gk非奇异时，由（6-3-14）式可解得使（6-3-10）式成

10、立的： （6-3-15）当充分小时，能使（6-3-10）式成立。但由于未知，故不能充分小，需不断迭代，直至充分小，其迭代公式为： （6-3-16）（6-3-16）式就是牛顿迭代的基本公式，迭代终止条件： （6-3-17）或 =0 （6-3-18）由于是一个绝对值较大的数，而的各元素的绝对值都很小，因此，由于计算机有效数字的限制，以（6-3-17）式作为迭代收敛条件比（6-3-18）式作为迭代收敛条件收敛要快一些。牛顿法的迭代步骤为：（1）选取初值，并令k=0。（2）按（6-3-11）式计算梯度方向，若=0则转至（7）。（3）计算Hessian矩阵。（4）解线性方程组（6-3-14）式，得。（5

11、）按（6-3-16）式计算新的近似值。（6）计算目标函数值，若则转至（2）继续迭代。（7）终止迭代，输出和，结束。例6-2 在例6-1中，仍设，用牛顿法求例6-1中非线性模型的非线性最小二乘平差值。解：由例6-1知P=I，故目标函数为： 将代入计算，G0后，按以上迭代程序迭代，结果列于表6-2。表6-2 牛顿法迭代计算k123456-1.2050249080.3991833382.028*0.0001691624122-3.949210-9 -2.401210-9 -17.15305037.0372427130.49484074240.0029380122642.403910-7 -1.556

12、910-7 5.3330132655.417198095.4227080035.4427445655.4227445935.422744582-0.2539145225-0.2542573375-025*-0.2556720853-0.2556720877-0.2556720866-40.21054702-40.58524686-40.63522342-40.63549278-40.63549281-40.63549281迭代6次后，有= -40.63549281，所以停止迭代，得X的非线性最小二乘解为 则， 由例6-1知，本迭代解与其真值的距离比线性近似解与其真值的距离要小一个数量级。当初值

13、取时，迭代发散，这说明牛顿法对初值很敏感。2信赖域法牛顿法具有很快的收敛速度，但它总是局部收敛的。因为牛顿法的基本思想是用二次函数 去逼近。只有当充分小时，才能很好地逼近。既然只有当充分小时，才能逼近，那么可以对dX加以限制，然后在限制条件下来寻求的极小值。这个思想相当于求解下列约束最优化问题：目标函数： （6-3-19）约束条件：式中：为一正数，它随迭代而变化。约束条件 限制了，使的长度不大于，这样总在一个给定的小区域中活动。这个区域是可信赖的，所以称该方法为信赖域法。常数取决于对的逼近程度。这个逼近程度可用下式来描述 （6-3-20）越接近于1，对的逼近程度越好，于是 （6-3-21）这样，可总结出信赖域法的迭代程序：（1）选取初值，。（2）按（6-3-11）式和（6-3-12）式计算梯度方向和矩阵，若=0则转至（7）。（3）按（6-3-15）式计算，并检查是否满足约束条件。若不满足，则采取适当方法对予以压缩。然后在区域内求使=min的。（4）计算的新的近似值。（5）按（6-3-20）式计算，并按（6-3-21）式确定。（6）检查是成立。若不成立，则转（2）继续迭代。（7）终止迭代，输出和，结束。例6-3，设，=0.08，用信赖域法求解例6-1中非线性模型的非线性最小二乘平差值。、g和Gk的表达式同例6-2，用信赖域法迭代计算的结果列于表6