c0d0.m

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

function [c0,d0]=c0d0(x,a,b,h,lb,g,l,m)
s=l-m+2;
n=l+m+2;
y(1:2:2*m)=-g(s+1:2:n);
y(2:2:2*m)=g(s+2:2:n);
if mod(s,2)==0
    y=-y;
end
p1=mfunc2(x(s+1:n),a(1:s-1),x(1:s-2));
p2=mfunc2(x(s+1:n),b(1:s-1),x(1:s-2));
for i=2:(2*m+1)
    r(i-1)=prod(x(s+i-1)-x(1:s-1));
end
deno=h(s+1:n)-p1;
c0(1)=1;
d0(1)=0;
c0(2:(2*m+1))=a(s)*r./deno;
d0(2:2*m+1)=(b(s)*r-c0(2:2*m+1).*(y-p2))./deno;