数值插值方法

数值插值方法

基本目标是根据函数的离散数据构造一个简单易于计算的函数P(x)来代替原有的复杂函数。应用十分广泛，比如根据离散数据绘制光滑曲线、图形放大算法等等。

拉格朗日插值

定义域D上的f(x) 对于 n+1个点$((x_0,f(x_0)),(x_1,f(x_1)),...,(x_n,f(x_n)))$的拉格朗日插值
$$L_n(x)=\sum _{i=0}^{n}f(x_i)\prod _{j\ne i,0\le j\le n}\frac{x-x_j}{x_i-x_j}$$
定义余项
$$R(x)=f(x)-L_n(x)$$
假设f(x) n+1 阶连续可导，则对于固定的一个x，存在ξ,使得
$$R(x)=\frac{1}{(n+1)!}{f}^{(n+1)}(\xi )\prod _{j=0}^n(x-x_j)$$

证明

任取取一个$x^{\prime }\in D,且x^{\prime }\ne x_i$
$$k=R(x^{\prime })/\prod _{j=1}^n(x^{\prime }-x_j)$$
定义函数
$$g(x)=R(x)-\frac{1}{(n+1)!}k\prod _{j=0}^n(x-x_j)$$
显然 g(x) 也n+1 阶连续可导， 因为g(x)有n+2个零点(包括x’)， 运用n+1次罗尔定理，则存在ξ,使得$g^{(n+1)}(\xi )=0$.而$R^{(n+1)}(x)=f^{(n+1)}(\xi )$, $({\prod }_{j=0}^n(x-x_j){)}^{(n+1)}=(n+1)!$, 所以(n+1)!k = f^{(n+1)}(ξ)。

所以对于固定的一个x，存在ξ,使得
$$R(x)=\frac{1}{(n+1)!}{f}^{(n+1)}(\xi )\prod _{j=0}^n(x-x_j)$$

牛顿插值

定义域D上的f(x) 对于 n+1个点$((x_0,f(x_0)),(x_1,f(x_1)),...,(x_n,f(x_n)))$的牛顿插值
$$N_n(x)=a_0+a_1(x-x_0)+...+a_n\prod _{i=0}^{n-1}(x-x_0)$$
定义差商函数族 $g(i_0),g(i_0,i_1),...,g(i_0,i_1,...,i_n)$
$$g(i_0)=f(x_{i_0})$$
$$g(i_0,i_1)=\frac{g(i_0)-g(i_1)}{x_{i_0}-x_{i_1}}$$
$$g(i_0,i_1,...,i_n)=\frac{g(i_0,i_1,...,i_{n-1})-g(i_1,i_2,...,i_n)}{x_{i_0}-x_{i_n}}$$
$$a_0=g(0),a_1=g(0,1),...,a_n=g(0,1,...,n)$$
牛顿插值本质上和拉格朗日插值一致，只是形式上不同， 方便后续加入点。

证明

引理1

$$g(i_0,i_1,...,i_m)=\sum _{p=0}^{m}f(x_{i_p})/\prod _{q\ne p,0\le q\le m}(x_{i_p}-x_{i_q})(m=0,...,n)$$

引理1证明

对m使用数学归纳法证明

1° 当m=0时，显然成立

2° 假设当m=k时成立。当m=k+1时

$$\begin{array}{rl}g(i_0,i_1,...,i_k)& =\sum _{p=0}^{k}f(x_{i_p})/\prod _{q\ne p,0\le q\le k}(x_{i_p}-x_{i_q})\\ & =f(x_{i_0})/\prod _{1\le q\le k}(x_{i_0}-x_{i_q})+\sum _{p=1}^{k}f(x_{i_p})/(x_{i_p}-x_{i_0})/\prod _{q\ne p,1\le q\le k}(x_{i_p}-x_{i_q})\end{array}$$
$$\begin{array}{rl}g(i_1,i_2,...,i_{k+1})& =\sum _{p=1}^{k+1}f(x_{i_p})/\prod _{q\ne p,1\le q\le k+1}(x_{i_p}-x_{i_q})\\ & =f(x_{i_{k+1}})/\prod _{1\le q\le k}(x_{i_{k+1}}-x_{i_q})+\sum _{p=1}^{k}f(x_{i_p})/(x_{i_p}-x_{i_{k+1}})/\prod _{q\ne p,1\le q\le k}(x_{i_p}-x_{i_q})\end{array}$$
$$\begin{array}{rl}g(i_0,i_1,...,i_{k+1})& =\frac{g(i_0,i_1,...,i_k)-g(i_1,i_2,...,i_{k+1})}{x_{i_0}-x_{i_{k+1}}}\\ & =f(x_{i_0})/\prod _{1\le q\le k+1}(x_{i_0}-x_{i_q})+f(x_{i_{k+1}})/\prod _{0\le q\le k}(x_{i_{k+1}}-x_{i_q})+\sum _{p=1}^{k}f(x_{i_p})/(x_{i_p}-x_{i_{k+1}})/\prod _{q\ne p,0\le q\le k+1}(x_{i_p}-x_{i_q})\\ & =p=0\end{array}$$