lagrange.m

来自「数值分析, 同济大学教材<<现代数值数学和计算>>数值分」· M 代码 · 共 12 行

function yh = lagrange(x, y, xh)
% xh - desired point array
n = length(x);
m = length(xh);
yh = zeros(1, m);
c1 = ones(n-1, 1);
c2 = ones(1, m);
for i=1 : n
  xp = x([1:i-1 i+1:n]);
  yh = yh + y(i) * prod((c1 * xh - xp' * c2) ./ (x(i) - xp' * c2));
end