复合梯形和复合辛普森MATLAB程序.docVIP

#### 实 验 报 告 课程名称 数值分析 实验项目名称 数值积分 实验类型 上机 实验学时 班级 学号 姓名 指导教师 实验室名称 实验时间 2014.11.19 实验成绩 预习部分 实验过程 表现 实验报告 部分 总成绩 教师签字 日期 实验三 数值积分 一．数值积分的基本思想 1.复合梯形公式：Tn=2； 2.复合辛普森公式：Sn=[f(a)+f(b)+2+4]； 以上两种算法都是将a-b之间分成多个小区间（n）,则h=(b-a)/n,xk=a+kh, xk+1/2=a+(k+1/2)h,利用梯形求积根据两公式便可。 3.龙贝格算法：在指定区间内将步长依次二分的过程中运用如下公式 （1）Sn=T2n-Tn （2）Cn=S2n-Sn （3）Rn=C2n-Cn4 T =T - T，k = 1，2，… 二、计算流程图 1、复合梯形和复合辛普森算法框图： 下图是龙贝格算法框图： 自适应辛普森积分算法流程框图： 二.实验题目及实验目的 实验题目：用不同数值方法计算积分 = -。 （1）取不同的步长h。分别用复合梯形及复合辛普森求积计算积分，给出误差中关于h的函数，并与积分精确值比较两个公式的精度，是否存在一个最小的h，使得精度不能再被改善？ （2）用龙贝格求积计算完成问题（1）。 （3）用自