用岭估计研究不全3 维趋势面在GPS 工程水准测量中的应用.docVIP

用岭估计研究不完全3维趋势面在GPS工程水准测量中的应用 李崇贵　牟玉香　赵书河(石油大学建筑工程系　257062) 　　【摘　要】　为减小GPS工程水准测量中拟合误差的影响，本文采用岭估计代替最小二乘估计，用考虑地形起伏影响的不完全3维趋势面代替2维趋势面，确定求测区(似)大地水准面和参考椭球面差距的数学模型。并通过实例论证了不完全3维趋势面的预报效果优于2维趋势面。　　一、引　言测绘信息网　　GPS网的高程是相对于参考椭球面的大地高，工程测量中要求的是相对于(似)大地水准面的正高或正常高。如图1所示，hP为任一点P的大地高，HP为正(常)高，NP为(似)大地水准面与参考椭球面的差距，它们之间的关系为 NP=hP-HP　　　　　　　　　　　　　　　　　　　　　　(1) 图1 设在地面上有若干点，它们既有GPS高，又有正(常)高，可求出这些点处参考椭球面和(似)大地水准面的差距。由这些差距采用适当的数学方法，找出测区计算N值的回归方程，就能由GPS高和N值求出正(常)高，无需进行常规水准测量。求N值的精度除与GPS测量的精度有关外，还与已知水准点的精度、数量、在网中的分布及测区地形起伏的变化程度有关。回归模型的不合理、已知水准点的数量及分布可能造成回归模型未知量间存在复共线性。为研究测区地形起伏、已知水准点的不同分布及复共线性对求N值回归方程的影响，本文采用岭估计及不完全3维趋势面确定求测区N值的回归方程。测绘信息网　　二、岭估计控制LS估计的条件　　设有线性回归模型 Y＝Xβ+e, E(e)=0,COV(e)=σ2I　　　　　　　　　　　　(2) 式中，X为n×m的设计阵，且R(X)＝m。　　β的LS估计为 LS=(XTX)-1XTY　　　　　　　　　　　　　　　　(3) LS为β的无偏估计，其均方误差为测绘信息网 MSE(LS)=σ2tr(XTX)-1　　　　　　　　　　　　　　　　　(4) 　　记λ1≥λ2≥…≥λm≥0为XTX的特征根，因XTX可逆，(XTX)-1的特征根为：λ-11,λ-12,…，λ-1m，(4)式可变为 　　　　　　　　　　　　　　　(5) 　　在上式中，若XTX至少有一个特征根很小，即很接近于0，MSE（LS）就很大。这时，尽管高斯-马尔可夫定理保证了σ2tr（XTX）-1在线性无偏估计类中很小，但它本身的值却很大，则LS不是β的良好估计。XTX至少有一特征根很小，称设计阵X呈病态，或存在复共线性。设计阵X是否存在复共线性，可通过计算矩阵XTX的条件数或方差扩大因子来确定［1～2］。　　线性模型(2)的岭估计为测绘信息网 (k)=(XTX+kI)-1XTY　　　　　　　　　　　　　　(6) 式中，k0，为岭参数。　　与LS估计相比，岭估计把XTX换成了XTX+kI。当XTX的特征根至少有一个很接近于0时，XTX+kI的特征根λ1+k,λ2+k,…，λm+k接近于0的程度就会得到改善，使MSE（(k)）MSE(LS)，即在均方误差意义下，岭估计优于LS估计。用岭估计代替LS估计，需求使MSE（(k)）达到最小的k值，已证明这样的k一定存在［1～2］，本文采用岭迹法求k值［1］。　　　　三、不完全3维趋势面模型测绘信息网　　定义完全3维趋势面模型为 　　Ni=β1xi+β2yi+β3zi+β4x2i+β5y2i+β6z2i+β7xiyi+β8xizi+β9yizi+β10　　　(7) 式中，xi,yi,zi分别为各控制点的平面坐标及GPS高程，以10 km为单位。　　上式模型中包括高程项，即考虑了地形起伏对N值的影响。对具体的GPS网，模型可能不尽合理，可用岭迹法及方差扩大因子法对其中的变量进行取舍，去掉不合理项，精化后的模型称为不完全3维趋势面。　　四、用岭估计及不完全三维趋势面确定测区N值回归方程的实例 　　表1为某控制网的GPS定位成果，包括15个控制点的平面坐标、GPS高及正常高。测绘信息网 表1　某控制网的GPS定位成果 序号 点名 x/10 km y/10 km GPS高h/m 正常高H/m （H-h）/m 1 2001 4.856 220 -1.578 891 10.485 8.313 -2.172 2 2002 4.826 826 -1.013 125 8.913 6.547 -2.366 3 2008 3.748 925 -0.745 264 8.616 6.107 -2.509 4 2010 3.764 185 -1.850 432 9.013 6.192 -2.839 5 3015 4.627 002 -1.624 061 7.792 5.586 -2.206 6 3018 4.433 584 -1.700 734 8.508