拉格朗日插值法是基于基函数的插值方法,插值多项式可以表示为:

其中称为 i 次基函数

Matlab中拉格朗日插值法函数为:Language

功能:求已知点数据点的拉格朗日多项式

调用格式:f=Lagrange(x,y) 或者 f ’=Lagrange(x,y,x0)

其中,x为已知数据点的x坐标向量

y为已知数据点的y坐标向量

x0为插值点的x坐标

f ’为求得的拉格朗日多项式或x0处的插值

Lagrange.m

function f = Lagrange(x,y,x0)
%求已知数据点的拉格朗日多项式
%x是已知数据点的x坐标向量
%y是已知数据点的y坐标向量
%x0为插值点的x坐标
%f为求得的拉格朗日多项式或x0出的插值
syms t;
if(length(x)==length(y))n=length(x);
elsedisp('x和y的维数不相等!');return;
end                        %检错f=0.0;
for (i=1:n)l=y(i);for (j=1:i-1)l=l*(t-x(j))/(x(i)-x(j));end;for(j=i+1:n)l=l*(t-x(j))/(x(i)-x(j));%计算拉格朗日基函数end;f=f+l;       %计算拉格朗日插值函数simplify(f);    %化简if(i==n)if(nargin==3)f=subs(f,'t',x0);   %计算插值点的函数值.subs是替换函数,吧x0用t替换elsef=collect(f);          %将插值多项式展开f=vpa(f,6);            %将插值多项式的系数化成6位精度的小数endend
end

LagrangeInsert.m

x=[1 1.2 1.8 2.5 4];
y=[0.8415 0.9320 0.9738 0.5985 -0.7568];
f=Lagrange(x,y);
f=Lagrange(x,y,1.6)x1=0:2*pi;
y1=sin(x1);
x2=0:0.5:2*pi;
y2=Lagrange(x1,y1,x2);
plot(x1,y1,'b:',x2,y2,'r')

拉格朗日插值法(Lagrange)相关推荐

  1. 线性插值、抛物插值、Lagrange插值 | Lagrange拉格朗日插值法(一)

    Lagrange(拉格朗日)插值法 Lagrange插值法是一种多项式插值方法. 1. 线性插值(两点插值或一次插值) 线性插值就是通过两个采样点 ( x 0 , y 0 ) (x_0,y_0) (x ...

  2. 拉格朗日插值法(Lagrange插值法)

    插值介绍: 在离散数据的基础上补插连续函数,使得这条连续曲线通过全部给定的离散数据点. 插值是离散函数逼近的重要方法,利用它可通过函数在有限个点处的取值状况,估算出函数在其他点处的近似值. 这是百度百 ...

  3. 拉格朗日(Lagrange)插值法(C语言)

    题目 拉格朗日插值法流程图 (1)二次插值: #include<stdio.h> float lagelangri(float x[],float y[],float xx,int n) ...

  4. 实验二:Lagrange拉格朗日插值法之C语言代码

    拉格朗日插值多项式的算法就比前面的简单些,30行代码左右可以搞定,不过为了通俗易懂,这里我写了比较多的注释.题目:已知下列函数表: x | 0.56160 | 0.56280 | 0.56401 | ...

  5. 2021-01-07 matlab数值分析  插值法 拉格朗日插值法 牛顿插值法

    matlab数值分析  插值法 1 拉格朗日插值法 function yh=lagrange(x,y,xh) n=length(x); m=length(xh); yh=zeros(1,m); for ...

  6. 拉格朗日c语言实验报告,拉格朗日插值法C语言的实现(实验报告)(9页)-原创力文档...

    韩山师范学院 C语言程序设计 实验名称: 实验一:拉格朗日插值法C语言的实现 日期: 2011-9-19 任课教师: 赖国明 专业: 数学与应用数学 学号: 2009111437 姓名: 郑永锋 [实 ...

  7. 缺失值处理 - 拉格朗日插值法 - Python代码

    目录 缺失值处理 拉格朗日差值法的理论基础 拉格朗日插值法代码实现 其他数据预处理方法 缺失值处理 处理缺失值常用的办法可分为三类:删除记录.数据插补.不处理. 其中常见的数据插补法有: 如果通过删除 ...

  8. 缺失值填充6——拉格朗日插值法填充序列缺失值

    拉格朗日插值法的直观表述可以参考:https://www.zhihu.com/question/58333118 拉格朗日插值法 在python中可以直接使用: from scipy.interpol ...

  9. python拉格朗日插值法_【统计学】拉格朗日插值法的一种python实现方式

    昨天的一篇文章中提到了数据清洗中涉及缺失值,可通过删除数据.填补空值以及无视等方式进行处理.在空值填补方面,可用平均值.众数.中位数.固定值或者临近值进行填补. 删除数据这种方式比较适用于缺失值较少的 ...

最新文章

  1. 对JDBC的一些理解
  2. 人工智能落地,商业化比技术更重要
  3. 简约而不简单|值得收藏的Numpy小抄表(含主要语法、代码)
  4. Spark常用RDD算子 - saveAsTextFile、saveAsObjectFile 可保存到本地文件或hdfs系统中
  5. html 拼接onmouseout,HTML onmouseout事件用法及代码示例
  6. JavaScript高级程序设计学习笔记--事件
  7. 令人蛋疼的错误提示 0xcdcdcdcd ,0xdddddddd ,0xfeeefeee ,0xcccccccc ,0xabababab
  8. 软件项目开发流程逻辑图
  9. 电脑安装Chrome OS
  10. 使用python压缩文件夹
  11. ajax到底怎么读呢
  12. js中Object.freeze()函数的作用
  13. 单词学习-长久的过程
  14. matlab绘制磁场图,基于Matlab的电磁场图示化教学
  15. cocos creator接入微信登陆sdk ios篇
  16. 注册表键值的导出与导入
  17. 安卓的绝对布局AbsoluteLayout
  18. 金融工资和计算机工资待遇,【成都京东金融工资】计算机待遇-看准网
  19. Redis学习之路(八)封装Redis工具类
  20. QT 5.15 Android Windows 10开发环境搭建

热门文章

  1. linux树莓派网易云音乐,基于树莓派的红外遥控版网易云音乐播放器
  2. Educational Codeforces Round 88 (Rated for Div. 2)(AB)
  3. mysql中存储日期的类型_选择合适的 MySQL 日期时间类型来存储你的时间
  4. linux音量模块,Linux下的音量控制器alsamixer
  5. unity实现图片轮播效果_Unity实现图片轮播组件
  6. 解决Ubuntu IDEA 不能输入中文
  7. 数学--数论--HDU1222 狼和兔子(最大公约数)
  8. 数学--数论--组合数(卢卡斯+扩展卢卡斯)模板
  9. 数学--数论--二次探测定理
  10. JMeter压测入门简单使用