华科计算方法实验报告.docVIP

课 程 实 验 报 告 课程名称： 计算方法 计算机科学与技术学院 目 录 实验一 Lagrange插值 1实验内容与要求 3 2 公式与算法描述 3 3 程序源码 4 4 实验结果 5 5 实验小结 5 实验二 牛顿插值 1实验内容与要求 5 2 公式与算法描述 6 3 程序源码 6 4 实验结果 7 5 实验小结 7 实验三 复化simpson公式 1实验内容与要求 8 2 公式与算法描述 8 3 程序源码 8 4 实验结果 10 5 实验小结 10 实验四 变步长梯形法 1实验内容与要求 10 2 公式与算法描述 11 3 程序源码 11 4 实验结果 13 5 实验小结 13 实验一、Lagrange插值 1.实验内容与要求 利用三次Lagrange插值多项式计算函数值，编写C程序实现其功能。 已知函数表如下： Xi 0.56160 0.56280 0.56401 0.56521 Yi 0.82741 0.82659 0.82577 0.82495 试用lagrange插值多项式求 X = 0.5635时的函数值。 公式与算法描述 Lagrange插值基函数 Lagrange插值公式 ↓ ↓ ↓ ↓ 3 程序源码 #include iostream #include conio.h #include malloc.h float lagrange(float *x,float *y,float xx,int n) /*拉格朗日插值算法*/ { int i,j; float *a,yy=0.0; /*a作为临时变量，记录拉格朗日插值多项式*/ 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; } int main() { int i; int n; float x[20],y[20],xx,yy; printf(Input n:); scanf(%d,n); if(n=20) { printf(Error!The value of n must in (0,20).); getch();return 1; } if(n=0) { printf(Error! The value of n must in (0,20).); getch(); return 1; } for(i=0;i=n-1;i++) { printf(x[%d]:,i); scanf(%f,x[i]); } printf(\n); for(i=0;i=n-1;i++) { printf(y[%d]:,i);scanf(%f,y[i]); } printf(\n); printf(Input xx:); scanf(%f,xx); yy=lagrange(x,y,xx,n); printf(x=%f,y=%f\n,xx,yy); getch(); } 实验结果 实验小结 本次实验，我用C语言实现了Langrange插值。本实验算法简单，只是有很多地方的细节需要注意，比如变量类型要用float，如果使用double型，则对于本实验的数据来说太大了。本次实验简单，也很快解决了