IGS精密星历内插GPS卫星位置研究.pdfVIP

第 3期 矿 山 测 量 No．3 2010年 6月 MINE SURVEYING J．n．2010 doi：10．3969／j．issn．1001—358X．2010．03．019 IGS精密星历 内插 GPS卫星位置研究 任 迎华 (无锡水文工程地质勘察院，江苏 无锡 214063) 摘要 ：阐述 了用拉格 朗日多项式 内插计算GPS卫星轨道位置的方法．并利用 IGS跟踪站给 出的精密 星历作为实例运用NATLAB进行编程计算．给出了拉格 朗 日内插法得到的卫星位置误差与多项式 阶数的关系。结果表明，用拉格朗日多项式内插法得到的卫星位置精度能够满足精密定位的要求。 关键词：IGS精密星历 ；MATLAB；拉格 朗El多项式；卫星轨道 中图分类号：P228．4 文献标识码 ：B 文章编号 ：1001—358X 2【010)03—0054—03 在 GPS定位 中，将 GPS卫星的位置作为 已知 1．1 拉格朗日多项式插值的数学原理 值，通过计算卫星至接收机 的距离来计算测站点的 以IGS跟踪站提供给用户的精密星历为基础，使用 坐标。因此 ，卫星的轨道信息是定位的基本要素，而 插值方法计算出所需任意点的轨道状态量，插值法中使 卫星的轨道信息是通过星历给出的，星历误差将成 用的节点数 目称为该插值方法的阶。普遍使用的内插和 为一种起算数据误差。这种误差直接影响着单点定 外推的方法是拉格朗日多项式插值法[1]。其定义如下： 位的精度 ，对绝对定位 也有一定的影响。根据 IGS 设)，= )是区间[n，b]上的一个实函数 ， (i 提供的数据．在 SA政策取消后，利用广播星历进行 ： 0，1…，n)是 [口，b]上 n+1个互异实数 ，且 y： 卫星导航和单点定位时，精度一般只能达到数米至 )在 的值为Y = )，则区间[n，b]上任意一 数十米…。一种有效的提高定位精度的方法是采用 点 的n阶拉格朗 日插值多项式的代数表达式为 n n IGS精密星历 ，但是精密星历是按照一定的时间间隔 ^． 一 (通常为 15分钟)给出的，要想获得观测瞬间卫星的 ，( ：∑yFI( ) (1) t u J u — i —J J≠ i 位置可 以采用内插法求得，其中拉格朗 日(Lagrange) 式 中点 (i：0，1，…，)称为插值节点 ，包含插 法被广泛的采用。因此 ，就涉及到插值的阶数选择 值节点的区间[n，b]称为插值 区间。 问题 ，如果选用的阶数太低 ，精度