数值分析实验报告包含源程序.docVIP

课程实验报告 课程名称 计算方法 班级 实验日期 姓名 学号 实验成绩 实验名称 实验一 插值多项式 实 验 目 的 及 要 求 掌握拉格郎日插值多项式的用法，适用范围及精确度。 掌握牛顿插值多项式的用法，适用范围及精确度。 实 验 环 境 Microsoft Visual C++ 6.0 实 验 内 容 1、Lagrange插值多项式：，其中 2、Newton插值多项式： 3、计算实例以2.0,2.1,2.2,2.3,2.4为节点,用牛顿插值求函数在2.51处的近似值； 取正弦函数 4、如果有错，修改直至运行成功，查看运行结果 算 法 描 述 及 实 验 步 骤 1、实验1实验步骤： （1）程序设计 （2）计算实例 （3）撰写实验报告 2、实验2实验步骤 （1）程序设计 （2）计算实例以2.0,2.1,2.2,2.3,2.4为节点,用牛顿插值求函数 在2.51处的近似值； （3）撰写实验报告 调 试 过 程 及 实 验 结 果 拉格朗日插值实验结果如图： 以0.32,0.34,0.36为节点,分别用线性插值和抛物插值求正弦函数在0.3367处的近似值 牛顿插值实验结果如图： 进行测试： 使用牛顿插值计算实例以2.0,2.1,2.2,2.3,2.4为节点,用牛顿插值求函数 在2.51处的近似值： 总 结 通过本次实验，使我进一步加深了对拉格朗日插值多项式的了解，掌握拉格朗日插值多项式的用法，能够正确的运用拉格朗日插值多项式进行运算。 通过本次实验，使我加深了对牛顿插值的认识，能够熟练的运用公式进行运算，掌握牛顿插值多项式的用法。 附 录 拉格朗日插值程序清单： #includestdio.h #includeconio.h double Lagrange(double x,int n,double a[],double b[]); void main() { double a[20],b[20],x,y; int i=0,n; printf(请输入数值个数：n=); scanf(%d,n); printf(请输入%d个不同的x的数值：\n,n); for(i=0;in;i++) { scanf(%lf,a[i]); } printf(请按顺序输入x对应的y的数值：\n); for(i=0;in;i++) { scanf(%lf,b[i]); } /* for(i=0;in;i++) { printf(%f ,a[i]); }*/ printf(请输入需计算的数：x=); scanf(%lf,x); y=Lagrange(x,n,a,b); printf(根据拉格朗日插值公式计算后的数为：y=%f\n,y); } double Lagrange(double x,int n,double a[],double b[]) { double y=0,L; int i,j; for(i=0;in;i++) { L=1; for(j=0;jn;j++) { if(j==i)continue; else { L=L*(x-a[j])/(a[i]-a[j]); } } y+=L*b[i]; } return y; } 牛顿插值程序清单： #includestdio.h #includestdlib.h #includeiostream.h double Newton(double x,int n,double a[],double b[]); double Chashang(int a1,int b1,double a[],double b[]); void main() { double a[20],b[20],x,y; int i=0,n; printf(请输入数值个数：n=); scanf(%d,n); printf(请输入%d个不同的x的数值：\n,n); for(i=0;in;i++) { scanf(%lf,a[i]); } printf(请按顺序输入x对应的y的数值：\n); for(i=0;in;i++) { scanf(%lf,b[i]); } printf(请输入需计算的数：x=); scanf(%lf,x); y=Newton(x,n,a,b); printf(根据牛顿插值公式计算后的数为：y=%f\n,y); } double Newton(double x,int n,double a[],double b[]) { int i; double t=1.0,zhi=b[0],p=0.0;