多项式插值在GPS卫星坐标与钟差计算中应用.docVIP

多项式插值在GPS卫星坐标与钟差计算中应用 　　 摘要：本文介绍了用多项式内插计算GPS卫星坐标和钟差的方法，并以拉格朗日多项式为例利用IGS跟踪站给出的精密星历作为实例运用MATLAB进行编程计算，分析了Runge现象对拉格朗日多项式内插精度的影响，给出了拉格朗日内插法得到的卫星坐标和钟差的误差与多项式阶数、插值间隔的关系。 关键词：拉格朗日多项式；MATLAB；Runge现象；卫星坐标；钟差 Abstract: This article describes an algorithm for Lagrange Polynomial interpolation of GPS satellite coordinates and clock error and its implementation in MATLAB, then analyzes the effect of Runge phenomenon on the Lagrange Polynomial interpolation and the relations of satellite coordinates and clock errors computing by Lagrange polynomial interpolation with different orders of polynomial and fit intervals. Key words: Lagrange Polynomial; MATLAB; Runge phenomenon; Satellite coordinates; clock error 中图分类号：P185.18 文献标识码：A 文章编号： 1 引言 在GPS定位中，将GPS卫星的位置作为已知值，通过计算观测瞬间卫星至接收机的距离来计算测站点的坐标。因此，卫星的轨道信息是定位的基本要素，而卫星的轨道信息是通过星历给出的，星历误差将成为一种起算数据误差。这种误差直接影响着单点定位的精度，对绝对定位也有一定的影响。根据文献[1]的分析，在SA政策取消后，利用广播星历进行卫星导航和单点定位时，精度一般只能达到数米至数十米。一种有效的提高定位精度的方法是采用IGS精密星历，精密星历是按照一定的时间间隔(通常为15分钟)给出的，要想获得观测瞬间卫星的位置可以采用内插法求得，其中主要的内插方法有多项式内插法和三角内插法。根据分析，采用多项式内插法和三角内插法得到的精度基本相当。对于多项式内插法来说，要想获得理想的精度就要合理的选择多项式的阶数和插值的时间间隔（步长），本文以GPS基线解算软件中常用的拉格朗日多项式为例，通过比较计算来确定具体参数的选择。 2拉格朗日多项式插值原理 对于多项式插值，用户可以根据需要来选择不同的插值阶数和步长，插值多项式的标准形式为 (1) 其中为插值函数，点为插值节点，节点数目称为该插值方法的阶。在众多的多项式插值方法中拉格朗日插值多项式算法简单并且容易编程实现，因此在很多方面得到了应用，其定义如下： 设是区间[a,b]上的一个实函数，是[a,b]上n+1个互异实数，且在的值为，则区间[a,b]上任意一点x的n阶拉格朗日插值多项式的代数表达式为[2] (2) 利用(2)式进行卫星星历的内插计算，为观测时刻，为精密星历给出的时刻卫星的坐标或者钟差。因此，要想获得任意时刻卫星星历的n阶拉格朗日内插值，必须至少有n+1个时刻的卫星位置和钟差，并且内插时刻要在已知时刻的区间内。如果待求点位于插值区间的中央，则可以保证最佳的内插效果，即对称内插的方法，如图1所示（横轴为插值时刻）。 图1：对称内插法示意图 3 程序实现 用精密星历进行卫星星历的内插时将涉及到大量的计算，可以用计算机编程实现插值计算。MATLAB具有强大的科学计算功能，它不仅可以实现简单的求和、正弦、余弦，也可以方便的进行矩阵求逆、傅里叶变换等 [3][4]，并且MATLAB独特的以数组存储变量的方式可以方便的将读取的卫星星历数据按照原星历格式存储于计算机内存中。 作者在MATLAB7.0环境下编写了精密星历读取及拉格朗日多项式算法程序，部分代码如下： %......读精密星历文件................% function sp3=readsp3(filename) fid=fopen(filename,’rt’); sp3=textread(filename,’%s’); fclose(fid); % ---- Lagrange 插值多项式 ---- %X为插值节点，Y为相应的函数值 %x为插值时