数值分析例题II.ppt

*/25 五、六章内容小结 典型例题分析 练习题分析 课外练习题 ? ? ? ? 数值分析典型例题 II 若插值结点 x0, x1,…,xn 是(n+1)个互异点,则满足 插值条件P(xk)= yk (k = 0,1,···,n)的n次插值多项式 P(x)=a0 + a1x +……+ anxn 存在而且惟一。 多项式插值存在唯一性定理 Laglarge插值公式 插值基 ( k = 0, 1, 2, ······ , n ) 插值误差(余项) 其中, 注记1: n 次多项式插值的插值结点最佳选择是(n+1)次切比雪夫多项式零点。 已知 x0, x1, ···, xn 处的值 f(x0), ···, f(xn). ( j = 0,1, ······, n-1 ) ( j = 0,1,…,n-2 ) 均差定义 牛顿插值公式 ( k=1,2,···,n ) 注记2:均差具有对称性： 牛顿插值余项 ( j = 0, 1 ) 三次Hermite插值 割线法收敛速度:设割线法收敛, 分析其收敛阶 解: 设 x* 是方程的 f(x)=0 解.割线法迭代格式 利用带余项的牛顿插值公式 其中 ? ? ? 设迭代法 p 阶收敛, 则 ? 割线法收敛速度为: 近似为 1.618 阶 给定[a , b ]的分划: a = x0 x1 … xn = b. 已知f(xj) = yj (j = 0,1,···,n), 如果 满足: (1) S(x)在 [xj，xj+ 1]上为三次多项式; (2) S”(x)在区间[a，b]上连续; (3) S(xj) = yj ( j = 0,1,···,n). 则称S(x)为三次样条插值函数. 三次样条的定义 拟合函数: ?(x)=a0 ?0(x) + a1 ?1(x) + ······ +an ?n(x) 数据拟合的线性模型 离散数据 x x1 x2 ·········· xm f(x) y1 y2 ·········· ym 超定方程组 超定方程组最小二乘解: 勒让德多项式正交性定理 ? 证: 由勒让德方程 两式相减 ? 积分 ? ? 勒让德多项式的模平方 证: 由递推公式 ? 递推公式? 例1.设x0，x1，……，xn 是互异的插值结点，l0(x) 为对应于x0的拉格朗日插值基函数，试证明 x x0 x1 ·········· xn f(x) 1 0 ·········· 0 证 由基函数插值条件计算差商 ········· ········· 代入牛顿插值公式,并注意插值误差为零,则有 例2.设x0, x1, x2, …, xn为互异的结点,求证 Lagrange 插值基函数满足下列恒等式 (1) (2) ( k = 1,···,n ) 证: (1)令 在插值结点处 Pn(xj) = 0 ( j = 0,1,2,···,n ) n 次多项式 Pn(x)有 n+1 个相异零点?Pn(x) = 0 ? 所以 将 f(x) = xk (k ≤n) 代入, 得 (k =0,1,2,······,n) 问题: f(x)是(n+1)次多项式且最高次项系数为1，取互异的插值结点x0，x1，……，xn，构造插值多项式Pn(x)，证明： f(x) = Pn(x) + (x – x0) (x – x1)……(x – xn) (2) 取 f(x) = xk ? f(n+1)(x)=0 ? Rn(x) =0 例3. 设 P(x) 是不超过 n 次的多项式,而 ?n+1(x) =(x – x0)(x – x1)······(x – xn) 证明存在常数Ak( k =0，1，…，n)使得 证 由n次多项式插值得 其中 ? 证明: F[ x0, x1,······, xn ] = 例4. 记 ?n+1(x) =(x – x0)(x – x1)······(x – xn) ( j = 1,2, ···, n ) 对比Lagrange插值和Newton插值中 xn 的系数, 得 F[ x0, x1,······, xn ] = 例5. 如果X*是方程组GTGX=GTb的解，则X*是超定方程组GX=b的最小二乘解 证 由题设,有GT(b