数值分析第五章.ppt

例1. 函数 取等距插值结点: -5, -4, -3, -2, -1, 0, 1, 2, 3, 4, 5 x∈[-5, 5] ?11(x)=(x+5)(x+4)(x+3)(x+2)(x+1)x(x-1)(x-2)(x-3)(x-4)(x-5) ?11(x) ? -4.9491 -4.5482 -3.7787 -2.7032 -1.4087 0.0000 1.4087 2.7032 3.7787 4.5482 4.9491 在[-5, 5]区间上,取11个切比雪夫结点 ( k=10, 9, 8, ···, 1, 0 ) ?11(x)=(x – x0)(x – x1)(x – x2)······(x – x10) ?11(x) ? 插值函数 H(x)= a0 + a1x + a2x2 + a3x3 插值条件 三次Hermite插值问题 m1 m0 H’(x) y1 y0 H(x) x1 x0 x 设 已知 例2. 插值条件: 求3次插值函数. 解: 设 积分,得 利用函数插值条件, 得 C2 = 0 , C1 = 6 ? H(x) = 3x2– 2x3 = (3 – 2x)x2 x 0 1 H(x) 0 1 H’(x) 0 0 利用基函数表示Hermite插值 x0 x1 1 0 0 0 0 1 0 0 x x0 x1 0 0 1 0 0 0 0 1 x 令 则 设 则 ? ? 0 1 1 0 0 0 0 1 0 0 设 ? x0 x1 0 0 1 0 0 0 0 1 x ? 设 两点Hermite插值的误差估计式 证明: 由插值条件知 构造辅助函数 利用 f(x) – H(x)=C(x)(x – x0)2(x – x1)2 取 x 异于 x0 和 x1, 有 ? ? 显然,F(t) 有三个零点 x0, x, x1, 由Roll定理知,存在 两个零点 t0, t1. 故 有四个相异零点 反复应用 Roll 定理, 得 F(4)(t) 知一个零点设为 * 拉格朗日插值公式 均差与牛顿插值公式 埃尔米特插值函数 样条插值方法 《数值分析》 Ch5 ? ? ? ? 引例.函数 sin x 的逼近 (1) 线性函数逼近 y0 = x (2)泰勒级数逼近 y1(x)= x – x3/3! + x5/5! (3)抛物线逼近(error=0.0559) y2=4x(π – x)/π2 (4)帕特逼近(error=0.0036) (1)复杂函数的计算; (2)函数表中非表格点计算 (3)光滑曲线的绘制; (4)提高照片分辩率算法 (5)定积分的离散化处理; (6)微分方程的离散化处理; (7)积分方程的离散化处理; 插值法的应用背景: 已知f(x)在点xi上的函数值 yi=f(xi), (i=0,1,2,···,n) 则称 P(x) 为 f(x) 的插值函数. 称 x0, x1, ······, xn为 插值结点; 称 f(x) 为被插值函数. 如果 P(x)=a0 + a1x +···+ anxn 满足: P(xk)= yk (k = 0,1,…,n) 设 f(x)∈C [a , b], 取点 a ≤x0＜x1＜···＜xn≤b 代数插值问题 插值条件 点,则满足插值条件 P(xk)= yk (k = 0,1,…,n) 的 n 次插值多项式 P(x)=a0 + a1x +……+ anxn 存在而且是唯一的。 证明: 由插值条件 P(x0)= y0 P(x1)=y1 ·············· P(xn)=yn 定理5.1 若插值结点x0, x1,···, xn 是 (n+1)个互异 ? 方程组系数矩阵取行列式 故方程组有唯一解. 从而插值多项式 P(x) 存在而且是唯一的. 例5.1 误差函数表可构造6次插值函数 x 0 0.5000 1.0000 1.5000 2.0000 2.5000 3.0000 y 0 0.5205