数值分析实验报告2.pdfVIP

数值分析实验报告2 实验报告 一、实验名称 复合梯形求积公式、复合辛普森求积公式、龙贝格求积公式及自适应辛普森积 分。 二、实验目的及要求 1. 掌握复合梯形求积计算积分、复合辛普森求积计算积分、龙贝格求积计算积 分和自适应辛普森积分的基本思路和步骤. 2. 培养 Matlab 编程与上机调试能力. 三、实验环境 计算机，MATLAB 软件 四、实验内容 4 1.用不同数值方法计算积分 1 xlnxdx 。 0 9 （1）取不同的步长h。分别用复合梯形及复合辛普森求积计算积分，给出误 差中关于 h 的函数，并与积分精确指比较两个公式的精度，是否存在一个最小的 h，使得精度不能再被改善。 （2）用龙贝格求积计算完成问题（1）。 -4 （3）用自适应辛普森积分，使其精度达到 10 。 五、算法描述及实验步骤 1.复合梯形公式 将区间 [a,b]划分为 n 等份，分点 x =a+ah,h=(b-a)/h,k=0,1,...,n ，在每个 k 子区间 [x,x +1](k=0,1,...,n-1)上采用梯形公式（1.1），得 k k ba bf (x)dx [f (a)f (b)] a 2 （1.1） hn1 h n1   T  [f (x )f (x )] [f (b)2 f (x )f (b)] n 2 k k1 2 k k0 k1 （1.2） ba 2   R (f ) h f ( ), (a,b) n 12 （1.3） 其中 Tn 称为复合梯形公式，Rn 为复合梯形公式的余项。 2.复合辛普森求积公式 2 第 页 3 第 页 4 第 页 七、总结 通过本次学习Matlab，掌握了复合梯形求积公式、复合辛普森求积公式、龙贝格 求积公式及自适应辛普森积分的程序和算法，为以后处理数据提供一种更加简便， 准确的方法。 八、附录（源程序清单） 复合梯形 1. function s=fuhetixing(f,a,b,n) %f 为被积分函数 %a，b 是积分上下限 %n 是子区间个数 %s 是积分值 h=(b-a)/n; s=0; for k=1:(n-1) x=a+h*k; s=s+feval( f,x); end format long s=h*(feval(f,a)+feval(f,b))/2+h*s; 复合辛普森 2. function S=Comsimpson(f,a,b,n) %f 为被积分函数 %a，b 是积分上下限 %n 是子区间个数 %s 是积分值 h=(b-a)/(