（全新）数值分析第五版上机实验答案实验一~实验六.docVIP

实验一 Lagrange插值算法 实验目的：掌握拉格朗日（Lagrange）插值算法的基本原理，理解插值基函数的性质，掌握基本的误差概念。学习用计算机语言编写程序实现算法。 [参考程序] #include stdio.h //定义插值节点及所求点数据，根据题目不同而修改 double x[]={0.32,0.34,0.36}; double y[]={0.314567,0.333487,0.352274}; double xx=0.3367; // Lagrange插值算法函数，利用循环计算具有对称性的基函数和最终结果 double Lagrange(double xxx,int n) { int i; double result=0,temp; for(i=0;i=n;i++) { temp=1; for(int j=0;j= n;j++) { if(j!=i) { temp=temp*(xxx-x[j])/(x[i]-x[j]); } } result=result+temp*y[i]; } return result; } void main() { int n; printf(Please input n:); scanf(%d,n); printf(Sin(%f) = %f \n,xx,Lagrange(xx,n)); } 实验二 Newton均差插值算法 实验目的：掌握Newton均差插值算法的基本原理，理解均差的概念，掌握均差表的计算方法。学习用计算机语言编写程序实现算法。 [参考程序] #include stdio.h #define N 10 double f[N][N]; //定义插值节点及所求点数据，根据题目不同而修改 double x[]={0.4,0.55,0.65,0.80,0.90,1.05}; double y[]={0.41075,0.57815,0.69675,0.88811,1.02652,1.25382}; double fx(int i,int j); double S(int start,int end,double xx); main() { int loopi,loopj,n; double result,xx; scanf(%d,n); scanf(%lf,xx); for(loopi=0;loopi=n;loopi++) { //零阶均差作为均差表二维数组的第0列 f[loopi][0]=y[loopi]; } //循环计算均差表中的所有数据 for(loopi=1;loopi=n;loopi++) { for(loopj=1;loopj=loopi;loopj++) { f[loopi][loopj]=fx(loopi,loopj); } } result=S(0,n,xx); printf(Result is: %f,result); return 1; } //求均差的函数 double fx(int i,int j) { if(j==0) { return f[i][j]; } else { //这种表示方法需要注意两个x的下标 return (fx(i,j-1)-fx(i-1,j-1))/(x[i]-x[i-j]); } } //用秦九韶算法计算插值多项式结果 double S(int start,int end,double xx) { if(start==end) { return f[end][end