lagrange.m

来自「计算插值与曲线拟合的MATLAB程序」· M 代码 · 共 24 行

%程序名称	lagrange.m
%调用格式	y=lagrange(x0,y0,x)
%程序功能	求插值点 处的函数值 
%输入变量	x0，y0为已知的离散数据
%			x为插值节点
%输入变量	y为插值节点对应的函数值
%程序
function y=lagrange(x0,y0,x)
n=length(x0);m=length(x);
for i=1:m
z=x(i);
s=0.0;
for k=1:n
 p=1.0;
    for j=1:n
    if j~=k
    p=p*(z-x0(j))/(x0(k)-x0(j));
    end
    end
 s=p*y0(k)+s;
end
y(i)=s;
end