fx_lagrange_interpolating_chebyshev.m

来自「对于切比雪夫点 利用 Matlab 进行lagrange插值」· M 代码 · 共 14 行

N=input('Enter the number of N:');
pointx=ones(1,N+1);
pointy=ones(1,N+1);
%*******************
for i=1:N+1
    pointx(i)=cos(((2.*i-1)./(2.*N+2)).*pi);
    pointy(i)=(1+25.*(pointx(i).^2)).^(-1);
end
%**********************
x=linspace(-1,1,1000);
y=(1+25.*(x.^2)).^(-1);
xi=linspace(-1,1,1000);
yi=lagrange_interpolating(xi,pointx,pointy);
plot(x,y,xi,yi,'--');