数学实验3-插值与值积分-安振华-2012011837.docxVIP

实验3：插值与数值积分 化学工程系 分2 安振华 2012011837 【实验目的】 1．掌握用MATLAB 计算拉格朗日、分段线性、三次样条三种插值的方法，改变节点的数目，以及能对三种插值结果进行初步分析； 2．熟练掌握用MATLAB 及梯形公式、辛普森公式计算数值积分； 3．通过实例学习用插值和数值积分解决实际问题。 【实验内容】 预备： 1编制计算拉格朗日插值的M文件(lagr.m) 2分段线性插值函数（interp.m） 【实验三：习题10】 表3.10给出的x,y数据位于机翼断面的轮廓线上，y1与y2分别对应轮廓的上下线，假设需要得到x坐标每改变0.1时的y坐标。试完成加工所需数据，画出曲线，求机翼剖面的面积。 表3.10 机翼剖面轮廓线数据 x 0 3 5 7 9 11 12 13 14 15 y1 0 1.8 2.2 2.7 3.0 3.1 2.9 2.5 2.0 1.6 y2 0 1.2 1.7 2.0 2.1 2.0 1.8 1.2 1.0 1.6 【分析】 由于给定的数据点是有限的，要想确定更多的有效数据，就要运用插值方法。而加工断面的面积，则应通过数值积分分别求出上下轮廓线对x 轴围成的面积，然后作差求得。 【解答】 运用学习的三种插值方法逐一计算，具体如下： 1、拉格朗日多项式插值 MATLAB程序: 运行结果： 给出机翼剖面面积 S=40.3044 2、分段线性插值 MATLAB程序: 运行结果： 给出机翼剖面面积 S=10.7500 3、三次样条插值 MATLAB程序: 运行结果： 给出机翼剖面面积 S=11.3444 整合程序，将三种差值方法集成在同一程序，可以直观地比较各自得出的结果。 整合后的MATLAB程序: 运行结果： 【结果分析】 由图可知拉格朗日插值法得出图形与实际不符，不予考虑，可见拉格朗日插值方法的收敛性无法得到保证；分段线性插值的大致形状与原图相近，但它得到的图线不平滑，应用到机翼这样的部件上也是无法使用的。而三次样条插值比较光滑，误差最小。 由梯形公式得到加工剖面面积为11.3444。根据梯形公式和辛普森公式的应用特点（已知被积函数的情况下应用辛普森函数更方便，而对于此处插值得到的数据则用梯形公式更方便），故此处应用梯形公式。又由于在插值之后，步长较短，因此推测梯形公式也较为准确。 综上，此问题的最佳解决方案是采用三次样条插值方法，机翼剖面的面积为11.3444。 【实验三：习题11】 图3.11是欧洲一个国家的地图，为算出它的国土面积，首先对该国地图做出如下测量：以由西向东方向为x轴，由南到北方向为y轴，选择方便的原点，并将从最西边界点到最东边界点在x轴上的区间适当的划分为若干段，在每个分点的y方向测出南边界点和北边界点的y坐标y1和y2，得到表3.11（单位：mm）。 根据地图的比例尺我们知道18mm相当于40km，试由测量数据计算该国国土的近似面积，与它的精确值41288km2作比较。 图3.11 表3.11 地图边界点数据 x 7.0 10.5 13.0 17.5 34.0 40.5 44.5 48.0 56.0 61.0 68.5 76.5 80.5 91.0 y1 44 45 47 50 50 38 30 30 34 36 34 41 45 46 y2 44 59 70 72 93 100 110 110 110 117 118 116 118 118 x 96.0 101.0 104.0 106.5 111.5 118.0 123.5 136.5 142.0 146.0 150.0 157.0 158.0 y1 43 37 33 28 32 65 55 54 52 50 66 66 68 y2 121 124 121 121 121 122 116 83 81 82 86 85 68 【分析】 由3.10解题思路知，可以先用三种插值计算方法分别进行插值，得到图像后分析选择与实际符合最好的方案进行计算。计算时可以选择辛普森积分或梯形积分方式进行积分运算。将结果与精确值进行比较和误差分析即可 【解答】 运用学习的三种插值方法分别计算，具体如下： MATLAB程序: 运行结果： 由以上图像分析可得，拉格朗日插值所得结果与实际偏差非常大，顾舍弃该方法。选择三次样条插值或分线段性插值法分析计算。 MATLAB程序: 运行结果： S1 = 42407 S2 = 42458 【结果分析】 这两种方法得到的最终结果都很相近，和精确值41288km2虽然也有一定误差，但并不是很大，分别为2.71%和2.83%，在误差可以接受的