get_lower_inverse
Funcția este folosită pentru a calcula inversa unei matrice triunghiulare inferioare L.
Implementare
Funcția va returna matricea P
, care este inversa matricei triunghiulare inferioare L
.
determin dimensiunea matricei
L
și se inițializez matriceaP
ca o matrice identitate de dimensiunen
xn
.n = size(L, 1); P = eye(n);calculez elementele diagonale ale matricei
P
prin inversarea elementelor diagonale ale matriceiL
.for j = 1:n P(j, j) = 1 / L(j, j);Calculez elementele sub diagonala 1 ale matricei
P
. Pentru fiecare elementP(i, j)
, se înmulțește secțiunea corespunzătoare a linieii
dinL
cu coloanaj
dinP
(până la elementuli-1
), se negativizează rezultatul și se împarte la elementul diagonalL(i, i)
.for i = j+1:n P(i, j) = -L(i, j:i-1) * P(j:i-1, j) / L(i, i);
Cod sursa
function P = get_lower_inverse(L)
n = size(L, 1);
P = eye(n);
for j = 1:n
P(j, j) = 1 / L(j, j);
for i = j+1:n
P(i, j) = -L(i, j:i-1) * P(j:i-1, j) / L(i, i);
end
end
end
Last modified: 30 April 2024