数值分析上机实验——数值积分.docVIP

实 验 报 告 课程名称 数值分析 实验项目名称 数值积分 实验类型 上机 实验学时 2 班级学号 2011113130 姓名 张振 指导教师 沈艳 实验室名称 理学楼407 实验时间 2013.11.18 实验成绩 预习部分 实验过程 表现 实验报告 部分 总成绩 教师签字 日期 哈尔滨工程大学教务处 制 实验三 数值积分 一．数值积分的基本思想 1.复合梯形公式：Tn=2； 2.复合辛普森公式：Sn=[f(a)+f(b)+2+4]； 以上两种算法都是将a-b之间分成多个小区间（n）T2n-Tn （2）Cn=S2n-Sn （3）Rn=C2n-Cn4 T =T - T，k = 1，2，… 二.实验题目及实验目的 (第4章计算实习题第1题)用不同数值方法计算积分 = -。 （1）取不同的步长h。分别用复合梯形及复合辛普森求积计算积分，给出误差中关于h的函数，并与积分精确值比较两个公式的精度，是否存在一个最小的h，使得精度不能再被改善？ （2）用龙贝格求积计算完成问题（1）。 （3）用自适应辛普森积分，使其精度达到10。 三．实验手段： 指操作环境和平台: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-1); S=S+sqrt(x).*log(x); end R(J+1,1)=R(J,1)/2+h*S; M=2*M; for k=1:J R(J+1,k+1)=R(J+1,k)+(R(J+1,k)-R(J,k))/(4^k-1); end err=abs(R(J+1,J)-R(J+1,J+1)); end q=R(J+1,J+1); 控制台输入代码： （1） a=0; b=1; h=0.1; t=TiXing_quad(a,b,h) s=Simpson_quad(a,b,h) h=0.01; t=TiXing_quad(a,b,h) s=Simpson_quad(a,b,h) h=0.001; t=TiXing_quad(a,b,h) s=Simpson_quad(a,b,h) （2） a=0; b=1; eps=10^-8; [quad,R]=Romberg(a,b,eps) （3） a=0; b=1; eps=10^-4; q=ZiShiYingSimpson(sqrt(x).*log(x),a,b,eps) 实验结果比较与分析 （1） h = 0.1时 h = 0.01时 h = 0.001时 由结果（1）可知对于同一步长h，复合辛普森法求积分精度明显比复合梯形法求积的精度要高，且当步长取不同值时即h越小时，积分精度越高。实验结果说明不存在一个最小的h,使得精度不能再被改善。又两个相应的关于h的误差（余项）Rn(f)=-h2f’’(η)；Rn(f)=-(h/2)4f(4)( η),其中η属于a到b。可知h愈小，余项愈小，从而积分精度越高。 （2） （注：看不清的话附有图片文件，可放大） 求的积分q = -0.444291362290625 （3） 求得积分q = -0.434745027462563 六.学习心得 对于同一步长h，复合辛普森法求积分精度明显比复合梯形法求积的精度要高，且当步长取不同值时即