基于拉格朗日与切比雪夫方法的精密星历插值研究
计算方法论文浅谈拉格朗日插值法
计算方法论文浅谈拉格朗日插值法拉格朗日插值法是一种常用的数值计算方法，用于构造一个多项式来逼近一些已知的离散数据点。
它被广泛应用于插值问题，如图像处理、物理实验数据处理、曲线拟合以及信号处理等领域。
本文将从原理、计算步骤以及优缺点三个方面，对拉格朗日插值法进行探讨。
拉格朗日插值法的基本原理是利用多项式的线性组合来逼近函数。
假设已知n+1个数据点：(x0, y0), (x1, y1), ... , (xn, yn)，其中x0, x1, ... , xn是互不相同的。
我们的目标是通过已知的数据点构造一个多项式P(x)，使得在这n+1个数据点上有P(xi) = yi。
根据插值定理，只要这些数据点满足一定的条件，存在唯一的插值多项式。
下面我们来具体讨论拉格朗日插值法的计算步骤。
首先，我们需要构造一个基于已知数据点的拉格朗日基函数。
对于每个数据点(xi, yi)，我们定义一个拉格朗日基函数Li(x)，它满足在xi处取值为1，而在其他数据点xj上取值为0。
拉格朗日基函数的定义如下：Li(x) = Π(j=0, j≠i, n)(x - xj) / Π(j=0, j≠i, n)(xi - xj)其中，Π表示一系列数的乘积符号。
接下来，我们需要将基函数与其对应的函数值进行线性组合，得到插值多项式P(x)。
插值多项式的表达式如下：P(x) = Σ(i=0, n)Li(x) * yi最后，我们可以利用插值多项式来计算任意点的函数值。
拉格朗日插值法的优点在于相对简单和容易理解，它能够精确地通过已知的n+1个数据点来构造一个次数不超过n的多项式，实现对函数的逼近。
然而，拉格朗日插值法也存在一些缺点。
首先，拉格朗日插值法对于数据点的选择非常敏感，如果数据点的密度不均匀或者存在较大误差，那么插值结果可能会出现较大的误差。
此外，拉格朗日插值法在计算多项式系数时需要进行大量的乘法和除法运算，这在数据规模较大时可能会导致计算效率降低。
拉格朗日插值法在数值分析中的应用研究
拉格朗日插值法在数值分析中的应用研究拉格朗日插值法是一种常用的数值分析方法，广泛应用于函数逼近、数据拟合、信号处理等领域。
本文将探讨拉格朗日插值法的原理、优缺点以及其在数值分析中的具体应用。
一、拉格朗日插值法原理拉格朗日插值法基于一个简单的思想：通过已知的离散数据点，构建一个多项式函数，该函数能够在给定的区间内，以已知数据点为插值节点，对未知数据进行逼近。
插值的多项式函数称为拉格朗日插值多项式。
设已知的离散数据为{(x₀, y₀), (x₁, y₁), ..., (xₙ, yₙ)}，其中xi为已知的节点，yi为相应数据点的函数值。
拉格朗日插值多项式L(x)可以表示为：L(x) = Σ(yᵢ * Li(x))其中Li(x)称为基函数，满足条件：Li(xi) = 1，Li(xj) = 0 (i ≠ j)。
二、拉格朗日插值法的优缺点拉格朗日插值法具有以下几个优点：1. 简单易懂：拉格朗日插值法的原理简单明了，易于理解和实现。
2. 精度较高：在节点较密集的情况下，拉格朗日插值多项式可以准确地逼近原始函数。
3. 适用范围广：拉格朗日插值法适用于各种类型的数据，包括等间隔数据和非等间隔数据。
然而，拉格朗日插值法也存在一些缺点：1. 多项式次数过高时，可能出现龙格现象：在某些情况下，拉格朗日插值多项式次数过高会引起振荡，降低插值的准确性。
2. 对于大规模数据的计算量较大：当节点数量较多时，计算拉格朗日插值多项式的复杂度较高。
三、拉格朗日插值法的应用拉格朗日插值法在数值分析中有着广泛的应用，以下是几个常见的应用场景：1. 数据拟合：给定一组离散数据点，我们可以使用拉格朗日插值法拟合出一个多项式函数，从而对未知的数据点进行估计。
这在科学实验中常用于实验数据处理和结果预测。
2. 函数逼近：对于已知的函数，我们可以通过设定一组插值节点，使用拉格朗日插值法将这个函数逼近为一个多项式函数。
这在数学建模和函数分析中非常有用。
GPS精密星历插值法与拟合法的精度分析
始 时间 ，△ 为拟合 区间的长 度 。 将变量 ｔ Ｅ。 ￡ 先 ∈ ｔ，
ｔ＋ Ａ ］变为 变量 ｒ ［ １ １ 。 ｔ ∈ 一 ，］
ｒ一 （ 一 ）一 １ （） ２
时刻 的的卫星 坐标 ， 因此需要 对精 密星历 进行插 值 或 拟合 以获得 采样 间隔 更 高 的卫 星 坐标 。常用 的 插 值法 和拟合 法有 拉 格 朗 日多项 式 插 值 和切 比雪 夫 多项 式拟合 。 ］
需 用精密 星历 。
假设 需要在 时间 间隔 Ｅ。ｔ＋ａ ］上用 阶切 ｔ， 。 ｔ
比雪 夫多项 式拟合 卫 星轨 道 ，其 中 ｔ 为 拟 合 的初 。
Ｉ Ｓ发布 的精密 星历是 以 １ ｎ为间 隔 的卫 Ｇ ５ｍｉ
星坐标口 ， ］ 而在 Ｇ Ｓ事后 数据 处理 中 ， 用 到任 意 Ｐ 要
中广播 星历 的精 度 偏 低 ， 为 ２ｍ［ ， 由最 终 精 约 １而 ］
象［ ， 以应尽 可能使 插值 点位 于插值 区间 的 中间 ５所 ］
部分 。
Ｉ ２ 切 比 雪 夫 多 项 式 拟 合 原 理 ． 。
密星 历 （ＧＳ 提 供 的精 密 星 历 精 度 较 高 ， 于 ５ Ｉ ） 优 ｃ 在对 精度 要求 较高的 定位 工 作 中 ， 据处 理 时 ｍ． 数
时刻 ｔ 可 以用 拉 格 朗 日多 项 式 来 计 算 卫 星 的 坐 ，
标Ｅ ：
ｒ卜 ｒ１ 井 １
．
ｚ￡一 ∑Ｃ， （ ） ｚ （ 丁
ｆ ０ ＝
式 中，
， ห้องสมุดไป่ตู้， Ｃ
为切 比雪夫多 项式 的系数 ， 是
厂￡一∑Ⅱ（ ｚ （ ） ÷）
收稿 日期 ｉ ０ ０ １ －６ １ — ２１ ２
一种新的精密星历内插方法
一种新的精密星历内插方法孙鹏;赵长胜;王仁【摘要】Precise Point Positioning(PPP) is based on GPS precise ephemeris interpola‐tion .This article uses Generalized Barycentric Rational Interpolation for precise ephemeris interpolation ,and verifies that Sliding Generalized Barycentric Rational Interpolation is also feasible .T he experiment show s that ,Generalized Barycentric Rational Interpolation is pre‐cise and highly stable and can meet the demands of PPP.%G PS精密星历插值是精密单点定位的基础，本文将广义重心有理插值应用于G PS精密星历插值，同时验证了滑动广义重心有理插值算法在精密星历插值中的可行性，实验表明，滑动广义重心有理插值具有较好的稳定性及很高的内插精度，可以满足精密定位的需要。
【期刊名称】《全球定位系统》【年(卷),期】2015(000)006【总页数】3页(P89-91)【关键词】精密星历;拉格朗日插值;重心拉格朗日插值;广义重心有理插值【作者】孙鹏;赵长胜;王仁【作者单位】江苏师范大学测绘学院，徐州221116;江苏师范大学测绘学院，徐州221116;江苏师范大学测绘学院，徐州221116【正文语种】中文【中图分类】P228.4获取GPS卫星的在轨位置是GPS定位的关键问题,获取卫星坐标有两种方式,一种是广播星历,一种是精密星历。
广播星历通过导航电文的数据块Ⅱ直接发送给用户接收机,精度较低,不能满足精密定位需要。
数值分析实验报告--实验2--插值法
1 / 21数值分析实验二：插值法1 多项式插值的震荡现象1.1 问题描述考虑一个固定的区间上用插值逼近一个函数。
显然拉格朗日插值中使用的节点越多，插值多项式的次数就越高。
我们自然关心插值多项式的次数增加时， 是否也更加靠近被逼近的函数。
龙格（Runge ）给出一个例子是极著名并富有启发性的。
设区间[-1,1]上函数21()125f x x=+ (1)考虑区间[-1,1]的一个等距划分，分点为n i nix i ,,2,1,0,21 =+-= 则拉格朗日插值多项式为201()()125nn ii iL x l x x ==+∑(2)其中的(),0,1,2,,i l x i n =是n 次拉格朗日插值基函数。
实验要求：(1) 选择不断增大的分点数目n=2, 3 …. ，画出原函数f(x)及插值多项式函数()n L x 在[-1，1]上的图像，比较并分析实验结果。
(2) 选择其他的函数，例如定义在区间[-5，5]上的函数x x g xxx h arctan )(,1)(4=+=重复上述的实验看其结果如何。
(3) 区间[a,b]上切比雪夫点的定义为 (21)cos ,1,2,,1222(1)k b a b ak x k n n π⎛⎫+--=+=+ ⎪+⎝⎭(3)以121,,n x x x +为插值节点构造上述各函数的拉格朗日插值多项式，比较其结果,试分析2 / 21原因。
1.2 算法设计使用Matlab 函数进行实验, 在理解了插值法的基础上，根据拉格朗日插值多项式编写Matlab 脚本，其中把拉格朗日插值部分单独编写为f_lagrange.m 函数，方便调用。
1.3 实验结果1.3.1 f(x)在[-1,1]上的拉格朗日插值函数依次取n=2、3、4、5、6、7、10、15、20，画出原函数和拉格朗日插值函数的图像，如图1所示。
Matlab 脚本文件为Experiment2_1_1fx.m 。
可以看出，当n 较小时，拉格朗日多项式插值的函数图像随着次数n 的增加而更加接近于f(x)，即插值效果越来越好。
基于拉格朗日的IGS精密星历和钟差插值分析
基于拉格朗日的IGS精密星历和钟差插值分析
段清超;许宝成
【期刊名称】《水利科技与经济》
【年(卷),期】2016(022)009
【摘要】在GPS精密单点定位计算中,需要高精度卫星轨道位置及高采样率的钟差产品.采用滑动式拉格朗日内插方法,利用Matlab编程,对国际地球动力学服务机构(IGS)提供的精密星历和钟差进行加密计算,通过改变插值阶次和采用不同的采样率钟差产品进行插值计算,得出一些有益的结论.
【总页数】4页(P50-53)
【作者】段清超;许宝成
【作者单位】华北水利水电大学资源与环境学院,郑州450000;郑州工商学院建筑与测绘工程系,郑州450000
【正文语种】中文
【中图分类】TV212
【相关文献】
1.基于IGS精密星历的卫星坐标和钟差插值 [J], 宫厚诚;李全海
2.基于拉格朗日插值方法的GPS IGS精密星历插值分析 [J], 何玉晶;杨力
3.广义延拓插值法在 IGS 精密钟差插值中的应用 [J], 焦宁;田龙华;孙秀宁;韦铖
4.IGS精密星历和钟差的算法比较研究 [J], 张养安;李俊锋;薛兆元;李飞
5.基于频谱分析的IGS精密星历卫星钟差精度分析研究 [J], 黄观文;张勤;许国昌;王利
因版权原因，仅展示原文概要，查看原文内容请购买。
计算方法拉格朗日插值
计算方法拉格朗日插值拉格朗日插值是一种用于在给定数据点间进行插值的方法，它基于拉格朗日多项式的性质来进行计算。
拉格朗日插值可以用于任何数量的数据点，无论是线性插值还是高阶插值。
拉格朗日插值的基本思想是，使用多个插值点的拉格朗日多项式来逼近给定数据点。
具体而言，对于给定的插值点(x0, y0)，(x1, y1)，...，(xn, yn)，我们需要找到一个多项式P(x)来满足以下条件：P(xi) = yi，其中 i = 0, 1, ..., n。
假设我们要计算的插值点为x，那么根据拉格朗日插值的公式，多项式P(x)可以写为：P(x) = Σyi * Li(x)，其中 i = 0, 1, ..., n。
在上述公式中，Li(x)是拉格朗日基函数，可以用以下公式表示：Li(x) = Π(x - xj) / Π(xi - xj)，其中j ≠ i，i, j = 0,1, ..., n。
现在我们可以根据上述公式进行计算，以下是拉格朗日插值的详细步骤：1. 输入数据点的坐标 (x0, y0)，(x1, y1)，...，(xn, yn) 和待插值点的坐标 x。
2. 对于每个插值点(xi, yi)，计算拉格朗日基函数Li(x)。
3. 对于每个插值点(xi, yi)，计算插值多项式中对应的项 yi *Li(x)。
4.将所有项相加，得到插值多项式P(x)。
5.根据插值多项式P(x)，计算插值点x的函数值，即P(x)=y。
拉格朗日插值的优点是简单易懂，计算过程相对简单，但它也存在一些缺点。
拉格朗日插值的计算复杂度为O(n^2)，这意味着当数据点的数量较多时，计算会变得非常耗时。
此外，拉格朗日插值在边界点附近的插值结果可能会出现较大的误差。
为了减小计算量和提高插值的准确性，还有其他更高效的插值方法，如牛顿插值和样条插值。
这些方法在实际应用中经常被使用，具有更好的性能和更准确的插值结果。
利用精密星历计算北斗卫星坐标的方法探讨
利用精密星历计算北斗卫星坐标的方法探讨苏怡婷;肖琴琴【摘要】本文以北斗卫星的精密星历作为实验数据,采用不同的拟合方法计算卫星的三维坐标,对3种常用的方法,即多项式拟合法、拉格朗日插值法和切比雪夫多项式拟合法,所计算的卫星坐标结果进行分析与比较,最后得到了3种拟合方法中卫星位置误差与阶数的关系,得出了适合计算任意时刻北斗卫星坐标的拟合方法和阶数,并获得了不同轨道卫星与拟合方法的关系,可为今后的研究提供一定的理论基础.【期刊名称】《湖南城市学院学报（自然科学版）》【年(卷),期】2018(027)006【总页数】6页(P57-62)【关键词】精密星历;卫星坐标;多项式拟合法;拉格朗日插值法;切比雪夫拟合法【作者】苏怡婷;肖琴琴【作者单位】天津城建大学地质与测绘学院,天津 300384;湖南城市学院市政与测绘工程学院,湖南益阳 413000【正文语种】中文【中图分类】P228.1北斗卫星导航系统(Beidou Navigation Satellite System, BDS)是我国自主研制的一种全球卫星导航系统，是继全球定位系统(Global Position System, GPS)、伽利略卫星导航系统(Galileo Satellite Navigation System, GLONASS)之后第三个成熟的卫星导航系统[1]，是军事、民用和科学研究中发挥着重要作用的基础设施﹒近几年随着计算机技术的改进更新，北斗系统发展迅速，至2018年年底，北斗卫星导航系统将服务“一带一路”沿线国家和地区[2]﹒北斗卫星导航系统的功能有很多，其中获取卫星的位置是处理问题的关键，由于子午仪系统存在的卫星较少、不能实时定位和有间隔时间等影响，所以通过对国际卫星导航服务(Internation GPS Service for Geodynamics, IGS)等国际组织发布的精密星历(SP3格式)[3]进行内插，获取任意时刻的卫星坐标，从而进行卫星定位﹒北斗卫星星座是由中地球轨道卫星(Middle Earth Orbit, MEO)、倾斜地球同步轨道卫星(Inclined Geosynchronous Satellite Orbit, IGSO)、地球静止轨道卫星(The Geostationary Orbit, GEO)组成的混合星座[4]，而卫星轨道的不同也是影响拟合方法适用性的一个重要因素﹒本文采用2018年间隔5 min的不同轨道北斗卫星精密星历作为实验基础数据，根据多项式拟合法、切比雪夫拟合法和拉格朗日插值法的拟合原理，结合Matlab 软件编程[5]，并通过计算得出拟合误差，分析探讨各方法的优势与劣势，进而得出一些理论为今后的深入研究提供参考﹒多项式函数[6]是形式比较简单的函数，即从本质上讲，多项式拟合也是一个线性模型，其数学表达式为其中，n为多项式最高次数；代表的是x的i次幂；是的系数﹒对于卫星坐标，用拉格朗日插值法分别对X坐标、Y坐标和Z坐标进行拟合，已知函数的n+1个节点，，，，及其相应函数，，，，，对于插值区间任一点t，都可用如下公式计算[7-8]﹒其中t属于[0, 135]，并且在选取区间内满足所选取的卫星时刻至少为n+1个时刻，这里n为插值阶数，n+1为插值节点个数﹒利用公式计算切比雪夫多项式系数，设t属于[0, 135]，作以下变换[9-11]：切比雪夫多项式的递推公式为卫星坐标可用多项式表示为采用MEO、IGSO、GEO卫星2018年2月11日、2018年4月14日及2018年8月1日的144组精密星历作为实验数据，编写多项式拟合程序进行计算，得出残差表如表1~表3所示﹒从多项式拟合结果来看，MEO卫星和IGSO卫星，虽然2个轨道卫星的误差有所不同，但是总体走向大体一致﹒随着阶数增加，拟合误差明显减小；但阶数继续增加时，误差明显增大；拟合中间节点到两端节点呈逐渐增大的趋势；端点坐标的误差最大﹒对于MEO卫星，9阶和10阶误差较小，10阶Z轴误差较大，接近1 cm；IGSO卫星多项式拟合9阶误差最小；其它阶数拟合的误差含有异常值，拟合精度没有9阶高﹒GEO卫星与其它2个轨道的拟合误差有较大区别﹒在X轴上所有阶数误差都较大，从起始节点开始误差呈增大的趋势，末端节点的误差最大；Z轴所有阶数拟合误差相似，除去端点误差，其余节点拟合误差在5 mm内，误差较小；Y轴拟合误差随阶数增加而减小，但阶数到10以上时，误差增大，其9阶拟合误差较小，在2 dm内﹒综合3个坐标来看，拟合效果无法满足高精度卫星的定位要求，所以，多项式拟合法不适合GEO卫星坐标的拟合﹒采用MEO、IGSO、GEO卫星2018年2月11日、2018年4月14日及2018年8月1日的144组精密星历作为实验数据，利用拉格朗日插值程序进行计算，得出残差表如表4~表6所示﹒结果显示，拉格朗日插值法9阶及以下插值因实验数据较少，故拟合精度不高，误差值达到米级﹒MEO卫星9阶端点拟合误差在分米级，10阶以上拟合效果较好，拟合误差在4 mm内；IGSO卫星10阶以上拟合误差在8 mm内，11阶拟合误差小于2 mm，除去边缘误差，其余节点误差不超过1 mm，拟合效果最好，13阶拟合误差明显增大﹒其它所有阶数拟合效果相同的地方就是在首尾两端误差明显增大，出现不收敛的现象，这种现象叫做“龙格现象”，消除“龙格现象”的方法是尽可能使插值点位于插值区间的中间部分﹒GEO卫星X轴7~9阶右端节点处的拟合误差最大，随阶数增大整体误差逐渐减小，10~13阶拟合误差相对较小；Y轴10阶以上的拟合误差均在1 mm内；Z轴拟合误差没有明显的规律，7阶拟合实验数据较少，且误差较大﹒综合分析可知，11阶拟合误差趋于0，整体拟合效果最好﹒采用MEO、IGSO、GEO卫星2018年2月11日、2018年4月14日及2018年8月1日的144组精密星历作为实验数据，利用切比雪夫拟合程序进行计算，得出残差表如表7~表9所示﹒综合来看，切比雪夫拟合与拉格朗日拟合的效果相似，9阶误差达到厘米级，其它阶数拟合误差都在7 mm内，较大误差集中在端点﹒IGSO卫星与MEO卫星11阶拟合的X轴与Z轴拟合误差在1 mm内，与拉格朗日插值法11阶拟合效果相同，相比于其它阶数拟合，效果较好﹒GEO卫星坐标拟合9阶及以上与同卫星的拉格朗日法拟合效果相同，X轴上整体拟合效果较好，所有阶数拟合误差在8 mm内，11阶相对于其它阶数拟合效果较好，11阶的Y轴拟合误差在1 mm内﹒综合以上分析，11阶拟合效果可以满足高精度卫星的定位要求﹒通过在官网上下载SP3精密星历文件，选取起始历元在2018年2月11日0时0分0秒、2018年4月14日6时0分0秒以及2018年8月1日0时0分0秒的SP3精密星历文件作为实验基本数据﹒选取不同轨道卫星不同拟合方法的最优拟合阶数，即9阶多项式拟合、11阶拉格朗日插值拟合和11阶切比雪夫拟合，对其误差数据进行分析比对，通过以上3种方法，还可以拟合出任意时刻的卫星坐标﹒为了检查拟合的质量，选取量测时间间隔为5 min的北斗卫星PC10、PC11和PC01坐标，根据残差对其进行评定(在本文中，所得到的残差为相应时刻的卫星坐标减去所拟合的坐标)，然后将3种方法所得之最优精度拟合阶数的误差结果进行比较，得出最优拟合的方法，结果分别如表10~表12所示﹒群体标准偏差公式：其中，是各个坐标的误差；是3个坐标误差的平均值，即；N为样本个数﹒从以上结果看来，本次实验的综合结论是采用11阶的拉格朗日插值方法和11阶切比雪夫拟合法可以满足卫星坐标拟合的高精度要求，除去端点误差，拟合误差均在5 mm内﹒1)实验结果证明了MEO卫星轨道计算方法适用于IGSO卫星，但不适合GEO卫星(其轨道接近于0°)﹒2)多项式拟合法9阶拟合效果最好，但不适合GEO卫星的拟合，11阶拉格朗日插值法与真实卫星坐标最为接近，但拉格朗日法最大的不足是，当需要增加新的节点时，原有的编程公式要重新建立，此外还需要考虑节点间的距离与插值次数，若结点间距离很小但插值次数较高时就容易出现“龙格现象”，消除“龙格现象”的方法是尽可能使插值点位于插值区间的中间部分﹒3)切比雪夫拟合法在利用精密星历计算北斗卫星坐标的实践中，10阶拟合精度也比较好，但其数学模型相对复杂，计算中涉及矩阵求逆等运算，运算量较大，拟合效果次于拉格朗日10阶的拟合效果﹒从误差比较和实际操作来看，拉格朗日插值法最适合此次实例研究，且原理简单，可为生产实践提供参考﹒【相关文献】[1]赖山东, 王有亮, 黄河清. 利用精密星历计算卫星位置方法的比较[J]. 江西测绘, 2009(2): 35-37.[2]刘季. 北斗GEO卫星位置计算方法探究[J]. 测绘地理信息, 2012, 37(5): 33-36.[3]洪樱, 欧吉坤, 彭碧波. GPS卫星精密星历和钟差三种内插方法的比较[J]. 武汉大学学报: 信息科学版, 2006, 31(6): 516-518.[4]陈端阳, 王忠军, 陈洪卿. 北斗系统时间(BDT)的认知与应用[J]. 数字通信世界, 2013(8): 44-48.[5]万亚豪, 张书毕, 侯东阳. GPS精密星历插值法与拟合法的精度分析[J]. 全球定位系统, 2011(2): 40-44.[6]魏二虎, 柴华. GPS精密星历插值方法的比较研究[J]. 全球定位系统, 2006(5): 13-15.[7]赵辉, 张书毕, 张秋昭, 等. 基于质心拉格朗日插值的GPS轨道标准化方法[J]. 全球定位系统, 2011(2): 15-18.[8]彭泽泉. GPS精密星历拟合方法的研究[J]. 测绘科学, 2010, 35(增1): 63-65.[9]余鹏, 孙学金, 赵世军. GPS定位中卫星坐标计算的切比雪夫多项式拟合法[J]. 气象科技, 2004, 32(3): 198-200.[10]杨学锋, 程鹏飞, 方爱平, 等. 利用切比雪夫多项式拟合卫星轨道坐标的研究[J]. 测绘通报, 2008(12): 1-3.[11]孔巧丽. 用切贝雪夫多项式拟合GPS卫星精密坐标[J]. 测绘通报, 2006(8): 1-3.。
高等数值分析拉格朗日插值多项式切比雪夫高斯龙格现象复合梯形辛普森求积公式
高等数值分析拉格朗日插值多项式切比雪夫高斯龙格现象复合梯形辛普森求积公式解答：1.拉格朗日插值函数：function y=lagrange (a,b,x)y=0;if length(a)==length(b)n=length(a);else disp('ERROR!length(a)!=length(b)')return;endfor i=1:nk=1;for j=1:nif j~=ik=k.*(x-a(j))/(a(i)-a(j));endendy=y+k*b(i);end2.问题（a）：function Q_am=100;n=10;x=-1:2/n:1;y=1./(1+9*x.^2);x0=-1:2/m:1;y0=lagrange(x,y,x0);y1=1./(1+9*x0.^2);plot(x0,y0,'--r');hold on;plot(x0,y1,'-b');end3.问题（b）：function Q_bm=100;n=10;x=zeros(1,n+1);for i=1:n+1x(i)=cos((2*i-1)*pi/(2*n+2)); endy=1./(1+9*x.^2);x0=-1:2/m:1;y0=lagrange(x,y,x0);y1=1./(1+9*x0.^2);plot(x0,y0,'--r');hold on;plot(x0,y1,'-b');end4.问题（c）：main.m（m文件）figure(1)Q_a()figure(2)Q_b()syms xy=1/(1+9*x^2);I0=int(y,-1,1);%准确值n=10;x=-1:2/n:1;y=1./(1+9*x.^2);I1=trapz(x,y);%复合梯形x0=zeros(1,n);for i=1:nx0(i)=(x(i)+x(i+1))/2;endy0=2/n*1./(1+9*x0.^2);I2=I1/3+2*sum(y0)/3;%复合辛普森x1=[-0.5384693101 0.5384693101 -0.9061798459 0.9061798459 0];y1=1./(1+9*x1.^2);A=[0.4786286705 0.4786286705 0.2369268851 0.2369268851 0.5688888889]; I3=y1*A'; %高斯5总结：(1).使用等距节点构造的高次拉格朗日插值多项式在正负1附件，插值值与真实值偏差非常大，存在较大的震荡。
切比雪夫多项式拟合GPS卫星星历精度分析
根据 最 小 二 乘 原 理 可 知 Ｃ：Ｎ Ｂ ， 其 中 Ｎ＝
Ｂ Ｐ Ｂ 。将 求得 的 ｃ代 回式 （ ２ ） ， 即 可求 出任 意 时 刻 的Ｇ Ｐ Ｓ卫 星 的坐标 。 根据 Ｃ ｈ ｅ ｂ ｙ ｓ ｈ ｅ ｖ多项 式 （ ｒ ） 结构 ， 可 以看 出在 区
２ ０ １ ５ 年 ８月
第 ４期 文章编号 ： １ ６ ７ ２ — ８ ２ ６ ２ （ ２ ０ １ ５ ） ０ ４ － ７ ４ － ０ ４
城
市
勘
测
Ａｕ ｇ ． ２ ０ １ ５
Ｎｏ ． ４
Ｕｒ ｂ ａ ｎ Ｇｅ ｏ ｔ ｅ ｃ ｈ ｎ ｉ ｃ ａ ｌ Ｉ ｎ ｖ ｅ ｓ ｔ ｉ ｇ ａ ｔ ｉ ｏ ｎ＆ Ｓ ｕ ｒ ｖ ｅ ｙ ｉ ｎ ｇ
３ ． 广 西空 间信息与测绘重点实验室 ， 广西 桂林
摘
５ ４ １ ０ ０ ４ ）
要： 采 用切 比 雪 夫 多项 式 拟 合 Ｇ Ｐ Ｓ精 密星 历 ， 节 点 的 选 择 对 精 度 产 生 一 定 的 影 响 。利 用 Ｉ Ｇ Ｓ提 供 的 ３种 精 密 星
历（ 最终 星历一 Ｉ Ｇ Ｆ 、 快速 星历 一 Ｉ Ｇ Ｒ、 超快速 星历一 Ｉ Ｇ Ｕ） 分析 随机 选择插值 节 点和选择 切 比雪夫点 拟合轨 道 的精 度差 异， 并比较 ３种星历拟舍轨道和 Ｉ Ｇ Ｓ提供 的轨道 差值 ， 分析 其精度 ， 结果表明采用切比 雪夫点对精密星历进行拟合 ， 最
中图分类号 ： Ｐ ２ ２ ８
文献 标 识 码 ： Ａ
切比雪夫多项式拟合 Ｇ Ｐ Ｓ卫星星历精度分析
杨兴 跃 ， 任 超 ， 吕东 ， 孙 建伟
５ ４ １ ０ ０ ４ ； （ １ ．南 阳市 城 乡 规 划 测 绘 院 ， 河南 南 阳 ４ ７ ３ ０ ０ ０； ２ ． 桂 林 理工 大学 ， 广西 桂林
三种精密星历内插方法的比较
收稿日期：2015-05-06作者简介：祝程程（），女，硕士，山东英才学院建筑工程学院教师。
三种精密星历内插方法的比较祝程程（山东英才学院建筑工程学院，山东济南250104）摘要：GPS （Global Pos itioning Sys tem ）主要包括三种精密星历的内插方法，分别是拉格朗日内插、牛顿内插、契比雪夫多项式拟合。
作者运用MATLAB 编程计算阐述了内插精度与内插阶数的关系，并通过实例计算得到了一些比较有意义的结论。
关键词：精密星历；拉格朗日多项式；牛顿多项式；契比雪夫多项式0引言在G PS 定位中，将GPS 卫星的位置作为已知值，通过计算卫星至接收机的距离来计算测站点的坐标。
因此，卫星的轨道信息是定位的基本要素，而卫星轨道信息是通过星历给出的，星历误差将成为一种起算数据误差[1]。
在GPS 数据后处理中，用户通常是需要根据已有的精密星历来选择合理的插值或拟合方法解算任意观测时刻的卫星位置。
本文就是使用IGS 精密星历，分别采用拉格朗日插值方法、牛顿插值方法以及契比雪夫多项式拟合方法进行精密星历内插解算，从而检验其是否满足进行精密定位的要求，并对三种方法进行了一定的比较。
1内插模型的计算原理1.1拉格朗日插值设在n+1个时间为121,,......,nt t t 插值节点上卫星的坐标分别是112211(),(),......,()n n x t x t x t ，那么在任一时刻t 卫星的坐标可以表示为1111()[]n n j ii j ijj i t t x t x tt （1）利用公式（1）在X 、Y 和Z 方向上分别进行插值处理,即可得到任意时刻卫星的位置。
拉格朗日多项式插值是代数多项式插值中较为简单，格式整齐、对称和规范，便于程序设计的一种形式。
其中公式（1）的：11()n j i j ijj i t t l t tt （2）称为拉格朗日插值基函数。
其实，拉格朗日插值函数就是拉格朗日插值基函数的一个线性组合。
两种IGS精密星历多项式插值方法的比较分析
７ （－ｏ － ，∈［０￡ ＝ ｔｔ） １￡ ￡， ＋△￡ ｏ ］
Ｕ
（） １
则 卫星坐 标 ，，ｚ分 量 可 用 如 下切 比雪 犬 多项 ｌ ， 式表 示 ： （ ） ｃ Ｔ（ ＝ ）
ｙ ＝ｉ （）
Ｇ ；
插值 、 拉格 朗 日多项式插值、 牛顿多项式插值 、 三角函 数 多项 式 插值 等 。 目前在 Ｇ Ｓ数 据 处理 中 ， Ｐ 应用 比较
．
在式 （ ） ， ２ 中 ｎ为切 比雪夫 多项式 的系数 ， ．Ｃ Ｃ ，
分别为 坐标分量 、 ｙ坐标分量 、 ｚ坐标 分量切 比雪 犬多
项式的系数。切比雪夫多项式 ７ 用以下递推公式确定 ：
， ｒ ＝１ ｒ（ ） ０ １
２ 切 比雪夫 多项式 插值 与拉格 朗 日多项式 插值
２ １ 切 比雪夫 多项式 插值 ．
（ ） ＝
Ｉ （） ３
（Ｉ＝２ 一 ） 一 ） Ｉ ｌ ， ≥２ ７ ） ｒ １ 一 ２ ； ≤１ ｎ （ （ Ｊ
切 比雪 夫多项 式拟 合就 是根据 给 定的数 据 拟合 出
一
根据精密星历 ， 设卫星的 坐标为观测值 ， 则误
（－ ７）
原理 利用 Ｖ Ｃ编 制了程序来实 现两种插值方法 。文 中选
， ｒ（ ） ｎ
（Ｉ ７） ： （ ｔ） 『
（ Ｔ（ １ ） １丁 ） ｒ （－ ７） ２
取 了 ＩＳ 布 的 ２０ Ｇ 发 ０９年 ３月 ２ 日的精密星历数据 作为 ９ 算例来 研 究不 同 阶数 多项 式 下 两 种插 值 方 法 的精 度 。 由于精密星历给 出的卫 星位置 的时 间间 隔为 １ ｎ 因 ５ｍｉ， 此本 文计算 中 以 ３ ｉ为时 间 间隔选 取 内插点 内插 计 ０ｍｎ 算 得到各个 内插 时间段 中间时刻 的卫 星位置 ， 以便将 两
数值分析常用的插值方法
数值分析常用的插值方法数值分析中常用的插值方法有线性插值、拉格朗日插值、分段线性插值、Newton插值、Hermite插值、样条插值等。
下面将对这些插值方法进行详细介绍。
一、线性插值（linear interpolation）线性插值是最简单的插值方法之一、假设已知函数在两个点上的函数值，通过这两个点之间的直线来估计中间点的函数值。
线性插值公式为：f(x)=f(x0)+(x-x0)*(f(x1)-f(x0))/(x1-x0)其中，f(x)表示要求的插值点的函数值，f(x0)和f(x1)是已知的两个点上的函数值，x0和x1是已知的两个点的横坐标。
二、拉格朗日插值（Lagrange interpolation）拉格朗日插值是一种基于多项式的插值方法。
它通过多个已知点的函数值构造一个多项式，并利用这个多项式来估计其他点的函数值。
拉格朗日插值多项式的一般形式为：f(x) = Σ[f(xi) * Li(x)] (i=0,1,2,...,n)其中，f(x)表示要求的插值点的函数值，f(xi)是已知的多个点的函数值，Li(x)是拉格朗日基函数。
拉格朗日基函数的表达式为：Li(x) = Π[(x-xj)/(xi-xj)] (i≠j，i,j=0,1,2,...,n)三、分段线性插值（piecewise linear interpolation）分段线性插值是一种逐段线性近似函数的方法。
通过将整个插值区间分成多个小段，在每个小段上使用线性插值来估计函数的值。
分段线性插值的过程分为两步：首先确定要插值的点所在的小段，在小段上进行线性插值来估计函数值。
四、Newton插值（Newton interpolation）Newton插值也是一种基于多项式的插值方法。
利用差商的概念来构造插值多项式。
Newton插值多项式的一般形式为：f(x)=f(x0)+(x-x0)*f[x0,x1]+(x-x0)*(x-x1)*f[x0,x1,x2]+...其中，f(x)表示要求的插值点的函数值，f(x0)是已知的一个点的函数值，f[xi,xi+1,...,xi+k]是k阶差商。
基于拉格朗日插值方法的GPS IGS精密星历插值分析
密， 因此 ， 要想 利用 某 一 时刻 的卫 星 位 置 ， ９部队 ， 黑龙江 哈 尔滨 １０ ３ ；２ 信息工程大学 测绘 学院， ５０９ ． 河南 郑 州 ４ ０ ５ ） ５０ ２ 摘 要 ： Ｓ数据处理 中 ， ＧＰ 为使用 Ｉ ＧＳ精密星历获取 ＧＰ Ｓ卫 星精确 的、 更高 采样 率的轨道位置 ， 必须要对精密星历进
第２ ０卷第 ５ 期
２１ 年 １ ０１ Ｏ月
测
绘
工
程
Ｖ０． Ｏ № ． １２ ５
０ｃ ． ２ １ ｔ，０ １
ＥＮＧＩ ＥＲＩ ＮＥ ＮＧ ＵＲＶＥＹＩ ＯＦ Ｓ ＮＧ ＡＮＤ ＡＰ Ｉ Ｍ Ｐ ＮＧ
基 于拉 格 朗 日插 值 方 法 的 ＧＰ ＧＳ精 密 星 历 插 值 分 析 ＳＩ
ｚｏ ５０ ２ ｈ ｕ ４ ０ ５ ，Ｃｈ ｎ ） ｉａ
Ａｂ ｔａ ｔ Ｉ ｓｒ ｃ ：ｎ ＧＰＳ ｄ ｔ ｒ ｃ ｓｉ ｇ ａ ａ ｐ ｏ ｅ ｓｎ ，ｗｅ ｍｕ ｔｉｔ ｒ ｏ ａｅ ｔ ｅ ｐ ｅ ｉｅ ｅ ｈ ｍｅ ｉ ｒ ｖｄ ｄ ｂ ＧＳ ｆ ｒｇ ｔｉｇ ｓ ｎ ｅ ｐ ｌｔ ｈ ｒｃｓ ｐ ｅ ｒｓ ｐ ｏ ｉｅ ｙ Ｉ ｏ ｅｔｎ ｐ ｅ ｉｅａ ｄｍ ｏ ｅ ｓ ｍｐ ｉｇ ｐ ｓｔ ｎ ｏ Ｓ ｓ ｔ ｌ ｔｓ Ｔｈ ａ ｅ ｉ ｌ ｎ ｒ ｄ ｃ ｓ ｔ ｅ ｂ ｓｃｐ ｉｃｐ ｅ ｏ ｒｃｓ ｎ ｒ ａ ｌ ｏ ｉ ｏ ｆＧＰ ａｅｌ ｅ ． ｎ ｉ ｉ ｅ ｐ ｐ ｒ ｓｍｐ ｙ ｉ ｔｏ ｕ ｅ ｈ ａ ｉ ｒｎ ｉｌ ｆ Ｌａ ｒ ｎ ｅｉ ｔｒ ｏａ ｉｎ ｎ ｌ ｓ ｓｔ ｅｉｔ ｒ ｏ ａｉｎ ｒ ｓ ｌｓｉ ｅａｌ ｎ ｌｏｇ ｔ ｏ ｅｈ ｌｆ ｌ ｏ ｃｕ ｉｎ ． ｇ ａ ｇ ｎ ｅ ｐ ｌｔ ；ａ ａｙ ｅ ｈ ｎ ｅ ｐ ｌｔ ｅ ｕ ｔ ｎｄ ｔ ｉ，ａ ｄａｓ ｅ ｓｓ ｍ ｅｐ ｕ ｎ ｌ ｓｏ ｓ ｏ ｏ ｃ
精密星历轨道内插方法比较
15/20
三、算例及结果比较
由图4可以看出： a. 插值曲线仍为U型，但是比 拉格朗日插值曲线要平滑。 b. 两端和中间插值的相差大部
分在0.1以内，插值过程那个
不会出现大的震荡和跳跃。 c. 9阶、10阶和11阶的切比雪
夫计算卫星位置不稳定，达到
12阶后，插值点位于节点中央 时，随着阶数增加，切比雪夫 图4. 切比雪夫多项式法内插卫星位置 多项式内插位置比拉格朗日插
其中 t 0 为起始时刻，t 为拟合时间区间的长度。首先将变量 t [t0 , t0 t ]
2 (t t0 ) 1, t [t0 , t0 t ] t
（3）
7/20
二、插值方法
则卫星坐标X、Y、Z的切比雪夫多项式为：
f x (t ) axiTi ( )
(5)
8/20
二、插值方法
1 根据最小二乘原理，求出 M ( BT B) 1 ( BT L)，即可求出在 t [t0 , t0 t ]
的 n 阶切比雪夫多项式系数。然后根据这些系数计算观测时刻卫星的坐标。
9/20
二、插值方法
2.3 Neville插值法
2.2 切比雪夫多项式拟合法
以切比雪夫多项式为基函数拟合构成的卫星位置函数为：
f x (t ) axiTi (t )
i 0
n
f y (t ) a yiTi (t )
i 0
n
f z (t ) aziTi (t )
i 0
n
（2）
假设需要在时间间隔 [t0 , t0 t ] ，计算n阶切比雪夫多项式系数。 变换为 [1,1] ，则
Neville插值是一种线性逐次插值，其基本思想是通过低一次多项式的组合 来获得高一次插值多项式。
基于精密星历的切比雪夫多项式卫星轨道坐标拟合研究
（ Ｔ（ ） 丁）
（ １ Ｔ（ ） ｒ） 。 。 （ｚ ） （ ） （ ） … （ｚ … ） （ ） … （ ｒ） （ ） Ｙ（３ ｎｒ ）
Ｃ０ ｘ
２ 切 比雪夫 多项式拟合原理和实现
２ １ 切 比雪夫 多项 式拟合 原理 ．
１
误 差方程 的矩 阵展开式 为
中图 分 类 号 ：２ ８ Ｐ ２
文 献标 识 码 ： Ａ
基于精密星历的切 比雪夫 多项 式卫星轨道坐标拟合研究
刘 刚
（ 武汉大学 测绘学 院， 湖北 武汉 摘 ４０７ ） ３０ ９ 要 ： 星位置的计算是 Ｇ Ｓ卫 星定位 中的一个关键环 节。本文介绍 了基于精密星历 的切比雪夫 多项 式轨 道拟合 卫 Ｐ
（ ） Ｔ（ 】 １『 ）
（ ） …
（ ．） ｒ
３ 计算实例及精度分析
作者选取 了 ＩＳ发布 的 ２０ Ｇ ０ ９年 ３月 ２ ９日的精 密 星历 数据文件 作 为算 例 ， 了分 析 切 比雪夫 多项 式插 为 值 的特 点 ， 别 按 照 ３ ｎ ４ ｉ，０ｍｎ 间间 隔 分 ０ｍｉ，５ｍｎ ６ ｉ时 选取节 点计算 了不 同阶次 的轨道 拟合 结 果 。当采用 ｎ 阶切 比雪夫 多项 式进 行轨 道拟合 时 ， 至少 需要 ｎ １ ＋ 个 节点 的坐标数 据 ， 了检验轨道 拟合 的质量 ， 取 ｎ １ 为 选 ＋ 个节点 中间 已知坐标 的 ，个点 作 为插 值 点 ， ｚ 计算 出这
【 （ ） Ｔ Ｂ （ 一 ２丁 ； ≤１ ｎ ２ 丁 ＝２ Ｔ— ） 一（ ） ｌ ｌ ｌ ， Ｉ ＞
（） ３
根 据精密 星历 ， 卫 星 的 坐 标 为观 测 值 ， 设 则误
差方程 为 ：
基于拉格朗日多项式的求解方法研究
基于拉格朗日多项式的求解方法研究拉格朗日多项式是一种用于插值函数的方法，该方法可以将给定的一组离散数据点在数据范围内插值得到一个连续的函数表达式。
在后面的数值分析中，拉格朗日多项式被广泛应用于数值积分和微分、拟合、求解差分方程等方面。
在实际应用中，拉格朗日多项式求解方法的应用价值是不可低估的。
首先，拉格朗日多项式是一种相对比较简单的插值方法，对于离散数据点进行插值的时候，求解过程十分简单，且适用范围广，可以适合各种数据类型、数据分布情况。
其次，拉格朗日多项式法具备高精度的特点，能够通过增加节点的数量进行优化，从而得到更高精度的计算结果。
最后，拉格朗日多项式是一种优美的方法，不仅简单易懂，而且可以直观地理解插值的过程。
下面我们将从各种角度详细阐述基于拉格朗日多项式的求解方法研究。
一、拉格朗日多项式的定义和性质拉格朗日多项式是一种多项式函数，常见形式是：$L(x) = \sum_{i=0}^{n}y_i\prod_{j=0,j\neq i}^{n}\frac{x-x_j}{x_i - x_j}$其中$x_0,x_1,\cdots,x_n$为离散数据点的横坐标，$y_0,y_1,\cdots,y_n$为对应的纵坐标。
这里$x$是一个自变量，是需要求解的新数据点。
将该自变量代入$L(x)$中，得到$L(x)$函数的返回值，即为求解结果。
拉格朗日多项式有两个重要的性质：1. 任意$n+1$个点都可用$n$次拉格朗日多项式进行插值2. 如果$n+1$个点均不同，$L(x)$函数在每个点$x_i(i=0,1,\cdots,n)$的值都是$1$，在其他点的值都是$0$。
二、基于拉格朗日多项式的插值方法不同的插值方法都有自己独特的步骤，基于拉格朗日多项式的插值方法步骤如下：1. 确定插值数据点集合$(x_0,y_0),(x_1,y_1),\cdots,(x_n,y_n)$2. 根据插值数据点集合，构建一个$x$的$n$次多项式$L(x)$3. 将$x$的值代入$L(x)$中，求解对应的$y$值下面我们将详细分析各个步骤的实现过程。
IGS精密星历文件的读取及内插方法研究
IGS精密星历文件的读取及内插方法研究柳迪;逢淑涛;董绪荣【期刊名称】《全球定位系统》【年(卷),期】2011(036)005【摘要】利用matlab软件编程实现了精密星历文件的读取及信息筛选。
在GPS 精密单点定位的过程中,需要比国际全球导航卫星系统服务（IGS）组织提供的15min更小的精密星历时间间隔的观测数据,介绍了拉格朗日和牛顿两种插值方法,对IGS精密星历进行了内插,并对插值的结果进行对比分析,得到了一些有益的结论。
%This article had completed the read and information calendar of precise ephemeris file by using Matlab.In the process of PPP,it needs to observe data of precise ephemeris whose time interval is less than 15 min,so two methods of interpolation are introduced in the article：Lagrange and Newton.Insert in the precise ephemeris,and comparative analysis on the results of interpolation,so that we can get some beneficial conclusions.【总页数】4页(P46-48,64)【作者】柳迪;逢淑涛;董绪荣【作者单位】装备指挥技术学院,北京101416;装备指挥技术学院,北京101416;装备指挥技术学院,北京101416【正文语种】中文【中图分类】P228.4【相关文献】1.基于IGS精密星历的卫星位置内插方法比较 [J], 马俊;唐诗华;黄鹰;王文杰2.IGS精密星历内插及外推方法比较分析 [J], 刘翔;时振伟3.基于滑动式算法的IGS精密星历内插与拟合精度 [J], 吴伟;任超;王文杰;黄征凯4.IGS精密星历内插方法的比较 [J], 张朋飞;谢劭峰;刘立龙5.IGS精密星历内插方法研究 [J], 高明超;徐泮林;谷彦斐因版权原因，仅展示原文概要，查看原文内容请购买。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
Ｋｅ ｙ ｗｏ ｒ ｄ ｓ ： Ｃ ｈ ｅ ｂ ｙ ｓ ｈ ｅ ｖ ；Ｌ ａ ｇ ｒ ａ ｎ ｇ ｅ ；ｉ ｎ ｔ ｅ ｒ ｐ ｏ ｌ ａ ｔ ｉ ｏ ｎ；ｆ ｉ ｔ ｉ ｔ ｎ ｇ
Ｏ 引 言
在 高精度 Ｇ Ｐ Ｓ精 密 测 量 中 ， 基 线 处 理 采 用 了 经 过后
ｎ阶 拉 格 朗 日插 值 多 项 式 的 代 数 表 达 式 点 （ ｉ ＝０ ， １ ， …，
合 法等 。本文 采用 了拉格 朗 日方 法和 切 比雪 夫方 法对 １ ５
轨道状态时， 需要选取一个 区间［ ｔ ， ｔ ］ ， 使得 ｔ ∈［ ｔ 。 ， ｔ ］ 。
Ｉ Ｇ Ｓ 提供 的精 密星历 文件 中 ， 是以 １ ５ ｍ ｉ ｎ 为 间隔提供 观测 卫星 的三维 坐标 。 以 ３ ０ ｍ ｉ ｎ 为间隔， 读 取数 据 ， 将 两 数据 点之 间 １ ５ ｍ ｉ ｎ时 刻 的 内插结 果 和 Ｓ Ｐ ３原 始 坐标 做 比较 ， 得 到 内插精 度 。想 要获 得任 意时刻 ｎ 阶拉 格朗 日内差值 ，
必须有 ｎ＋１ 个 时刻 的卫 星位置 。其 中 ，ｉ ｎ的精 密 星 历 加 密 并 与 Ｉ Ｇ Ｓ提 供 原 始 数 据 进 行 对 比。 详 细分析 了精 密 星历加 密 的 两种 方 法在 不 同 阶数 时 的插
值精 度 。
求 值 内插 函数为 ：
ｓ ｉ ｏ ｎ ｓ ａ ｔ ｄ ｉ ｆ ｅ ｒ ｅ ｎ ｔ ｏ ｒ ｄ ｅ ｒ ｓ ｆ ｏ ｒ ｂ ｏ ｔ ｈ ｍｅ ｔ ｈ ｏ ｄ ｓ ａ ｃ ｃ ｏ ｒ ｄ ｉ ｎ ｇ ｔ ｏ ｗｈ ｉ ｃ ｈ ｗ ｅ ｃ ａ ｎ ｃ ｈ ｏ ｏ ｓ ｅ ｔ ｈ ｅ ｍｏ ｓ ｔ ａ ｐ ｐ ｒ ｏ ｐ ｉａ ｒ ｔ ｅ ｏ ｒ ｄ ｅ ｒ ．
ｍｅ ｔ ｈ ｏ ｄ ａ ｎ ｄ ｔ ｈ ｅ ｎ ｒ ｅ ｓ ｅｃ ｐ ｔ ｉ ｖ ｅ ｌ ｙ ｃ ｏ ｍｐ ａ ｒ ｅ ｓ ｔ ｈ ｅ ｏ ｕ ｔ ｃ ｏ ｍｅ ｓ ｗｉ ｔ ｈ ｏ ｉ ｒ ｇ ｉ ｎ ａ ｌ ｅ ｐ ｈ ｅ ｍｅ ｉ ｒ ｓ ｇ ｅ ｔ ｔ ｉ ｎ ｇ ｐ ｒ ｅ ｃ ｉ ｓ ｉ ｏ ｎ ｏ ｆ ｔ ｈ ｉ ｓ ｔ ｗｏ ａ ｌ ｇ ｏ ｉ ｒ ｔ ｈ ｍｓ ．Ｗ ｅ ａ ｎ ｌｙ ａ ｚ ｅ ｐ ｒ ｅ ｅ ｉ —
ｎ ） 称为插值节点 ， 包含插值节点 的区间［ ａ ， ｂ ］ 称 为插值
区 间。
处理的 Ｉ Ｇ Ｓ 精密星历 ， 它在 Ｇ Ｐ Ｓ 精密定位 、 精密单点定位
（ Ｐ Ｐ Ｐ ） 和低 轨 卫 星精 密 定 轨 中 发挥 了 至关 重 要 的作 用 。
一
Ｈ【 】 （ ０ ． Ⅶ ） （ １ ）
将 拉 格 朗 日多项 式 应 用 到求 解任 意 时 刻 ｔ 某 卫 星 的
般接 收机 获取 的广 播 星历 时 间间 隔 为 ２ ｈ ， Ｉ Ｇ Ｓ网 站上
公 布 的精密卫 星 星历 的时间 间隔 为 １ ５ ｍ ｉ ｎ ， 而 实 际情 况下 的单 点定 位 中 ， 用户所 设 置 的 接收 信 号 间 隔一 般 为 ３ ０ Ｓ ， １ ５ Ｓ 或 者更 小 间隔 ， 因此需 要运 用一 定 的算 法模 型对 卫 星 星历进行 加 密 。方 法 包 括 拉格 朗 日内插 法 ， 切 比雪 夫 拟
孙腾 科
（ 河海大学 ， 江苏 南京 ２ １ ０ ０ ９ ８ ）
摘 要 ： 采用拉 格朗 日方法和切比雪夫方法对间隔为 １ ５ ｍ ｉ ｎ的精 密星历进行 内插或拟合 ， 将 结果和原始星历进行
比较 ， 得 到 内插 拟合 的精度 。详细分析 了在不 同阶数时的两种方法各 自的精度 ， 比较 了它们的优 劣。
Ｂａ ｓ ｅ ｄ ｏ ｎ Ｌａ ｇ ｒ ａ ｎ ｇ ｅ ａ ｎ ｄ Ｃｈｅ ｂ ｙ ｓ ｈ ｅ ｖ
Ｓ ＵＮ Ｔ ｅ ｎ ｇ—ｋ ｅ
（ Ｈ ｏ ｈ ａ ｉ Ｕ ｎ ｉ ｖ ｅ ｒ ｓ ｉ ｔ ｙ ， Ｎａ ｎ ｊ ｉ ｎ ｇ ２ １ ０ ０ ９ ８ ，Ｃ ｈ ｉ ｎ ａ ）
关键词 ： 拉 格 朗 日方 法 ； 切 比雪夫方法 ； 内插 ； 拟合
中图分 类号 ： Ｐ ２ ２ ８
文献标识码 ： Ｂ
文章编号 ： １ ６ ７ ２— ５ ８ ６ ７ （ ２ ０ １ ４ ） ０ ２— ０ ０ ３ ３— ０ ５
Ｔｈ ｅ Ｒｅ ｓ ｅ ａ ｒ ｃ ｈ ｏ ｆ Ｇ ＰＳ Ｐｒ ｅ ｃ ｉ ｓ ｅ Ｉ ｎ ｔ ｅ ｒ ｐ ｏ ｌ ａ ｔ ｉ ｏ ｎ Ｍｅ ｔ ｈ ｏ ｄ ｓ
＝ ＝
】 ］ （ ２ ） ． ｉ ｝ 【 ］
第３ ７卷 第 ２期
２ ０１ ４年 ２ 月
测绘 与 空 间地 理 信 息
Ｇ ＥＯ ＭＡＴ Ｉ Ｃ Ｓ＆ Ｓ Ｐ ＡＴ Ｉ Ａ Ｌ ｌ ＮＦ ｏ ＲＭＡ Ｔ ｌ ｏＮ Ｔ ＥＣ ＨＮｏ Ｌ ＯＧ Ｙ
Ｖｏ １ ． ３ ７， Ｎｏ ． ２ Ｆ ｅ ｂ ． ，２ ０ １ ４
基 于拉 格 朗 日与切 比雪夫 方 法 的精 密 星历 插 值 研 究