插值法

目录

• 插值法
  • 插值
  • 拉格朗日插值
    • 拉格朗日多项式
    • 定义
    • 应用
    • 公式
    • Code
    • 例子 && 证明：存在性与唯一性 && 几何性质
    • 优点与缺点
    • 重心拉格朗日插值法
  • 牛顿插值法（格雷戈里-牛顿公式）
    • 定义
    • Code
    • 优点

---

插值法

2018.1.9 学习笔记 （就当是个汇总吧）
Tags: 数学，数论

---

自己以前写的好naive啊QAQ丢个链接走人：
https://www.cnblogs.com/zwfymqz/p/10063039.html
https://blog.csdn.net/qq_35649707/article/details/78018944
不要看这篇了...

---

Another place.

插值

引用青峰碧陋室
  在工程应用和科学研究中，经常要研究变量之间的关系(y=f(x))。但对于函数(f(x))，常常得不到一个具体的解析表达式，它可能是通过观测或实验得到的一组数据((x,f(x)))，(x)为一向量;或则是解析表达式非常复杂，不便于计算和使用。因此我们需要寻找一个计算比较简单的函数(S(x))近似代替(f(x))，并使得(S(x)=f(x))，这种方法就称为插值法。

常见的插值法有：
　　一维插值法：拉格朗日插值、牛顿插值、分段低次插值、埃尔米特差值、样条插值。
　　二维插值法：双线性插值、双二次插值。

拉格朗日插值

摘自Angel_Kitty、attack's blog

拉格朗日多项式

  许多实际问题中都用函数来表示某种内在联系或规律，而不少函数都只能通过实验和观测来了解。如对实践中某个物理量进行观测，在若干个不同的地方得到相应的观测值，拉格朗日插值法可以找到一个多项式，其恰好在各个观测的点取到观测到的值。这样的多项式称为拉格朗日（插值）多项式。
　　对于给定的(n+1)个点((x_0,y_0),(x_1,y_1),cdots,(x_n,y_n))，对应于它们的次数不超过n的拉格朗日多项式L只有一个。如果计入次数更高的多项式，则有无穷个，因为所有与L相差 (lambda(x-x_0)(x-x_1)cdots(x-x_n))的多项式都满足条件。

定义

　　对某个多项式函数，已知有给定的(k+1)个取值点：((x_0,y_0),(x_1,y_1),cdots,(x_k,y_k))，其中(x_j)对应着自变量的位置，而(y_j)对应着函数在这个位置的取值。
　　假设任意两个不同的(x_j)都互不相同，那么应用拉格朗日插值公式所得到的拉格朗日插值多项式为：

[L(x)=sum_{j=0}^{k}y_jl_j(x) ]

其中每个(l_j(x))为拉格朗日基本多项式（或称插值基函数），其表达式为：

[l_j(x)=prod_{i=0,i eq j}^{k}frac{x-x_i}{x_j-x_i}=frac{x-x_0}{x_j-x_0}cdotsfrac{x-x_{j-1}}{x_j-x_{j-1}}frac{x-x_{j+1}}{x_j-x_{j+1}}cdotsfrac{x-x_k}{x_j-x_k} ]

　　拉格朗日基本多项式(l_j(x))的特点是在(x_j)上取值为1，在其它的点(x_i,i eq j)上取值为0。

应用

　　若给定(n)个互不相同的点((x,y))，则该(n-1)次多项式被唯一确定。利用拉格朗日插值，可以在(O(n))的时间内计算出某一项的值，还可以在(O(n^2))的时间内计算出给定的x所对应的y。

公式

　　若对于(n-1)次多项式，给定(n)个互不相同的点((x,y))，那么对于给定的(x)，
　　第(i)项的值为$$l(i)=y_iprod_{j=1,j eq i}^{n}frac{x-x_j}{x_i-x_j}$$
　　所对应的(y)为$$egin{aligned}y&=sum_{i=1}^{n}l(i)&=sum_{i=1}^{n}y_iprod_{j=1,j eq i}{n}frac{x-x_j}{x_i-x_j}end{aligned}$$

Code

double Lagrange(double *x,double *y,double X)//x和y分别为x和f(x)的取值，X为要求的点，返回f(X) 
{
	double res=0,tmp;
	for(int i=1;i<=n;++i)
	{
		tmp=1;
		for(int j=1;j<=n;++j)
			if(i!=j) tmp=tmp*(X-x[i])/(x[i]-x[j]);
		res+=y[i]*tmp;
	}
	return res;
}

例子 && 证明：存在性与唯一性 && 几何性质

　　见blog。

优点与缺点

　　拉格朗日插值法的公式结构整齐紧凑，在理论分析中十分方便，然而在计算中，当插值点增加或减少一个时，所对应的基本多项式就需要全部重新计算，于是整个公式都会变化，非常繁琐。这时可用重心拉格朗日插值法或牛顿插值法来代替。此外，当插值点比较多的时候，拉格朗日插值多项式的次数可能会很高，因此具有数值不稳定的特点，也就是说尽管在已知的几个点取到给定数值，但在附近却会和“实际上”的值之间有很大偏差，这类现象也称为龙格现象，解决方法是分段用较低次数的插值多项式。

重心拉格朗日插值法

　　见blog...已弃疗

---

牛顿插值法（格雷戈里-牛顿公式）

牛顿插值原理、利用差分的牛顿插值、利用均差的牛顿插值。弃疗*2

定义

　　已知函数(f(x))的(n+1)个互异节点(x_i)处的函数值(f(x_i))，则其牛顿插值多项式可以写为：

[N_n(x)=a_0+a_1(x-x_0)+cdots+a_n(x-x_0)(x-x_1)cdots(x-x_{n-1}) ]

　　其中，(a_k)为(f(x))的k阶差商（也叫均差），可以表示如下：

[a_k=f[x_0,x_1,cdots,x_{k-1},x_k] ]

　　也可以由函数值(f(x_i))线性表示为：

[a_k=sum_{i=0}^{n}frac{f(x_i)}{(x_i-x_0)cdots(x_i-x_{i-1})(x_i-x_{i+1})cdots(x_i-x_k)} ]

Code

double DifferenceQuotient(double *x,double *y,int k)//计算差商 
{
	double res=0,tmp;
	for(int i=1;i<=k;++i)
	{
		tmp=1;
		for(int j=1;j<=k;++j)
			if(i!=j) tmp/=(x[i]-x[j]);
		res+=y[i]*tmp;
	}
	return res;
}
double Newton(double *x,double *y,double X)
{
	double res=y[1],tmp=1;
	for(int i=2;i<=n;++i)
		tmp*=(X-x[i-1]), res+=DifferenceQuotient(x,y,i)*tmp;
	return res;
}

优点

　　当插值节点增加时，之前已计算的结果仍然能用，每增加一个节点，只要再增加一项即可，从而避免了重复计算。