Romberg数值积分.docxVIP

Romberg数值积分实验报告?一、实验目的1、加深外推法的原理理解, 掌握Romberg外推法的计算方法。2、用matlab软件实现Romberg数值积分来计算题目的运算。二、基本理论及背景1、理论推导：对区间[a， b]，令h=b-a构造梯形值序列{T2K}。 　　T1=h[f(a)+f(b)]/2 　　把区间二等分，每个小区间长度为 h/2=(b-a)/2，于是 　　T2 =T1/2+[h/2]f(a+h/2） 　　把区间四(2)等分，每个小区间长度为h/2 =（b-a）/4，于是 　　T4 =T2/2+[h/2][f（a+h/4）+f（a+3h/4)..................... 　　把[a，b] 2等分，分点xi=a+（b-a）/ 2 ·i （i =0，1，2 · · · 2k）每个小区间长度为（b-a）/ 2 .2、参考Romberg数值积分，实现积分的数值求解，完成下列题目：??三、算法设计及实现1、算法设计(a)function y=fun1(x)y=sqrt(4-(sin(x)).^2)；(b)function y=fun2 (x)y=log(1+2*x)/3;(c)function y=fun3(x)y=exp(-x).*sin(x.^2)四、实验步骤1、打开matlab软件，新建Romberg.m文件,在窗口中编辑Romber数值积分函数程序代码,并保存在指定的文件夹下,在Current Directory窗口右边点击《Browse For Folder》按钮指向Romberg.m文件；在Command Window中编辑相应要计算的题目的数值函数及相应的题目的表达式。2、输出结果和初步分析说明（见附件一）。五、使用说明实验结果分析1、在Command Window窗口中编辑要调用的函数名与指定的函数名字不同导致出现错误，通过改正与函数名相同即可。2、Romberg方法也称为逐次分半加速法。它是在梯形公式、辛卜生公式和柯特斯公式之间的关系的基础上，构造出一种加速计算积分的方法。作为一种外推算法，它在不增加计算量的前提下提高了误差的精度在等距基点的情况下，用计算机计算积分值通常都采用把区间逐次分半的方法进行,这样，前一次分割得到的函数值在分半以后仍可被利用，且易于编程。六、算法的改进和实验总结1、算法进一步的发展：在实际计算中，为防止假收敛，可设置最小二分次数k1，当kk1时，跳过精度判别而继续运算。2、用数值计算的方法来实现数值积分的求解，被调用的函数需要清楚的描述所要执行的问题的求解过程，在matlab函数调用时，执行函数一定要与函数同名。七、源程序（见附件二）附件一[val,M]=Romberg(fun1,[0 pi/2],1e-10)val = 2.9349M = 2.9311 0 0 0 0 0 0 2.9349 2.9362 0 0 0 0 0 2.9349 2.9349 2.9348 0 0 0 0 2.9349 2.9349 2.9349 2.9349 0 0 0 2.9349 2.9349 2.9349 2.9349 2.9349 0 0 2.9349 2.9349 2.9349 2.9349 2.9349 2.9349 0 2.9349 2.9349 2.9349 2.9349 2.9349 2.9349 2.9349[val,M]=Romberg(fun2,[1 3],1e-10)y = 0.3662 0.6486y = 0.5365y = 0.4621 0.5973y = 0.4176 0.5014 0.5682 0.6239y = 0.3929 0.4406 0.4823 0.5194 0.5527 0.5831 0.6109 0.6365y = 0.3798 0.4055 0.4293 0.4515 0.4724 0.4920 0.5105 0.5280 0.5447 0.5606 0.5757