chebeshev.m

来自「切比雪夫插值（利用切比雪夫点进行多项式插值）」· M 代码 · 共 20 行

close all;clear all;clc;
syms L x;
figure;
ezplot(sqrt(x),[0,1]);
set(findobj('Type','line'),'Color','r');
for n=3:2:7
    f=1;L=0;
    for i=1:n
        t(i)=cos((2*i-1)*pi/4/n);
        f=f*(x-t(i));
    end
    for i=1:n
        y=f/(x-t(i));
        yi=subs(y,x,t(i));
        L=L+sqrt(t(i))*y/yi;
    end
    %hold on;
    figure;ezplot(L,[0,1]);
end
set(findobj('Type','line'),'LineWidth',3);