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

具体介绍请参考http://www.cnblogs.com/yabin/p/6426849.html

#include 
#include 
#include “Spline.h“

using namespace std;
using namespace SplineSpace;

int main（void）
{
	//异常抛出测试
		
	//double x0[2]={12};		//已知的数据点
	//double y0[2]={34};
	//try
	//{
	//	SplineInterface* sp = new Spline（x0y02）;
	//}
	//catch（SplineFailure sf）
	//{
	//	cout<	//}
	//getchar（）;		//程序暂停

	//单点插值测试

	//double x0[5]={12456};		//已知的数据点
	//double y0[5]={13425};
	//try
	//{
	//	//Spline sp（x0y05GivenSecondOrder00）;
	//	SplineInterface* sp = new Spline（x0y05）;	//使用接口，且使用默认边界条件
	//	double x=4.5;
	//	double y;
	//	sp->SinglePointInterp（xy）;	//求x的插值结果y
	//	cout<<“x=“<	//}
	//catch（SplineFailure sf）
	//{
	//	cout<	//}
	//getchar（）;	//程序暂停

	//多点插值测试

	double x0[5]={12456};		//已知的数据点
	double y0[5]={13425};

	double x[4] = {1.52.53.54.5};	//插值点
	double y[4];
	double leftBound=0RightBound=0;	//边界导数

	try
	{
		Spline sp（x0y05GivenSecondOrderleftBoundRightBound）;
		sp.MultiPointInterp（x4y）;			//求x的插值结果y
		for（int i = 0;i < 4;i++）
		{
			cout<<“x=“<		}
	}
	catch（SplineFailure sf）
	{
		cout<	}
	getchar（）;	//程序暂停

	// 自动插值测试

	// double x0[5]={12456};		//已知的数据点
	// double y0[5]={13425};

	// double x[10];	//插值点
	// double y[10];

	// try
	// {
		// SplineInterface* sp = new Spline（x0y05）;	//使用接口，且使用默认边界条件
		// sp->AutoInterp（10xy）;			//求x的插值结果y

		// for（int i = 0;i < 10;i++）
			// cout<		// cout<		// for（int i = 0;i < 10;i++）
			// cout<	// }
	// catch（SplineFailure sf）
	// {
		// cout<	// }
	// getchar（）;	//程序暂停
}

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2017-02-24 13:21  Spline\
     文件        1915  2017-02-24 13:18  Spline\main.cpp
     文件         160  2017-02-24 13:11  Spline\make.bat
     文件         218  2017-02-24 12:41  Spline\makefile1
     文件         159  2017-02-24 13:21  Spline\readme.txt
     文件        5399  2017-02-24 12:57  Spline\Spline.cpp
     文件        3405  2017-02-24 12:49  Spline\Spline.h