龙贝格积分算法实验.docVIP

实验题目2 Romberg积分法 摘要 考虑积分 欲求其近似值，可以采用如下公式： （复化）梯形公式 （复化）辛卜生公式 （复化）柯特斯公式 这里，梯形公式显得算法简单，具有如下递推关系 因此，很容易实现从低阶的计算结果推算出高阶的近似值，而只需要花费较少的附加函数计算。但是，由于梯形公式收敛阶较低，收敛速度缓慢。所以，如何提高收敛速度，自然是人们极为关心的课题。为此，记为将区间进行等份的复化梯形积分结果，为将区间进行等份的复化辛卜生积分结果，为将区间进行等份的复化柯特斯积分结果。根据李查逊（Richardson）外推加速方法，可得到 可以证明，如果充分光滑，则有 （固定） 这是一个收敛速度更快的一个数值求积公式，我们称为龙贝格积分法。 该方法的计算可按下表进行 … … … … … 很明显，龙贝格计算过程在计算机上实现时，只需开辟一个一维数组，即每次计算的结果，可存放在位置上，其最终结果是存放在位置上。 前言 利用龙贝格积分法计算积分  程序设计流程： 1．准备初值，计算 且（为等份次数） 2．按梯形公式的递推关系，计算 3．按龙贝格公式计算加速值 4．精度控制。对给定的精度，若 则终止计算，并取作为所求结果；否则，重复2~4步，直到满足精度为止。  问题1 程序运行如下： I = Romberginterg(inline(x.^2.*exp(x)),0,1,25,1e-6) I = 0.7183 程序运行如下： I = Romberginterg(inline(exp(x).*sin(x)),1,3,25,1e-6) I = 10.9502 程序运行如下： I = Romberginterg(inline(4./(1+x.^2)),0,1,25,1e-6) I = 3.1416 程序运行如下： I = Romberginterg(inline(1./(1+x)),0,1,25,1e-6) I = 0.6931 问题2 程序运行如下： I = Romberginterg(inline(cos(x)./sqrt(1-x)),0,1,25,1e-6) I = NaN 因为函数在x=0处出现了0/0的情况，极限为1，所以Matlab的结果显示为NaN非数，解决方法是把下限0改为一个小数eps。 I = Romberginterg(inline(sin(x)./x),eps,1,25,1e-6) I = 0.9461 程序运行如下： I = Romberginterg(inline(cos(x)./sqrt(1-x)),0,1,25,1e-6) I = NaN 与（1）的原因相同，函数在x=1处出现了0/0的情况，结果为无穷，此时应选择变换 QUOTE ，将积分变为 QUOTE ，再进行变换 QUOTE ，将积分变为 QUOTE ，变换后输入命令： I = RombergInterg(inline(2*sin(x).*cos((sin(x)).^2)),0,pi/2,25,1e-6) I = 1.4996 程序运行如下： I = Romberginterg(inline(cos(x)./sqrt(x)),0,1,25,1e-6) I = NaN 函数在x=0处出现了1/0的情况，结果为无穷。先将积分变为 QUOTE ，再做变换 QUOTE ， QUOTE I = 2*Romberginterg(inline(cos(x.^2)),0,1,25,1e-6) I = 1.8090 程序运行如下： I = Romberginterg(inline(x.*sin(x)./(1-x.^2)),-1,1,25,1e-6) I = NaN 函数在x=1,-1处出现了sin1/0的情况，结果为无穷。被积函数为偶函数，做变换 QUOTE ，积分变为 QUOTE I = 2*Romberginterg(inline(sin(x).*sin(sin(x))),0,pi/2,25,1e-6) I = 1.382