计算方法上机题探讨.docVIP

第八题： 一，题目 从函数表 ≈。 选取距点x最近的三个结点，，进行二次插值，即取f(x)≈。 n=5, 三，计算结果 (1) x 0.15 0.31 0.47 f(x) 0.394039 0.380067 0.356927 (2) x 0.15 0.31 0.47 f(x) 0.394482 0.380225 0.361399 (3) x 0.15 0.31 0.47 f(x) 0.394473 0.380219 0.357222 四，结果分析 编程基本按照课本套公式，方法从（1）至（3），计算结果越来越精确。 五，源程序 （1） #includestdio.h int main(void) { int i,j; double a[6],b[6],x,y; printf(Enter x:); for(i=0;i6;i++) scanf(%lf,a[i]); printf(Enter f(x):); for(i=0;i6;i++) scanf(%lf,b[i]); for(j=0;j3;j++){ printf(Enter x:); scanf(%lf,x); for(i=0;xa[i];i++); y=b[i-1]*(x-a[i])/(a[i-1]-a[i])+b[i]*(x-a[i-1])/(a[i]-a[i-1]); printf(x=%lf , ,x); printf(y=%lf\n,y); } return 0; } （2） #includestdio.h int main(void) { int i,j; double a[6],b[6],x,y; printf(Enter x:); for(i=0;i6;i++) scanf(%lf,a[i]); printf(Enter f(x):); for(i=0;i6;i++) scanf(%lf,b[i]); for(j=0;j3;j++){ printf(Enter x:); scanf(%lf,x); for(i=0;xa[i];i++); if((a[i+1]-x)(x-a[i-2])) { y=b[i-1]*(x-a[i])/(a[i-1]-a[i])*(x-a[i+1])/(a[i-1]-a[i+1]); y=y+b[i]*(x-a[i-1])/(a[i]-a[i-1])*(x-a[i+1])/(a[i]-a[i+1]); y=y+b[i+1]*(x-a[i-1])/(a[i+1]-a[i-1])*(x-a[i])/(a[i+1]-a[i]); printf(x=%lf , ,x); printf(y=%lf\n,y); } else{ y=b[i-2]*(x-a[i-1])/(a[i-2]-a[i-1])*(x-a[i])/(a[i-2]-a[i]); y=y+b[i-1]*(x-a[i-2])/(a[i-1]-a[i-2])*(x-a[i])/(a[i-1]-a[i]); y=y+b[i]*(x-a[i-2])/(a[i]-a[i-2])*(x-a[i-1])/(a[i]-a[i-1]); printf(x=%lf , ,x); printf(y=%lf\n,y); } } return 0; } （3） #includestdio.h int main(void) { int i,j,k; double a[6],b[6],x,y,p; printf(Enter x:); for(i=0;i6;i++) scanf(%lf,a[i]); printf(Enter f(x):); for(i=0;i6;i++) scanf(%lf,b[i]); for(j=0;j3;j++) { y=0; printf(Enter x:); scanf(%lf,x); for(i=0;i6;i++){ p=1; for(k=0;k6;k++){ if(k!=i) p=p*(x-a[k])/(a[i]-a[k]);} y=y+p*b[i]; } printf(x=%lf , ,x); printf(y=%lf\n,y); } return 0; } 六，附图 （1） （2） （3） 第九题： 一，题目 已知直升机旋转机翼外形曲线轮廓线上的某些型值点及端点处的一阶导数值 y`(x0)=