Jacobi迭代法和Gauss-Seidel迭代法Matlab程序-豆柴文库

您所在位置：网站首页 / Jacobi迭代法和Gauss-Seidel迭代法Matlab程序.doc / 文档详情

在线预览结束，喜欢就下载吧，查找使用更方便

10 金币

下载文档

/ 4

下载提示

如果您无法下载资料，请参考说明：

1、部分资料下载需要金币，请确保您的账户上有足够的金币

2、已购买过的文档，再次下载不重复扣费

3、资料包下载后请先用软件解压，在使用对应软件打开

文本预览

．分别Jacobi迭代法和Gauss-Seidel迭代法求解线性方程组

迭代法计算停止的条件为：．

解（1）：
采用Jacobi迭代法时，Matlab计算程序为：
clear
clc
i=1;
a=[621-2;250-2;-2085;1327];
d=diag(diag(a));
l=d-tril(a);
u=d-triu(a);
d0=inv(d);
b=[4;7;-1;0];
x0=zeros(4,1);
B=d0*(l+u);
f=d0*b;
x=B*x0+f;
whilenorm(x-x0,inf)>=1e-6
x0=x;
x=B*x0+f;
i=i+1;
end
x
i
计算结果为：x=（0.0520；1.1509；0.2446；-0.5706）迭代次数i=15






采用Gauss-Seidel迭代法计算时，Matlab计算程序为：
clear
clc
i=1;
a=[621-2;250-2;-2085;1327];
d=diag(diag(a));
l=d-tril(a);
u=d-triu(a);
b=[4;7;-1;0];
x0=zeros(4,1);
B=inv(d-l)*u;
f=inv(d-l)*b;
x=B*x0+f;
whilenorm(x-x0,inf)>=1e-6
x0=x;
x=B*x0+f;
i=i+1;
end
x
i
计算结果为：x=（0.0520；1.1509；0.2446；-0.5706）迭代次数i=8