C语言实现三次样条插值

《数值分析》实验习题1

已知函数值如下表：

试用三次样条插值求f(5.472)的近似值。

/*****************************************************************
《数值分析》实验习题1

机械工程 
现代制造技术教育部重点实验室
*****************************************************************/

#include
#include
double * cal_matrix(int , double *, double *, double *);
int cpm_x(double , double *, int);

int main(){ 
	int n = 10, i;
//	printf("请输入数据个数：\n");
//	scanf("%d", &n);
	double h[n], u[n-1], d[n+2], *p, test = 5.472, result, M[n+1];
	double x[n] = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10};
	double fx[n+2] = {0, 0.69314718, 1.0986123, 1.3862944, 1.6094378, 1.7917595, 1.9459101, 2.079445, 2.1972246, 2.3025851, 1, 0.1};
//	printf("请输入x的值：\n");
//	for(i=0; i0; i--)
		x[i] = (z[i] - c[i] * x[i-1]) / r[i];
	return x; 
}

int cpm_x(double x0, double *x, int n){
	int i;
	for(i=0; i





/*****************************************************************
《数值分析》实验习题2
机械工程 
现代制造技术教育部重点实验室
*****************************************************************/

#include
#include
#include
#define PI 3.14159265

typedef struct node  
{
	int num;
    double T;
    double S;
    double C;
    double R;
    struct node *next; 
}ParaList;

double fun(double );
double calT(ParaList* ,double ,double );
double calS(ParaList*);
double calC(ParaList*);
double calR(ParaList*);

int main(){
	double a = 0, b = 3, eps = 0.00001;
	
	//第二题参数 
//	double a = 1, b = 3, eps = 0.00001;

//	printf("请输入积分下限\n");
//	scanf("%lf", &a);
//	printf("请输入积分上限\n");
//	scanf("%lf", &b);
//	printf("请输入eps");
//	scanf("%lf", &eps);
	//算前四行的T，S，C的值并存入参数链表中 
	ParaList *phead, *p;
	int i;
	phead = (ParaList*)malloc(sizeof(ParaList));
	phead->num = 1;
	phead->T = (b - a) / 2 * (fun(a) + fun(b));
	p= phead;
	for(i=2; ;i++){
		p->next = (ParaList*)malloc(sizeof(ParaList));
		p->next->num = i;
		p->next->T = calT(p, a, b);
		p->next->S = calS(p);
		p->next->C = calC(p);
		p->next->R = calR(p);
		if(i>=5 && (fabs(p->next->R-p->R)next->R, p->next->num);
			break;
		}
		p = p->next;
	} 
}

double fun(double x){
	//第一题函数 
	double fresult;
	fresult = x * pow(1 + pow(x, 2), 0.5);
	return fresult;
}

//double fun(double x){
//	//第二题函数 
//	double fresult;
//	fresult = pow(3, x)*pow(x, 1.4)*(5*x+7)*sin(pow(x, 2));
//	return fresult;
//}

double calT(ParaList* s, double a, double b){
	//传入某一行的指针和积分上下限，用于计算下一行的T值 
	int i;
	double sum = 0, tresult;
	for(i=0; inum -1); i++)
		sum += fun(a + (i + 0.5) * (b - a) / pow(2, s->num -1));
	tresult = 0.5 * (s->T + (b - a) / pow(2, s->num -1) * sum);
	return tresult;
}

double calS(ParaList* t){
	//传入某一行的指针，用于计算下一行的S值 
	double sresult;
	sresult =(4 * t->next->T - t->T) / 3;
	return sresult; 
} 

double calC(ParaList* s){
	//传入某一行的指针，用于计算下一行的C值 
	double cresult;
	cresult = (16 * s->next->S - s->S) / 15;
	return cresult;
}

double calR(ParaList* c){
	//传入某一行的指针，用于计算下一行的R值 
	double rresult;
	rresult = (64 * c->next->C - c->C) / 63;
	return rresult;
}