Romberg积分法计算定积分(Matlab程序)

原创 2012年03月30日 18:48:05

%Romberg积分法计算定积分
%参考教材:《数值分析》李乃成,梅立泉,科学出版社
%《计算方法教程》第二版 凌永祥,陈明逵

clear;clc;close all;
format long
% %被积函数为f(x)=4/(1+x^2);积分区间为[0,1]
% b=1;a=0;h=b-a;eps=10^(-5);
%误差界eps%被积函数为f(x)=(x^3+sin(x))/x;积分区间为[0.3,0.8]
%误差界eps=10^(-5)
b=0.8;a=0.3;h=b-a;eps=10^(-5);%误差界eps
kmax=10;%最大递推次数
T1=h/2*((a^3+sin(a))/a+(b^3+sin(b))/b)
S1=0;C1=0;C2=0;R1=0;R2=0;
for k=1:kmax
    h=(b-a)/2^k;
    i=1:2^(k-1);
    x=a+(2*i-1)*h;
    fx=sum((x.^3+sin(x))./x);
    T2=T1/2+fx*h
    S2=T2+(T2-T1)/3
    if(k<3)
        if k==2
            C2=S2+(S2-S1)/15
        end
    else
        C2=S2+(S2-S1)/15
        R2=C2+(C2-C1)/63
        if abs(R2-R1)<eps
            break;
        end
        R1=R2;
    end
    T1=T2;S1=S2;C1=C2;
end
fprintf('所求积分I=%9.8f\n',R2);

相关文章推荐

Romberg法求定积分

1 实验题目:2 算法组织:2.1 算法原理:由梯形算法可以推得:,这里h为2n个区间的步长,即h = (b - a) / 2n;如果有2^n个相同的子区间,则上式变为:,这里hn为2^n个区间的步长...

Romberg求积分算法

// Integral-romberg方法求积分.cpp : 定义控制台应用程序的入口点。 // /* romberg方法求积分 方法也称为逐次分半加速法。它是在梯形公式,simpson公式和newt...

Delphi7高级应用开发随书源码

  • 2003年04月30日 00:00
  • 676KB
  • 下载

Delphi7高级应用开发随书源码

  • 2003年04月30日 00:00
  • 676KB
  • 下载

Matlab计算运行时间方法

(其一) MATLAB实际单位时间计时函数的具体应用,在编写程序时,经常需要获知代码的执行实际时间,这就需要在程序中用到计时函数,matlab中提供了以下三种方法: 1.cputime(单位不...

Delphi7高级应用开发随书源码

  • 2003年04月30日 00:00
  • 676KB
  • 下载

详解协方差与协方差矩阵

协方差的定义   对于一般的分布,直接代入E(X)之类的就可以计算出来了,但真给你一个具体数值的分布,要计算协方差矩阵,根据这个公式来计算,还真不容易反应过来。网上值得参考的资料也不多,这里用...

kirchhoff积分法偏移的程序

  • 2012年08月28日 10:10
  • 12KB
  • 下载

龙贝格积分MATLAB程序

  • 2009年12月26日 15:13
  • 29KB
  • 下载

第4.4节 换元积分法

4.4 INTEGRATION  BY  CHANGE  OF  VARIABLES    We haveseen that the sum, constant, and power rule...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:Romberg积分法计算定积分(Matlab程序)
举报原因:
原因补充:

(最多只允许输入30个字)