测量平差最小二乘法与数学模型.doc

测量平差由于测量仪器的精度不完善和人为因素及外界条件的影响，测量误差总是不可避免的。为了提高成果的质量，处理好这些测量中存在的误差问题，观测值的个数往往要多于确定未知量所必须观测的个数，也就是要进行多余观测。有了多余观测，势必在观测结果之间产生矛盾，测量平差的目的就在于消除这些矛盾而求得观测量的最可靠结果并评定测量成果的精度。测量平差采用的原理就是“最小二乘法”。测量平差是德国数学家高斯于1821～1823年在汉诺威弧度测量的三角网平差中首次应用，以后经过许多科学家的不断完善，得到发展，测量平差已成为测绘学中很重要的、内容丰富的基础理论与数据处理技术之一。最小二乘法 最小二乘法 在我们研究两个变量x，y）之间的相互关系时，通常可以得到一系列成对的数据x1，y1、x2y2……xm，ym）；将这些数据描绘在x-y直角座标系中如图1若发现这些点在一条直线附近，可以令这条直线方程如式1-1）。 Y计=a0+a1X1-1） 其中：a0、a1是任意实数 为建立这直线方程就要确定a0和a1，应用《最小二乘法原理》，将实测值Yi与利用式1-1）计算值Y计=a0+a1X的离差Yi-Y计的平方和∑Yi-Y计2最小为优化判据。 令φ=∑（Yi-Y计2 （1-2） 把式1-1）代入式1-2）中得 φ=∑（Yi-a0-a1Xi）2 （1-3） 当∑Yi-Y计平方最小时，可用函数φ对a0、a1求偏导数，令这两个偏导数等于零。 1-4） （1-5） 亦即：ma0+∑Xi）a1=∑Yi （1-6） （∑Xi）a0+（∑Xi2）a1=∑（Xi，Yi） （1-7） 得到两个关于a0、a1为未知数的两个方程组，解这两个方程组得出： a0=∑Yi）/m-a1（∑Xi）/m （1-8） a1=［∑XiYi-（∑Xi∑Yi）/m］/［∑Xi2-（∑Xi）2/m］ （1-9） 这时把a0、a1代入式1-1）中此时的式1-1）就是我们回归的元线性方程即：数学模型。 在回归过程中，回归的关联式是不可能全部通过每个回归数据点x1，y1、x2，y2……xm，ym）的，为了判断关联式的好坏可借助相关系数R”，统计量F”，剩余标准偏差S”进行判断；R”越趋近于1越好；F”的绝对值越大越好；S”越趋近于0越好。 R = [∑XiYi - m (∑Xi / m)(∑Yi / m)]/ SQR{[∑Xi2 - m (∑Xi / m)2][∑Yi2 - m (∑Yi / m)2]}1-10） 在式1-10）中，m为样本容量，即实验次数；Xi、Yi分别任意一组实验X、Y的数值。 最小三乘法 当研究实际中两个变量x，y）之间的相互关系时，也可得到一系列成对的数据x1，y1、x2，y2……xm，ym）；将这些数据描绘在x-y直角座标系如图2中发现这些点在一条曲线附近，假设这条曲线的一元非线性方程如式2-1。 Y计=a0+a1Xk 2-1） 其中：a0、a1、k是任意实数 为建立曲线方程，就要确定a0、a1和k值，应用《最小二乘法》同样的方法将实测值Yi与计算值Y计Y计=a0+a1Xik的离差Yi-Y计的平方和i-Y计）2为依据： 令φ=∑（Yi-Y计2 （2-2） 把式2-1）代入式2-2）中得： φ=∑Yi-a0-a1Xik）2 （2-3） 用函数φ分别对a0、a1和k求偏导数，令这三个偏导数等于零即： （2-4） （2-5） （2-6） 得到三个关于a0、a1和k为未知数的三元方程组，解方程组即可得到数学模型。 判断数学模型的好坏，同样可借助相关系数R”，统计量F”，剩余标准偏差S”进行判断；R”越趋近于1越好；F”的绝对值越大越好；S”越趋近于0越好。这样的验证很好时，有的模型计算误差还是很大，为了更进一步的验证数学