数值分析大作业三.docVIP

《数值分析》大作业三 一、算法设计方案 1.分段低次插值 分段低次插值方法有Lagrange法和Newton法等，由于满足插值条件的n次多项式（相应的节点数为n+1）唯一，采用Lagrange法和Newton法所得到的结果是一致的。这里采用的方法是Newton插值法。 2.样条插值 采用三弯矩法求解三次样条插值。 二、全部源程序 #includestdio.h #includemath.h double A[2][13]={ {0, 4.74, 9.5, 19, 38, 57, 76, 95, 114, 133, 152, 171,190}, {0, 5.32, 8.1, 11.97,16.15,17.1,16.34,14.63,12.16,9.69,7.03,3.99,0} }; double result_of_jie1[12]={0}; double result_of_jie2[11]={0}; double result_of_jie3[10]={0}; double B[13][13]={0}; /*存储三次样条差值的系数矩阵*/ double b[13]={0}; /*存储三次样条差值的结果向量系数矩阵*/ double x=0,f1[191]={0},f2[191]={0},f3[191]={0},f4[191]={0}; FILE *pfile=NULL; char *filename=F:\\result.txt; void jie_1( ); void jie_2( ); void jie_3( ); void xian_sheng_cha_zhi( ); void er_ci_cha_zhi( ); void san_ci_cha_zhi( ); void san_ci_yangtiao_cha_zhi( ); void Init_B( ); void Init_b( ); void Gauss( ); /***************主函数***************/ main(void) { pfile=fopen(filename,w); if(pfile==NULL) { printf(Error opening %sfor writting,filename); } xian_sheng_cha_zhi( ); er_ci_cha_zhi( ); san_ci_cha_zhi( ); san_ci_yangtiao_cha_zhi( ); fclose(pfile); } /********************** 函数名称：xian_sheng_cha_zhi( ) 函数功能：用线性插值法细化数据并输出 *********************/ void xian_sheng_cha_zhi( ) { int i=0,j=0; jie_1(); for(j=0;j12;j++) for(;x=A[0][j+1];x+=1,i++) f1[i]=A[1][j]+result_of_jie1[j]*(x-A[0][j]); printf(分段线性插值法细化数据输出\n); for(j=0;j=190;j++) { printf(%f ,f1[j]); fprintf(pfile,%f ,,f1[j]); } printf(\n); fprintf(pfile,\n); } /********************** 函数名称：er_ci_cha_zhi( ) 函数功能：用二阶多项式插值法细化数据并输出 *********************/ void er_ci_cha_zhi( ) { int i=0,j=0; jie_2(); x=A[0][0]; for(j=0;j10;j++) for(;x(A[0][j+1]+A[0][j+2])/2;x+=1,i++) f2[i]=A[1][j]+result_of_jie1[j]*(x-A[0][j])+result_of_jie2[j]*(x-A[0][j])*(x-A[0][j+1]); for(;x=190;x++,i++) f2[i]=A[1][