《计算方法-第2章-插值法拉格朗日插值.ppt》由会员分享,可在线阅读,更多相关《计算方法-第2章-插值法拉格朗日插值.ppt(36页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第二章 插值法12/6/20221第二章 插值法 2.1 引言引言 2.2 拉格朗日插值拉格朗日插值 2.3 均差与牛顿插值公式均差与牛顿插值公式 2.4 埃尔米特插值埃尔米特插值 2.5 分段低次插值分段低次插值12/6/20222本章要点用简单的函数(如多项式函数)作为一个复杂函数的近似,最简单实用的方法就是插值.本章主要介绍有关插值法的一些基本概念,及多项式插值的基础理论和几个常用的插值方法:拉格朗日插值、分段线性插值、牛顿插值、埃尔米特插值。12/6/20223能否存在一个性能优良、便于计算的函数一、插值问题12/6/20224这就是插值问题,上式为插值条件其插值函数的图象如下图12/
2、6/2022512/6/20226二、插值法的类型且满足其中 为实数,就称P(x)为插值多项式,相应的插值法称为多项式插值;若P(x)为分段的多项式,就称为分段插值;若P(x)为三角多项式,就称为三角插值。本章只讨论多项式插值与分段插值12/6/20227此插值问题可表述为如下:问题问题 求作次数 多项式 ,使满足条件这就是所谓的拉格朗日(拉格朗日(Lagrange)插值)插值。12/6/20228问题问题 求作一次一次式 ,使满足条件 从几何图形上看,表示过两点 的直线,因此可表示为如下点斜式:12/6/20229从几何图形上看,表示过两点的直线,因此也可表示为如下对称形式:其中,显然,12
3、/6/202210线性插值举例线性插值举例例1:已知 ,求代入点斜式插值多项式得 y=10.71428精确值为 10.723805,故这个结果有3位有效数字。12/6/20221112/6/202212 问题问题 求作二次二次式 ,使满足条件二次插值的几何解释是用通过三个点 的抛物线来近似考察曲线,故称为拋物插值。类似于线性插值,构造基函数,要求满足下式:12/6/20221312/6/202214(x0 x1)(x0 x2)(xx1)(xx2)f(x0)+(x1x0)(x1x2)(xx0)(xx2)f(x1)+(x2x0)(x2x1)(xx0)(xx1)f(x2)L2(115)=x0=100
4、,x1=121,x2=144f(x0)=10,f(x1)=11,f(x2)=12(100121)(100144)(115121)(115144)*10+(121100)(121144)(115100)(115144)*11+(144100)(144121)(115100)(115121)*12=10.7228抛物插值举例抛物插值举例例例2 2:L L2 2(x)=(x)=和用线性插值相比,有效数字增加一位12/6/202215为了构造 ,我们先定义n次插值基函数。定义:若n次多项式在n+1个节点上满足条件12/6/202216n+1次多项式对n=1及n=2时的情况前面已经讨论,用类似的推导方法
5、,可得到n次插值基函数为:12/6/202217且从而12/6/202218其中总总结结称为y=f(x)的拉格朗日插值多项式称为n次拉格朗日插值基函数12/6/202219例3:求过点(2,0)(4,3)(6,5)(8,4)(10,1)的 拉格朗日插值多项式。12/6/20222012/6/20222112/6/20222212/6/202223拉格朗日插值多项式的缺点:(1)插值基函数计算复杂(2)高次插值的精度不一定高12/6/202224一、插值余项满足不会完全成立因此,插值多项式存在着截断误差,那么我们怎样估计这个截断误差呢?12/6/20222512/6/202226令设其中证明:证
6、明:假设在区间a,b上f(x)的插值多项式为12/6/202227若引入辅助函数12/6/202228根据罗尔定理,再由罗尔定理,依此类推由于12/6/202229所以因此12/6/202230则注意(1)余项表达式只有在f(x)的高阶导数存在时 才能应用。(2)在内的具体位置通常不可能给出,所以,设12/6/202231例1:解:12/6/20223212/6/202233例2.并作图比较.解:12/6/202234不同次数的拉格朗日插值多项式的比较图Runge现象12/6/202235结果表明,并不是插值多项式的次数越高,插值效果越好,精度也不一定是随次数的提高而升高,这种现象在上个世纪初由Runge发现,故称为Runge现象.P48 2、3、4本章作业12/6/202236