工程算法插值法作业.docVIP

插值：三次样条插值与分段插值(20分) 已知飞机下轮廓线数据如下： x 0 3 5 7 9 11 12 13 14 15 y 0 1.2 1.7 2.0 2.1 2.0 1.8 1.2 1.0 1.6 飞机下轮廓线形状大致如下图所示： 要求分别用拉格朗日插值法、Newton插值法、分段线性插值法和三次样条插值法计算x每改变0.5时y的值，即x 取 0.5, 1, 1.5, … , 14.5 时对应的y值。比较采用不同方法的计算工作量、计算结果和优缺点。 《工程常用算法》综合实践作业一 完成日期： 13年 4月 13日 班级 学号 姓名 主要工作说明 自评成绩 711 2010071112 关雨 编写程序 A 711 2010071116 黄诗博 画程序框图，总结 711 2010071126 王馨晨 改错，讨论 分段差值 #include stdio.h main() { int i,n=10; float x[10]={0,3,5,7,9,11,12,13,14,15},y[10]={0,1.2,1.7,2.0,2.1,2.0,1.8,1.2,1.0,1.6},xx=0,yy,a,b,c,d; while(xx=15) { for(i=0;in;i++) if(x[i]xxxxx[i+1]) { a=x[i]; b=x[i+1]; c=y[i]; d=y[i+1]; } yy=c*(xx-b)/(a-b)+d*(xx-a)/(b-a); printf(x=%f,y=%f\n,xx,yy); xx=xx+0.5; }}  拉格朗日 #include stdio.h main() { float x[10]={0,3,5,7,9,11,12,13,14,15},y[10]={0,1.2,1.7,2.0,2.1,2.0,1.8,1.2,1.0,1.6},l[10],xx=0,yy,m,k; int i,j; while(xx=15)//拉格朗日函数 { for(i=0;i=9;i++) { m=1; for(j=0;j=9;j++) {if(j!=i) m=m*(xx-x[j]);} k=1; for(j=0;j=9;j++) {if(j!=i) k=k*(x[i]-x[j]); } l[i]=m/k; } yy=0; for(i=0;i=9;i++) {yy=l[i]*y[i]+yy; } xx=xx+0.5; printf(x=%f,y=%f\n,xx-0.5,yy); }}  牛顿差值 #include stdio.h #include conio.h #include malloc.h main() { float x[10]={0,3,5,7,9,11,12,13,14,15},y[10]={0,1.2,1.7,2.0,2.1,2.0,1.8,1.2,1.0,1.6} ,m[10],s[10],xx=0,yy;//输入初值 int i,j,k,h,p; for(i=0;i=9;i++) //构造差商 { for(j=0;j=i;j++) {m[j]=1; for(k=0;k=i;k++) { if(k!=j) m[j]=m[j]*(x[j]-x[k]); } } s[i]=0; for(h=0;h=i;h++) {s[i]=s[i]+y[h]/m[h]; } } for(i=0;i=9;i++) do//牛顿差值函数 {xx+=0.5;yy=0; for(i=0;i=9;i++) {p=1; for(j=0;ji;j++) {if(j!=i) p=p*(xx-x[j]); } yy=yy+p*s[i]; } printf(x=%f,y=%f\n,xx,yy); }while (xx15);} } 三次样条 #includestdio.h main() { int i; float x[10]={0,3,5,7,9,11,12,13,14,15},y[10]={0,