实验2计算方法..docVIP

《计算方法》实验报告 学 院： 计算机学院 专 业： 计算机科学与技术 指导教师： 爨莹 班级学号： 201207010229 姓 名： 图尔荪托合提  实验二 插值法 一、实验目的： 1、掌握直接利用拉格郎日插值多项式计算函数在已知点的函数值；观察拉格郎日插值的龙格现象。 2、了解Hermite插值法原理，结合计算公式，确定函数值。 实验要求: 认真分析题目的条件和要求，复习相关的理论知识，选择适当的解决方案和算法； 2、编写上机实验程序，作好上机前的准备工作； 3、上机调试程序，并试算各种方案，记录计算的结果（包括必要的中间结果）； 4、分析和解释计算结果； 5、按照要求书写实验报告； 三、实验内容： 用拉格郎日插值公式、牛顿插值法确定函数值；对函数f(x)进行拉格郎日插值，并对f(x)与插值多项式的曲线作比较。 已知函数表：（0.56160,0.82741）、（0.56280,0.82659）、（0.56401,0.82577）、（0.56521,0.82495）用三次拉格朗日插值多项式、牛顿插值法求x=0.5635时函数近似值。 四、题目、原理、程序及结论： Lagrange插值 原理：1．Lagrange插值公式 输入，，，取次多项式作为插值函数，构造一组插值基函数 它满足，并称为拉格朗日插值基函数。 插值多项式为，其满足插值条件 ，并称此多项式为拉格朗日插值多项式。 特别地，当时称为线性插值，其插值多项式为 满足，从几何上看，为过两点的直线。 当时，称为抛物线插值，其插值多项式为 满足。从几何上看为过点和的一条抛物线。 程序： #includestdio.h #includeconio.h #includemalloc.h float Lagrange(float *x,float *y,float xx,int n) { int i,j; float *a,yy=0.0; a=(float *)malloc(n*sizeof(float)); for(i=0;i=n-1;i++) { a[i]=y[i]; for(j=0;j=n-1;j++) { if(j!=i) a[i]*=(xx-x[j])/(x[i]-x[j]); } yy+=a[i]; } free(a); return yy; } void main() { float x[4]={0.56160,0.56280,0.56401,0.56521}; float y[4]={0.82741,0.82659,0.82577,0.82495}; float xx=0.5635,yy; float Lagrange(float *,float *,float ,int); yy=Lagrange(x,y,xx,4); printf(x=%f,y=%f\n,xx,yy); getch(); } 运行结果： x=0.563500,y=0.826116 Newton插值 原理：Newton插值公式 输入，，，对，计算函数的各阶差商， 计算Newton插值函数值 程序： #includestdio.h void main() { float static x[4]={0.56160,0.56280,0.56401,0.56521}; float static y[4]={0.82741,0.82659,0.82577,0.82495}; int i,k; double c,p; for (k=1;k=3;k++) { printf(\n%dth is:,k); for(i=3;i=k;i--) { y[i]=(y[i]-y[i-1])/(x[i]-x[i-k]); printf(\n%8.6