【12月学习进度3/31—计算机图形学期末准备01】拉格朗日插值 + 三次Hermite插值

相关概念

• 型值点：事先给定的离散点
• 插值：得到的曲线通过所有的型值点
• 逼近：不要求通过给定的所有型值点，用给定型值点控制曲线形状
  

Hermite多项式

拉格朗日插值

拉格朗日插值是当 $j=0$ 时，即只给出函数值的Hermite多项式插值问题。

参考通俗易懂讲解视频

拉格朗日插值多项式：

其中 $l_i(x)$ 是插值基函数，它类似开关。
因为当且仅当 $x=x_i$ 时，$l_i(x_i)$为1，否则为0；
相当于$x=x_i$时，开关 $i$ 打开，其余开关均闭合，$f(x_i)=y_i$

这就保证了 给定 $n+1$ 个点，拉格朗日插值法可以构造出 $n$ 阶多项式，恰好穿过这 $n+1$ 个点
$f(x_i)=y_i$
$(i=0,1,2,...,n)$

Q：n+1个点为什么构造n阶插值多项式？
n阶多项式： $f(x)=a_0+a_1\ast x+a_2\ast x^2+a_3\ast x^3+...+a_n\ast x^n$
n阶多项式有 $n+1$ 个系数（$a_0,a_1,...a_n$）
n+1个点 <==> n+1个条件（方程）<=> 可解n+1个变量 <=> n阶多项式可确定

对于一组给定的点，拉格朗日插值法总能给出一条最低次数函数穿过这些点

但当一组连续均匀的点中出现一个“叛徒”（异常点），则拉格朗日插值就会想办法通过这个点，则就会出现龙格现象(Runge Phenomenon)。

【龙格现象】参考视频演示
利用多项式对某一函数进行逼近时，多项式次数越高，插值结果反而越偏离原函数的现象。（次数比较高时，产生剧烈的振荡）

例题

已知函数 $f(x)$，$f(144)=12,f(169)=13,f(225)=15$.
求 $f(x)$ 的二次拉格朗日插值多项式.

解：

1. 题目给定离散点：$(x_0,y_0),(x_1,y_1),(x_2,y_2)$
   其中 $x_0=144,y_0=12,x_1=169,y_1=13,x_2=225,y_2=15$ ;
2. 构造插值基函数 $l_i(x)$
   $l_0(x)=(x-x_1)(x-x_2)/(x_0-x_1)(x_0-x_2)$
   $l_1(x)=(x-x_0)(x-x_2)/(x_1-x_0)(x_1-x_2)$
   $l_2(x)=(x-x_0)(x-x_1)/(x_2-x_0)(x_2-x_1)$
3. 二次拉格朗日插值多项式 $L(x)$
   $L(x)=y_0\ast l_0(x)+y_1\ast l_1(x)+y_2\ast l_2(x)$
   （代入上述值即可求出）

三次Hermite插值

参考博文
（以下图片来自于参考博文）

三次Hermite插值要求在离散点处：

• 插值多项式的值 $P(x)=f(x)$
• 插值多项式的导数值 $P^{\prime }(x)=f^{\prime }(x)$

借鉴拉格朗日插值法构造插值基函数的思想

两点三次Hermite插值推导过程：
参考讲解视频：【《数值分析》| 华科 | 研究生基础课-哔哩哔哩】

【列出插值需满足的条件】

【构造基函数 ${\alpha }_0(x),{\alpha }_1(x),{\beta }_0(x),{\beta }_1(x)$】

以 ${\alpha }_0(x)$ 为例：

需满足条件：${\alpha }_0(x_0)=1,{\alpha }_0(x_1)=0,{\alpha }_0^{\prime }(x_0)=0,{\alpha }_0^{\prime }(x_1)=0$
则 $x_1$ 是 ${\alpha }_0(x)$ 的二重零点 → ${\alpha }_0(x)$ 有一因式为 $(x-x_1{)}^2$

借鉴拉格朗日插值基函数 $l_i(x)=(x-x_j)/(x_i-x_j)$ 在 $x_i$ 取1 ；
使用 $l_i(x)=(x-x_j{)}^2/(x_i-x_j{)}^2$

所以 ${\alpha }_0(x)=(x-x_1{)}^2/(x_0-x_1{)}^2\ast (a_{0}x+b_0)$
(最后乘上一次式 $a_{0}x+b_0$ 是为了构造 ${\alpha }_0(x)$ 为三次式）

同理构造三次式 ${\beta }_0(x),{\beta }_1(x)$ 乘以 $c_i$ 零次项系数

【求解待定系数 $a_i,b_i,c_i$】

将上述条件代入到构造的插值基函数式子中，解出待定系数。

三次Hermite插值多项式：