matlab拉格朗日插值代码

以下是MATLAB中的拉格朗日插值代码：

% x 和 y 分别为已知数据点的横纵坐标
% xi 为待求点的横坐标
function yi = lagrange_interp(x, y, xi)

n = length(x);
yi = 0;

for i = 1:n
    % 计算拉格朗日基函数
    L = 1;
    for j = 1:n
        if j ~= i
            L = L * (xi - x(j)) / (x(i) - x(j));
        end
    end
    
    % 利用拉格朗日基函数计算插值结果
    yi = yi + y(i) * L;
end

end

使用方法：

将上述代码保存为 `lagrange_interp.m` 文件，放在MATLAB工作路径下。

在命令窗口中输入以下代码即可进行插值计算：

% 定义已知点的横纵坐标
x = [1, 2, 3];
y = [4, 5, 6];

% 定义待求点的横坐标
xi = 2.5;

% 进行拉格朗日插值计算
yi = lagrange_interp(x, y, xi);

相关问题

matlab 拉格朗日插值代码

拉格朗日插值是一种常用的数值分析方法，其主要思想是通过已知的数据点来求解一个多项式，从而实现在这些数据点之间的近似值。在 MATLAB 中，可以使用 polyfit 和 polyval 函数来实现拉格朗日插值。

具体操作步骤如下：

1. 假设有 n 个数据点 (x1, y1), (x2, y2), ... , (xn, yn)，其中 x1, x2, ... , xn 是已知的数据点的横坐标，y1, y2, ... , yn 是对应的纵坐标。

2. 定义一个函数 f(x)，它可以通过这些已知数据点进行插值。

3. 构造一个 n-1 阶多项式 P(x)，使得 P(xi) = yi (i=1,2,...,n)。

4.＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊

matlab拉格朗日插值法代码

以下是使用拉格朗日插值法进行插值的 MATLAB 代码示例：

function [y_interp] = lagrange_interp(x, y, x_interp)
% Lagrange Interpolation
%
% Inputs:
% x: vector of known x values
% y: vector of known y values
% x_interp: vector of x values to be interpolated
%
% Output:
% y_interp: vector of interpolated y values at x_interp

n = length(x);
m = length(x_interp);
y_interp = zeros(1, m);

for i = 1:m
    % Compute the value of the interpolated point using Lagrange formula
    L = ones(n, 1);
    for j = 1:n
        for k = 1:n
            if j ~= k
                L(j) = L(j) * (x_interp(i) - x(k)) / (x(j) - x(k));
            end
        end
    end
    y_interp(i) = sum(y .* L);
end
end

其中，`x` 和 `y` 分别是已知的数据点的 x 和 y 坐标，`x_interp` 是要进行插值的 x 坐标，返回 `y_interp` 是对应的插值结果。