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

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

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

龙贝格积分MATLAB编程实验报告

一、问题描述：

编写函数实现龙贝格积分。

二、实验步骤（过程）：

（一）龙贝格积分——以书本例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）;

end

forj=2:

n

n=n-1;

fori=1:

n

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

end

end

T=vpa（T,7）

T=eval（T）;

n=8;

w=ones（1,7）;

forj=2:

n

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

break;

end

end

a=T（1,j）；

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

三、结论：

（一）龙贝格积分