matlab向前插值代码

根据您提供的代码，这是一个用MATLAB编写的牛顿插值的向前插值代码：

function yi = forwardInterpolation(x, y, xi)
    m = length(x);
    n = length(y);
    if m ~= n
        error('向量x与y的长度必须一致');
    end
    A = zeros(n);
    A(:, 1) = y;
    for j = 2:n
        for i = 1:(n-j+1)
            A(i, j) = (A(i+1, j-1) - A(i, j-1)) / (x(i+j-1) - x(i));
        end
    end
    N(1) = A(1, 1);
    for j = 2:n
        T = 1;
        for i = 1:j-1
            T = T * (xi - x(i));
        end
        yi = N(j-1) + A(1, j) * T;
    end
end