基于IGS精密星历的卫星位置内插方法比较.pdfVIP

2011年 1O月 城 市 勘 测 0et．2011 第5期 UrbanGeotechnicalInvestigation ＆ Surveying No．5 文章编号：1672—8262(2011)05—89-05 中图分类号：P228 文献标识码：A 基于 IGS精密星历的卫星位置内插方法比较 马俊 ，唐诗华，黄鹰，王文杰 (桂林理工大学 土木与建筑工程学院，广西 桂林 541004) 摘 要：在高精度的GPS测量数据处理过程中，获取高精度卫星轨道是重要的环节，内插是获取任意历元的精密轨 道信息的重要手段 ，本文利用IGS提供的精密星历用MATLAB作为工具编程计算，用实例分析和比较 了拉格 朗日 插值法、牛顿插值法、内维尔逐次插值法以及最小二乘拟合在求卫星坐标当中的应用效果。前三种方法在 16阶时 内插精度最高，内维尔插值和牛顿插值有很多相似的特性，方便灵活，龙格现象没有拉格朗 日插值那么剧烈，将会 得到广泛应用。 关键词：GPS精密星历；拉格朗日插值 ；牛顿插值；内维尔插值；最小二乘拟合 1 引 言 … ，求得在这些节点上n+1个 11,次拉格朗日基函 在GPS定位中，将 GPS卫星的位置作为已知值， 数2。()，f，()，z：()，…f()。 通过计算卫星至接收机的距离来计算测站点的坐标。 去 因此，卫星的轨道信息是定位的基本要素，而卫星的轨 道信息是通过星历给出的，星历误差将成为一种起算 = 兀一 (k=0，1，…，n) 数据误差。这种误差直接影响着单点定位的精度，对 k J 于是所求n次拉格朗日插值多项式为： 绝对定位也有一定的影响。为了降低星历误差对单点 定位和相对定位的影响，提高 目标在地心坐标系的绝 L()=Lykfk() … 对定位精度，IGS(InternationalGeodeticService)免费给 拉格朗日多项式函数模型简单，插值效率高，收敛 出了若干精度等级和不同时间延迟的精密星历。精密 速度快，是经典的插值方法，但是要增加或者删除一个 星历是按一定时间间隔(通常为 15min)给出相应历 插值节点需要重新构造多项式。 元时刻GPS卫星在地心坐标系的三维坐标、三维速度 2．2 牛顿多项式插值 及卫星钟改正数等信息。 定义 ： GPS精密数据处理需要得到间隔 1s的卫星坐标。 i^] 为 )在 i及 j两点的一 GPS接收机的采样率一般为 30s或者 l5s甚至更密， 因此，要想得到某瞬时的卫星坐标就必须对精密星历 阶差商。 和进行高精度、快速的内插或拟合。其中较常用的插 f~Xi,Xj,Xk]_j 为厂()在三点 i ， 值方法是拉格朗 日多项式插值法，和牛顿多项式插值 及 的二阶差商。 法。近来也有一