第二章分段插值PPT课件

上传人:仙*** 文档编号:96009229 上传时间:2022-05-25 格式:PPT 页数:12 大小:119KB
收藏 版权申诉 举报 下载
第二章分段插值PPT课件_第1页
第1页 / 共12页
第二章分段插值PPT课件_第2页
第2页 / 共12页
第二章分段插值PPT课件_第3页
第3页 / 共12页
资源描述:

《第二章分段插值PPT课件》由会员分享,可在线阅读,更多相关《第二章分段插值PPT课件(12页珍藏版)》请在装配图网上搜索。

1、引言我们已经知道插值有多种方法:Lagrange 插值、 Newton插值、Hermit 插值等多种方式。插值的目的就是数值逼近的一种手段,而数值逼近,为得是得到一个数学问题的精确解或足够精确的解。那么,是否插值多项式的次数越高,越能够达到这个目的呢?现在,我们来讨论一下这个问题。 我们已经知道:f(x)在n+1个节点xi(i=0,1,2,n) 上的n次插值多项式Pn (x) 的余项 设想当节点数增多时会出现什么情况。由插值余项可知,当f(x)充分光滑时时,余项随n增大而趋于0的,这说明可用增加节点的方法达到这个目的,那么实际是这样吗? niinnxxnfxPxfxR011)()!()()()

2、()()(1901年龙格(Runge) 给出一个例子: 定义在区间-1,1上,这是一个光滑函数,它的任意阶导数都存在,对它在-1,1上作等距节点插值时,插值多项式情况,见图:-1-0.500.51-1.00-0.500.000.501.00f(x)n=8n=4从图中,可见,在靠近-1或1时,余项会随n值增大而增大,如P12(0.96)=36!但f(0.96)=0.25 22511xxf )(从图中,还可看见,在0附近插值效果是好的,即余项较小,另一种现象是插值多项式随节点增多而振动更多。 这种插值多项式当节点增加时反而不能更好地接近被插之数的现象,称为龙格现象龙格现象。这个任意阶可导的光滑函数

3、之所以出现这种现象,跟它在复平面上有x=1/5是奇点有关。俄罗斯数学家伯恩斯坦在1916年还给出如下定理:定理1:函数f(x)=|x|在-1,1上取n+1个等距节点x0=-1, xn=1,构造n次插值多项式Pn (x),当n增大时,除了-1,0,1,三点外,在-1,1中任何点处Pn(x)都不收敛于|x|。 上述现象和定理,告诉我们用高次插值多项式是不妥当的,从数值计算上可解释为高次插值多项式的计算会带来舍入误差的增大,从而引起计算失真。因此,实践上作插值时一般只用一次、二次最多用三次插值多项式。 那么如何提高插值精度呢?采用分段插值是一种办法。设f(x)是定义在a,b上的函数,在a,b上节点

4、a= x0 x1x2xn-1xn=b,的函数值为 y0 , y1 ,y2 ,yn-1 ,yn ,若函数(x)满足条件 (1) (x)在区间a , b上连续; (2) (x)在每个子区间xi , xi+1(i=0,1,2,n-1)上是次数为m的多项式; 则称(x)是f(x)在a ,b上的分段分段m m次插值多项式。次插值多项式。 m=1称为分段线性插值 m=2称为分段抛物线插值定义:分段线性插值的构造: 由定义, (x)在每个子区间xi , xi+1(i=0,1,2,n-1)上是一次插值多项式;分段线性插值的余项:定理:设f(x)在a,b上有二阶连续导数f(x) ,且| f(x)| m2,记:

5、h = max |xi+1-xi|,就有估计: |f(x)- (x) |=|R(x)| m2h2/8, xa, b。注意到h随分段增多而减少,因此用分段法提高精度是很好的途径.证明:由Lagrange 余项公式,当xxi, xi+1时 |f(x)- (x) |=|R(x)| = |f()(x-xi)(x- xi+1 )|/2! m2max |(x-xi)(x- xi+1 )|/ 2m2h2/8,上式右端与小区间的位置无关,证毕。11111 iiiiiiiiiixxxxxxxyxxxxyx,)(分段线性插值曲线图:y=f(x)x0 x1x2xnXY例:设 -1 x 1(1)将-1,1 10 等份

6、,用分段线性插值近似计算f(-0.96)。(2)将-1,1 n 等份,用分段线性插值近似计算,问如何选择步长h可使近似计算误差R10-4?解:(1)插值节点为xi=-1+ i/5 (i=0,1,10),h=1/5因为 -0.96-1,-0.8,取此区间为线性插值区间,其上的插值函数为所以f(-0.96) (-0.96)=0.0425322511xxf )(80112941080192302018020801.)(.).(.).(.)()( xxxxfxfx(2)插值节点为xi=-1+ ih (i=0,1,n),h=(b-a)/2=2/n由分段线性插值的余项估计: |f(x)- (x) |=|R

7、(x)| m2h2/8002801051210025117550251504211232222.)()(max|)(|)(|)()( hhxRxfmxxxfxxxfx分段二次插值分段二次插值即:选取跟节点x最近的三个节点xi-1,xi, xi+1进行二次插值,即在区间xi-1, xi+1,取:这种分段的低次插值叫分段二次插值,在几何上就是用分段抛物线代替y=f(x),故分段二次插值又和分段抛物插值。 11112iikikjijikjikxxxxyxLxf)()()()(实际上,上面介绍的分段低次插值,虽然具有计实际上,上面介绍的分段低次插值,虽然具有计算简便,收敛性有保证,数值稳定性又好且易在算简便,收敛性有保证,数值稳定性又好且易在计算机上实现等优点,但它却不能保证整条曲线计算机上实现等优点,但它却不能保证整条曲线的光滑性,从而不能满足某些工程技术上的要求,的光滑性,从而不能满足某些工程技术上的要求,从六十年代开始,首先由于航空、造船等工程设从六十年代开始,首先由于航空、造船等工程设计的需要而发展起来的样条插值(计的需要而发展起来的样条插值(spline)方法,方法,既保留了分段低次插值的各种优点,又提高了插既保留了分段低次插值的各种优点,又提高了插值函数的光滑性,在许多领域显得越来越广泛的值函数的光滑性,在许多领域显得越来越广泛的应用。应用。

展开阅读全文
温馨提示:
1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2: 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
3.本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 装配图网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

copyright@ 2023-2025  zhuangpeitu.com 装配图网版权所有   联系电话:18123376007

备案号:ICP2024067431-1 川公网安备51140202000466号


本站为文档C2C交易模式,即用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。装配图网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知装配图网,我们立即给予删除!