matlab拉格朗日插值拟合

Matlab中可以使用拉格朗日插值法进行数据拟合。拉格朗日插值的函数可以通过以下代码实现：

function y = lagrange(x0, y0, x)
    n = length(x0);
    m = length(x);
    for i = 1:m
        z = x(i);
        s = 0.0;
        for k = 1:n
            p = 1.0;
            for j = 1:n
                if j ~= k
                    p = p * (z - x0(j)) / (x0(k) - x0(j));
                end
            end
            s = s + p * y0(k);
        end
        y(i) = s;
    end
end

要进行拉格朗日插值拟合，需要提供已知数据点的横坐标和纵坐标，分别存储在数组`x0`和`y0`中。然后可以通过调用`lagrange`函数来得到拟合后的结果。例如，对于给定的数据点：

x0 = [3, 5, 7, 9, 11, 12, 13, 14, 15];
y0 = [1.2, 1.7, 2.0, 2.1, 2.0, 1.8, 1.2, 1.0, 1.6];

可以使用以下代码进行拟合：

x = 1:0.1:20;
z = lagrange(x0, y0, x);
plot(x, z, ':', x0, y0, 'ko');

这将会得到拉格朗日插值拟合后的曲线图像，其中`x`为插值点，`z`为插值结果。<span class="em">1</span><span class="em">2</span><span class="em">3</span>
#### 引用[.reference_title]
- *1* *2* *3* [Matlab实现插值与拟合](https://blog.csdn.net/qq_34374664/article/details/79167690)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v92^chatsearchT3_1"}}] [.reference_item style="max-width: 100%"]
[ .reference_list ]