《非线性最小二乘平差》由会员分享，可在线阅读，更多相关《非线性最小二乘平差（26页珍藏版）》请在金锄头文库上搜索。
1、非线性最小二乘平差6-1问题的提出经典平差是基于线性模型的平差方法。然而在现实世界中，严格的线性模型并不多见。测量上大量 的数学模型也是非线性模型。传统的线性模型平差中的很多理论在非线性模型平差中就不一定适用；线 性模型平差中的很多结论在非线性模型平差中就不一定成立;线性模型平差中的很多优良统计性质在非 线性模型平差中就不一定存在。例如，在线性模型平差中，当随机误差服从正态分布时，未知参数X 的最小二乘估计占必具有一致无偏性和方差最小性。但在非线性模型平差中，即使随机误差严格服从 正态分布，未知参数X的非线性最小二乘估计乂芯也是有偏的。其方差一般都不能达到最小值。对于测量中大量的非线性模型，在
2、经典平差中总是进行线性近似（经典的测量平差中称之为线性 化），即将其展开为台劳级数，并取至一次项，略去二次以上各项。如此线性近似，必然会引起模型误 差。过去由于测量精度不高，线性近似所引起的模型误差往往小于观测误差，故可忽略不计。随着科学 技术的不断发展，现在的观测精度已大大提高，致使因线性近似所产生的模型误差与观测误差相当，有 些甚至还会大于观测误差。例如，GPS载波相位观测值的精度很高，往往小于因线性近似所产生的模型 误差。因此，用近似的理论、模型、方法去处理具有很高精度的观测结果，从而导致精度的损失，这显 然是不合理的。现代科学技术要求估计结果的精度尽可能高。这样，传统线性近似的方法就不
3、一定能满 足当今科学技术的要求。另外，有些非线性模型对参数的近似值十分敏感，若近似值精度较差，则线性 化会产生较大的模型误差。由于线性近似后，没有顾及因线性近似所引起的模型误差，而用线性模型的 精度评定理论去评定估计结果的精度，从而得到一些虚假的优良统计性质，人为地拔高了估计结果的精 度。鉴于上述各种原因，对非线性模型平差进行深入的研究是很有必要的。非线性模型的平差和精度 估计以及相应的误差理论研究也是当前国内外测绘界研究的前沿课题之一。一、非线性误差方程测量中大量的观测方程是非线性方程。比如导线测量中，以待定点坐标为未知参数的角度观测方程和边长观测方程分别为：思=如源一 arctg + A用
4、心-七 皿一七 ？禺二十A# 一（6-2-1）T- 】An WA.-BM式中：矿 乒为待定点坐标的真值，片 斗分别为角度观测值X和边长观测值的真误差。角度观测值和边长观测值的观测方程（6-2-1）式是待定点坐标真值（两凹）的非线性函数。又如在GPS伪距测量中，第j颗卫星至测站k的几何距离的观测方程为：底二/& _ C 十厂疗十（气_事沪十成十乱也是测站点k的待定坐标真值（% .气）的非线性函数。一般地，用L表示活xl的观测向量，用X表示京1的未知参数向量的真值，用表示於幻的真误差向量，则非线性观测方程可写为:（6-2-2）.T式中：/二皿为0）-），是由n个顼的非线性函数组成的E的向量;。（6
5、-2-2）式就是我们所要讨论的一般的非线性模型。在一般的非线性模型（6-2-2）式中，用未知参数向量和真误差向量的估计值代替其真值，得非 线性误差方程如下:（6-2-3）式中：V为观测值的改正数向量（残差向量）；乂为参数向量的估值。由非线性误差方程（6-2-3）式知，非线性误差方程（6-2-3）式中仅有n个方程，而有n + t个 未知数（n个观测值的改正数和t个参数）。因此非线性误差方程（6-2-3）式是非线性不定方程组，有 无穷组解。在这无穷组解中，必然有一组解能使尸 PF= S 刈-iyp(ftx)-l) = min(6-2-4)6 2 4)i.我们将满足（6-2-4）式的一组解作为最优解
6、，并称（6-2-4）式所确定的预为席的一个非线性最小二乘23 估计。本书中将求解非线性最小二乘估计的过程称为非线性模型平差。可见，非线性模型平差与线性模型平差的是完全一致的。（6-2-4）式的几何意义就是观测空间至解空间的距离最短，或者说是解轨迹n上离观测值L最近的点（见图6-1）。1到丸的距离就是|VII。图6-1在非线性模型（6-2-3）式中，若了（力存在一阶连续偏导数，且顶的非线性最小二乘估计量京存在，则残差向量V在侦处垂直于切空间T （见图6-1） 24。一、非线性最小二乘估计的近似解当非线性模型（6-2-3）式的非线性强度2较弱时，可以将非线性模型在o处线性近似，并用线性模型的求解理
7、论和方法来近似地求解非线性模型（6-2-3）式。这也就是我们大家所熟悉的传统方法线性化方法，即将非线性模型（6-2-3）式在处用台劳级数展开，取至一次项，得：SX-L-f（Xx=x（6-3-1）g）断*）g）l%唳X）死（顶）您药S_凯工=耳瓦壶项=与B =（6-3-2）（6-3-3）则（6-3-1）式可写为:（6-3-4）VB6X-l（6-3-4）式就是我们熟悉的间接平差的误差方程。由间接平差知，根据最小二乘原理可解得:（6-3-5）于是参数X的非线性平差结果为:（6-3-6）I.=心产X X.例6-1 （本例取自参考文献24）已知非线性模型为1。其中参数1和的真值为X （5.4201361
8、S7 0.25436189）r人古住-中*混古住 v 笛，曰、4口l 人n业士。的5个真值（用参数的真值X算得）和相应的5个同精度独立观测值列于表6-1。表6-14的真值和相应的观测值i12345真值4.2028343.2589242.5270061.9594691.519394观测值4.203.252.521.951.51观测方程为:乌二互舟+A3取参数X的近似值为繇一伉1得误差方程：ri%） 一（54 -0.3）。将观测方程在X口处线性近似，V = BSxl =0.74084-000431996、0.54885.92720.28640.40666.58640.32450.30126.505
9、80.3236k 0.22316.0245；03051)由（6-3-5）式得:I -0.005858021 3X =0.049953787 J于是，由(6-3-6)式得参数X的平差值为:X = Xq+3X =5394141979、-0.250246213；参数估值京的真误差为:-0.0259942000.004315630其范数为:宜= 0.02631801二、非线性最小二乘平差的迭代解当非线性模型的非线性强度很强时，线性近似可能产生大于观测误差的模型误差，所以对于非线性模型，一般采用迭代的方法求解。求解非线性误差方程（6-2-3）式的最小二乘平差值，就是求参数X的估值，使VTPVX)=fX-
10、L)TPfX)-(6-3-7)(6-3-8)=2fr (X) PL + rL = min由于PL是一常量，所以（6-3-7）式等价于目标函数为心)=f (艾)Pf(X)-2f(X)PL = min的非线性无约束最优化问题。因为是芳的非线性函数，所以对（6-3-8）式求一阶偏导数，并令其为零，得不至虐的显：4：表达式。故求不出甘的解析解。因此，我们只能设法寻找某一近似解X，使（6-3-9）1*4*成立。寻找使（6-3-9）式成立的近似解X，一般只有采用迭代的方法。为此，下面介绍几种常用的 迭代方法。1. 牛顿法设（*）的极小值X*的-个近似值为X*），在W近将A（）展为台劳级数，取至二 次项得：
11、二氏（工因+衣=氏吠 ） + ）赍3） + （疚林fG* 二 min（6-3-10）式中：=0劳中派 dRdRI， ， ， ， ， ，L鬼 t3(6-3-11)G?R d2Rd2R 蚌 弘如如如QC 瑟瓦=戏I II I I I E卷寸E寸RX = X闻鬼 吹既秘(6-3-12)*称为X）处的Hessian矩阵。*）=XX（6-3-13）W是配在彳何处的梯度方向。由于X）是X*的一个已知的近似值，故（6-3-10）式只是dX的函数，为了求得使（6-3-10）式成立的，将（6-3-10）式对求偏导，并令其为零，得：&凶+（疚）以=0移项后两边转置，顾及（6-3-12）式，得冬衣何二 一（/）（6
12、-3-14）当GkaE奇异时，由（6-3-14）式可解得使（6-3-10）式成立的宓=-（驴丁（6-3-15）当充分小时，能使（6-3-10）式成立。但由于X未知，故不能充分小，需不断迭代，直至充分小，其迭代公式为:事上+D =X“）二X时-G：遗f（6-3-16）（6-3-16）式就是牛顿迭代的基本公式，迭代终止条件:R（X（S）二 ）（6-3-17）=0（6-3-18）由于R（x）是一个绝对值较大的数，而#幻的各元素的绝对值都很小，因此，由于计算机有效数字的限制，以（6-3-17）式作为迭代收敛条件比（6-3-18）式作为迭代收敛条件收敛要快一些。牛顿法的迭代步骤为:(1)(2)并令k=0
13、。5J 5）(3)计算Hessian矩阵(4)解线性方程组（6-3-14）式，得(5)按（6-3-16）式计算新的近似值(6)计算目标函数值汽，若汽顷当次若）则转至继续迭代。(7)终止迭代，输出 和汽（*） ，结束。例6-2在例6-1中，仍设W 罗晋）*4，用牛顿法求例6-1中非线按（6-3-11）式计算梯度方向占，若占 =0则转至（7）。性模型的非线性最小二乘平差值。解：由例6-1知P=I，故目标函数为:将工代入计算 V，G0后，按以上迭代程序迭代，结果列于表6-2。表6-2牛顿法迭代计算k123456-1.2050249080.39918333820.028893980180.0001691624122-3.9492X10-9-2.4012X10-9-17.15305037.0372427130.49484074240.0029380122642.4039X10-7-1.5569X 10-75.3330132655.417198095.4227080035.4427445655.4227445935.422744582-0.2539145225-0.254257
