matlab for求积分,MATLAB数值分析之数值积分（一）

【实验原理】

一、复合辛普森方法

1)复合辛普森方法基本思路

把积分区间分成若干子区间，再在每个子区间上用低阶辛普森求积公式。

2)复合辛普森算法描述

二、龙贝格方法

1)基本思路

实质上是梯形公式的递推化，是在梯形公式、辛普森公式和柯特斯公式关系的基础上构造出的加速计算的方法，是一种外推算法。

2)算法描述

【实验内容】

一、回答下面的问题

1.什么是复合求积公式?写出事前误差估计和事后误差估计的公式。

2.当积分区间是奇数等份时，如何搭配使用n=2和n=3的辛普森公式，举例说明。

3.试说明高斯-勒让德积分、牛顿-柯特斯的区别和联系。

4.时间测试下，请对自己编写的龙贝格算法的执行效率进行分析，说明对分次数和计算执行时间的关系？

二、编程计算课后习题中规定题目，交回实验报告与计算

1.用不同数值方法计算积分

(1)取不同的步长h，分别用复合梯形及复合辛普森求积计算积分，给出误差中关于h的函数，并与积分精确值比较两个公式的精度，是否存在