matlab 矩阵jocobi迭代_雅可比迭代法的MATLAB程序

雅可比迭代法的

MA

TLAB

程序：

Function[x,k,index]=Jacobi(A,b,ep,it-max)

%

求线性方程组的雅可比法

;

% A

为方程组的系数矩阵

;

% b

为方程组的右端项

;

% x

为方程组的解

;

% ep

为精度要求

,

缺省值为

le-5;

% it_max

为最大迭代次数

,

缺省值为

100;

% k

为迭代次数

;

% index

为指标变量

,index=0

表示计算失败

,index=1

表示计算成功

;

if nargin<4

it_max=100;

end

if nargin<3

ep=le-5;

end

n=length(A);k=0;x=zeros(n,1);y=zeros(n,1);index=1;

while k<=it_max

for i=1:n

if abs (A(i,i))

index=0;

return;

end

y(i)=(b(i)-A(i,1:n)*x(1:n)+A(i,i)*x(i))/A(i,i);

end

if norm(y-x,inf)

break;

end

k=k+1;

x=y;

end

高斯

-

赛德尔迭代的

MATLAB

程序

Function[x,k,index]=Gau-seidel(A,b,ep,it-max)

%

求线性方程组的高斯

-

赛德尔迭代法

;

% A

为方程组的系数矩阵

;

% b

为方程组的右端项

;

% x

为方程组的解

;

% ep

为精度要求

,

缺省值为

le-5;

% it_max

为最大迭代次数

,

缺省值为

100;

% k

为迭代次数

;

% index

为指标变量

,index=0

表示计算失败

,index=1

表示计算成功

;

if nargin<4

it_max=100;