Matlab 中的拉格朗日插值问题

Question

我在完成这个拉格朗日函数时遇到问题，我不确定哪里出错了。

我一直收到 'Index exceeds matrix dimensions' 错误。我认为这是因为我没有正确移动数组。我尝试了多种方法，但似乎继续回到这个错误。

function y = lagrange(X, Y, x)
n = length(X);
if n ~= length(Y)
 error('X and Y must have the same length.');
end
y = zeros(size(x)); 
for i = 1:n 
  L = ones(size(x)); 
    for j = [1:i-1 i+1:n] 
      if (i~=j)
       L(1:i-1, i+1:n) = L(1:i-1, i+1:n).*(x-X(j))/(X(i)-X(j));
      end
    end
  y = y+Y(i)*L(1:i-1);
end

Answer 1

我想下面的功能就是你需要的。

function y = lagrangian(X, Y, x)
if length(X) ~= length(Y); error('X and Y must have the same length.'); end
y = zeros(size(x));
for i = 1:length(X)
   L = ones(size(x));
   for j = [1:i-1 i+1:length(X)]
       L = L.*(x-X(j))/(X(i)-X(j));
   end
   y = y+Y(i)*L;
end

Matlab 中的拉格朗日插值问题

Issue with Lagrange interpolation in Matlab

matlab

interpolation

polynomial-math