GPS精密星历外推插值方法研究.pdfVIP

146 童瞧晨舛技 2011年第1期 GPS精密星历外推插值方法研究 陈 卫，刘焕新，孙清奎 (徐州长城基础工程有限公司，江苏 徐州 221006) 摘 要 使用Lagrange插值、Newton插值和Neville插值三种算法对IGS精密星历进行外推插值计算，通过实验得到，当插值方法为Neville法且 插值阶数为11阶时运算精度和效率最优。 关键词 GPS IGS精密星历 外推插值 中图分类号P228．4 文献标识码 A StudyonpredictedinterpolatingofpreciseGPSephemeris ChenWei，LiuHuanxin，SunQingkui (xuzhouGreatwallFoundationEngineeringCO．，hdNO．87HuangheRoadXuZhou，221116) Abstract TheprincipleofLagrangeinterpolation，Ne~oninterpolationandNevilleinterpolationwereanalysizedtopredictthelacationofGPSsatellite withIGSpreciseephemeris．BasedOntheexperiment，theconclusionisobtained，whentheinterpolationalgorithm isNevilleandtheorderofinterpolation is11．theprecisionandefficiencyareoptimum． Keywords GPS 1GSpreciseephemeris predictinterpolatde ille内插。本文在比较以上三种插值方法基础上 ，选择 适用于 IGS精密星历外推插值的方法。 1 引言 2 插值数学模型原理 GPS定位是利用 GPS卫星在空中瞬时位置来计算 地面接收机三维坐标 ，所以精确确定 GPS卫星的位置 2．1 Lagrange插值 在GPS绝对定位和相对定位中有着至关重要的作用。 设有rt+1个节点时刻，设 f0()，llx)，…，z()都 目前求取卫星瞬间位置方法有两种 ：广播星历和精密 是 rt次多项式，且满足 星历。其中广播星历精度为 10m左右 ，GPS定位有一 Ii( o12，… (1) 定的误差；精密星历由IGS(国际GPS服务中心)提供， 精度优于5cm，可以满足大部分用户使用。 贝Ⅱ有 ()=yolo()+Y。Z。()+…+ Z() (2) 由于 IGS提供的精密星历历元间隔为 15rain，而接 为所求的Lagrange插值多项式。 收机接受的时间间隔为 1min左右，所以有必要利用合 为求得rt次多项式 z()，根据条件 (1)可知应该 理的数学模型计算所需历元时刻的卫星瞬时位置。目 有如下形式 ： 前利用 IGS精密星历进行内插求取卫星任意时刻位置 Z()=c(x—0…(--Xl—1)(一 +1)…(一 ) (3) 的内插方法主要有 Lagrange内插 、Newton插值和 Nev 其中c为待定系数，再根据条件 (1)的fl()=1 可知常数c为： $收稿 日期 ：2010—07—12