实验1拉格日插值与牛顿插值.docVIP

西华数学与计算机学院上机实践报告 课程：级：：名称： 上机实践日期： 上机实践编号： 一、目的 二、内容与y=f(x)的数据表如下，求t=0.63处的函数值z=f(t) 。 i xi yi i xi yi 1 0.10 0.904837 6 0.57 0.565525 2 0.15 0.860708 7 0.70 0.496585 3 0.25 0.778801 8 0.85 0.427415 4 0.40 0.670320 9 0.93 0.394554 5 0.50 0.606531 10 1.00 0.367879 三、用环境 四#includestdio.h main() { float x[8],y[8]; float t,z=0.0,s; int n=0,i=0,tp=0,start,end,j; printf(请输入节点控制数n); printf(\n); scanf(%d,n); printf(请输入节点t的值); printf(\n); scanf(%f,t); //节点赋值 for(i=0;in;i++) { printf(请输入插点序列x%d与y%d的值\n,i,i); scanf(%f %f,x[i],y[i]); } //记录t节点的位置 while(x[tp]t tpn) { tp++; } //确定起点下标和终点下标位置//用8个节点就选取值离t最近的8个节点 start=tp-4; end=tp+3; if(start0) { start=0; end=start+7; } if(endn-1) { end=n-1; start=end-7; } //用所选插值点做拉格朗日插值多项式运算 for(i=start;iend;i++) { s=1.0;//插值基函数初值为1 for(j=start;jend;j++) { if(j!=i) { s=s*(t-x[j])/(x[i]-x[j]); } } z=z+s*y[i]; } printf(拉格朗日插值结果为%f,z); printf(\n); } 2.牛顿插值 #includestdio.h #define max 20 main() { double x[max],y[max];//插值点xi以及对应函数yi的函数值 double t,fw;//t为待求的参数点fw为对应要求的值 double ft[max][max];//表示各阶差商第二维为阶数 int n=0,i=0,j=0; double b[max]; /*printf(请输入节点控制数n); printf(\n); scanf(%d,n); printf(请输入节点t的值); printf(\n); scanf(%f,t); //节点赋值 for(i=0;in;i++) { printf(请输入插点序列x%d与y%d的值,i,i); printf(\n); scanf(%f%f,x[i],y[i]); } printf(\n);*/ n=10; t=0.63; x[0]=0.10;y[0]=0.904837; x[1]=0.15;y[1]=0.860708; x[2]=0.25;y[2]=0.778801; x[3]=0.40;y[3]=0.670320; x[4]=0.50;y[4]=0.606531; x[5]=0.57;y[5]=0.565525; x[6]=0.70;y[6]=0.496585; x[7]=0.85;y[7]=0.427415; x[8]=0.93;y[8]=0.394554; x[9]=1.00;y[9]=0.367879; for(i=0;in;i++) { printf(%f\t%f\n,x[i],y[i]); } //计算差商 for(i=0;in;i++) { ft[i][0]=y[i];//计算0阶差商 } for(i=1;in;i++) { for(j=1;j=i;j++) { ft[i][j]=(ft[i][j-1]-ft[i-1][j-1])/(x[i]-x[i-j]); } } //累加计算牛顿多项式 b[n-1]=ft[n-1][n-1]; for(i=n-1;i=1;i--) { b[i-1]=ft[i-1][i-1]+b[i]*(t-x[i-1]); } fw=b[0]; printf(牛顿插值的结果为%f\n,fw); } 五、总结 西华大学计算机系学生上机报告 5 页 共 5 页 第 1 页 共