fx_lagrange_interpolating.m

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

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

%********************************
%subplot(1,2,1);
%plot(x,y);
%subplot(1,2,2);
%plot(xi,yi);