拉格朗日多项式插值.doc
《拉格朗日多项式插值.doc》由会员分享,可在线阅读,更多相关《拉格朗日多项式插值.doc(9页珍藏版)》请在咨信网上搜索。
拉格朗日多项式插值法浅析 摘要 拉格朗日插值多项式是一种最常见的多项式插值法,也是一种最常用的逼近工具。“学以致用 ”是每一门学科都致力追求的境界,数学自然也不例外。下面,探讨拉格朗日插值法的基本原理、如何构造拉格朗日多项式、拉格朗日多项式的误差界,并用 MATLAB程序来实现这一数学算法的自动化,为复杂的分析研究提供了一条数学算法的捷径。 【关键词】:拉格朗日多项式 算法实现 MATLAB 在科学研究和实际的工程设计中,几乎所有的问题都可以用来表示其某种内在规律的数量关系。但理想化的函数关系在实际工程应用中是很难寻找 的,对于那些没有明显解析式的函数关系表达式则只能通过实验观察的数据,利用多项式对某一函数的进行逼近,使得这个逼近函数能够反映的特性,而且利用多项式就可以简便的计算相应的函数值。例如我们不知道气温随日期变化的具体函数关系,但是我们可以测量一些孤立的日期的气温值,并假定此气温随日期变化的函数满足某一多项式。这样,利用已经测的数据,应用待定系数法便可以求得一个多项式函数f(x)。应用此函数就可以计算或者说预测其他日期的气温值。一般情况下,多项式的次数越多,需要的数据就越多,而预测也就越 准确。当然,构造组合多项式方法比较多,如线性方程求解、拉格朗日系数多项式以及构造牛顿多项式的分段差分和系数表等等,这里只对拉格朗日多项式插值法进行深入探讨。 一、拉格朗日多项式插值算法基本原理 函数在区间[a,b]上有定义,在是[ a,b]上取定的 N + 1个互异节点, 且在这些点处的函数值 , ,…,为已知, 即 yi =f (xi ) , (),若存在一个和近似的函数,满足 () (1) 则称 φ(x) 为 f (x) 的一个插值函数, 点为插值节点,(1)称为插值条件, 区间[a,b]称为插值区间, 而误差函数称为插值余项。即是求一个不超过N次多项式 () 满足 () 则成为的N次拉格朗日插值多项式。 二、拉格朗日插值多项式的构造 1、线性插值 当 n = 1时即为线性插值, 这也是代数插值最简单的形式。 根据给定函数在两个互异节点、的值、,用线性函数来近似代替。 由点斜式直线方程可得: (2) 公式(1)可整理写成: (3) 式(2)的右端的每一项都包含了一个线性因子,记 (4) 很容易看出来,,,因此式(3)中的多项式也给定两个定点: (5) 式(3)中的项和称为基于节点和的拉格朗日系数多项式(线性插值基函数)。利用这种记法,式(2)可以记为和式: (6) 也可以写成如下的矩阵: (7) 2、二次插值 当 n = 1时即为线性插值, 这也是常用代数插值。 根据给定函数在两个互异节点、、的值、、,构造次数不超过二次的多项式 来近似代替。使满足二次插值条件()。的参数直接由插值条件决定,并满足下面方程组: (6) 仿线性插值,用基函数的方法求解方程组。求二次式,,,因、是的两个零点,因此设,又,确定系数c=,从而导出: (7) 同理,构造出条件满足,,的插值多项式 (8) 构造出条件满足,,的插值多项式 (9) 式(7)(8)(9)中的项、和称为基于节点、和的拉格朗日系数多项式(二次插值基函数)。利用这种记法,相应的有: (10) 也可以写成如下的矩阵: 3、N次插值 当插值点增加到 N+ 1个时, 就可以通过 N+ 1个不同的已知点() 来构造一个次数为n的代数多项式 P (x)。类似二次插值, 先构造一个特殊的 n 次多项式,使其各点满足,,,因、…是的N个零点,因此设,又,确定系数,从而导出: (12) 相应的有: (13) 也可以写成如下的矩阵: 4、Lagrange插值余项 设,且,,,...,[a,b]为N+1个节点。如果x[a,b],则 (14) 其中是可以用来逼近的多项式: (15) 误差项形如 (16) C为区间内的某个值。 三、拉格朗日多项式插值实现流程 1、根据初始数据X的取值求出相应的Y值; 2、建立W*W的矩阵; 3、利用卷积公式计算基于节点的Lagrange系数矩阵; 4、求 四、MATLAB程序代码 Lagrange多项式逼近程序 function [C,L]=lagran(X,Y) w=length(X); n=w-1; L=zeros(w,w); for k=1: n+1 V=1; for j=1: n+1 if k~=j V=conv(V,poly(X(j)))/(X(k)-X(j)); end end L(k,:)=V; end C=Y*L; 五、实验结果 考虑[0.0,1.2]上的曲线。 (1)利用节点=0.0和=1.2构造线性插值多项式; (2)利用节点=0.0,=0.8和=1.8构造线性插值多项式; (3)利用节点=0.0,=0.4,=0.8和=1.2构造线性插值多项式。 解答: (1)输入 X=[0.0,1.2]; Y=cos(X); [C,L]=lagran(X,Y) 输出 C = -0.5314 1.0000 L = -0.8333 1.0000 0.8333 0 则一次逼近函数为 误差函数为 函数图像和误差函数图像 (2) 输入 X=[0.0, 0.6,1.2]; Y=cos(X); [C,L]=lagran(X,Y) 输出 C = -0.4004 -0.0508 1.0000 L = 1.3889 -2.5000 1.0000 -2.7778 3.3333 0 1.3889 -0.8333 0 则一次逼近函数为 误差函数为 函数图像和误差函数图像 (3) 输入 X=[0.0,0.4,0.8 1.2]; Y=cos(X); [C,L]=lagran(X,Y) 输出 C = 0.0922 -0.5651 0.0139 1.0000 L = -2.6042 6.2500 -4.5833 1.0000 7.8125 -15.6250 7.5000 0 -7.8125 12.5000 -3.7500 0 2.6042 -3.1250 0.8333 0 则一次逼近函数为 误差函数为 函数图像和误差函数图像 六、实验分析 拉格朗日多项式插值模型简单,结构紧凑,是经典的插值法。这种算法模型在科学的各个领域都有良好的应用。但是由于拉格朗日的插值多项式和每个节点 都有关,当改变节点个数时,需要重新计算。 一般情况下,多项式的次数越多,需要的数据就越多,误差就越小,从而预测也就越准确。例外发生了,龙格在研究多项式插值的时候,发现有的情况下,并非取节点越多多项式就越精确。著名的例子是,它的插值 函数在两个端点处发生剧烈的波动,造成较大的误差。研究发现,是舍入误差造成的。 的多项式逼近,基于[-1,1]的等距离节点 七、总结 本学期学习了数值方法这门学科,对我来说是非常欣喜的。它让我知道了高等代数和数学分析不光是纯理论,是可以用于实践生活中的。也让我感受到了数学的魅力,算法的强大。这门课程是为数不多的用理论解决实际问题的课程,这让我在枯燥的数学理论学习中,看到了数学应用的曙光,也感受到了数学的前景。在这门课的学习中,我始终是兴奋的。因为这门课很多的知识都是在数学分析中学过的,这让我非常有成就感。当然,老师在课堂中,时不时的给我们注入数学思想,提高我们的科研能力,对我们在以后的学习和工作中是有极大的帮助的,在这里,请允许我真诚的说句感谢! 这本书的第一章主要讲了函数的分析性质和二进制的基础,这里介绍了用于多项式计算的霍纳方法,也就是嵌套乘法。 第二章主要讲方程的解法,主要介绍了不动点迭代法、波尔查诺二分法和牛顿-拉夫森割线法。其实这三个方法在数学分析中是有接触的。不动点迭代法和牛顿-拉夫森割线法在证明递推函数的收敛性经常会用到,二分法也就是零点定理。所以,这章学习起来也是非常容易的。迭代法必须知道迭代公式和给出初始值,再进行逐步迭代,在做一些函数时,是非常慢的。二分法是需要在某个区间,端点函数值异号时可用,但不能求重根而且收敛速度慢。割线法收敛速度快,但在的导数为0时,则可能存在被零除错误。所以每种方法都有它的优缺点,在应用时依情况而定。 第四章主要讲了插值与多项式逼近,讲了泰勒级数逼近、拉格朗日多项式逼近和牛顿多项式逼近。泰勒级数要求函数存在n阶导数,从而寻找解析函数在区间的精确逼近多项式,当远离时,误差会很多。拉格朗日多项式简单易懂, 但是由于拉格朗日的插值多项式和每个节点都有关,当改变节点个数时,需要重新计算,并且在次数较高时,在端点处会剧烈的震荡,也就是我们常说的龙格现象。牛顿多项式把与通过递推关系联系起来了,不需要每个节点都重新计算,并且在进行人工计算是也是有很好的操作性的。 第五章讲曲线拟合,主要讲了最小二乘法和样条插值函数。最小二乘法在高等代数的选学部分有一定得接触,哪里讲的是线性方程的最小二乘法,在这一章线性方程的最小二乘法更深入的讲解,并且也介绍了型的幂函数的最小二乘法。线性二乘法拟合主要是根据已知条件,确定最佳的系数进行拟合。样条插值函数主要是利用分段函数在端点处的4个条件来进行拟合。 总之,通过这一学期的学习,对“函数”这个词有了不同的见解,函数是从初中到现在一直接触的,在研究函数分析性质时,总是机械的分析,从来没有深入的去想为什么要去研究函数,研究函数到底有什么用处,更别说去想在生活中怎么应用函数了。特别在样条插值函数拟合时就有很大的感触,通过4个分析性质就可以解决这么大的难题,我想我刚学习的时候,只能用目瞪口呆来形容。通过数值方法我似乎有了些其他的想法。- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 拉格朗日 多项式
咨信网温馨提示:
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,个别因单元格分列造成显示页码不一将协商解决,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【xrp****65】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【xrp****65】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,个别因单元格分列造成显示页码不一将协商解决,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【xrp****65】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【xrp****65】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。
关于本文