拉格朗日插值

function s=larg1(x,y,xi)
m=length(x);          %%求出插值节点向量长度
n=length(y);
if m ~= n , error('向量x与y的长度必须一致');
    %%这里肯定一致，只是为了消除直接选取x,y数据的时候出错而设定的
end;
s=0;
for i=1 : n
    z=ones(1,length(xi));%%这里定义一个长度为计算节点个数的零向量，可以加快运行速度
    for j=1 : n
        if j ~= i
           if x(j)==0.5  %%用于增加误差的时调用
               y(j)=y(j)+0.05;
           end
            z=z.*(xi-x(j))./(x(i)-x(j));%%拉格朗日插值基函数
        end
    end
    s=s+z*y(i);
end

fprintf('\t\t\t插值点X\t\t插值点Y\n');

for i=1:length(x)
    fprintf('I(%2.0f)\t\t%5.3f\t\t%5.3f\n',i,x(i),y(i));
end
fprintf('\t\t\t取值点X\t\t取值点点Y\n');
for i=1:length(s)
    fprintf('I(%2.0f)\t\t%5.3f\t\t%5.3f\n',i,xi(i),s(i));
    
end
hold on

%以下是设置绘图区域的相关指令
plot(x,y,'r*');
ezplot('x./(1+x.^4)',[-5,5]);
plot(xi,s,'k-');
plot(xi,s,'go');

title('拉格朗日插值');
legend('插值点','原始函数','拉格朗日插值函数','计算点');
xlabel('x的取值');
ylabel('y的取值');
axis([-5,5,min([min(y) min(s)])-0.2,max([max(y) max(s)])+0.2]);
set(gca,'xtick',[-5:0.25:5]);
set(gca,'ytick',[-10:0.1:10]);
grid on
hold off

拉格朗日插值