插值法实验报告.docVIP

PAGE  PAGE 4 实验2 插值法 学院专 业日期学号姓 名成绩实验题目插值法指导老师 一．实验目的 掌握拉格朗日插值法、牛顿插值法。 记录运行结果，回答问题，完成实验报告。 二．实验环境 实验环境：C++ 实验内容 思考问题：插值多项式是否阶次越高越好?数值积分与插值的关系是什么? 1．用拉格朗日插值法求2的平方根。? 提示：可以用抛物线插值，f(1.69)=1.3，f(1.96)=1.4，f(2.25)=1.5。? 2．用牛顿插值法求2的平方根。?? 提示：可以用抛物线插值，f(1.69)=1.3，f(1.96)=1.4，f(2.25)=1.5。? 实验原理 一般，若有个数据点，上述近似标准要求 即：，则称近似函数为满足插值条件(3.2)的插值函数，而点称为插值节点。若满足插值条件(3.2)的函数的基函数取 ，则是多项式 称满足插值条件(3.2) 的多项式为插值多项式。 源程序 拉格朗日插值法 #include stdio.h #include stdlib.h #include math.h #define MAXSIZE 50 void input (double x[MAXSIZE],double y[MAXSIZE],long n); int main(void) { double x[MAXSIZE],y[MAXSIZE],_x,_y,t; long n,i,j; printf(\n 请输入插值节点的个数:); scanf(%ld,n); input(x,y,n); printf(\n 请输入插值点:); scanf(%lf,_x); _y=0; for(i=0;i=n-1;i++) { t=1; for(j=0;j=n-1;j++) if(j!=i) t*=(_x-x[j])/(x[i]-x[j]); _y+=t*y[i]; } printf(\n 插值点 (x,y)=(%lf,%lf).,_x,_y); system (pause); } void input(double x[MAXSIZE],double y[MAXSIZE],long n) { long i; for(i=0;i=n-1;i++) { printf(\n 请输入插值节点 x[%ld],y[%ld]:,i,i); scanf(%lf,%lf,x[i],y[i]); } } 牛顿插值法 #includestdio.h #includestdlib.h #includemath.h #define MAXSIZE 50 void input(double x[MAXSIZE],double f[MAXSIZE],long n); int main(void) { double x[MAXSIZE],f[MAXSIZE],_x,_y; long n,i,j; printf(\n请输入插值节点的个数：); scanf(%ld,n); input(x,f,n); printf(\n请输入插值点：); scanf(%lf,_x); for(j=1;j=n-1;j++) for(i=n-1;i=j;i--) f[i]=(f[i]-f[i-1])/(x[i]-x[i-j]); _y=f[n-1]; for(i=n-2;i=0;i--) _y=f[i]+(_x-x[i])*_y; printf(\n插值点(x,y)=(%lf,%lf).,_x,_y); system(pause); } void input(double x[MAXSIZE],double f[MAXSIZE],long n) { long i; for(i=0;i=n-1;i++) { printf(\n请输入插值节点x[%ld],y[%ld]:,i,i); scanf(%lf%lf,x[i],f[i]); } } 六．运行结果及问题答案 截图 问题答案:多项式是阶次越高越好. 七.实验结论 通过这次实验，让我了解到掌握拉格朗日插值法、掌握数值积分常用算法：逐 次分半梯形公式求积。通过运用C语言编译运行程序，发现其中有很