数值积分上机实验报告.docxVIP

数值积分上机实验报告 141110038 桂贤进题一：数学上已经证明了成立，所以可以通过数值积分来求的近似值。1.分别使用复合梯形、复合Simpson求积公式计算的近似值。选择不同的h,对每种求积公式，试将误差刻画为h的函数，并比较两方法的精度。是否存在某个值，当低于这个值之后，再继续减小h的值，计算精度不再有所改进，为什么？2.实现Romberg求积方法，并重复上面的计算；3.实现自适应积分方法，并重复上面的计算。解：1.1公式分析：(a) 对于复合梯形公式离散误差为：(b) 对于复合Simpson公式离散误差为： 1.2实现算法结果：分别利用梯形公式和Simpson公式计算结果如下：（下表中 此处为MATLAB中的数,可以认为具有足够大的精度）iT(f)S(f)113.0000000000000000.141621/23.1000000000000000.04163.1333333333333330.008341/43.1311764705882350.010423.1415686274509802.4026e-0561/63.1369630664712640.004633.1415917809360438.7265e-0781/83.1389884944910890.002603.1415925024587071.5113e-07101/103.1399259889071590.001673.1415926139392153.9651e-08121/123.1404352468468510.001163.1415926403053801.3284e-08201/203.1411759869541294.1667e-043.1415926529697856.2001e-10301/303.1414074684073301.8519e-043.1415926535353595.4434e-11401/403.1414884869236121.0417e-043.1415926535801059.6878e-12501/503.1415259869232546.6667e-053.1415926535872532.5402e-121001/1003.1415759869231291.6667e-053.1415926535897533.9968e-142001/2003.1415884869231304.1667e-063.1415926535897930从上表中可以看出：复合Simpson公式比复合梯形公式精度高，误差收敛的速度快不少。1.3误差下降速度对比：从上图可以看出，复合Simpson公式误差的收敛速度比复合梯形公式的误差的收敛速度快不少，下面验证收敛阶。1.4验证收敛阶：本实验的实际误差主要由离散误差和计算过程中的舍入误差组成，这里离散误差起主导作用，故理论上实际误差的收敛阶应该与离散误差的收敛阶相同。下面利用如下公式来计算实际的收敛阶，并与理论分析所得出的离散误差的收敛阶作比较。对上表格中所列的区间长度值，逐次利用相邻两个小区间长通过上述公式来计算收敛阶,并绘制成图形。得到图形如下：(a) 对复合梯形公式：由上面公式(2)可知，离散误差关于h为二阶收敛，同时由上图可知实验结果的收敛阶将近为2，故与理论分析相符。(b) 对复合Simpson公式：这里却有些奇怪，由上面公式(4)可知，离散误差理论上为4阶收敛，可实验结果却是将近6阶收敛。下面将进一步深入探究。探究如下：考虑这是由于被积函数f(x) 的特殊性导致，而不是由于Simpson公式离散误差真的能达到6阶收敛。由误差的余项公式考虑到f(x)=1.5计算过程中舍入误差的影响从上表格可以看出：当区间数n从1到200时，随着h的减小，实际误差在减小。考虑如下问题：若不断减小h的值，即不断增加区间数n，是否实际误差会一直减小?1.5.1理论分析：我们知道影响实验结果的精确度的因素主要有离散误差和舍入误差 ，而离散误差的大小可以通过离散误差的余项来体现。由公式（2）和（4），可以看出当不断增加区间数，即区间长度h不断减小时，离散误差会越来越小。但相应地由于计算精度的限制，当h不断减小时，舍入误差却会变大。故理论上会存在一个阈值H。当h大于H时，由于离散误差起主导作用，随着区间长度h的减小，实际误差会变小；而当h小于H 时，此时舍入误差将在计算结果的精确度起主导作用，再减小h,会导致计算结果的精确度基本保持不变甚至可能会有减小。1.5.2实验检验：（a）对于复合梯形公式：注意到误差收敛的速度较小，故我们首先选取区间数进行分析，得到下面图形。从图中可以看出，复合梯形公式的阈值H在区间数n为10^6到10^8之间取到。下面对n处于区间10^6到10^8进行分析。由上图可