cholesky

Implementare

Funcția va returna o matrice triunghiulară inferioară L caracteristica descompunerii Cholesky

1. determin dimensiunea matricei A și initializez matricea L ca o matrice pătratică de dimensiune n x n

   [n, m] = size(A); L = zeros(n);

2. calculez elementele diagonale ale matricei L folosind rădăcina pătrată a diferenței dintre elementul corespunzător din A și suma pătratelor elementelor anterioare din aceeași linie.

   for i = 1:n L(i, i) = sqrt(A(i, i) - dot(L(i, 1:i-1), L(i, 1:i-1)));

3. calculez elementele non-diagonale ale matricei L prin scăderea produsului scalar al elementelor anterioare din linia i și coloana j și împărțirea rezultatului la elementul diagonal curent din L

   for j = i+1:n L(j, i) = (A(j, i) - dot(L(i, 1:i-1), L(j, 1:i-1))) / L(i, i);

Cod sursa