代数多项式插值问题-Read.PPTVIP

一般插值基函数的构造 思路： 分析（一） 分析（二） 简略记号 系数用结点处的二阶导数表示的三次样条函数 补充条件 牛顿差商型n次插值多项式(n次牛顿差商插值多项式) 计算牛顿差商插值多项式的步骤： （1）作差商表 （2）根据公式计算牛顿型插值多项式（表中对角线上各差商值就是Pn(x)的各项系数）。 余项公式的牛顿形式： 例4.5 已知函数y=f(x)的观测数据如例4.4，试用全部基点构造牛顿差商插值多项式，并用二次插值求f(3)的近似值。 解 用全部基点时，n=4，先作差商表,见例4.4。 P4(x)=f(0)+f[0,2](x- 0)+f[0,2,4](x- 0)(x- 2)+ f[0,2,4,5](x- 0)(x- 2)(x- 4)+f[0,2,4,5,6](x- 0)(x- 2)(x- 4)(x- 5) xi f(xi ) 1阶 2阶 3阶 4阶 0 2 4 5 6 2 2 -13 17 1 5 9 -4 13 0 -5 15 -1 5 1 =1+2x- x(x- 2)(x- 4)+x (x- 2)(x- 4)(x- 5) 用二次插值求f(3)时n=2，x=3，作内插取 x0=2, x1=4, x2=5 f(3)≈P2(3)=f(2)+ f[2,4](3 – 2)+ f[2,4,5](3- 2)(3- 4) =7- 5(3- 2)(3- 4) = 12 若插值基点等距分布，牛顿型插值多项式还可以利用差分得到简化,。 差分 设基点等距， xi=x0+ th, i=0,1,…,n. h称为步长。 一阶差分： 一阶向前差分、一阶前差 一阶向后差分、一阶后差 二阶差分： 二阶前差 二阶后差 k阶前差和k阶后差 零阶差分规定为 简记为 例4.6 设f(x)=x2+ x ,取步长 h=0.5,计算 解 利用差分可以简化等距基点的牛顿型插值公式，称为前差（或后差）插值公式 。 用数学归纳法可以证明前差、后差、差商之间的下述关系： △f(0)=f(0.5)- f(0) =0.75 , △f(0.5)=f(1)- f(0.5) =1.25 △2f(0)=△(△f(0))=△f(0.5)- △f(0) =1.25 – 0.75 = 0.5 令 x0=0, 则 x1=0.5, x2=1 4.4 三次样条插值 4.4.1 三次样条插值函数的概念 4.4.2 三次样条插值函数的求法 4.4.1 三次样条插值函数的概念 样条插值的思想：逐段选取适当的低次多项式，按一定的光滑性要求连接起来构成插值函数。 定义 设给定区间[a,b]上n+1个点 a=x0x1……xn=b 以及相应的函数值yi=f(xi ), i=0,1,…,n. 如果函数S(x)满足: (1)在每个子区间 [xk , xk+1](k=0,1,…,n-1)上， S(x)是不超过三次的多项式，且S(xi )=yi ,i=0,1,…,n. 在[a,b]上连续, 则称S(x)是f(x)在基点x0 , x1 , x2 ,…, xn上的三次样条插值函数. 称xoy平面上的点(xi , yi )(i=0,1,…,n)为样点。 例 给定区间[0,3]上3个点的函数值f(0)=0,f(1)=2,f(3)=4, 试求数 a,b,c,d,使函数 S(x)为给定点上的三次样条插值函数。 解 设 根据定义，由 得 d=0，故 则 由 得 由 得 由 得 由 得 求得 给定n+1个样点(xi , yi )(i=0,1,…,n) ，确定一个三次样条插值函数需要4n个独立条件。在定义中，已指定了4n-2个条件，即 所以，一般需补充指定2个边界条件。 * 第4章 插值法和曲线拟合 4.1 插值法的基本理论 4.2 拉格朗日插值多项式 4.3 牛顿差商插值多项式 4.4 三次样条插值 4.5 曲线拟合的最小二乘法 } 4.1.1 插值问题及代数多项式插值 函数y=f(x)给出一组函数值 x: x0 x1 x2 …… xn y: y0 y1 y2 …… yn 其中x0 ,x1,x2 ,…,xn是区间[a,b]上的互异点，要构造一个简单的函数φ(x) 作为f(x)的近似表达式，使满足 (插值原则、插值条件 ) 这类问题称为插值问题。 -----f(x)的插值函数, f(x) -----被插值函数 x0 ,x1,x2 ,…,xn -----插值节点 求插值函数的方法称为插值法。 若x∈[a,b]，需要计算f(x) 的近似值φ(x),则称x为插值点