数学实验——插与数值积分.docxVIP

实验3 插值与数值积分 分1 黄浩 2011011743 实验目的 掌握用matlab计算Lagrange、分段线性、分段三次和三次样条插值这四种插值方法，并通过改变节点的数目，分析插值结果 掌握用matlab及梯形公式、辛普森公式计算数值积分 通过实例学习用插值和数值积分解决实际问题 实验内容 1.（1）题： 问题叙述： 考虑函数：fx=11+x2 方案实施： 由于在matlab中，没有预设的lagrange插值法函数，因此我首先编制了Lagrange（程序见四.1）和分段线性的插值函数（程序见四.2），并进行了y=x2函数的插值检验（程序见四.3） 上图中，红、绿两线分别是自设Lagrange和分段线性的插值图，蓝色十字点是matlab自带的分段线性函数的插值图，蓝绿完全吻合、红绿基本吻合，证明自设函数是正确的。 然后进行本题目要求的插值操作：以（x0,y0）为节点，节点在区间内均匀分布且将节点数分别调为6、11、21；x为插值点，绘图时使用101个插值点（始终大于节点数），制表时使用的插值点数目依节点数调节；y为真值，y1、y2、y3、y4分别为使用Lagrange、分段线性、分段三次和三次样条插值的方法所得的插值结果。 节点数为6时：（即节点为-5,-3,…,3,5） x y y1 y2 y3 y4 0 1.0000 0.5673 0.5000 0.5000 0.5685 0.5000 0.8000 0.5501 0.5000 0.5000 0.5513 1.0000 0.5000 0.5000 0.5000 0.5000 0.5000 1.5000 0.3077 0.4213 0.4000 0.4425 0.4167 2.0000 0.2000 0.3212 0.3000 0.3133 0.3131 2.5000 0.1379 0.2097 0.2000 0.1775 0.2029 3.0000 0.1000 0.1000 0.1000 0.1000 0.1000 3.5000 0.0755 0.0078 0.0846 0.0754 0.0183 4.0000 0.0588 -0.0481 0.0692 0.0559 -0.0285 4.5000 0.0471 -0.0460 0.0538 0.0431 -0.0263 5.0000 0.0385 0.0385 0.0385 0.0385 0.0385 图像为： 由上图可见，当节点数为6时，插值效果皆不理想，在各节点之间的插值都有很大的误差，在x=0附近误差最大 节点数为11时：（即节点为-5,-4,…,4,5） x y y1 y2 y3 y4 0 1.0000 1.0000 1.0000 1.0000 1.0000 0.5000 0.8000 0.8434 0.7500 0.7969 0.8205 1.0000 0.5000 0.5000 0.5000 0.5000 0.5000 1.5000 0.3077 0.2353 0.3500 0.3219 0.2973 2.0000 0.2000 0.2000 0.2000 0.2000 0.2000 2.5000 0.1379 0.2538 0.1500 0.1385 0.1401 3.0000 0.1000 0.1000 0.1000 0.1000 0.1000 3.5000 0.0755 -0.2262 0.0794 0.0755 0.0745 4.0000 0.0588 0.0588 0.0588 0.0588 0.0588 4.5000 0.0471 1.5787 0.0486 0.0465 0.0484 5.0000 0.0385 0.0385 0.0385 0.0385 0.0385 图像为： 由上图可见，当节点数为11时，分段三次与三次样条插值的效果较好，其次是分段线性插值，而Lagrange插值在[-2,2]区间内精度较高，但在|x|3时，就会发生Runge振荡 节点数为21时（程序见四.4）：（即节点为-5,-4.5…,4.5,5） x y y1 y2 y3 y4 0 1.0000 1.0000 1.0000 1.0000 1.0000 0.1000 0.9901 0.9904 0.9600 0.9869 0.9891 0.2000 0.9615 0.9626 0.9200 0.9526 0.9594 0.3000 0.9174 0.9189 0.8800 0.9050 0.9152 0.4000 0.8621 0.8632 0.8400 0.8515 0.8606 0.5000 0.8000 0.8000 0.8000 0.8000 0.8000 0.6000 0.7353 0.733