matlab切比雪夫拟合

matlab中没有切比雪夫拟合的现成算法,这里把我程序中的这部分抽出来,说一下。

1、首先是切比雪夫计算式

function [ res ] = ChebyShev(num,i)

if i==0 
    res=1;
elseif i==1
    res=num;
else
    res=2*num*ChebyShev(num,i-1)-ChebyShev(num,i-2);
end
   
end

2、计算拟合系数

function [ coeff ] = ChebyCoeff( xfile,timefile )

x=load(xfile);
time=load(timefile);
B=zeros(length(time),11); %B
L=x;%L
sec_start=time(1);
sec_end=time(end);
for i=1:length(time)
     chang_sec=(2*time(i)-(sec_start+sec_end))/(sec_end-sec_start);
    for j=1:11
        B(i,j)=ChebyShev(chang_sec,j);
    end
end
coeff=(B'*B)\(B'*L);
end

3、根据系数计算拟合值

function [ res ] = ChebyValue( timefile,coeff,num )

res=0;
time=load(timefile);
change_time=(2*num-(time(1)+time(end)))/(time(end)-time(1));
for i=1:length(coeff)
    res=res+coeff(i)*ChebyShev(change_time,i);
end

end

 

posted @ 2016-11-29 21:28  啊王会  阅读(3611)  评论(0编辑  收藏  举报