《计算方法》第四章 插值方法.ppt

计 算 方 法;§4 插值方法;§4.0 引言 ;§4.0 引言 ;§4.1 多项式插值问题的一般提法 ; 插值函数 p (x) 作为 f (x) 的近似，可以选自不同类型的 函数, 如 p (x) 为代数多项式、三角多项式、有理分式； 其函数性态可以是光滑的、亦可以是分段光滑的。其 中，代数多项式类的插值函数占有重要地位： ;x0 , x1, … , xn 插值节点, ;例题:; 插值的几何意义;插值方法的研究问题;求 n 次多项式 使得:;注意到插值节点;n = 1;线性插值与其基函数示意图;显然, 是过 、 、 三点的一条抛物线。;显然, 是过 、 、 三点的一条抛物线。;抛物线插值基函数;希望找到 li (x)，i = 0, …, n 使得 li (xj) = ?ij ；然后令;设 函数表 则满足插值条件的多项式;以下的问题：如何分析插值的余项？; x -1 0 1 2 f (x) -2 -2 1 2 ;解：利用Lagrange插值法有 ; ，且 f 满足条件 ,;? Lagrange插值法的插值余项;? Lagrange插值法的插值余项;任意固定 x ? xi (i = 0, …, n), 考察;由于 是不能确定，因此我们并不能确定误差的大小 但如能求出 ，那么用 逼近 的截断误差限是： 当 时， 当 时;当 f (x) 为任一个次数? n 的多项式时， , 可知， 即插值多项式对于次数? n 的多项式是精确的。; 给定 xi = i +1, i = 0, 1, 2, 3, 4, 5. 下面哪个是 l2(x) 的图像？;算例1;算例1;其截断误差为： 其中, 因为 可取 于是： ;用抛物线插值时，取所有节点，得到 ;算例2;§4.3 差商与差分 ;设插值多项式P(x)具有如下形式：;§4.3.1 差商的概念; 一般地， 关于 的 n 阶差商: ;差商的基本性质;差商的基本性质;差商的计算-差商表;已知;§4.3.2 差分;差分的定义;高阶差分;不变算子 I、移位算子 E;差分的性质;性质3：差商与差分有如下关系： 性质4：差分与导数有如下关系：;差分的计算;4.4 牛顿插值公式;4.4 牛顿插值公式;其中 显然 满足插值条件，且次数不超过 ,它就 是插值多项式，其系数为： 我们称 为牛顿插值多项式.; 已知 的函数表，求4 次牛顿插值多项式, 并求 ;从表中可以看到4 阶差商几乎为0，故取4次插值多项式即可， 于是：;0.40; 和 均是 n 次多项式,且均满足插值条件: 由多项式的唯一性, ,因而,两个公式 的余项是相等的,即 当插值多项式从 n-1 次增加到 n 次时, 拉格朗日型插值必须重新计算所有的基本插值多项式; 而对于牛顿型插值,只需用表格再计算一个 n 阶差商, 然后加上一项即可。;4.5 分段插值公式; 设函数 ， 在该区间 上取 个等距节点, 构造 的 次 拉格朗日插值多项式为;% lagrange.m function y=lagrange (x0,y0,x) n=length(x0);m=length(x); for i=1:m z=x(i);s=0; for k=1:n L=1; for j=1:n if j~=k L=L*(z-x0(j))/(x0(k)-x0(j)); end end s=s+L*y0(k); end y(i)