拉格朗日插值法的matlab实现

一、基本原理

比如有如下这些点

x1 x2 x3 x4
y1 y2 y3 y4

那么在拉个朗日原理中可以把过这些点的曲线表示为:

其g(x)y叫做一个插值基函数 (开关),当x=x1时,g1(x)=1,而当x=x2,x3,x4时,g1(x)都为0,其他的也是如此,这样当x=x1时,f(x)中就只保留了g1(x)y1项,其他项都为0,就保证了过点(x1,y1);其他点以此类推;

原理讲解得很简单,具体可以参考其他文献资料。

二、Matlab代码

matlab 复制代码
clear
clc

x = [466 741 950 1422 1634];
y = [7.04 4.28 3.40 2.54 2.13];

x1 = min(x):0.5:max(x);
y1 = lagrange(x,y,x1);
plot(x,y,'ko',x1,y1,'r-')

function fx = lagrange(x0,y0,x)
    m = length(x);
    n = length(x0);
    for i = 1:m
        z = x(i);
        s = 0;
        for k = 1:n
            p = 1;
            for j = 1:n
                if j~=k
                    p = p * (z-x0(j))/(x0(k)-x0(j));
                end
            end
            s = p*y0(k)+s;
        end
        fx(i) = s;
    end
end
相关推荐
简简单单做算法3 分钟前
基于LSTM深度学习网络的视频类型分类算法matlab仿真
深度学习·matlab·分类·lstm·视频类型分类
悦悦子a啊34 分钟前
Python之--集合
开发语言·python·编程
运维帮手大橙子43 分钟前
字符串缓冲区和正则表达式
java·开发语言
慢慢沉2 小时前
Lua(数据库访问)
开发语言·数据库·lua
GISer_Jing2 小时前
50道JavaScript基础面试题:从基础到进阶
开发语言·javascript·ecmascript
Python涛哥2 小时前
PHP框架之Laravel框架教程:1. laravel搭建
开发语言·php·laravel
一百天成为python专家3 小时前
数据可视化
开发语言·人工智能·python·机器学习·信息可视化·numpy
武子康6 小时前
Java-82 深入浅出 MySQL 内部架构:服务层、存储引擎与文件系统全覆盖
java·开发语言·数据库·学习·mysql·spring·微服务
惜.己6 小时前
pytest中使用skip跳过某个函数
开发语言·python·测试工具·pytest
姜暮儿6 小时前
C++ 性能优化
开发语言·c++