matlab 高斯迭代的不同写法

来源:互联网 发布:淘宝天机平台不让申请 编辑:程序博客网 时间:2024/06/04 18:06
%高斯回代的高富帅写法
function X=GSHD_gaofushuai(A,B)
N=size(A);
n=N(1);
for i=n:-1:1
    if(i<n)
        deta=A(1,(i 1):n)*X((i 1):n,1);         
    else
        deta=0;
    X(i,1)=(B(i,1) deta)/A(i,i);
    end
end
end

%高斯回代的普通写法
function X=GSHD_putong(A,B)
N=size(A);
n=N(1);
X(n,1)=B(n,1)/A(n,n);
for i=n-1:-1:1
    deta=A(1,(i 1):n)*X((i 1):n,1);     
    X(i,1)=(B(i,1) deta)/A(i,i);
end
end

%高斯回代的屌丝写法
function X=GSHD_diaosi(A,B)
N=size(A);
n=N(1);
X(n,1)=B(n,1)/A(n,n);
for i=n-1:-1:1
    deta=0;
    for j=i 1:n        
        deta deta   A(i,j)*X(j,1);
        disp(deta);        
    end    
    X(i,1)=(B(i,1) deta)/A(i,i);
end
end
原创粉丝点击