数值分析完整版.doc

PAGE PAGE 1 实验一 Lagrange插值和Newton插值 一、实验目的及要求 1．掌握利用Lagrange插值法及Newton插值法求函数值并编程实现。 2．程序具有一定的通用性，程序运行时先输入节点的个数n，然后输入各节点的值（），最后输入要求的自变量x的值，输出对应的函数值。 二、实验设备和实验环境 安装有C、C++或MATLAB的计算机。 三、算法描述 插值的基本原理（求解插值问题的基本思路） 构造一个函数y=f(x)通过全部节点，即 (i=0、1、… n) 再用f(x)计算插值，即 拉格朗日（Lagrange）多项式插值 Lagrange插值多项式： 3．牛顿（Newton）插值公式 五、实验内容 1．给定，构造Lagrange插值函数计算。 2．已知4个点的函数值如下表，用Newton插值法求x=0.596时的函数值。 i 0 1 2 3 xi 0.40 0.55 0.65 0.80 yi 0.41075 0.57815 0.69675 0.88811 六、实验源程序 Lagrange插值源程序： #includeiostream using namespace std; double Lagrange(double x[20],double y[20],int n,double x1) { double y1=0; for(int k=0;kn;k++) { double l=1; for(int j=0;jn;j++) { if(k!=j) l=((x1-x[j])/(x[k]-x[j]))*l; } y1=y[k]*l+y1; } return y1; } int main() { int m; double x; cout=====拉格朗日插值=====endl; while(1) { cout请输入要插值的数据数量:; cinm; if (m=1) { cout输入数值必须大于1!\n请重新输入!endl; continue; } double a[20]; double b[20]; for (int i=0;im;i++) { cout请输入第i+1个x点:; cina[i]; cout请输入第i+1个y点:; cinb[i]; } cout请输入插值点:x=; cinx; cout结果输出插值点y=Lagrange(a,b,m,x)endl; } return 0; } 牛顿插值源程序： #includeiostream using namespace std; double Newton(double x[20],double y[20],int n,double x1) { double y1; double s[20]; for(int i=0;in;i++) s[i]=y[i]; for(int k=0;kn-1;k++) { for(int j=n-1;jk;j--) { s[j]=(s[j]-s[j-1])/(x[j]-x[j-1-k]); } double p=1; y1=s[0]; for(int h=0;hn-1;h++) { p=p*(x1-x[h]); y1=y1+p*s[h+1]; } } return y1; } int main() { int m; double x; cout=====牛顿插值=====endl; while(1) { cout请输入要插值的数据数量:; cinm; if (m=1) { cout输入数值必须大于1!\n请重新输入!endl; continue; } double a[20]; double b[20]; for (int i=0;im;i++) { cout请输入第i+1个x点:; cina[i]; cout请输入第i+1个y点:; cinb[i]; } cout请输入插值点:x=; cinx; cout结果输出插值点y=Newton(a,b,m,x)endl; return 0; } } 七、计算结果 1．给定，构造Lagrange插值函数计算。 程序运行结果图： 2．已知4个点的函数值如下表，用Newton插值法求x=0.596时的函数值。 i 0 1 2 3 xi 0.40 0.55 0.65 0.80 yi 0.41075 0.57815 0.69675 0.88811 程序运行结果图： 八、实