lagranline.asv

数值分析常用几种方法

function lagranline(x,y,xi)
%x为向量，全部的插值节点；
%y为向量，插值节点处的函数值；
%xi为标量，自变量x;
%yi为xi处的函数估计值
n=length(x);m=lenth(y);
if n~m
    error('the length of x and y must be equal');
    return;
end
for k=1:n-1
    if x(k)<=xi&xi<=x(k+1)
        yi=(xi-x(k+1))/(x(k)-x(k+1))*y(k)+(xi-x(k))/(x(k+1)-x(k))*y(k+1);
        return ;
    end
end
a=-5;b=5;n=8;h=(b-a)/n;
x=a: h: b; y =1./(1+x.^2);
x1=a:0.01:b; y1=1./(1+x1.^2);
m=length(x1); z=zeros(1,m);
for i=1:m
z(i)=lagranline(x,y,x1(i));
end
plot(x,y,’o’,x1,y1.’k:’,x1,z,’k-’);