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

课程名称 课程名称 数值分析 实验项目名称 数值积分 #### 实 验 报 告 实验类型 班级姓名 上机 实验学时 学 号 指导教师 实验室名称 实验时间 2014.11.19 实验成绩 预习部分 实验过程表现 实验报告部分 总成绩 教师签字 日期 实验三 数值积分 一．数值积分的基本思想 1.复合梯形公式：Tn= h [ f (a) ? f (b) ? 2 ?n?1 2  f (xk )] ； k ?1 h Sn= [f(a)+f(b)+2 2. 2. ?n?1复合辛普森公式： f (xk )] +4 ?n?1 f (x ? 1/ 2) ]； k ?1 k ?0 以上两种算法都是将 a-b 之间分成多个小区间（n）,则 h=(b-a)/n,x =a+kh, k x =a+(k+1/2)h,利用梯形求积根据两公式便可。 k+1/2 3.龙贝格算法：在指定区间内将步长依次二分的过程中运用如下公式 （1） 4 1 Sn= T2n- Tn 3 3 （2）Cn= 16 S2n- 1 Sn 15 15 （3）Rn= 64 C2n- 1 Cn4 T ( k ) 63 63 4m 1 = T ( k ?1) - T ( k ) ，k = 1 ，2，… m 4m ? 1 m ?1 4m ?1 m?1 二、计算流程图 1、复合梯形和复合辛普森算法框图： 韬人a, b.. H开始韬入a, 韬人 a, b.. H 开始 韬入 a, b R h＝店－屯）／N ＄二E Ca)＋f 心）J x二 邑 h= (b- 之)／N S=f (a)tf Ch), x=a x立 t h., x立 t h.,s二St 2f Gr) n=n.+1 S=Sh/2 ::-::=x＋应 2J S =S+4f 旬 ；x；；对应 2$ 5 +2 f 饮］ n=t1.+1 S=S凶6 辐出＄豁出 辐出＄ 豁出5 巴 三 茎台梯秏流程图 莹合辛普森流程图 下图是龙贝格算法框图： 开始 开始 读入 a,b,c h=b-a,T1=h[f(a)+f(b)/2,k=1 S=0,x=a+h/2 S=S+f(x) x=x+h Y xb ? N T2=T1/2+hS/2 S2=T2+(T2-T1)/3 Y k=1? k=k+1,h=h/2 T2=T1,S2=S1 N C2=S2+(S2-S1)/15 Y k=2 ? C2=C1 N R2=C2+(C2-C1)/63 Y k=3 ? R2=R1 N N │R2-R1│ε Y 输出R2 结束 自适应辛普森积分算法流程框图： 二.实验题目及实验目的 x实验题目：用不同数值方法计算积分? 1 x 0  ln xdx  = - 4 。 9 开始定义函数F输入A,B,NH=（B-A） 开始 定义函数F 输入A,B,N H=（B-A）/（2*N） S=0.5*（F（A）- F（B），调用函数 F I=1,N I ? ? S=S+2*F[A+（2*I-1*H）]+（F（A+2*I*H），调用函数F S=（B-A）/（3*N）S 输出 S 结束 用龙贝格求积计算完成问题（1）。 用自适应辛普森积分，使其精度达到10 4 。 实验目的： 1、了解并掌握 matlab 软件的基本编程、操作方法； 2、初步了解 matlab 中的部分函数，熟悉循环语句的使用； 3、通过上机进一步领悟用复合梯形、复合辛普森公式，以及用龙贝格求积方法计算积分的原理。 三．实验手段： 指操作环境和平台:win7系统下 MATLAB R2009a 程序语言：一种类似 C 语言的程序语言，但比 C 语言要宽松得多，非常方便。 四.程序 ①复合梯形求积程序 function t=TiXing_quad(a,b,.h) format long x=a:h:b; y=sqrt(x).*log(x); y(1)=0; t=0; for k=1:(b-a)/h, t= t+y(k)+y(k+1); end t=t*h/2; ②复合辛普森求积程序 function s=Simpson_quad(a,b,h) format long x=a:h:b; y=sqrt(x).*log(x); z=sqrt(x+h/2).*log(x+h/2); y(1)=0; s=0; for k=1:(b-a)/h, s= s+y(k)+y(k+1)+4*z(k); end s=s*h./6; ③龙贝格求积程序 function [q,R]=Romberg(a,b,eps) h=b-a; R(1,1)=h*(0+sqrt(b).*log(b))/2; M=1; J=0; err=1; while erreps J=J+1; h=h/2; S=0; for p=1:M x=a+h*(2*p-