三次样条插值算法C++实现

三次样条插值算法

1 总体说明

三次样条插值算法是一种计算量和效果都比较理想的插值算法。关于三次样条插值算法的原理这里不做过多的解释，下面的代码是我在网上收集了两种C++实现版本的基础上自己整合的一个版本。由于本人刚接触C++不久，水平有限。没有使用模板机制将代码做的更通用。关于算法实现有下面几点说明。

1. 所有有关的类都被包含到SplineSpace命名空间中。
2. SplineSpace中一个有三个类分别是异常类(SplineFailure),接口类(SplineInterface)和实现类(Spline)。有一个枚举类型说明边界条件(BoundaryCondition),取值为：GivenFirstOrder和GivenSecondOrder。分别对应I型边界条件和II型边界条件。
   
3. 接口类定义了Spline在实现的过程中必须要有的三个方法：单点插值、多点插值和自动生成插值序列。
   
4. 异常类是可能被实现类抛出的类，如果在实现类的运行过程中出现了已知数据过少构造失败、使用了外插值、设定输出点数过少等行为会抛出该类。因此应该将插值的过程用try...catch(SplineFailure sf)包裹起来。如：

double x0[2]={1,2};
double y0[2]={3,4};

try
{
  SplineInterface* sp = new Spline(x0,y0,2);
  //...
}
catch(SplineFailure sf)
{
  cout<<sf.GetMessage()<<endl;
}

上面代码就会抛出异常并显示“构造失败,已