半参数模型在GPS高程拟合中应用研究.docVIP

半参数模型在GPS高程拟合中应用研究 　　[摘要]介绍了几种GPS高程拟合的方法，通过对几种方法进行了比较以及试验分析，提出了将半参数模型应用到GPS高程拟合的实际应用中，通过对几种拟合模型进行分析，得出半参数模型在GPS高程拟合应用中精度更高，对建立GPS高程模型具有一定的参考意义。 　　[关键词]GPS高程拟合 高程异常 半参数模型 精度可靠性分析 　　[中图分类号] P228.4 [文献码] B [文章编号] 1000-405X（2014）-7-139-3 　　1引言 　　GPS测量所得的高程不具有物理意义。我国的高程系统采用正常高系统，因此，高程是相对于WGS-84椭球的大地高，它是一个几何量，仅具有几何意义，要用GPS定位技术确定观测点的正常高程，就必须将大地高转换成正常高，如图1所示。 　　获取精密的高程资料的主要工作在于通过GPS大地高的资料推算出点的正常高。利用GPS观测数据、重力数据和地球重力场模型实现GPS大地高向正常高的转换的是目前来说专家学者们认为的最好算法。目前我国似大地水准面（CQG-80）精度平均为2.7m，新的似大地水准面的精度在经济发达地区也仅有0.5m左右，难以满足很多测量工程的精度要求；除此之外，由于地面重力资料不齐全，在海洋地区、高山地区和边远地区等地区的重力资料相对较少，因此在实际运用中造难以操作。在布设GPS网时，通过水准点或者部分水准点对一些点进行统一测量，则这些点的正常高已知。如果能够将GPS测量资料和水准测量资料结合，精确地求定GPS点的正常高，必将严重影响到GPS在高程方面的应用。另外，也可以为应用到精密确定区域性似大地水准面的高程。本文介绍了几种GPS高程拟合的方法，提出引入一种半参数模型应用于GPS高程拟合中，通过实际比较分析得出了半参数模型应用于GPS高程拟合中精度可靠性。 　　2几种GPS拟合方法原理 　　2.1曲线拟合法 　　当GPS点布设成测线时，可通过曲线内插法求出选定点的正常高。 　　2.1.1多项式曲线拟合 　　设点的 和x（或y或拟合坐标）存在函数关系（i=0，1，...，n）可以用下面m（mn）次多项式来拟合： 　　ξ（x）=α0+α1x+α2x2+…+αmxm （1） 　　若mn，设Ri=ξ（xi）- ξi，在■ Ri2=min的条件下，通过求多项式系数αi和ξ，而获得该点的正常高。 　　2.1.2三次样条曲线拟合法 　　当布设的GPS网中的测线相对很长，结点ξ变化相对较大，插值点相对较多时，必须采用分段式计算来进行多项式拟合，在这种情况下，导致结点处的曲线不平滑，不连续，影响ξ的插值准确度，因此三样条曲线拟合法应运而生。 　　2.2曲面拟合法 　　在步设网状的GPS点时，正常高可以通过曲面拟合的方法求出，曲面拟合的基本原理是：根据GPS已知点的平面坐标和高程异常值，拟合出极大地水准面，通过内插法出待定点的高程异常就可以求出待定点的正常高。 　　2.2.1平面拟合法 　　在很小范围内，可以假设平面和局部似大地水准面相似。平面方法一般可表示为： 　　ξ（x，y）=α0+α1x+α2y （2） 　　式中，ξ（x，y）为与坐标值（x，y）相应的高程异常函数，αi（i=1，2，3）为模型参数。如果已知三个GPS点上的高程异常值，那么： 　　ξ=A?α （3） 　　从上式可得： 　　α= A -1ξ （4） 　　根据算出的模型参数，应用（4）可以计算出测区内任意点的高程异常。 　　如果已知高程异常值点的个数大于3时，通过最小二乘原理可知道： 　　α= （A TPA）-1ATPξ （5） 　　式中P为单位权。 　　2.2.2多项式曲面拟合 　　设点的ξ与平面坐标（x，y）有下列的关系式： 　　ξ=f（x，y）+ε （6） 　　式中f（x，y）为ξ的渐近值，ε为误差。下面假设： 　　f（x，y）= α0+α1x+α2y+α3x2+α4y2+α5xy+… （7） 　　变换成矩阵形式： 　　ξ=A?α+ε （8） 　　其中： 　　当m=n，有逆矩阵，如是α=A-1ξ 　　当mn，有无限解，通过最小二乘原理求得： 　　α= （ATPA）-1ATPξ （10） 　　2.2.3曲面样条拟合法 　　曲面样条拟合法是基于无限平板小扰度方程的数学模型，这点的异常值 和高与点的坐标x，y存在下列函数关系： 　　Xi，yi为已知点的坐标，x，y为待定点的坐标，Ai，Bi为待定系数。Pi为点的负载，D为刚度。对于么一个公共点都可以列出一个ξ（x，y）方程，对于n个公共点列出n+3个方程，求解出n+3个未知参数，α0， α1， α2，F1，F2，…，Fn求解方程组时至少应有三个公共点。 　　2.3多面函