拉格朗日插值的优缺点_拉格朗日插值法(图文详解)
对某个多项式函数,已知有给定的k + 1个取值点:
其中
对应着自变量的位置,而
对应着函数在这个位置的取值。
假设任意两个不同的xj都互不相同,那么应用拉格朗日插值公式所得到的拉格朗日插值多项式为:
其中每个
为拉格朗日基本多项式(或称插值基函数),其表达式为:
拉格朗日基本多项式
的特点是在
上取值为1,在其它的点
上取值为0。
范例
假设有某个二次多项式函数
,已知它在三个点上的取值为:
要求
的值。
首先写出每个拉格朗日基本多项式:
然后应用拉格朗日插值法,就可以得到
的表达式(
为函数
的插值函数):
此时代入数值
就可以求出所需之值:
。
证明
存在性
对于给定的k+1个点:
,拉格朗日插值法的思路是找到一个在一点
取值为1,而在其他点取值都是0的多项式
。这样,多项式
在点
取值为
,而在其他点取值都是0。而多项式
就可以满足
在其它点取值为0的多项式容易找到,例如:
它在点
取值为:
。由于已经假定
两两互不相同,因此上面的取值不等于0。于是,将多项式除以这个取值,就得到一个满足“在
取值为1,而在其他点取值都是0的多项式”:
这就是拉格朗日基本多项式。
唯一性
次数不超过k的拉格朗日多项式至多只有一个,因为对任意两个次数不超过k的拉格朗日多项式:
和
,它们的差
在所有k+1个点上取值都是0,因此必然是多项式
的倍数。因此,如果这个差
不等于0,次数就一定不小于k+1。但是
是两个次数不超过k的多项式之差,它的次数也不超过k。所以
,也就是说
。这样就证明了唯一性
几何性质
拉格朗日插值法中用到的拉格朗日基本多项式
(由某一组
确定)可以看做是由次数不超过n的多项式所组成的线性空间:
的一组基底。首先,如果存在一组系数:
使得,
,
那么,一方面多项式P是满足
的拉格朗日插值多项式,另一方面P是零多项式,所以取值永远是0。所以
。
这证明了
是线性无关的。同时它一共包含n+1个多项式,恰好等于
的维数。所以
构成了
的一组基底。
拉格朗日基本多项式作为基底的好处是所有的多项式都是齐次的(都是n次多项式)。
优点与缺点
拉格朗日插值法的公式结构整齐紧凑,在理论分析中十分方便,然而在计算中,当插值点增加或减少一个时,所对应的基本多项式就需要全部重新计算,于是整个公式都会变化,非常繁琐牛顿插值法来代替。此外,当插值点比较多的时候,拉格朗日插值多项式的次数可能会很高,因此具有数值不稳定的特点,也就是说尽管在已知的几个点取到给定的数值,但在附近却会和“实际上”的值之间有很大的偏差(如右下图)龙格现象,解决的办法是分段用较低次数的插值多项式。
重心拉格朗日插值法
重心拉格朗日插值法是拉格朗日插值法的一种改进。在拉格朗日插值法中,运用多项式
拉格朗日插值法的数值稳定性:如图,用于模拟一个十分平稳的函数时,插值多项式的取值可能会突然出现一个大的偏差(图中的14至15中间)
可以将拉格朗日基本多项式重新写为:
上面的表达式可以简化为:
于是拉格朗日插值多项式变为:
即所谓的重心拉格朗日插值公式(第一型)或改进拉格朗日插值公式。它的优点是当插值点的个数增加一个时,将每个
都除以
,就可以得到新的重心权
,计算复杂度为
,比重新计算每个基本多项式所需要的复杂度
降了一个量级。
将以上的拉格朗日插值多项式用来对函数
插值,可以得到:
因为
是一个多项式。
因此,将
除以
后可得到:
这个公式被称为重心拉格朗日插值公式(第二型)或真正的重心拉格朗日插值公式。它继承了(1)式容易计算的特点,并且在代入x值计算
的时候不必计算多项式
切比雪夫节点进行插值的话,可以很好地模拟给定的函数,使得插值点个数趋于无穷时,最大偏差趋于零切比雪夫节点进行插值可以达到极佳的数值稳定性。第一型拉格朗日插值是向后稳定的,而第二型拉格朗日插值是向前稳定的,并且勒贝格常数很小
参考来源
69: 59–67.
(英文)E. Meijering. A chronology of interpolation: From ancient astronomy to modern signal and image processing,. Proceedings of the IEEE: 323.
(英文)Julius Orion Smith III. Lagrange_Interpolation. Center for Computer Research in Music and Acoustics (CCRMA), Stanford University.
22 (5): 447–453.
The numerical stability of barycentric Lagrange Interpolation. IMA Journal of Numerical Analysis. 2004, 24 (4): 547–556.
(中文)李庆扬,王能超,易大义. 《数值分析》第4版. 清华大学出版社. 2001. ISBN 7-302-04561-5.
(中文)冯有前. 《数值分析》. 清华大学出版社. 2001. ISBN 7-810-82495-3.
(中文)拉格朗日插值多项式. 太原理工大学.
拉格朗日插值的优缺点_拉格朗日插值法(图文详解)相关推荐
- 拉格朗日插值的优缺点_拉格朗日插值法与牛顿插值法的比较
第 1 页 共 7 页 拉格朗日插值法与牛顿插值法的比较 [ 摘 要 ] 在生产和科研中出现的函数是多样的.对于一些函数很难找出其解析表达式.即使在某些情况下,可以写出函 数的解析表达式,但由于解析表 ...
- 拉格朗日插值的优缺点_拉格朗日插值的优缺点_拉格朗日插值法与牛顿插值法的比较...
第 1 页 共 7 页 拉格朗日插值法与牛顿插值法的比较 [ 摘 要 ] 在生产和科研中出现的函数是多样的.对于一些函数很难找出其解析表达式.即使在某些情况下,可以写出函 数的解析表达式,但由于解析表 ...
- 拉格朗日插值的优缺点_拉格朗日插值方法的理论基础。
只要做数据处理,一步不可避免的工作就是插值.而插值里面比较常用的方法之一就是拉格朗日插值法,这篇文章就跟大家一起来白活白活讲讲拉格朗日插值的理论基础. 为什么需要进行插值 我们进行数据处理的理想,当然 ...
- 拉格朗日插值的优缺点_拉格朗日与牛顿插值法的比较
第 1 页 共 7 页 拉格朗日插值法与牛顿插值法的比较 一. 背景 在工程和科学研究中出现的函数是多种多样的.常常会遇到这样的情况:在某个实际 问题中,虽然可以断定所考虑的函数 ) ( x f 在区 ...
- 结构专业规范大全_钢结构大全图文详解,手把手教你钢结构识图,反正我收藏了!...
钢结构大全图文详解,手把手教你钢结构识图,反正我收藏了! 钢结构是由钢制材料组成的结构,是主要的建筑结构类型之一.结构主要由型钢和钢板等制成的梁钢.钢柱.钢桁架等构件组成,并采用硅烷化.纯锰磷化.水洗 ...
- 拉格朗日插值的优缺点_浅谈拉格朗日插值
浅谈拉格朗日插值 在数值分析中,拉格朗日插值法是以法国十八世纪数学家约瑟夫·拉格朗日命名的一种多项式插值方法.许多实际问题中都用函数来表示某种内在联系或规律,而不少函数都只能通过实验和观测来了解.拉格 ...
- 拉格朗日插值的优缺点_「笔记」拉格朗日插值
简介 对于 \(k\) 次多项式函数 \(F(x)\). 若已知 \(k+1\) 个点值,则可构造出多项式. 有: \[F(x) = \sum_{i=1}^{k+1}y_i\prod_{i\not = ...
- 拉格朗日插值法(图文详解)
参考:https://www.cnblogs.com/ECJTUACM-873284962/p/6833391.html 在数值分析中,拉格朗日插值法是以法国十八世纪数学家约瑟夫·拉格朗日命名的一种多 ...
- pcb地线应该不应该做成环路_干货|图文详解EMC的PCB设计技术
除了元器件的选择和电路设计之外,良好的印制电路板(PCB)设计在电磁兼容性中也是一个非常重要的因素.PCB EMC设计的关键,是尽可能减小回流面积,让回流路径按照设计的方向流动.最常见返回电流问题来自 ...
最新文章
- 时间序列(五)股票分析
- Xcode的一些按钮使用
- python海贼王logo_Python实现的下载op海贼王网的图片
- linux vi模式替换,linux基础命令之:vi模式下查找和替换
- 这是啥?也太秀了吧?
- 如何批量在文件夹中建立php,怎么批量创建文件夹_一次性创建多个文件夹方法_一聚教程网...
- Linux shell脚本调试
- sql server根据表中数据生成insert语句
- python基础之psutil模块和发邮件(smtplib和yagmail)
- arcpy 查看arcgis的产品类型、产品许可与安装信息
- 计算机思维导论在线作业答案,大学计算机计算思维导论第讲习题及解析.pdf
- my security / csdn / 2000w / chakaifang.info
- 关于51仿真器驱动安装时进度卡在10%检测不到PZ51Tracker,以及额外下载的driver无法更新的解决方案
- ico 图标 生成 工具 网站
- 椭球面的c语言方程,椭球面
- 关于Google大陆手机号不能验证的问题
- 在线英语广播==好好学习,天天向上
- 信息安全系统所需要遵循的基本原则有哪些?
- APK 本地化/修改利器 —— 《android-apktool》
- 7-13 喊山 (30 分)
热门文章
- 分享微信好友代码php,js实现微信分享给好友功能
- 佛吉尼亚大学计算机世界排名,弗吉尼亚大学计算机科学专业介绍_计算机科学专业排名及就业方向和前景-小站留学...
- 【学科介绍】计算机科学
- 免费快递查询、天气查询、IP查询等接口
- 2017-2018 ACM-ICPC, Asia Daejeon Regional Contest 【11/12】
- Bartender2021实现安全远程标签打印,全新发布
- 个人微信订阅号(公众号)的注册流程
- ExtJs4 Ext.tab.Panel 选项卡
- SpringBoot 获取客户端 ip
- 八位七段数码管——学习笔记3