一．前言

数学是研究客观事物的数量关系和空间形式的科学。数学方法是对研究对象加以定量化，然后进行定量处理的方法。在测绘科学中，数学方法应用的范围和深度，已成为测绘科学发展的重要标志之一。

在全球定位系统（GPS）定位技术中，由于其测量定位技术的物理机制，其平面位置的精度可以达到很高的水准，而其高程精度较其平面精度约低2～5倍。尤其是在WGS84坐标系向1954年北京坐标系（1980年西安坐标系），1956年黄海高程系（或1985年黄海高程系）的转换过程中，由于WGS84的大地高仅有几何意义而无物理内涵，而黄海高程系统的正常高既有几何意义，又有地球内部质量密度分布不均匀这样一个物理现象。因此，目前此种转换是国际国内测绘学界积极探索的热门课题之一。

本文是在GPS技术测定高程的研究基础上，充分运用数学方法及数学函数所表达的抽象内涵，试图在GPS高程拟合这一课题中取得有益的研究成果。

二．数学方法的应用

1．基本概念

在GPS高程拟合中运用数学方法的关键，是针对GPS高程拟合的实际，根据似大地水准面的变化状况及水准点的分布情况，提炼一个合适的数学模型。

数学模型，是对于某个特定对象或一定问题，采用形式化数学语言来描述其特征及数量相依关系的一种数学结构，它是一组数学关系式或一套具体的数学算法。

数学模型方法，是通过建立和研究客观对象的数学模型来揭示对象本质特征和变化规律的方法。

要构建GPS高程拟合的数学模型，必须具备两个条件：（1）对GPS高程拟合机制的认识必须成熟到可以直接数量化的程度；（2）能找到适于刻划GPS高程拟合机制的数学手段。以上这两点已基本具备。

2．构造数学模型的基本步骤

（1）提出问题

GPS所测的高程是沿法线方向到WGS84椭球面的高度，即以简单的数学曲面为基准面，它具有明确的几何意义而缺乏物理意义，而工程测量中要求的正常高是沿垂线到似大地水准面的高度，即以不规则的有起伏的重力等位面为基准面，具有严格的物理意义，这两种基准面是不一致的，它们之间的差距称为高程异常。有关系式：ζ＝H－h

式中ζ为高程异常，表示似大地水准面至参考椭球面的距离。H为大地高，h为正常高。

在进行GPS测量后，由GPS三维平差可得到各点的大地高，若网中有部分GPS点是水准联测点，则这些点的正常高H是已知的，即可求得这些点的高程异常。在一定范围内高程异常不为常数但可以认为在此范围内变化平缓，可用一些数学函数来拟合，求得能反映GPS网控制范围中高程异常变化的函数，然后通过内插求得网中其它各类的高程异常。

（2）构造模型

在GPS高程拟合中，根据问题的实际，可采用曲线拟合法，曲面拟合法等确定性数学模型，而不采用随机性、模糊性、突变性数学模型。若测区承带状分布，则可采用曲线拟合，依据GPS点至起始重合点距离在首末重合点连线上的投影值，用多项式曲线函数拟合测段似大地水准面剖面曲线。

若测区似大地水准面起伏较大，且重合点数较多时，可采用曲面拟合来逼近似大地水准面。
为了说明问题，以曲面拟合二次多项式来表述，

（3）模型的计算和分析

式中A阵是曲面中二次型部分ζ*（x,y）的系数矩阵，它是二阶实对称阵，总存在正交变换，将其化为标准型：

若λ1λ2同号，则曲面为椭圆抛物面；若λ1λ2异号，则曲面为双曲抛物面。二次多项式表征的曲面分为且仅分为这两种类型。至于a0,a1x,a2y通过坐标变换都会吸收进平方项中，它们不会改变曲面的类型和形状。如果ζ*（x,y）=0，则曲面退化为平面。