龙贝格积分MATLAB编程实验报告_精品文档Word文件下载.doc

龙贝格积分MATLAB编程实验报告_精品文档Word文件下载.doc

《龙贝格积分MATLAB编程实验报告_精品文档Word文件下载.doc》由会员分享，可在线阅读，更多相关《龙贝格积分MATLAB编程实验报告_精品文档Word文件下载.doc（2页珍藏版）》请在冰豆网上搜索。

二、实验步骤（过程）：

（一）龙贝格积分——以书本例4.6来说明

clc,clear;

symsx;

n=8;

a=0;

b=1;

R=0.5*10^（-6）;

%取足够大的0矩阵来构建Romberg积分表

T=zeros（n,n）;

y=sin（x）/x;

p=subs（y,x,a）;

q=subs（y,x,b）;

ifa==0

p=1;

end

T（1,1）=（b-a）/2*（p+q）;

fori=2:

n

f=0;

forj=1:

2^（i-2）

t=a+（（2*j-1）/2^（i-1））*（b-a）;

z=subs（y,x,t）;

f=f+z;

end

T（i,1）=0.5*T（i-1,1）+f*（b-a）/2^（i-1）;

forj=2:

n=n-1;

fori=1:

T（i,j）=（4^（j-1）*T（i+1,j-1）-T（i,j-1））/（4^（j-1）-1）;

T=vpa（T,7）

T=eval（T）;

w=ones（1,7）;

forj=2:

ifabs（T（1,j）-T（1,j-1））<

=R

break;

end

a=T（1,j）；

a=vpa（a,7）%a即为所求答案

三、结论：

（一）龙贝格积分