数学建模社区-数学中国»论坛 › 算法论坛 › 数值分析算法 › 关于数值积分的程序

查看: 4837|回复: 1

[问题求助] 关于数值积分的程序

[复制链接]

字体大小: 正常放大

753085848

8 主题	3 听众	448 积分

升级 49.33%

TA的每日心情

	开心 2011-10-7 20:45

签到天数: 3 天

[LV.2]偶尔看看I

自我介绍: 目前正在数模培训中

电梯直达

1^#

发表于 2010-11-22 21:08 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

function [I,n]=TX(f,a,b,eps)
h=b-a;
n=1;
eps=1e-3;
T1=0;
T2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
while abs(T2-T1)>eps
h=h/2;
T1=T2;
x=a+h;
s=0;
while x<b
s=s+subs(sym(f),findsym(sym(f)),x);
x=x+h;
end
T2=1/2*(T1+h*s);
n=n+1;
end
I=T2+1/3*(T2-T1);

这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！

zan

转播0 淘帖0 分享0 收藏0 支持0 反对0 微信

使用道具举报

753085848

8 主题	3 听众	448 积分

升级 49.33%

TA的每日心情

	开心 2011-10-7 20:45

签到天数: 3 天

[LV.2]偶尔看看I

自我介绍: 目前正在数模培训中

2^#

发表于 2010-11-23 17:43 |只看该作者 |招呼Ta 关注Ta

自己修改过的，是正确的！
function [I,n]=FHTX(f,a,b,ep)
a=a+eps;
h=b-a;
n=1;
ep=1e-6;
T1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
tol=1;
while tol>ep
u=h/2;
x=a+u;
H=0;
while x<b
H=H+subs(sym(f),findsym(sym(f)),x);
x=x+h;
end
T2=1/2*(T1+h*H);
tol=abs(T2-T1);
h=u;
T1=T2;
n=n+1;
end
I=T2+1/3*(T2-T1);