《数值分析实验报告(插值法).docx》由会员分享，可在线阅读，更多相关《数值分析实验报告(插值法).docx（9页珍藏版）》请在三一办公上搜索。
1、学生学号实验课成绩武汉理工大学学生实验报告书实验课程名称数值分析开课学院计算机科学与技术学院指导老师姓名学生姓名 学生专业班级2010-2010学年第一学期实验课程名称： 数值分析实验项目名称实验成绩实验者专业班级组别同组者实验日期年 月曰第一部分：实验分析与设计(可加页)一、实验内容描述(问题域描述)1、分别画出Lagrange插值公式、Newton插值公式、分段插值公式和Hermite插值公式的算 法流程图2、分别用Lagrange插值公式和Newton插值公式通过编程计算函数f(x)的近似值已知对于f(x)-ex,有数据表如下：x00.51.02.0f(x.)1.000001.64872
2、2.718287.38906(1) 对x0-0,X-0.5利用线性插值计算f(0.25)的近似值；对x0-0.5,X-1利用线性插值计算f(0.75)的近似值；1(2) 对x0-0,xi-0.5,x2-2利用二次插值计算f(0.25)和f(0.75)的近似值(3) 对 x0-0,xi-0.5,x2-2 求 f(x)的 Hermite 插值多项式 H5 (x)；(4) 分析和比较各插值算法的精度差异3、通过编程计算函数f(x)的近似值。已知对于f(x)-，有数据表如下：x2.02.12.22.32.4f(x)1.4142141.4491381.4833401.5165751.549193(1)
3、计算各阶插值多项式在不同点的值：f(2.05),f(2.15),f(2.45)；(2) 利用分段线性插直和分段抛物插值计算(1)中的函数值；(3) 分析和比较算法的效率差异和精度差异(同时注意插值点的位置与精度之间的关系)。4、用不同万式万法编程给出计算Langrange插值和Newton插值的算法，分析和比较两种算法 的编程难易以及算法的效率差异总计算量之间的关系。5、写出实习报告二、实验基本原理与设计(包括实验方案设计，实验手段的确定，试验步骤等，用硬件逻辑或者算法描述)【拉格朗日插值法算法流程图】【拉格朗日插值法源程序】#includeusing namespace std;int ma
4、in()cout请输入坐标点个数：count;double point1002;int count1=0;cout请输入坐标：endl;while(count1pointcount10pointcount11;count1+;cout计算 f(x)请输入 x： x;double f=0,Lu=1;for(int i=0;icount;i+)for(int j=0;jcount;j+)if(j=i)continue;Lu=Lu*(x-pointj0)/(pointi0-pointj0);Lu=Lu*pointi1;f=f+Lu;Lu=1;coutf(x)的值为：fendl;return 0;【牛
5、顿插值法源程序】#includeusing namespace std;int main()cout请输入坐标点个数：count;double point1002;int count1=0;cout请输入坐标：endl;while(count1pointcount10pointcount11;count1+;cout计算 f(x)请输入 x： x;double d100;for(int i=0;icount;i+)di=pointi1;for(int j=1;j=j;i)di = (di-di-1)/(pointi0-pointi-j0);double f=d0,Lu=1,L;for(i=1;
6、icount;i+)Lu=Lu*(x-pointi-10);L=Lu*di;f=f+L;coutf(x)的值为：fendl;return 0;【埃米尔特插值法源程序】#includeusing namespace std;struct pointdouble x;double y;double d;point100;int main()cout请输入坐标点个数：count;int count1=0;cout请依次输入坐标的横纵坐标和对应的导数：endl;while(count1pointcount1.xpointcount1.ypointcount1.d;count1+;cout计算 f(x)
7、请输入 x： x;double f=0,Lu=1,Laa=0,La,Lb;for(int i=0;icount;i+)for(int j=0;jcount;j+)if(j=i)continue;Lu=Lu*(x-pointj.x)/(pointi.x-pointj.x);Laa=Laa+1.0/(pointi.x-pointj.x);La=pointi.y*(1-2*(x-pointi.x)*Laa)*Lu*Lu;Lb=pointi.d*(x-pointi.x)*Lu*Lu;f=f+La+Lb;Lu=1;Laa=0;coutf(x)的值为：fendl;return 0;三、主要仪器设备及耗材1
8、. PC 机2. 开发环境(比如：VC, Eclipse)第二部分：实验调试与结果分析请输入坐标点个数:拉格朗日插值法J请输入坐标：01.000000.51.64872肝算*0请输人*.25 的值为:1-32436请输入坐标点个数：肩输入坐标：.51.64872.02.71828计算请输入X：0.75“X)的值为:2.1835八上巳T.ILE牛顿插值法H育输入坐标：6.51.64872阡算x请输又c：.250的值为:1.32436兰三 T.JAV情输入坐标：0.51.648721.02.71828计算f奴请输入X：0.75的值为:2.1835(可加页)-、调试过程(包括调试方法描述、实验数据记
9、录，实验现象记录，实验过程发现的问题等)(1)用拉格朗日插值法计算时，输入及运行结果如下：(2)利用二次插值计算时，输入及运行结果如下:拉格朗日插值法日请输入坐标点个数：底输入坐标：1.648727.38906计算“X请输入X：.253的值为:1-24531请输入坐标点个数：备输入坐标：010.51.6487227.38906计算Ft X)请输入x：0.75F3的值为：牛顿插值法请输入坐标点个数：蓄输入坐标：010.51.6487227.38906计算fX请输入X：0.25FB的值为：1.24531请输入坐标点个数：蓄输入坐标：0.51.6487227.38906计算f tx请输入X：0.75
10、的值为:2.21022(3) 用艾尔米特插值法计算时,f(x)的插值多项式H5(x) = (1+4*x)*(x-0.5)*(x-0.5)*(x-2)*(x-2) + (3.90807-6.03838*x)*(x-2)*(x-2)*x*x+ (2.34573-4.16674*x)*x*x*(x-0.5)*(x-0.5)(4) 各插值算法的精度差异比较经过比较，拉格朗日插值法要比牛顿插值法算法的计算量多一些，拉格朗日插值法后一次计算时用到 了前一次计算的结果，提高了运算的效率，但拉格朗日插值法在构造艾尔米特插值法时很方便，将坐标点 和对应的导数结合起来的精度比线性插值的精度又要高一些。但从实验数据
11、来看，在坐标不是很多的情况下，已知的点越多精度也就相对较高。对于实验要求的第二组数据用拉格朗日插值法（或者牛顿插值法）实验结果如下: 一下分别是二阶、请输入坐标点个数:情输入坐标：R1.414214E.11.449138仕算g请输入x：的值为：1.43168“X请输入X：备X的值为：1.4666请输入x：的值为：1.57137三阶、四阶、五阶插值得到的结果的值为：1-56569有输入坐标：21.4142142.1 1.4491382.2 1.483340计算代X请输入X：2.05F3的值为：1.43177计算ftx请输入X：2.15“X的值为：1-46633摩没篇入X：2.45以上只是实验结果
12、的一部分，改变插值的位置时，得到的实验结果精度分析可知，插值次数并不是越多越好，多了反而会让结果更加偏离真实结果，这充分说廿在“病态性质”，在已知点很多的情况下应该采用分段低次插值，将拉格朗日插值法和牛顿插值法运用到 分段低次插值法当中，这样得到的结果可能胡更加精确。对于分段低次插值本实验没有给出实验结果，但从实践上来看，分段低次插值的精度要比线性插值精 度高，但当插值阶数比较少的时候没有必要采用分段低次插值。果 高次插值存二、实验小结、建议及体会各种插值法都有自己的利与弊，拉格朗日插值法运算过程相对复杂，但当和导数结合起来，组成抛物 插值的时候，精度就可以提高很多。牛顿插值法、拉格朗日插值法等线性插值法只能适合在已知点不多的 情况下使用，当已知的坐标点很多时候应该将区间分成小段进行分段线性插值或者分段抛物插值。采用组合的思想是数值分析常用的思想和技巧之一，单个的方法得到的结果虽然不是很理想，但将多 个方法按照某种方式结合在一起就能改进实验方法，我们应该触类旁通，在以后的学习中学会使用这种思 想。
