数值计算方法第四章—插值法

插值法

本文参考书为马东升著《数值计算方法》

代数插值

• 插值

  已知某些（有限）点的函数值求其余点的函数值

• 代数多项式插值

  $n+1$ 个互异节点处满足插值条件 $P(x_i)=y_i,i=0,1,\cdots ,n$ 的 $n$ 次代数多项式 $P(x)=a_0+a_{1}x+a_2{x}^2+\cdots +a_n{x}^n$ 是唯一的

  ---

  当 $f(x)$ 是次数不超过 $n$ 的多项式时，其 $n$ 次插值多项式就是 $f(x)$ 本身

  例： 已知函数 $f(x)=56x^3+24x^2+5$ 在点 $2^0,2^1,2^5,2^7$ 的函数值，求其三次插值多项式

  解： 因为 $f(x)$ 次数为 $3$ 不超过三次，且其为多项式，故其三次插值多项式 $P(x)=f(x)=56x^3+24x^2+5$

• 抛物线（二次）插值

  已知 $f(x)$ 在三个互异点 $x_0,x_1,x_2$ 的函数值 $y_0,y_1,y2$ 构造一个次数不超过二次的多项式
  $$P(x)=a_0+a_x+a_2{x}^2$$
  使其满足插值条件
  $$P(x_i)=y_i,(i=0,1,2)$$

---

拉格朗日插值

• 线性插值

  已知 $f(x_0)=y_0,f(x_1)=y_1,x_0\ne x_1$ 要构造线性函数 $P(x)=a_0+a_{1}x$ ，使其满足插值条件 $P(x_0)=y_0,P(x_1)=y_1$
  $$\frac{y-y_0}{x-x_0}=\frac{y_1-y_0}{x_1-x_0}P(x)=y_0+\frac{y_1-y_0}{x_1-x_0}(x-x_0)$$
  表达式为
  $$P(x)=(\frac{x-x_1}{x_0-x_1})y_0+(\frac{x-x_0}{x_1-x_0})y_1=l_0(x)y_0+l_1(x)y_1$$
  公式的结构，它是两个一次函数的线性组合

  线性插值基函数
  $$l_0(x)=\frac{x-x_1}{x_0-x_1},l_1(x)=\frac{x-x_0}{x_1-x_0}$$

  例： 给定平方表

  x	100	121
  y	10	11

  试用线性插值求 $115$ 的平方根

  解： $x_0=100,x_1=121,x=115$

  故$P(x)=(\frac{x-x_1}{x_0-x_1})y_0+(\frac{x-x_0}{x_1-x_0})y_1$

  $\sqrt{115}\approx P(115)=10.714$

  ---

• 抛物线插值

  已知 $f(x)$ 在三个互异点 $x_0,x_1,x_2$ 的函数值 $y_0,y_1,y2$ 构造一个次数不超过二次的多项式
  $$P(x)=a_0+a_x+a_2{x}^2$$
  使其满足插值条件
  $$P(x_i)=y_i,(i=0,1,2)$$
  表达式为
  $$P(x)=l_0(x)y_0+l_1(x)y_1+l_2(x)y_2$$
  其中插值基函数为
  $$l_0=\frac{(x-x_1)(x-x_2)}{(x_0-x_1)(x_0-x_2)}{l}_1=\frac{(x-x_0)(x-x_2)}{(x_1-x_0)(x_1-x_2)}{l}_2=\frac{(x-x_0)(x-x_1)}{(x_2-x_0)(x_2-x_1)}$$

  例： 已知 $100,121,144$ 的平方根，计算 $115$ 的平方根的近似值

  x	100	121	144
  y	10	11	12

  解： $P(x)=\frac{(x-121)(x-144)}{(100-121)(100-144)}\times 10+\frac{(x-100)(x-144)}{(121-100)(121-144)}\times 11+\frac{(x-100)(x-121)}{(144-100)(144-121)}\times 12$

  $f(115)\approx P(115)=10.7228$

  ---

• 拉格朗日插值多项式

  • 若 l k ( x ) = ∏ i = 0   , i ≠ k n