diedai.m

自己写的matlab程序 数值分析课程里的程序

clc
A=[8 2 3 -1;-2 7 3 1;3 0 6 -2;3 -2 2 8];
b=[3 0 -5 4]';
e=10^(-7);
X0=[1 0 2 0.5]';
n=length(A);i=1;
% 求迭代矩阵B以及g
for i=1:n
    d(i)=A(i,i);
end
D=diag(d);
B=eye(n)-inv(D)*A
for i=1:n
   g(i)=b(i)/A(i,i); 
end   
% while(i<=n&&A(i,i))
%    for j=1:n
%        B(i,j)=-A(i,j)/A(i,i);
%    end
%    g(i)=b(i)/A(i,i);   
%    i=i+1;
% end
% for i=1:n
%     B(i,i)=0;
% end
for j=1:n
    L(j+1:n,j)=B(j+1:n,j);
end
for i=1:n
    U(i,i+1:n)=B(i,i+1:n);
end
a=inv(eye(n)-L);B1=a*U;g1=a*g';
q=norm(B1)
X1=B1*X0+g1;
while norm(X1-X0)>q*e/(1-q)
    X0=X1;
    X1=B1*X0+g1;
end
X1