matlab ¾ØÕójocobiµü´ú_MATLABʵÏÖÑſ˱ȵü´ú

%

Ñſ˱ȵü´ú

clear

all

;

close

all

;

clc;

tic;

format

longe

disp(

'

ÇëÊäÈë²ÎÊý

'

);

K=input(

'

άÊý

K='

);

tic

A=100*rand(K);

% A

ÔªËØÊÇ

0-100

for

i=1:K

A(i,i)=sum(abs(A(i,:)))+25*rand(1);

%

¶Ô½ÇÕ¼ÓŵÄÁ¿Îª

0~25

end

b=zeros(K,1);

for

i=1:K;

x=0;

for

r=1:K;

x=x+A(i,r);

end

b(i,1)=x;

end

%

²úÉú

b

¾ØÕó，

b

ÖеÄÔªËØΪ

A

ÖжÔÓ¦ÐеĺÍ，Ä¿µÄÊÇʹ·½³Ì½âȫΪ

1

jd=input(

'

¿ØÖƾ«¶È

jd='

);

I=eye(K);

%

µ¥Î»Õó

L=-tril(A,-1);

%

ÏÂÈý½Ç

U=-triu(A,1);

%

ÉÏÈý½Ç

D=diag(diag(A));

%

¶Ô½Ç¾ØÕó

B=D\(L+U);

%

µü´ú¾ØÕó

»ò

B=I-D\A;

f=D\b;

%f

x0=zeros(K,1);

%

³õʼµü´ú¾ØÕó

y=B*x0+f;

%

µü´ú¹«Ê½

s=1;

%

µü´ú´ÎÊý

while

norm(y-x0)>=jd&&s<3000;

%

µü´úÌõ¼þ

x0=y;

y=B*x0+f;

%

µü´ú¹«Ê½

s=s+1;

end

y

s

toc

t=1:K;

yy=abs(y'-1)/1;