牛顿插值法计算 y = cos(x) 在 x = 0.05 时的值

要使用牛顿插值法计算 y = cos(x) 在 x = 0.05 时的值，我们需要选择一组已知数据点来构建插值多项式。

在这个例子中，我们可以选择一组数据点来自函数 y = cos(x) 的离散值。让我们选择以下数据点：

x = [0, 0.1, 0.2, 0.3, 0.4] y = [1, 0.995, 0.980, 0.955, 0.921]

我们将使用这些已知数据点来进行牛顿插值，并计算 x = 0.05 时的插值值。

以下是使用 MATLAB 编写的示例程序：

function y_interp = newton_interpolation(x, y, x_interp)
    n = length(x);
    F = zeros(n, n);  % 差商表

    % 计算差商表
    F(:, 1) = y;
    for j = 2:n
        for i = j:n
            F(i, j) = (F(i, j-1) - F(i-1, j-1)) / (x(i) - x(i-j+1));
        end
    end

    % 构建插值多项式并计算插值值
    y_interp = F(1, 1);
    for j = 2:n
        term = F(j, j);
        for i = 1:j-1
            term = term * (x_interp - x(i));
        end
        y_interp = y_interp + term;
    end
end

% 已知数据点
x = [0, 0.1, 0.2, 0.3, 0.4];
y = [1, 0.995, 0.980, 0.955, 0.921];

% 需要插值的点的x坐标
x_interp = 0.05;
y_interp = newton_interpolation(x, y, x_interp);  % 使用牛顿插值法进行插值

% 显示插值结果
disp(y_interp);

在上述示例中，我们选择了一组已知数据点，然后使用牛顿插值法计算 x = 0.05 时的插值值。运行上述代码将会输出插值结果。

希望这个示例对您有所帮助。如果您有其他问题，请随时提问。