实验报告课程名称数值计算方法题目名称用三次样条插值设计山地.docVIP

 实 验 报 告 课程名称 数值计算方法 题目名称 用三次样条插值设计山地公路平面曲线 学生学院 应用数学学院 专业班级 05信计4班 学 号 3105009571 学生姓名 庞小波 指导教师 徐圣兵(老师) 2007年12月17日 应用数学 学院 信息与计算科学 专业4班、学号3105009571 姓名 庞小波 合作者______________ 教师评定 实验题目 用三次样条插值设计山地公路平面曲线 实验目的与要求 目的： １．熟悉插值法的有关理论和方法； ２．会综合运用数学建模知识和数值计算方法方面的知识； ３．会用三次插值解决山地公路平面曲线设计实际问题。 要求与步骤： 了解问题的背景。 建立数学模型。 运用计算方法编出可行程序。 写出运行结果并进行数据分析。 实验方案 求出M的值，列出方程在各个小区间上表达式。 单位为（1单位绳长） 程序如下： #includestdio.h #includemath.h float h=0.5; //给出h的值 int i; int j; float u=0.5/(0.5+0.5); //给出u的值 float r=0.5/(0.5+0.5); //给出r的值 double x[10]={1.5,2.0,2.5,3.0,3.5,4.0,4.5,5.0,5.5,6.0}; double y[10]={0.90,0.52,0.26,0.08,0.0,0.08,0.21,0.40,0.72,0.91}; float f[10],F[9],G[8]; void d() //求d的函数，结果为d=G[i],表示从0到7的d的值 { printf(求出的d值如下：\n); for(i=0;i=9;i++) f[i]=y[i]; for(i=0;i=8;i++) F[i]=(f[i+1]-f[i])/(x[i+1]-x[i]); for(i=0;i=7;i++) G[i]=6*(F[i+1]-F[i])/(x[i+2]-x[i]); for(i=0;i=7;i++) printf(d[%d]=%f\n,i,G[i]); } void gg() //输入方程 { int i,j; double c[8][9]; printf(生成的增广矩阵如下:\n); c[0][0]=2;c[0][1]=0.5;c[0][8]=G[0]; //第一行的数 for(j=2;j=7;j++) c[0][j]=0; c[1][0]=u;c[1][1]=2;c[1][2]=0.5;c[1][8]=G[1];//第二行 for(j=3;j=7;j++) c[1][j]=0; c[2][0]=0;c[2][1]=u;c[2][2]=2;c[2][3]=r;c[2][8]=G[2];//第三行 for(j=4;j=7;j++) c[2][j]=0; c[3][0]=0;c[3][1]=0;c[3][2]=u;c[3][3]=2;c[3][4]=r;c[3][8]=G[3];//第四行 for(j=5;j=7;j++) c[3][j]=0; c[4][0]=0;c[4][1]=0;c[4][2]=0;c[4][3]=u;c[4][4]=2;c[4][5]=r;c[4][8]=G[4];//第五行 for(j=6;j=7;j++) c[4][j]=0; c[5][4]=u;c[5][5]=2;c[5][6]=r;c[5][7]=0;c[5][8]=G[5];//第六行 for(j=0;j=3;j++) c[5][j]=0; c[6][5]=u;c[6][6]=2;c[6][7]=r;c[6][8]=G[6];//第七行 for(j=0;j=4;j++) c[6][j]=0; c[7][6]=u;c[7][7]=2;c[7][8]=G[7];//第八行 for(j=0;j=5;j++) c[7][j]=0; for(i=0;i=7;i++) { for(j=0;j=8;j++) printf(%9f,c[i][j]); printf(\n); } //实现追赶法 double A[8]; for(i=0;i=7;i++) A[i]=2; for(i=1;i=7;i++ )