function?[C1C2C3SPRNmcc]?=?chebshev2
%?n介切比雪夫多項式，由擬合階段的長度c決定，
clc;
clear;
n=12;???????????%擬合階數n=10+（終止時間-初始時間）
fprintf（‘\nStart?processing:\n‘）;

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

t0=SP3time;
x0=SP3X;????????%x0為以分鐘為單位的時間序列，y0為坐標
y0=SP3Y;
z0=SP3Z;

t=min（t0）:15:max（t0）;???????%按15分鐘進行插值計算，時間間隔可調，單位分鐘
nn=length（t0）;???%歷元數
m=length（t）;
cc=t0（end）-t0（1）;???%c:區間擬合長度，x0（1）起始歷元時刻

r0=2*（t0-t0（1））/cc-1;?????????%變量變為（-1，1）區間
T=zeros（nnn）;
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;?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1196??2017-12-18?09:18??切比雪夫軌道擬合\chebshev2.asv

?????文件????????896??2017-12-18?09:19??切比雪夫軌道擬合\chebshev2.m

?????文件??????87377??2017-12-16?22:23??切比雪夫軌道擬合\data_doucment.mat

?????文件??????70823??2017-12-16?20:56??切比雪夫軌道擬合\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

?????文件?????245454??2013-04-11?09:00??切比雪夫軌道擬合\精密星歷文件15分鐘\igs16600.sp3\igs16600.sp3

?????文件?????124971??2013-04-11?09:00??切比雪夫軌道擬合\精密星歷文件30分鐘\igs16600.sp3?-?副本\igs16600.sp3

?????目錄??????????0??2017-12-18?09:22??切比雪夫軌道擬合\精密星歷文件15分鐘\igs16600.sp3

?????目錄??????????0??2017-12-18?09:22??切比雪夫軌道擬合\精密星歷文件30分鐘\igs16600.sp3?-?副本

?????目錄??????????0??2017-12-18?09:22??切比雪夫軌道擬合\精密星歷文件15分鐘

?????目錄??????????0??2017-12-18?09:22??切比雪夫軌道擬合\精密星歷文件30分鐘

?????目錄??????????0??2017-12-18?09:22??切比雪夫軌道擬合

-----------?---------??----------?-----??----

???????????????538481????????????????????14