function?[ddrhour1hour2hour3]?=?chebshev2
%?n介切比雪夫多項(xiàng)式，由擬合階段的長(zhǎng)度c決定，
clc;
clear;
n=34;???????????%擬合階數(shù)n=10+（終止時(shí)間-初始時(shí)間）
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為以分鐘為單位的時(shí)間序列，y0為坐標(biāo)
y0=SP3Y;
z0=SP3Z;

t=min（t0）:15:max（t0）;???????%按15分鐘進(jìn)行插值計(jì)算，時(shí)間間隔可調(diào)，單位分鐘
nn=length（t0）;???%歷元數(shù)
m=length（t）;
c=t0（end）-t0（1）;???%c:區(qū)間擬合長(zhǎng)度，x0（1）起始?xì)v元時(shí)刻

r0=2*（t0-t0（1））/c-1;?????????%變量變?yōu)椋?1，1）區(qū)間
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;???????%最小二乘求切比雪夫多項(xiàng)式系數(shù)
C2=inv（T‘*T）*T‘*y0;?
C3=inv（T‘*T）*T‘*z0;?
%*————————求未知?dú)v元擬合坐標(biāo)————

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----
?????目錄???????????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