рдпрд╣ рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рдХреНрдпрд╛ рдХрд░рддрд╛ рд╣реИ
рдпрд╣ рдЯреВрд▓ рдХрд┐рд╕реА рд╡рд░реНрдЧрд╛рдХрд╛рд░ n├Чn рд╡рд╛рд╕реНрддрд╡рд┐рдХ рдореИрдЯреНрд░рд┐рдХреНрд╕ A рдХрд╛ рд╡реНрдпреБрддреНрдХреНрд░рдо (inverse) рдирд┐рдХрд╛рд▓рддрд╛ рд╣реИред рд╡реНрдпреБрддреНрдХреНрд░рдо \(A^{-1}\) рд╡рд╣ рдЕрдиреЛрдЦрд╛ рдореИрдЯреНрд░рд┐рдХреНрд╕ рд╣реЛрддрд╛ рд╣реИ рдЬреЛ \(A \times A^{-1} = I\) рддрддреНрд╕рдордХ рдореИрдЯреНрд░рд┐рдХреНрд╕ (identity matrix) I рдХреЛ рд╕рдВрддреБрд╖реНрдЯ рдХрд░рддрд╛ рд╣реИред рдзреНрдпрд╛рди рд░рдЦреЗрдВ рдХрд┐ рд╣рд░ рдореИрдЯреНрд░рд┐рдХреНрд╕ рдХрд╛ рд╡реНрдпреБрддреНрдХреНрд░рдо рдирд╣реАрдВ рд╣реЛрддрд╛: рдХреЗрд╡рд▓ рд╡рд╣реА рд╡рд░реНрдЧрд╛рдХрд╛рд░ рдореИрдЯреНрд░рд┐рдХреНрд╕ рд╡реНрдпреБрддреНрдХреНрд░рдордгреАрдп рд╣реЛрддреЗ рд╣реИрдВ рдЬрд┐рдирдХрд╛ рд╕рд╛рд░рдгрд┐рдХ (determinant) рд╢реВрдиреНрдп рди рд╣реЛ, рдпрд╛рдиреА рдЬреЛ non-singular рд╣реЛрдВред рдЕрдЧрд░ рдЖрдкрдХрд╛ рдореИрдЯреНрд░рд┐рдХреНрд╕ рд╕рд┐рдВрдЧреБрд▓рд░ рд╣реИ, рддреЛ рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рдЖрдкрдХреЛ рдмреЗрдХрд╛рд░ рдХреЗ рдЕрдВрдХ рджреЗрдиреЗ рдХреЗ рдмрдЬрд╛рдп рд╕рд╛рдлрд╝ рдмрддрд╛ рджреЗрддрд╛ рд╣реИред
рдЗрд╕реЗ рдХреИрд╕реЗ рдЗрд╕реНрддреЗрдорд╛рд▓ рдХрд░реЗрдВ
рдбреНрд░реЙрдкрдбрд╛рдЙрди рд╕реЗ рдореИрдЯреНрд░рд┐рдХреНрд╕ рдХрд╛ рдЖрдХрд╛рд░ n рдЪреБрдиреЗрдВ, рдЬрд┐рд╕рд╕реЗ рдЗрдирдкреБрдЯ рдЧреНрд░рд┐рдб рдЕрдкрдиреЗ рдЖрдк n рдкрдВрдХреНрддрд┐рдпреЛрдВ рдФрд░ n рд╕реНрддрдВрднреЛрдВ рдореЗрдВ рдмрджрд▓ рдЬрд╛рддрд╛ рд╣реИред рдореИрдЯреНрд░рд┐рдХреНрд╕ A рдХреЗ рд╣рд░ рдЦрд╛рдиреЗ рдореЗрдВ рдПрдХ рд╡рд╛рд╕реНрддрд╡рд┐рдХ рд╕рдВрдЦреНрдпрд╛ рднрд░реЗрдВред рдлрд┐рд░ рддрдп рдХрд░реЗрдВ рдХрд┐ рдкрд░рд┐рдгрд╛рдо рдореЗрдВ рдХрд┐рддрдиреЗ рд╕рд╛рд░реНрдердХ рдЕрдВрдХ (significant digits) рджрд┐рдЦрд╛рдиреЗ рд╣реИрдВ, рдФрд░ рд╡реНрдпреБрддреНрдХреНрд░рдо рдореИрдЯреНрд░рд┐рдХреНрд╕, рдЙрд╕рдХрд╛ рд╕рд╛рд░рдгрд┐рдХ рддрдерд╛ рдкреНрд░рдпреБрдХреНрдд рд╡рд┐рдзрд┐ рджреЗрдЦреЗрдВред рдпрд╛рдж рд░рдЦреЗрдВ тАФ рд╕рд╛рд░реНрдердХ рдЕрдВрдХреЛрдВ рдХреА рд╕реЗрдЯрд┐рдВрдЧ рд╕рд┐рд░реНрдлрд╝ рджрд┐рдЦрд╛рдП рдЬрд╛рдиреЗ рд╡рд╛рд▓реЗ рдкрд░рд┐рдгрд╛рдо рдХреА рд░рд╛рдЙрдВрдбрд┐рдВрдЧ рдХреЛ рдкреНрд░рднрд╛рд╡рд┐рдд рдХрд░рддреА рд╣реИ, рднреАрддрд░реА рдЧрдгрдирд╛ рдХреЛ рдирд╣реАрдВ; рдЖрдВрддрд░рд┐рдХ рдЧрдгрдирд╛ рд╣рдореЗрд╢рд╛ рдкреВрд░реА рдбрдмрд▓ рдкреНрд░рд┐рд╕реАрдЬрд╝рди рдореЗрдВ рд╣реЛрддреА рд╣реИред
рд╡рд┐рдзрд┐ рдХреА рд╡реНрдпрд╛рдЦреНрдпрд╛
рдпрд╣ рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рдЖрдВрд╢рд┐рдХ рдкрд┐рд╡рдЯрд┐рдВрдЧ (partial pivoting) рдХреЗ рд╕рд╛рде LU рд╡рд┐рдШрдЯрди рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рддрд╛ рд╣реИред рд╕рдмрд╕реЗ рдкрд╣рд▓реЗ рдпрд╣ рдореИрдЯреНрд░рд┐рдХреНрд╕ рдХреЛ \(PA = LU\) рдХреЗ рд░реВрдк рдореЗрдВ рдЧреБрдгрдирдЦрдВрдбрд┐рдд рдХрд░рддрд╛ рд╣реИ, рдЬрд╣рд╛рдБ P рдПрдХ рдХреНрд░рдордЪрдп рдореИрдЯреНрд░рд┐рдХреНрд╕ (permutation matrix) рд╣реИ рдЬреЛ рдкрдВрдХреНрддрд┐рдпреЛрдВ рдХреА рдЕрджрд▓рд╛-рдмрджрд▓реА рдХрд░рдХреЗ рд╕рдмрд╕реЗ рдмрдбрд╝реЗ рдЙрдкрд▓рдмреНрдз рдкрд┐рд╡рдЯ рдХреЛ рд╡рд┐рдХрд░реНрдг рдкрд░ рд░рдЦрддрд╛ рд╣реИ (рдЗрд╕рд╕реЗ рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рд╕реНрдерд┐рд░рддрд╛ рдмреЗрд╣рддрд░ рд╣реЛрддреА рд╣реИ рдФрд░ рдмрд╣реБрдд рдЫреЛрдЯреА рд╕рдВрдЦреНрдпрд╛рдУрдВ рд╕реЗ рднрд╛рдЧ рджреЗрдиреЗ рд╕реЗ рдмрдЪрд╛рд╡ рд╣реЛрддрд╛ рд╣реИ), L рдПрдХ рдЗрдХрд╛рдИ рдирд┐рдореНрди-рддреНрд░рд┐рдХреЛрдгреАрдп (unit lower-triangular) рдореИрдЯреНрд░рд┐рдХреНрд╕ рд╣реИ, рдФрд░ U рдКрдкрд░реА-рддреНрд░рд┐рдХреЛрдгреАрдп (upper-triangular) рд╣реИред
$$PA = LU \quad\Longrightarrow\quad A^{-1} = U^{-1} L^{-1} P$$
рдЗрд╕рдХреЗ рдмрд╛рдж, рддрддреНрд╕рдордХ рдореИрдЯреНрд░рд┐рдХреНрд╕ рдХреЗ рд╣рд░ рд╕реНрддрдВрдн \(e_k\) рдХреЗ рд▓рд┐рдП, рдпрд╣ рдЖрдЧреЗ рдкреНрд░рддрд┐рд╕реНрдерд╛рдкрди (forward substitution) рд╕реЗ \(L y = P e_k\) рдФрд░ рдкреАрдЫреЗ рдкреНрд░рддрд┐рд╕реНрдерд╛рдкрди (back substitution) рд╕реЗ \(U x = y\) рд╣рд▓ рдХрд░рддрд╛ рд╣реИ; рдЗрд╕ рдкреНрд░рдХрд╛рд░ рдкреНрд░рд╛рдкреНрдд x рд╣реА \(A^{-1}\) рдХрд╛ k-рд╡рд╛рдБ рд╕реНрддрдВрдн рдмрдирддрд╛ рд╣реИред рд╕рд╛рд░рдгрд┐рдХ рдХрд╛ рдорд╛рди U рдХреЗ рд╡рд┐рдХрд░реНрдг рдкрджреЛрдВ рдХреЗ рдЧреБрдгрдирдлрд▓ рдХреЛ рдкрдВрдХреНрддрд┐-рдХреНрд░рдордЪрдп рдХреЗ рдЪрд┐рд╣реНрди рд╕реЗ рдЧреБрдгрд╛ рдХрд░рдХреЗ рдирд┐рдХрд╛рд▓рд╛ рдЬрд╛рддрд╛ рд╣реИред
$$\det(A) = (-1)^{s}\prod_{k=1}^{n} U_{kk}$$
рд╣рд▓ рдХрд┐рдпрд╛ рд╣реБрдЖ рдЙрджрд╛рд╣рд░рдг
рдорд╛рди рд▓реАрдЬрд┐рдП \(A = [[4, 3], [6, 3]]\)ред рдЗрд╕рдХрд╛ рд╕рд╛рд░рдгрд┐рдХ рд╣реИ $$4 \times 3 - 3 \times 6 = 12 - 18 = -6,$$ рдЬреЛ рд╢реВрдиреНрдп рдирд╣реАрдВ рд╣реИ, рдЗрд╕рд▓рд┐рдП A рд╡реНрдпреБрддреНрдХреНрд░рдордгреАрдп рд╣реИред 2├Ч2 рдореИрдЯреНрд░рд┐рдХреНрд╕ рдХреЗ рд▓рд┐рдП рд╕реАрдзреЗ рд╕реВрддреНрд░ рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рддреЗ рд╣реБрдП, $$A^{-1} = \frac{1}{\det} \times [[3, -3], [-6, 4]] = [[-0.5, 0.5], [1, -0.6666666667]].$$ рдЕрдм A рдХреЛ рдЗрд╕ рд╡реНрдпреБрддреНрдХреНрд░рдо рд╕реЗ рдЧреБрдгрд╛ рдХрд░рдиреЗ рдкрд░ рддрддреНрд╕рдордХ рдореИрдЯреНрд░рд┐рдХреНрд╕ рд╡рд╛рдкрд╕ рдорд┐рд▓рддрд╛ рд╣реИ, рдЬреЛ рд╣рдорд╛рд░реЗ рдЙрддреНрддрд░ рдХреА рдкреБрд╖реНрдЯрд┐ рдХрд░рддрд╛ рд╣реИред
рдЕрдХреНрд╕рд░ рдкреВрдЫреЗ рдЬрд╛рдиреЗ рд╡рд╛рд▓реЗ рд╕рд╡рд╛рд▓
рдЕрдЧрд░ рдореЗрд░рд╛ рдореИрдЯреНрд░рд┐рдХреНрд╕ рд╕рд┐рдВрдЧреБрд▓рд░ рд╣реЛ рддреЛ рдХреНрдпрд╛ рд╣реЛрдЧрд╛? рдЕрдЧрд░ рд╕рд╛рд░рдгрд┐рдХ рд╢реВрдиреНрдп рд╣реЛ (рдХрд┐рд╕реА рдмрд╣реБрдд рдЫреЛрдЯреА рд╕рд╣рдирд╢реАрд▓рддрд╛ рд╕реАрдорд╛ рдХреЗ рднреАрддрд░), рддреЛ рдореИрдЯреНрд░рд┐рдХреНрд╕ рдХрд╛ рдХреЛрдИ рд╡реНрдпреБрддреНрдХреНрд░рдо рдирд╣реАрдВ рд╣реЛрддрд╛ рдФрд░ рдХреИрд▓рдХреБрд▓реЗрдЯрд░ "рд╕рд┐рдВрдЧреБрд▓рд░ / рд╡реНрдпреБрддреНрдХреНрд░рдордгреАрдп рдирд╣реАрдВ" рдХрд╛ рд╕рдВрджреЗрд╢ рджрд┐рдЦрд╛рддрд╛ рд╣реИред
рдПрдбрдЬреБрдЧреЗрдЯ рд╕реВрддреНрд░ рдХреА рдЬрдЧрд╣ LU рд╡рд┐рдШрдЯрди рдХреНрдпреЛрдВ рдЗрд╕реНрддреЗрдорд╛рд▓ рдХрд░реЗрдВ? рдмрдбрд╝реЗ рдореИрдЯреНрд░рд┐рдХреНрд╕ рдХреЗ рд▓рд┐рдП рдЖрдВрд╢рд┐рдХ рдкрд┐рд╡рдЯрд┐рдВрдЧ рд╡рд╛рд▓рд╛ LU рд╡рд┐рдШрдЯрди рдХреЛрдлрд╝реИрдХреНрдЯрд░ рд╡рд┐рд╕реНрддрд╛рд░ (cofactor expansion) рдХреА рддреБрд▓рдирд╛ рдореЗрдВ рдХрд╣реАрдВ рдЕрдзрд┐рдХ рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рд░реВрдк рд╕реЗ рд╕реНрдерд┐рд░ рдФрд░ рдХреБрд╢рд▓ рд╣реИ, рдЬрдмрдХрд┐ рдХреЛрдлрд╝реИрдХреНрдЯрд░ рд╡рд┐рдзрд┐ рдХреА рд▓рд╛рдЧрдд рдлрд╝реИрдХреНрдЯреЛрд░рд┐рдпрд▓ рджрд░ рд╕реЗ рдмрдврд╝рддреА рдЬрд╛рддреА рд╣реИред
рдХреНрдпрд╛ рд╕рд╛рд░реНрдердХ рдЕрдВрдХреЛрдВ рдХрд╛ рдЪрдпрди рдЧрдгрд┐рдд рдХреЛ рдмрджрд▓ рджреЗрддрд╛ рд╣реИ? рдирд╣реАрдВред рдЧрдгрдирд╛ рд╣рдореЗрд╢рд╛ рдкреВрд░реА рдкреНрд░рд┐рд╕реАрдЬрд╝рди рдореЗрдВ рд╣реЛрддреА рд╣реИ; рдпрд╣ рд╕реЗрдЯрд┐рдВрдЧ рдХреЗрд╡рд▓ рдпрд╣ рддрдп рдХрд░рддреА рд╣реИ рдХрд┐ рдкрд░рд┐рдгрд╛рдо рдореЗрдВ рдХрд┐рддрдиреЗ рд╕рд╛рд░реНрдердХ рдЕрдВрдХ рджрд┐рдЦрд╛рдП рдЬрд╛рдПрдБред