前要
在数模比赛中,通常需要根据已知的数据点进行数据、模型的处理和分析,但有时候现有的数据是极少的,不足以支撑分析的进行,在这个时候我们可以通过 “模拟产生”一些新的但是又比较靠谱的值来满足需求 ,这就是插值的作用
一 一维插值
假设已经有 n + 1 n+1 n+1个结点 ( x i , y i ) , ( i = 0 , 1 , . . . , n ) (x_{i},y_{i}),(i=0,1,...,n) (xi,yi),(i=0,1,...,n),其中x_{i}互不相同,不妨假设 a = x 0 < x 1 < . . . < x n = b a=x_{0}<x_{1}<...<x_{n}=b a=x0<x1<...<xn=b,求任一插值点 x ∗ ( ≠ x i ) x^*(\neq x_{i}) x∗(=xi)处的插值 y ∗ y^* y∗
红点为插值数据:
思路:通过构造函数 y = f ( x ) y=f(x) y=f(x),使得 f ( x ) f(x) f(x)过所有节点,求 f ( x ∗ ) f(x^*) f(x∗)即可得到 y ∗ y^* y∗,即把函数构造出来以后,可以求出任意一 x x x值对应的 y y y值1
1.1 插值法的分类
1.2 插值法的概念
- 若 P ( x ) P(x) P(x)是次数不超过 n n n的代数多项式,即 P ( x ) = a 0 + a 1 x + . . . + a n x n , P(x) = a_{0}+a_{1}x+...+a_{n}x^n, P(x)=a0+a1x+...+anxn,
- 若 P ( x ) P(x) P(x)为分段多项式,则称为分段插值
- 若 P ( x ) P(x) P(x)为三角多项式,则成为三角插值
这里只讨论多项式插值和分段插值2
1.3 一般插值多项式原理
1.3.1 插值法原理
定理: 假设有 n + 1 n+1 n+1个互不相同的结点 ( x i , y i ) , ( i = 0 , 1 , 2 , . . . , n ) (x_{i},y_{i}),(i=0,1,2,...,n) (xi,yi),(i=0,1,2,...,n),则存在唯一的多项式:
L n ( x ) = a 0 + a 1 x + a 2 x 2 + . . . + a n x n 使 L n ( x j ) = y i ( j = 0 , 1 , 2 , . . . , n ) L_{n}(x)=a_{0}+a_{1}x+a_{2}x^2+...+a_{n}x^n \\ 使 \qquad L_{n}(x_{j})=y_{i} \quad(j=0,1,2,...,n) Ln(x)=a0+a1x+a2x2+...+anxn使Ln(xj)=yi(j=0,1,2,...,n)
可证:
构造方程组为
f ( x ) = { a 0 + a 1 x 0 + a 2 x 0 2 + . . . + a n x 0 n = y 0 a 0 + a 1 x 1 + a 2 x 1 2 + . . . + a n x 1 n = y 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . a 0 + a 1 x n + a 2 x n 2 + . . . + a n x n n = y n f(x)=\left\{ \begin{aligned} a_{0}+a_{1}x_{0}+a_{2}x_{0}^2+...+a_{n}x_{0}^n=y_{0}\\ a_{0}+a_{1}x_{1}+a_{2}x_{1}^2+...+a_{n}x_{1}^n=y_{1}\\ ....................................................\\ a_{0}+a_{1}x_{n}+a_{2}x_{n}^2+...+a_{n}x_{n}^n=y_{n}\\ \end{aligned} \right. f(x)=⎩⎪⎪⎪⎪⎨⎪⎪⎪⎪⎧a0+a1x0+a2x02+...+anx0n=y0a0+a1x1+