数值分析小白学习之路(二) 插值法

最新推荐文章于 2022-05-24 22:34:08 发布

leon_pythonic

最新推荐文章于 2022-05-24 22:34:08 发布

阅读量2k

点赞数 2

分类专栏：数值分析文章标签：数值分析插值法

本文链接：https://blog.csdn.net/Major_out/article/details/78681234

版权

数值分析专栏收录该内容

1 篇文章 0 订阅

订阅专栏

插值法

设函数 $y = f(x), 对于 x \in [a,b], a \le x_0 < x_1 < \dots < x_n \le b 的函数值 y_1,y_2, \dots, y_n,$ , 存在函数 $P(x)$ , 使得 $P(x_i) = y_i, i = 0,1,\dots,n$ 成立, 则称 $P(x)$ 为 $f(x)$ 的插值函数. $x_i, i = 0,1,\dots,n$ 为插值节点, $[a,b]$ 为插值区间, 求插值函数 $P(x)$ 的方法称为插值法.
插值函数的作用就是用于函数的拟合. 通常插值函数的应用就是通过函数表来构造一个足够简单, 并且能够反映函数 $f(x)$ 的特征的函数.
一般来说, 我们可以将每个函数都看做是一个多项式, 对于一个不超过 n 次的多项式函数, 我们给出其定义 $P(x)= a_0 + a_1x + \dots + a_nx^n$ .
根据 $P(x)$ 的类型, 我们对 $P(x)$ 有不同的称呼. 如果插值函数是一个多项式, 那么我们就称 $P(x)$ 为插值多项式, 相应的插值法称之为多项式插值. 若 $P(x)$ 为分段的多项式, 则称之为分段插值. 若 $P(x)$ 为三角多项式, 就称之为三角插值.

插值函数的唯一性

对于区间 $[a,b]$ , 我们给出这个区间内的 n + 1个点 $a \le x_0 < x_1 < \dots < x_n \le b$ 和其对应的函数值 $y_i = f(x_i) (i =0,1\dots,n)$ , 对于给定的这个条件, 我们可以得到一个简单函数并且使其满足 $P(x_i) = y_i, i = 0,1,\dots,n$ .
将 n+1个点依次代入函数 $P(x)= a_0 + a_1x + \dots + a_nx^n$ 中, 我们就可以得到一个关于系数 $a_0, a_1, \dots,a_n$ 的 n+1元线性方程组.

$⎧ ⎩ ⎨ ⎪ ⎪ a 0 + a 1 x 0 + \dots + a n x n 0 = y 0, a 0 + a 1 x 1 + \dots + a n x n 1 = y 1, a 0 + a 1 x n + \dots + a n x n n = y n$ $\begin{equation} \left\{ \begin{array}{} a_0 + a_1x_0 + \dots + a_nx_0^n = y_0, &\\ a_0 + a_1x_1 + \dots + a_nx_1^n = y_1, \\ a_0 + a_1x_n + \dots + a_nx_n^n = y_n \end{array} \right. \end{equation}$
系数矩阵

$A = ⎡ ⎣ ⎢ ⎢ ⎢ ⎢ 11 ⋮ 1 x 0 x 1 \dots x n \dots \dots \dots x n 0 x n 1 ⋮ x n n ⎤ ⎦ ⎥ ⎥ ⎥ ⎥$ $\begin{equation} A= { \left[ \begin{array}{ccc} & 1 & x_0 & \cdots& x_0^n & \\ & 1 & x_1 & \cdots&x_1^n\\ & \vdots & \cdots & & \vdots \\ & 1 & x_n & \cdots&x_n^n\\ \end{array} \right ] } \end{equation}$
对于这个行列式, 我们可以通过范德蒙德公式来求得其行列式的值.
$det A = \prod _{i,j=0, i>j}^{n-1} (x_i - x_j \not = 0)$
$det A \not = 0$ , 则根据行列式的性质, 可以知道A 是可逆的, 即 AX = 0 只存在零解, AX = b只存在唯一解.
对应的参数确定下来, 则由参数唯一确定的插值多项式的也随之确定下来了, 故我们得到的插值函数 $P(x)$ 也是唯一的.