三次样条插值算法
1 总体说明
三次样条插值算法是一种计算量和效果都比较理想的插值算法。关于三次样条插值算法的原理这里不做过多的解释,下面的代码是我在网上收集了两种C++实现版本的基础上自己整合的一个版本。由于本人刚接触C++不久,水平有限。没有使用模板机制将代码做的更通用。关于算法实现有下面几点说明。
- 所有有关的类都被包含到SplineSpace命名空间中。
- SplineSpace中一个有三个类分别是异常类(SplineFailure),接口类(SplineInterface)和实现类(Spline)。有一个枚举类型说明边界条件(BoundaryCondition),取值为:GivenFirstOrder和GivenSecondOrder。分别对应I型边界条件和II型边界条件。
- 接口类定义了Spline在实现的过程中必须要有的三个方法:单点插值、多点插值和自动生成插值序列。
- 异常类是可能被实现类抛出的类,如果在实现类的运行过程中出现了已知数据过少构造失败、使用了外插值、设定输出点数过少等行为会抛出该类。因此应该将插值的过程用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;
}
上面代码就会抛出异常并显示“构造失败,已