多种曲线插值方法

一、线性插值

       线性插值是获得数据点之间位置值的最简单方法。只需用直线段将各点连接起来即可。每个线段（以两个数据点为界）都可以独立进行插值。参数 mu 定义了在内插线上估计值的位置，在第一点为 0，在第二点为 1。对于两点之间的内插值，mu 的取值范围在 0 和 1 之间。以下所有方法都遵循这一惯例。与后续一样，下面是一段纯 C 代码来描述数学。

double LinearInterpolate(
   double y1,double y2,
   double mu)
{
   return(y1*(1-mu)+y2*mu);
}

一、余弦插值

        线性插值会导致每个点的不连续性。所以需要更平滑的插值函数，也许是最简单的是余弦插值。余弦的合适定向部分功能用于在相邻之间提供平滑过渡段。

double CosineInterpolate(
   double y1,double y2,
   double mu)
{
   double mu2;

   mu2 = (1-cos(mu*PI))/2;
   return(y1*(1-mu2)+y2*mu2);
}

 二、三次插值

        三次插值法是最简单的方法，它能提供线段之间真正的连续性。因此，它不仅需要线段的两个端点，还需要其两侧的两个点。因此，在下面的代码中，该函数需要 4 个点，分别标记为 y0、y1、y2 和 y3。这确实会引起如何在第一段和最后一段之间进行插值的问题。在这里的示例中，我并没有考虑这个问题。常见的解决方法是在序列的起点和终点多创建两个点，新点的斜率等于起点或终点的斜率。

double CubicInterpolate(
   double y0,double y1,
   double y2,double y3,
   double mu)
{
   double a0,a1,a2,a3,mu2;

   mu2 = mu*mu;
   a0 = y3 - y2 - y0 + y1;
   a1 = y0 - y1 - a0;
   a2 = y2 - y0;
   a3 = y1;

   return(a0*mu*mu2+a1*mu2+a2*mu+a3);
}

Catmull-Rom 样条曲线

        Paul Breeuwsma 提出了以下系数用于更平滑的插值曲线，它使用上一点与下一点之间的斜率作为当前点的导数。这就是通常所说的 Catmull-Rom 样条曲线。

   a0 = -0.5*y0 + 1.5*y1 - 1.5*y2 + 0.5*y3;
   a1 = y0 - 2.5*y1 + 2*y2 - 0.5*y3;
   a2 = -0.5*y0 + 0.5*y2;
   a3 = y1;

Hermite 插值法与三次插值法一样，需要 4 个点，因此可以实现更高的连续性。此外，它还有很好的张力和偏置控制。张力可用于收紧已知点的曲率。偏置则用于扭曲已知点的曲线。此处显示的示例默认张力和偏置值均为 0，可自行探索不同的张力和偏置值。 

/*
   Tension: 1 is high, 0 normal, -1 is low
   Bias: 0 is even,
         positive is towards first segment,
         negative towards the other
*/
double HermiteInterpolate(
   double y0,double y1,
   double y2,double y3,
   double mu,
   double tension,
   double bias)
{
   double m0,m1,mu2,mu3;
   double a0,a1,a2,a3;

	mu2 = mu * mu;
	mu3 = mu2 * mu;
   m0  = (y1-y0)*(1+bias)*(1-tension)/2;
   m0 += (y2-y1)*(1-bias)*(1-tension)/2;
   m1  = (y2-y1)*(1+bias)*(1-tension)/2;
   m1 += (y3-y2)*(1-bias)*(1-tension)/2;
   a0 =  2*mu3 - 3*mu2 + 1;
   a1 =    mu3 - 2*mu2 + mu;
   a2 =    mu3 -   mu2;
   a3 = -2*mu3 + 3*mu2;

   return(a0*y1+a1*m0+a2*m1+a3*y2);
}