计算方法与实习大作业.docVIP

  1. 1、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
计算方法与实习大作业,数值计算方法大作业,计算方法与实习,计算方法上机作业,计算方法与实习第五版,大学生计算机实习报告,计算传热学的大作业,湖南工大计算机作业,计算流体力学大作业,大乐透计算方法

计算方法与实习大作业 班级:021011班 学号 姓名:师梦艳 插值法:拉格朗日插值多项式 在数值分析中,拉格朗日插值法是以法国十八世纪数学家约瑟夫·路易斯·拉格朗日命名的一种多项式插值方法。许多实际问题中都用函数来表示某种内在联系或规律,而不少函数都只能通过实验和观测来了解。如对实践中的某个物理量进行观测,在若干个不同的地方得到相应的观测值,拉格朗日插值法可以找到一个多项式,其恰好在各个观测的点取到观测到的值。这样的多项式称为拉格朗日(插值)多项式。数学上来说,拉格朗日插值法可以给出一个恰好穿过二维平面上若干个已知点的多项式函数。拉格朗日插值法最早被英国数学家爱德华·华林于1779年发现,不久后(1783年)由莱昂哈德·欧拉再次发现。1795年,拉格朗日在其著作《师范学校数学基础教程》中发表了这个插值方法,从此他的名字就和这个方法联系在一起。 数据建模有两大方法:一类是插值方法,另一类是拟合函数一般的说,插值法比较适合数据准确或数据量小的情形。然而Lagrange插值有很多种,1阶,2阶,…n阶。我们可以利用拉格朗日插值求方程,根据它的程序求原方程的图像。下面我具体介绍分析一下拉格朗日插值的算法设计及应用。 具体算法 1、基本概念 已知函数y=f(x)在若干点的函数值=(i=0,1,,n)一个差值问题就是求一“简单”的函数p(x):p()=,i=0,1,,n, (1) 则p(x)为f(x)的插值函数,而f(x)为被插值函数会插值原函数,,,,...,为插值节点,式(1)为插值条件,如果对固定点求f()数值解,我们称为一个插值节点,f()p()称为点的插值,当[min(,,,...,),max(,,,...,)]时,称为内插,否则称为外插式外推,特别地,当p(x)为不超过n次多项式时称为n阶Lagrange插值。 Lagrange插值公式 (1)线性插值 设已知 , 及=f() ,=f(),为不超过一次多项式且满足=,=,几何上,为过(,),(,)的直线,从而得到 =+(x-). (2) 为了推广到高阶问题,我们将式(2)变成对称式 =(x)+(x). 其中, (x)=,(x)=。均为1次多项式且满足 (x)=1且(x)=0。或(x)=0且(x)=1。 两关系式可统一写成= 。 (3) (2)n阶Lagrange插值 设已知,,,...,及=f()(i=0,1,.....,n),为不超过n次多项式且满足(i=0,1,...n). 易知=(x)+....+. 其中,均为n次多项式且满足式(3)(i,j=0,1,...,n),再由(ji)为n次多项式的n个根知=c.最后,由 c=,i=0,1,...,n. 总之,=,=式为n阶Lagrange插值公式,其中,(i=0,1,...n)称为n阶Lagrange插值的基函数。 3,Lagrange插值余项 设,,,...,[a,b],f(x)在[a,b]上有连续的n+1阶导数,为f(x)关于节点,,,...,的n阶Lagrange插值多项式,则对任意x[a,b], 其中,位于,,,...,及x之间(依赖于x),(x)=,此时只有在f(x)的(n+1) Xi 0.30 0.42 0.50 0.58 0.66 0.72 Yi 1.04403 1.08462 1.11803 1.15603 1.19817 1.23223 做五次拉格朗日插值多项式,求出x=0.46时的函数近似值 C语言下程序如下: #includestdio.h float Lagrange(float x[],float y[],float xx,int n) {int i,j; float *a,yy=0; a=new float[n]; for(i=0;j=n-1;i++) { a[i]=y[i]; for(j=0;j=n-1;j++) if(j!=i)a[i]*=(xx-x[j])/(x[i]-x[j]); yy+=a[i]; } delete a; return yy; } void main() { float x[6]={0.30,0.42,0.50,0.58,0.66,0.72}; float y[6]={1.04403,1.08462,1.11803,1.15603,1.19817,1.23223}; float xx=0.46,yy; yy=Lagrange(x,y,xx,6); printf(x=%f,y=%f\n,xx,yy); } 运行结果:

文档评论(0)

tianma2015 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档