线性插值

概述

最近，在一个项目中使用到了线性插值相关知识。首先我们从最简单的线性插值说起。在理解线性插值的基础上，我们可以去构建更高维度的插值算法。

原理

假设我们已知坐标 $(x_0,y_0)$ 与 $(x_1,y_1)$ ，要得到 $[x_0,x_1]$ 区间内某一位置 $x$ 在直线上的值。根据图中所示，我们得到

$$\frac{y-y_0}{x-x_0}=\frac{y_1-y_0}{x_1-x_0}$$
由于$x$值已知，所以可以从公式得到$y$ 的值
$$y=y_0+(x-x_0)\frac{y_1-y_0}{x_1-x_0}=y_0+\frac{(x-x_0)y_1-(x-x_0)y_0}{x_1-x_0}$$
已知$y$ 求 $x$的过程与以上过程相同，只是$x$ 与$y$要进行交换。

1. 1维双线性插值

void OneDimansionalInterpolation(int point_n,double x,double* Xn,double *Fn,double output)
{
	int i;
	for(i=0;i<point_n;i++)
	{
		if(x < Xn[i])
		{
			if(i == 0)
			{
				*output = ((Xn[i+1] - x)*Fn(i) - (Xn[i] - x)*Fn[i+1])/(Xn[i+1] - Xn[i]);
			}
			else
			{
				*output = ((x - Xn[i-1])*Fn(i) - (x - Xn[i])*Fn[i-1])/(Xn[i] - Xn[i-1]);
			}
		}
	}
}

目前先写到这，后续会添加2维插值，3维插值，4维插值相关代码。