sympos1.m

线性方程组的直接发的源代码

function [x,L]= SymPos1 (A,b)
N = size(A);
n = N(1);
L(1,1) = sqrt(A(1,1));
L(2:n,1) = A(2:n,1)/L(1,1);   %L的第一列

for k=2:n
    L(k,k) = sqrt(A(k,k)-L(k,1:(k-1))*transpose(L(k,1:(k-1)))); 
    %L的对角元素
    for i=(k+1):n
        L(i,k) = (A(i,k)-L(i,1:(k-1))*transpose(L(k,1:(k-1))))/L(k,k);
        %L的第k列
    end
end

y = SolveDownTriangle(L,b);
x = SolveUpTriangle(transpose(L),y);   %求解方程