콜레스테롤 기능을 어떻게 비교합니까?

Question

Cholesky 메서드의 솔루션 단계에서 대체 알고리즘을 앞뒤로 추론하려면 어떻게합니까?

choleskiSol의 기능을 어떻게 비교합니까?

여기에 MATLAB에서 choleskisol

function x = choleskiSol(L,b) 
% Solves [L][L’]{x} = {b} 
% USAGE: x = choleskiSol(L,b) 

n = length(b); 

if size(b,2) > 1 
    b = b’; 
end % {b} must be column vector 

for k = 1:n % Solution of [L]{y} = {b} 
    b(k) = (b(k) - dot(L(k,1:k-1),b(1:k-1)’))/L(k,k); 
end 

for k = n:-1:1 % Solution of {L}’{x} = {y} 
    b(k) = (b(k) - dot(L(k+1:n,k),b(k+1:n)))/L(k,k); 
end 
x = b; 

Answer 1

표준 콜레 분해 (CHOL (A))에 대한 상위 삼각형 형태로 대칭 (포지티브 명확한) 매트릭스,을 분해 내 코드입니다. 선형 방정식 시스템을 풀려면 상 삼각형 형태를 취하고 역 대체를 통해이를 해결해야합니다. 그러면 시스템의 변수 값이 산출됩니다.

L = chol(A); % A must be sym and det(A) > 0 
x = (L \ (L' \ b)); % L' is lower-triangular 

:

은/파라미터 행렬 A 및 출력 B 벡터 w 매트랩 용액을 완성하려면