一. Lagrange插值

最新推荐文章于 2023-10-17 21:39:46 发布

XIAOHE1004

最新推荐文章于 2023-10-17 21:39:46 发布

阅读量224

点赞数

文章标签： matlab 算法开发语言

原文链接：https://blog.csdn.net/forest_LL/article/details/124359478?ops_request_misc=%257B%2522request%255Fid%2522%253A%2522167040752716782425662201%2522%252C%2522scm%2522%253A%252220140713.130102334..%2522%257D&request_id=167040752716782425662201&biz_id=0&utm_medi

版权

一. Lagrange插值

1.1 数学解释
一. Lagrange插值
1.1 数学解释
对给定的n个插值点(x_1,y_1),(x_2,y_2),\cdots,(x_n,y_n)，可以构造n-1次Lagrange插值多项式。对插值区间内的任意x，对应的y值可由如下公式计算：

在这里插入图片描述

1.2 MATLAB实现
MATLAB中没有lagrange函数，需要提前自己构造。

构造函数代码如下：

function y=lagrange(x0,y0,x)
ii=1:length(x0);
y=zeros(size(x));
for i=ii
    ij=find(ii~=i);
    y1=1;
    for j=1:length(ij),y1=y1.*(x-x0(ij(j)));
    end
    y=y+y1*y0(i)/prod(x0(i)-x0(ij));
end
end

例题1
在这里插入图片描述

clc;clear;
x=[0.4:0.1:0.8];
y=[-0.916291,-0.693147,-0.510826,-0.356675,-0.223144];
lagrange(x,y,0.54)
 
%函数部分
function y=lagrange(x0,y0,x)
ii=1:length(x0);
y=zeros(size(x));
for i=ii
    ij=find(ii~=i);
    y1=1;
    for j=1:length(ij),y1=y1.*(x-x0(ij(j)));
    end
    y=y+y1*y0(i)/prod(x0(i)-x0(ij));
end
end

运行结果如下：
在这里插入图片描述