高程拟合的matlab代码,利用MATLAB结合C语言实现GPS高程二次曲面拟合

熊志强+肖腾飞

摘 要：对GPS高程曲面拟合的研究有很多，通常人们实现GPS高程曲面拟合都是用单一的编程语言。而文章则介绍结合两种编程语言更加方便地实现GPS高程二次曲面拟合。

关键词：GPS高程；二次曲面拟合；MATLAB；c语言

1 概述

随着GNSS技术的发展，GPS高程测量也更加受到重视，但GPS直接测出的为地面点的大地高，与水准高之间存在高程异常，因此，若想直接利用GPS测定水准高，需要拟合出当地的似大地水准面，通常采用二次曲面拟合，用编程实现。而编程语言较多，利用MATLAB和c语言各自的优势相结合实现二次曲面拟合会更加便捷。

2 GPS高程及二次曲面拟合原理

在实际应用中，地面点的高程采用正常高系统。地面点的正常高H_r是地面点沿铅垂线至似大地水准面的距离。这种高程是通过水准测量来确定的。地面点的大地高H是地面点沿铅垂线至参考椭球面的距离。这就有必要找出GPS点H与Hr的关系，并用一定的方法将H转换为Hr。似大地水准面至椭球面间的高差ζ，叫做高程异常。如果知道了各GPS点的高程异常ζ值，则可由各GPS点的大地高H求得各点的正常高Hr。

当GPS点布设成一定区域面时，可以应用数学曲面拟合法求待定点的正常高。其原理是，根据测区中已知点的平面坐标x，y和ζ值，拟合出测区似大地水准面，再内插出待求点的ζ，从