用Romberg方法求解积分Word文档格式.docx
- 文档编号:21344839
- 上传时间:2023-01-29
- 格式:DOCX
- 页数:10
- 大小:26.21KB
用Romberg方法求解积分Word文档格式.docx
《用Romberg方法求解积分Word文档格式.docx》由会员分享,可在线阅读,更多相关《用Romberg方法求解积分Word文档格式.docx(10页珍藏版)》请在冰豆网上搜索。
S=zeros(1,1)
S(1,1)=(H/2)*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b))
whileEr>
EPS
k=k+1
f1=0
H=H/2
fori=1:
m
x=a+H*(2*i-1)
f1=f1+subs(sym(f),findsym(sym(f)),x)
end
S(k+1,1)=S(k,1)/2+H*f1
m=2*m
forn=1:
k
S(k+1,n+1)=S(k+1,n)+(S(k+1,n)-S(k,n))/(4^n-1)
Er=abs(S(k+1,n+1)-S(k,n))
end
I=S(k+1,k+1)
step=k
命令:
clear
clc
formatshort
a=0;
b=0.8;
EPS=1e-2;
[I,step]=Romberg('
x^(1/2)'
a,b,EPS)
计算结果:
m=
1
k=
0
Er=
H=
S=
f1=
x=
2
0.35780
0.46030
4
0.357800
0.43190.45660
3
0.470900
8
0.47090.47450
0.3578000
0.43190.456600
0.46030.46980.47070
I=
step=
2.设方程组
试用Jacobi迭代法求解此方程,
,当
时终止迭代。
Jacobi.m文件:
functionJacobi(A,b,max,eps)%max为最大迭代次数,eps为容许误差
n=length(A);
x=zeros(n,1);
x1=zeros(n,1);
k=0;
while1
x1
(1)=(b
(1)-A(1,2:
n)*x(2:
n,1))/A(1,1)
fori=2:
n-1
x1(i)=(b(i)-A(i,1:
i-1)*x(1:
i-1,1)-A(i,i+1:
n)*x(i+1:
n,1))/A(i,i)
x1(n)=(b(n)-A(n,1:
n-1)*x(1:
n-1,1))/A(n,n)
ifsum(abs(x1-x))<
eps
fprintf('
number=%d\n'
k)
break
ifk>
=max
TheMethodisdisconvergent\n'
)
x=x1
ifk<
max
n
fprintf('
x[%d]=%f\n'
i,x1(i))
A=[521;
-142;
2-310];
b=[-12203]'
;
max=100;
eps=1e-5
Jacobi(A,b,max,eps)
i=
A=
521
-142
2-310
b=
-12
20
D=
500
040
0010
L=
000
100
-230
U=
0-2-1
00-2
D0=
0.200000
00.25000
x0=
B=
-0.20000.30000
f=
5
6
7
9
10
11
12
13
14
15
16
17
18
19
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- Romberg 方法 求解 积分