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