chebshev GPS精密星历卫星轨道插值运算（采用切比雪夫拟合）matlab程序。直接运行chebshev2.mat文件即可

GPS精密星历卫星轨道插值运算（采用切比雪夫拟合）matlab程序。直接运行chebshev2.mat文件即可，结果显示在doucument2.mat文件中

function [ddrhour1hour2hour3] = chebshev2
% n介切比雪夫多项式，由拟合阶段的长度c决定，
clc;
clear;
n=34;           %拟合阶数n=10+（终止时间-初始时间）
fprintf（‘\nStart processing:\n‘）;

[GCminGPS_TimeS_PRNdata_SP3Xdata_SP3Ydata_SP3ZSP3clk]=refrinex;
SP3time=GCmin;        SP3X=data_SP3X;       SP3Y=data_SP3Y;       SP3Z=data_SP3Z; 

t0=SP3time;
x0=SP3X;        %x0为以分钟为单位的时间序列，y0为坐标
y0=SP3Y;
z0=SP3Z;

t=min（t0）:15:max（t0）;       %按15分钟进行插值计算，时间间隔可调，单位分钟
nn=length（t0）;   %历元数
m=length（t）;
c=t0（end）-t0（1）;   %c:区间拟合长度，x0（1）起始历元时刻

r0=2*（t0-t0（1））/c-1;         %变量变为（-1，1）区间
T=zeros（nn）;
for i=1:nn
	T（i1）=1;
	T（i2）=r0（i）;
	for j=3:n
		T（ij）=2*r0（i）*T（ij-1）-T（ij-2）;
	end
end

C1=inv（T‘*T）*T‘*x0;       %最小二乘求切比雪夫多项式系数
C2=inv（T‘*T）*T‘*y0; 
C3=inv（T‘*T）*T‘*z0; 
%*————————求未知历元拟合坐标————

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2013-04-12 14:12  切比雪夫轨道拟合\
     文件        1185  2013-04-12 14:08  切比雪夫轨道拟合\chebshev2.m
     文件       44601  2013-04-12 14:06  切比雪夫轨道拟合\data_doucment.mat
     文件       67874  2013-04-12 14:06  切比雪夫轨道拟合\doucment2.mat
     文件        5260  2013-04-09 16:32  切比雪夫轨道拟合\Read_SP3.m
     文件        2282  2013-04-12 10:06  切比雪夫轨道拟合\refrinex.m
     文件         222  2013-04-11 08:48  切比雪夫轨道拟合\TimetoJD.m
     目录           0  2013-04-12 14:12  切比雪夫轨道拟合\精密星历文件15分钟\
     目录           0  2013-04-12 14:11  切比雪夫轨道拟合\精密星历文件15分钟\igs16600.sp3\
     文件      245454  2013-04-11 09:00  切比雪夫轨道拟合\精密星历文件15分钟\igs16600.sp3\igs16600.sp3
     目录           0  2013-04-12 14:12  切比雪夫轨道拟合\精密星历文件30分钟\
     目录           0  2013-04-12 14:11  切比雪夫轨道拟合\精密星历文件30分钟\igs16600.sp3 - 副本\
     文件      124971  2013-04-11 09:00  切比雪夫轨道拟合\精密星历文件30分钟\igs16600.sp3 - 副本\igs16600.sp3