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算法计算积分 matlab程序

  • 2011年03月26日 13:31
  • 582B
  • 下载

matlab实现Romberg算法

  • 2010年01月02日 12:56
  • 582B
  • 下载

Romberg法求定积分

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

matlab下利用romberg求积分近似值

  • 2011年12月24日 23:26
  • 742B
  • 下载

Richardson外推加速技术(含Romberg详细分析)的Matlab实现

Richardson外推加速技术可以非常快地收敛到指定的精度,对于提高精度和效率都有非常显著的作用。...
  • lol_IP
  • lol_IP
  • 2017年11月03日 18:55
  • 414

Romberg(龙贝格)求积公式求解数值积分时的注意事项

《数值分析》第5版(李庆扬编著)的第四章课后习题第8-(2)题中,要求使用Romberg(龙贝格)求积公式求解f(x)=xsinx在区间[0,2pi]上的积分,要求误差小于10^(-5)。 针对此问...
  • lbsgrowup
  • lbsgrowup
  • 2015年11月22日 15:05
  • 2474

Romberg求积分算法

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

数值积分-龙贝格(Romberg)积分

  数值积分在工程上是个比较有用的数学工具。在工程上有很多数学问题,看似简单,计算所用的数学公式不算复杂,但是求解起来却很困难,很难获得解析解的公式,这个时候就需要用到数值求解的办法来获取满足工程需要...
  • bslzyp
  • bslzyp
  • 2011年03月03日 10:44
  • 1584

数值积分之Romberg求积法

//Romberg求积法 #include #include #include using namespace std; class romberg { private:  in...
  • u011378809
  • u011378809
  • 2013年08月08日 18:39
  • 706

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

%Romberg积分法计算定积分 %参考教材:《数值分析》李乃成,梅立泉,科学出版社 %《计算方法教程》第二版 凌永祥,陈明逵 clear;clc;close all; format long...
  • zhangchao3322218
  • zhangchao3322218
  • 2012年03月30日 18:48
  • 6945
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:Romberg积分法计算定积分(Matlab程序)
举报原因:
原因补充:

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