计算方法实验六数值积分1.doc

《计算方法实验六数值积分1.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、实验六数值积分（1）复化梯形积分1、实验程序实现复化梯形积分的MATLAB函数文件agui_trapz.m在MATLAB命令窗口输入及实验结果及操作界面（2）复化辛普生求积公式1、实验程序实现复化辛普生求积公式的MATLAB函数文件agui_simpson.m在MATLAB命令窗口输入及实验结果及操作界面（3）龙贝格积分1、实验程序实现龙贝格积分的MATLAB函数文件agui_rbg.m在MATLAB命令窗口输入及实验结果及操作界面结果分析从上面三种方法，我们可以得出如下表格，方便进行比较：复化梯形复化辛普生龙贝格数值精确值计算结果7.389056127230227.389056126

2、214687.389056098931697.389056098930650误差范围2~3E-82~3E-81~2E-12等分数或二分次数n=7019n=24k=4复化梯形复化辛普生龙贝格数值精确值计算结果3.141592641093333.141592653588753.141592653638243.141592653589793误差范围1~2E-81~2E-124~5E-11等分数或二分次数n=3652n=29k=5在求和时，复化梯形公式、复化辛普生公式和龙贝格数值积分公式都有着较高的精度，其中龙贝格数值积分公式精度基本上是最高的。而在对积分区间作同样的分割的条件下，复合辛普生求

3、积公式比复合梯形公式的计算精度高。我们已知，数值求积的误差除了与被积函数有关之外，还与积分区间的长度（b-a）有关，积分区间越小，则求积公式的截断误差也越小，因此在求积分时，常把积分区间分成若干小区间，在每个区间上采用次数不高的求积公式，这就是基本的原理。在计算速度方面，从表中可看出，复化梯形公式的等分数要比其它两个大得多，且从计算结果上很容易知道复化辛普生公式也比复化梯形公式的收敛速度快得多。而龙贝格数值积分公式的计算量是最少的。由上可知，龙贝格数值积分公式在精度和计算速度上都是最高的。而使用龙贝格公式通过对梯形值进行外推加速的处理，能使精度快速提高。