插值法拉格朗日插值PPT学习教案

《插值法拉格朗日插值PPT学习教案》由会员分享，可在线阅读，更多相关《插值法拉格朗日插值PPT学习教案（17页珍藏版）》请在装配图网上搜索。

1、会计学1插值法拉格朗日插值插值法拉格朗日插值1问题的提出函数y = f(x)1）解析式未知；2）虽有解析式但表达式较复杂，通过实验计算得到的一组数据，即在某个区间a,b上给出一系列点的函数值yi=f(xi)，xx0 x1x2 xny=f(x)y0y1y2yn3）列表函数问题：无法求出不在表中的点的函数值，也不能进一步研究函数的其他性质，如函数的积分和导数等。因此需寻找y = f(x)的近似函数p(x),但要求p(xi) = f(xi) 。插值问题插值问题第1页/共17页已知精确函数已知精确函数 y = f(x) 在一系列节点在一系列节点 x0 xn 处测得函数值处测得函数值 y0 = f(x0

2、), yn = f(xn)，由，由此构造一个简单易算的近似函数此构造一个简单易算的近似函数 p(x) f(x)，满足条件，满足条件p(xi) = f(xi) (i = 0, n)。这里。这里的的 p(x) 称为称为f(x) 的的插值函数插值函数。最常用的插值。最常用的插值函数是函数是 ?多项式多项式x0 x1x2x3x4xp(x) f(x)第2页/共17页Taylor插值函数y = f(x)在点x0处展开有Taylor 多项式:200000)(! 2)( )()()(xxxfxxxfxfxpnnnxxnxf)(!)(.00)(可见: Pn(k)(x0)= f (k)(x0) k=0,1,n因此

3、, Pn(x)在点x0邻近会很好的逼近f(x). Taylor展开方法就是一种插值方法.泰勒插值要求提供 f(x) 在点x0处的各阶导数,这仅仅适用于 f(x) 相当简单的情况.第3页/共17页 设函数y = f(x)在区间a,b上有定义，且给出一系列点上的函数值yi=f(xi) (i=0,1,2,n)，求作n次多项式pn(x) 使得 pn (xi)= yi (i=0,1,2,n) 函数pn (x)为f(x)的插值函数；称x0,x1, xn称为插值节点或简称节点。插值节点所界的区间a,b称为插值区间。pn (xi)= yi 称为插值条件。 构造的n次多项式可表示为: Pn(x)= a0 + a

4、1x + a2x2+ anxn1.2 Lagrange插值第4页/共17页定理定理 (插值多项式的插值多项式的存在唯一性存在唯一性) 满足满足 的的 n 阶插值多项式是唯一存在的。阶插值多项式是唯一存在的。niyxPii,., 0,)( 证明：证明： ( 利用利用Vandermonde 行列式行列式论证论证)nnnnnnnnnyxaxaayxaxaayxaxaa.101111000010这是一个关于这是一个关于a0 , a1 , an 的的n+1元线性方程组元线性方程组,其系数行列式其系数行列式:10110)(),.,(ijjininnxxxxxV由于由于i j时时, xi xj ,因此因此

5、,即方程组有唯一解即方程组有唯一解. . 0),.,(10nnxxxV第5页/共17页2 拉格朗日插值公式拉格朗日插值公式 niyxPiin,.,0,)( 求求 n 次多项式次多项式 使得使得nnnxaxaaxP 10)(条件：条件：无重合节点，即无重合节点，即jixx ji n = 1已知已知 x0 , x1 ; y0 , y1 ，求，求xaaxP101)( 使得使得111001)(,)(yxPyxP 可见可见 P1(x) 是过是过 ( x0 , y0 ) 和和 ( x1, y1 ) 两点的直线。两点的直线。)()(0010101xxxxyyyxP 101xxxx 010 xxxx = y0

6、 + y1l0(x)l1(x) 10)(iiiyxl称为称为拉氏基函数拉氏基函数第6页/共17页直线方程的两点式：线性插值线性插值101001011)(yxxxxyxxxxxLl0(x)l1(x) 10)(iiiyxlL1(x)第7页/共17页抛物插值抛物插值2120210121012002010212)()()()()()()(yxxxxxxxxyxxxxxxxxyxxxxxxxxxLl0(x)l1(x)l2(x)第8页/共17页n 1li(x)每个每个 li 有有 n 个根个根 x0 xi xn njj i jiniiixxCxxxxxxCxl00)().().()( j i jiiiix

7、xCxl)(11)( njijjijixxxxxl0)()()( niiinyxlxL0)()(N次拉格朗日插值多项式次拉格朗日插值多项式与与 有关，而与有关，而与 无关无关节点节点f希望找到希望找到li(x)，i = 0, , n 使得使得 li(xj)= ；然后令；然后令 niiinyxlxP0)()(，则显然有，则显然有Pn(xi) = yi 。01ji ji n次多项式次多项式第9页/共17页 插值余项插值余项 /* Remainder */设节点设节点)1( nf在在a , b内存在内存在, 考察截断误差考察截断误差)()()(xLxfxRnn , baCfn bxxxan 10，且

8、，且 f 满足条件满足条件 ,用简单的插值函数用简单的插值函数L n(x)代替原复杂函数代替原复杂函数f(x),其精度取决于截断误差其精度取决于截断误差,即插值余项即插值余项.)()()()!1()()(210) 1(nnnxxxxxxxxnfxR即niinnxxnfxR0) 1()(! ) 1()()(,ba其中拉格朗日余项定理拉格朗日余项定理第10页/共17页注：注： 通常不能确定通常不能确定 , 而是估计而是估计 , x (a,b) 将将 作为误差估计上限。作为误差估计上限。1)1()( nnMxf niinxxnM01|)!1(当当 f(x) 为任一个次数为任一个次数 n 的的多项式多

9、项式时，时， , 可知可知 ，即插值多项式对于次数，即插值多项式对于次数 n 的的多项式是多项式是精确精确的。的。0)()1( xfn0)( xRn第11页/共17页例：例：已知已知233sin,214sin,216sin 分别利用分别利用 sin x 的的1次、次、2次次 Lagrange 插值计算插值计算 sin 50 并估计误差。并估计误差。 解：解：0 x1x2x185500 n = 1分别利用分别利用x0, x1 以及以及 x1, x2 计算计算4,610 xx利用利用216/4/6/214/6/4/)(1 xxxL这里这里)3,6(,sin)(,sin)()2( xxxfxxf而而

10、)4)(6(!2)()(,23sin21)2(1 xxfxRxx00762. 0)185(01319. 01 Rsin 50 = 0.7660444)185(50sin10 L0.77614外推外推 /* extrapolation */ 的实际误差的实际误差 3,421 xx利用利用sin 50 0.76008, 00660. 018500538. 01 R内插内插 /* interpolation */ 的实际误差的实际误差 内插通常优于外推。选择要计算的内插通常优于外推。选择要计算的 x 所在的区间的端点，插值效果较好。所在的区间的端点，插值效果较好。第12页/共17页n = 223)(

11、)(21)()(21)()()(4363463464363646342 xxxxxxxL)185(50sin20 L0.7654323cos21;)3)(4)(6(!3cos)(2 xxxxxxR 00077. 018500044. 02 Rsin 50 = 0.76604442次插值的实际误差次插值的实际误差 高次插值通常优于低次插值高次插值通常优于低次插值但绝对不是次数越高就越好，嘿嘿但绝对不是次数越高就越好，嘿嘿第13页/共17页ininjijjijnyxxxxxL)()()(00 拉格朗日插值多项式编程容易，只需双重循环拉格朗日插值多项式编程容易，只需双重循环 如果发现当前的插值方法不

12、够精确，就要增如果发现当前的插值方法不够精确，就要增加插值点的个数，则拉格朗日插值基函数加插值点的个数，则拉格朗日插值基函数 li(x)都将重新计算。都将重新计算。 牛顿插值法将讨论该问题牛顿插值法将讨论该问题。第14页/共17页例：已知数据表 xk10111213f(xk)2.302 62.397 92.484 92.564 9试用二次插值计算f(11.75)(计算过程保留4位小数) 解：因为11.75更接近12，故应取11，12，13三点作二次插值先作插值基函数已知x0=11,y0=2.397 9，x1=12,y0=2.484 9 ，x2=13,y2=2.564 9 2)13)(12()(

13、)()(2010210 xxxxxxxxxxxl1)13)(11()()()(2101201xxxxxxxxxxxl2)12)(11()()()(1202102xxxxxxxxxxxl 2(x)=9484. 21)13)(11(9397. 22)13)(12(xxxx9564. 22)12)(11(xx f(11.75)2(11.75)= 9484. 21)135 .11)(1175.11(9397. 22)1375.11)(1275.11(8 2.4639564. 22)1275.11)(1175.11(第15页/共17页例 已知x=1,4,9的平方根值，用拉格朗日插值公式求71/2解：x0=1, x1=4, x2=9f(x0)=1, f(x1)=2, f(x2)=3 L2(7) =(14)(19)(74)(79)* 1+(41)(49)(71)(79)* 2+(91)(94)(71)(74)* 3= 2.7(x0 x1)(x0 x2)(xx1)(xx2)f(x0)+(x1x0)(x1x2)(xx0)(xx2)f(x1)+(x2x0)(x2x1)(xx0)(xx1)f(x2)L2(x) =第16页/共17页