matlab内维尔差值,基于IGS精密星历的卫星位置内插方法比较

2011 年 10 月第 5 期 城 市 勘 测 Urban Geotechnical Investigation ＆ Surveying Oct． 2011 No． 5 文章编号:1672－8262( 2011) 05－89－05 中图分类号: P228 文献标识码: A 基于 IGS 精密星历的卫星位置内插方法比较马俊* ，唐诗华，黄鹰，王文杰 * 收稿日期:2011—04—01 作者简介: 马俊( 1985—) ，男，硕士研究生，主要研究方向: GPS 单点定位技术。 ( 桂林理工大学 土木与建筑工程学院，广西 桂林 541004) 摘 要: 在高精度的 GPS 测量数据处理过程中，获取高精度卫星轨道是重要的环节，内插是获取任意历元的精密轨道信息的重要手段，本文利用 IGS 提供的精密星历用 MATLAB 作为工具编程计算，用实例分析和比较了拉格朗日插值法、牛顿插值法、内维尔逐次插值法以及最小二乘拟合在求卫星坐标当中的应用效果。前三种方法在 16 阶时内插精度最高，内维尔插值和牛顿插值有很多相似的特性，方便灵活，龙格现象没有拉格朗日插值那么剧烈，将会 得到广泛应用。 关键词: GPS 精密星历; 拉格朗日插值; 牛顿插值; 内维尔插值; 最小二乘拟合 1 引 言 在 GPS 定位中，将 GPS 卫星的位置作为已知值，通过计算卫星至接收机的距离来计算测站点的坐标。因此，卫星的轨道信息是定位的基本要素，而卫星的轨道信息是通过星历给出的，星历误差将成为一种起算数据误差。这种误差直接影响着单点定位的精度，对绝对定位也有一定的影响。为了降低星历误差对单点定位和相对定位的影响，提高目标在地心坐标系的绝对定位精度，IGS( International Geodetic Service) 免费给出了若干精度等级和不同时间延迟的精密星历。精密星历是按一定时间间隔( 通常为 15 min) 给出相应历元时刻 GPS 卫星在地心坐标系的三维坐标、三维速度及卫星钟改正数等信息。 GPS 精密数据处理需要得到间隔 1 s的卫星坐标。GPS 接收机的采样率一般为 30 s或者 15 s甚至更密，因此，要想得到某瞬时的卫星坐标就必须对精密星历和进行高精度、快速的内插或拟合。其中较常用的插值方法是拉格朗日多项式插值法，和牛顿多项式插值法。近来也有一些学者和论述中提到内维尔逐次线性插值和切比雪夫多项式拟合法。本文着重分析了拉格朗日多项式插值、牛顿多项式插值、内维尔逐次插值以及最小二乘拟合算法，将其在精密星历卫星位置内插应用方面做了详细的比较，得出了一些有益的结论。 2 插值数学模型 2. 1 拉格朗日多项式插值 一般情况，对于给定的 n+1 个插值节点，x0 ＜x1 ＜x2 ＜…＜xn，求得在这些节点上 n+1 个 n 次拉格朗日基函 数 l0( x) ，l 1( x) ，l 2( x) ，…ln( x) 。 lk( x) = ( x－x0) …( x－xk－1) ( x－xk+1) …( x－xn) ( xk－x0) …( xk－xk－1) ( xk－xk+1) …( xk－xn) = ∏ n j=0 j≠k x－xj xk－xj ( k=0，1，…，n) 于是所求 n 次拉格朗日插值多项式为: Ln( x) = ∑ n k=0 yklk( x) 拉格朗日多项式函数模型简单，插值效率高，收敛速度快，是经典的插值方法，但是要增加或者删除一个插值节点需要重新构造多项式。 2. 2 牛顿多项式插值定义: f［xi，xj］= f［xj］－f［xi］xj－xi 为 f( x) 在