切比雪夫轨道拟合MATLAB代码

利用切比雪夫进行轨道拟合详细代码，MATLAB代码，利用30s的精密星历来拟合5s的星历卫星坐标

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