agsdl.m

电力自动化专业领域的一个重要内容

A=[40,2,6,-3,7,10,-5,6;
    7,32,1,2,-2,9,-8,1;
    9,3,55,-2,15,7,6,-9;
    1,3,2,40,13,-4,10,4;
    -2,-2,11,5,45,-8,7,6;
    -7,8,-9,6,4,51,12,-4;
    4,-14,8,-7,9,10,60,7;
    5,8,-15,10,6,9,-10,65;]
B=[55;63;90;88;101;77;149;123]
X=zeros(8,1);
f='fail';
L=100;
P=100;
EPS=0.000001
while (P>=EPS)&(L~=0)
    P=0.0;
    L=L-1;
for i=1:8
    T=X(i);
    S=0.0;
    for j=1:8
        if j~=i,S=S+A(i,j)*X(j);
        end
    end
    X(i)=(B(i)-S)/A(i,i);
    Q=abs(X(i)-T)/(1+abs(X(i)));
    if Q>P,P=Q;
    end 
end
if P<=EPS
    X
    break
end 
if L==0
    f
    break
end 
end