多项式插值——拉格朗日多项式

拉格朗日多项式

如果将线性插值多项式看成由直线相连的两个点的加权平均，那么它的公式可写为：

其中L是加权系数：

同理，二次拉格朗日插值多项式可以写为：

推广到更高次的多项式：

其中：

函数实现：

function yint = Lagrange(x, y, xx)
n = length(x);
s = 0;
for i = 1:n
    product = y(i);
    for j = 1:n
        if i ~=j
            product = product*(xx-x(j)) / (x(i) - x(j));
        end
    end
    s = s+product;
end
yint = s;
end