对应教材 5.2、定理 5.2、例 3。老师课堂未讲,按教材补齐,标“教材有、课堂未讲”。
与普通插值的区别:不仅要求函数值相等,还要求若干阶导数相等。
给定互异节点 (x_0,\dots,x_n),以及
f(xi), f′(xi), …, f(mi)(xi),i=0,…,n.
求次数不超过
N=i=0∑n(mi+1)−1
的多项式 (H_N(x)),使
HN(r)(xi)=f(r)(xi),r=0,…,mi.
这就是 Hermite 插值问题;(H_N) 称为 Hermite 插值多项式。
定理 5.2(存在唯一与余项):上述 Hermite 插值多项式存在且唯一。若 (f) 足够光滑,则
f(x)−HN(x)=(N+1)!f(N+1)(ξ)i=0∏n(x−xi)mi+1,
其中 (\xi) 在包含 (x) 与所有节点的区间内。
证明要点:把条件写成关于 (N+1) 个系数的线性方程组;系数矩阵是“广义范德蒙德”矩阵,由节点互异及导数条件可证非奇异。余项用 Rolle 定理反复求零(或构造辅助函数),与 Lagrange 余项同法。■
方法一(待定系数):设 (H_N(x)=\sum_{k=0}^N a_kx^k),代入所有条件,解线性方程组。
方法二(Newton 型):把重节点看成“节点重合”,用带重节点的差商表;重节点差商用导数代替:
f[xi,xi]=f′(xi),f[xi,xi,xi]=2!f′′(xi), …
方法三(降阶法):先作满足函数值的 (p_n(x)),再令
HN(x)=pn(x)+q(x)i=0∏n(x−xi)mi+1,
用导数条件定 (q)。例 3 就用这个办法。
给定
xf(x)−110014116215
并要求 (H_5’(-1)=1,\ H_5’(1)=0.1)。
先作满足函数值的三次插值 (p_3),再令
H5(x)=p3(x)+(ax+b)(x+1)x(x−1)(x−2).
由两个导数条件解出
a=36059,b=360161.
于是
H5(x)=p3(x)+36059x+161(x+1)x(x−1)(x−2).
这个例子说明:Hermite 插值不是“另起炉灶”,而是在普通插值上加一个在节点处为零的修正项,再用导数条件定系数。
覆盖核对
本节覆盖教材 5.2、定理 5.2、例 3;老师课堂未讲。下一节:5.3 三次样条。