cd1.m

用有理多项式做Chebyshev逼近的一种直接方法的matlab程序

function [c,d,c1,d1]=cd1(x,c0,d0,l,m)
n=2*m;
cc=zeros(n,n-1);
cc(1:n,1)=((x(l-m+3:l+m+2)-x(l-m+2:l+m+1))./(diff(c0)))';
dd(1:n,1)=cc(1:n,1)./(diff(c0))'.*(-diff(d0))';
if m>1 
    cc(2:n,2)=(x(l-m+4:l+m+2)-x(l-m+2:l+m))'./diff(cc(1:n,1))+c0(2:n)';
    dd(2:n,2)=cc(2:n,2)./diff(cc(1:n,1)).*(-diff(dd(1:n,1)))+d0(2:n)';
    if m>2
        for i=3:n-1
            cc(i:n,i)=(x(l-m+1+i+1:l+m+2)-x(l-m+2:l+m+2-i))'./diff(cc(i-1:n,i-1))+cc(i-1:n-1,i-2);%error
            dd(i:n,i)=cc(i:n,i)./diff(cc(i-1:n,i-1)).*(-diff(dd(i-1:n,i-1)))+dd(i-1:n-1,i-2); %error
        end
    end
end
c(1)=1;
c(2)=cc(1);
d(2)=dd(1);
if n>2
c(2:n)=diag(cc)';
d(2:n)=diag(dd)';
end
d(1)=0;
c1=cc(n,n-1);
d1=dd(n,n-1);

%c=1;d=1;c1=1;d1=1;