拉格朗日插值法(Lagrange)
拉格朗日插值法是基于基函数的插值方法,插值多项式可以表示为:
其中称为 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)相关推荐
- 线性插值、抛物插值、Lagrange插值 | Lagrange拉格朗日插值法(一)
Lagrange(拉格朗日)插值法 Lagrange插值法是一种多项式插值方法. 1. 线性插值(两点插值或一次插值) 线性插值就是通过两个采样点 ( x 0 , y 0 ) (x_0,y_0) (x ...
- 拉格朗日插值法(Lagrange插值法)
插值介绍: 在离散数据的基础上补插连续函数,使得这条连续曲线通过全部给定的离散数据点. 插值是离散函数逼近的重要方法,利用它可通过函数在有限个点处的取值状况,估算出函数在其他点处的近似值. 这是百度百 ...
- 拉格朗日(Lagrange)插值法(C语言)
题目 拉格朗日插值法流程图 (1)二次插值: #include<stdio.h> float lagelangri(float x[],float y[],float xx,int n) ...
- 实验二:Lagrange拉格朗日插值法之C语言代码
拉格朗日插值多项式的算法就比前面的简单些,30行代码左右可以搞定,不过为了通俗易懂,这里我写了比较多的注释.题目:已知下列函数表: x | 0.56160 | 0.56280 | 0.56401 | ...
- 2021-01-07 matlab数值分析 插值法 拉格朗日插值法 牛顿插值法
matlab数值分析 插值法 1 拉格朗日插值法 function yh=lagrange(x,y,xh) n=length(x); m=length(xh); yh=zeros(1,m); for ...
- 拉格朗日c语言实验报告,拉格朗日插值法C语言的实现(实验报告)(9页)-原创力文档...
韩山师范学院 C语言程序设计 实验名称: 实验一:拉格朗日插值法C语言的实现 日期: 2011-9-19 任课教师: 赖国明 专业: 数学与应用数学 学号: 2009111437 姓名: 郑永锋 [实 ...
- 缺失值处理 - 拉格朗日插值法 - Python代码
目录 缺失值处理 拉格朗日差值法的理论基础 拉格朗日插值法代码实现 其他数据预处理方法 缺失值处理 处理缺失值常用的办法可分为三类:删除记录.数据插补.不处理. 其中常见的数据插补法有: 如果通过删除 ...
- 缺失值填充6——拉格朗日插值法填充序列缺失值
拉格朗日插值法的直观表述可以参考:https://www.zhihu.com/question/58333118 拉格朗日插值法 在python中可以直接使用: from scipy.interpol ...
- python拉格朗日插值法_【统计学】拉格朗日插值法的一种python实现方式
昨天的一篇文章中提到了数据清洗中涉及缺失值,可通过删除数据.填补空值以及无视等方式进行处理.在空值填补方面,可用平均值.众数.中位数.固定值或者临近值进行填补. 删除数据这种方式比较适用于缺失值较少的 ...
最新文章
- 对JDBC的一些理解
- 人工智能落地,商业化比技术更重要
- 简约而不简单|值得收藏的Numpy小抄表(含主要语法、代码)
- Spark常用RDD算子 - saveAsTextFile、saveAsObjectFile 可保存到本地文件或hdfs系统中
- html 拼接onmouseout,HTML onmouseout事件用法及代码示例
- JavaScript高级程序设计学习笔记--事件
- 令人蛋疼的错误提示 0xcdcdcdcd ,0xdddddddd ,0xfeeefeee ,0xcccccccc ,0xabababab
- 软件项目开发流程逻辑图
- 电脑安装Chrome OS
- 使用python压缩文件夹
- ajax到底怎么读呢
- js中Object.freeze()函数的作用
- 单词学习-长久的过程
- matlab绘制磁场图,基于Matlab的电磁场图示化教学
- cocos creator接入微信登陆sdk ios篇
- 注册表键值的导出与导入
- 安卓的绝对布局AbsoluteLayout
- 金融工资和计算机工资待遇,【成都京东金融工资】计算机待遇-看准网
- Redis学习之路(八)封装Redis工具类
- QT 5.15 Android Windows 10开发环境搭建
热门文章
- linux树莓派网易云音乐,基于树莓派的红外遥控版网易云音乐播放器
- Educational Codeforces Round 88 (Rated for Div. 2)(AB)
- mysql中存储日期的类型_选择合适的 MySQL 日期时间类型来存储你的时间
- linux音量模块,Linux下的音量控制器alsamixer
- unity实现图片轮播效果_Unity实现图片轮播组件
- 解决Ubuntu IDEA 不能输入中文
- 数学--数论--HDU1222 狼和兔子(最大公约数)
- 数学--数论--组合数(卢卡斯+扩展卢卡斯)模板
- 数学--数论--二次探测定理
- JMeter压测入门简单使用