fdiedai.m

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

function X1=fdiedai(A,b,X0,e)
% 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;
if (length(b)-n)||(length(X0)-n)
    disp('Input Error');
    return
end
% 求迭代矩阵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   
% 求迭代矩阵B1以及g1
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';
% 计算线性方程组AX=b的解X1
q=norm(B1);
X1=B1*X0+g1;
while norm(X1-X0)>q*e/(1-q)
    X0=X1;
    X1=B1*X0+g1;
end