数学建模社区-数学中国
标题:
关于数值积分的程序
[打印本页]
作者:
753085848
时间:
2010-11-22 21:08
标题:
关于数值积分的程序
function [I,n]=TX(f,a,b,eps)
4 r% j' [7 x2 z9 L1 v' W
h=b-a;
6 n+ y) s# T" p
n=1;
& t1 Z v5 z& Y7 n
eps=1e-3;
; L# j1 e/ D3 D
T1=0;
! k6 e9 a% Z1 ?0 D& i# L
T2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
# N/ J. W- l% S9 |, a" ~, Z. g
while abs(T2-T1)>eps
' t1 m4 g, ?' ]" d2 H! ~) k
h=h/2;
# H" @6 }1 h8 s `
T1=T2;
7 Z& J+ p5 f4 e8 N& s, V
x=a+h;
: U) b0 U* N. c
s=0;
" L6 w( v; S. V e& G: Y
while x<b
0 m* ^% W$ L+ [' B
s=s+subs(sym(f),findsym(sym(f)),x);
1 e7 a5 z, K- b
x=x+h;
8 r' D5 ~* ?; O( [2 n7 u/ D
end
; L8 x. i' O6 M: C# I
T2=1/2*(T1+h*s);
- _" F7 E2 O* I. |$ T0 ~) d
n=n+1;
! K% T1 W) {. T: R( b) F2 a
end
5 G. E4 x, E- S; a% F6 u6 O1 J
I=T2+1/3*(T2-T1);
" G, g( ?0 ]+ H; O2 c, ^
这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！
作者:
753085848
时间:
2010-11-23 17:43
自己修改过的，是正确的！
; a Q0 H ?- M) U5 q; s9 l2 t& f8 ~
function [I,n]=FHTX(f,a,b,ep)
5 f" i. o3 c0 \+ f
a=a+eps;
, }( I* x8 K: E. _! B- m
h=b-a;
) H- h; t6 ]) l. s
n=1;
; f, a" I" T% z0 t4 ?* `5 V" G
ep=1e-6;
N! g9 v3 ?4 E! X7 \* p
T1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
0 ?9 z# {8 C& P' Y0 |- L' [
tol=1;
4 d, O3 d. ^& ]
while tol>ep
5 ^1 W4 I5 z' L
u=h/2;
; D; E! n8 k4 _' s! o* l. ]
x=a+u;
) O4 ]% r$ A3 q/ t, y9 x
H=0;
1 m: ~# p5 a0 v+ N1 D
while x<b
B: g4 |" B1 N2 Y* L; r( q
H=H+subs(sym(f),findsym(sym(f)),x);
9 H% t8 ^3 X7 B# @ p# g1 F7 Y) H
x=x+h;
/ Q' ~+ {" Z& o/ [! p2 R
end
) O- j3 [3 [4 e3 e Y* g# B6 q3 `
T2=1/2*(T1+h*H);
5 v6 a: E9 E: T
tol=abs(T2-T1);
L/ r. Q1 v+ @9 Z$ v
h=u;
1 F3 c! w. h9 g9 f
T1=T2;
/ R/ K3 D z8 q% N% ?- k; D
n=n+1;
& N) T4 g d j; _3 i5 g; ^
end
7 Y, A# ]; l" V. V; c- T
I=T2+1/3*(T2-T1);
2 K* L* B' q+ S1 Z) Z8 }2 Y# g* x
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
这是一个关于变步长梯形积分的程序，但是怎么调试就是得不出正确的结果，请各位帮忙看看。帮忙修改一下，顺便指出我的错误，感激不尽！