рдпрд╣ рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рдХреНрдпрд╛ рдХрд░рддрд╛ рд╣реИ
рдпрд╣ рдЯреВрд▓ n рдЕрдЬреНрдЮрд╛рддреЛрдВ рд╡рд╛рд▓реЗ n рд░реИрдЦрд┐рдХ рд╕рдореАрдХрд░рдгреЛрдВ рдХреЗ рдирд┐рдХрд╛рдп рдХреЛ рд╣рд▓ рдХрд░рддрд╛ рд╣реИ, рдЬрд┐рд╕реЗ рд╕рдВрдХреНрд╖реЗрдк рдореЗрдВ \(\mathbf{A}\cdot\mathbf{x} = \mathbf{b}\) рдХреЗ рд░реВрдк рдореЗрдВ рд▓рд┐рдЦрд╛ рдЬрд╛рддрд╛ рд╣реИред рдпрд╣рд╛рдБ A рдПрдХ n\(\times\)n рдЧреБрдгрд╛рдВрдХ рдореИрдЯреНрд░рд┐рдХреНрд╕ рд╣реИ, x рдЕрдЬреНрдЮрд╛рдд рд╡реЗрдХреНрдЯрд░ рд╣реИ, рдФрд░ b рд╕реНрдерд┐рд░рд╛рдВрдХ рд╡реЗрдХреНрдЯрд░ рд╣реИред рдпрд╣ рдЕрджреНрд╡рд┐рддреАрдп рд╣рд▓ рд╡реЗрдХреНрдЯрд░ x рдХреЗ рд╕рд╛рде-рд╕рд╛рде A рдХрд╛ рд╕рд╛рд░рдгрд┐рдХ (determinant) рднреА рджреЗрддрд╛ рд╣реИред рдпрд╣ рд╡рд┐рдзрд┐ рд╢реБрджреНрдз рд░реИрдЦрд┐рдХ рдмреАрдЬрдЧрдгрд┐рдд рдкрд░ рдЖрдзрд╛рд░рд┐рдд рд╣реИ рдФрд░ рд╣рд░ рдЬрдЧрд╣ рдПрдХ рдЬреИрд╕реА рдХрд╛рдо рдХрд░рддреА рд╣реИ тАФ рдЗрд╕рдореЗрдВ рдХрд┐рд╕реА рджреЗрд╢ рдпрд╛ рдЗрдХрд╛рдИ рдХреА рдХреЛрдИ рдирд┐рд░реНрднрд░рддрд╛ рдирд╣реАрдВ рд╣реИ; рд╣рд░ рдкреНрд░рд╡рд┐рд╖реНрдЯрд┐ рдмрд╕ рдПрдХ рд╡рд╛рд╕реНрддрд╡рд┐рдХ рд╕рдВрдЦреНрдпрд╛ рд╣реИред
рдЗрд╕рдХрд╛ рдЙрдкрдпреЛрдЧ рдХреИрд╕реЗ рдХрд░реЗрдВ
рд╕рдмрд╕реЗ рдкрд╣рд▓реЗ n рд╕реЗрдЯ рдХрд░реЗрдВ (рд╕рдореАрдХрд░рдгреЛрдВ рдФрд░ рдЕрдЬреНрдЮрд╛рддреЛрдВ рдХреА рд╕рдВрдЦреНрдпрд╛)ред рдЧреБрдгрд╛рдВрдХ рдореИрдЯреНрд░рд┐рдХреНрд╕ A рдХреЛ рдкреНрд░рддрд┐ рдкрдВрдХреНрддрд┐ рдПрдХ рд╣реА рд▓рд╛рдЗрди рдореЗрдВ рдЯрд╛рдЗрдк рдХрд░реЗрдВ, рдЬрд╣рд╛рдБ рд╕рдВрдЦреНрдпрд╛рдПрдБ рд╕реНрдкреЗрд╕ рдпрд╛ рдХреЙрдорд╛ рд╕реЗ рдЕрд▓рдЧ рд╣реЛрдВред рдлрд┐рд░ рд╕реНрдерд┐рд░рд╛рдВрдХ рд╡реЗрдХреНрдЯрд░ b рдХреЛ n рд▓рдВрдмрд╛рдИ рдХреА рд╕реВрдЪреА рдХреЗ рд░реВрдк рдореЗрдВ рджрд░реНрдЬ рдХрд░реЗрдВред рдкреНрд░рджрд░реНрд╢рди рдкрд░рд┐рд╢реБрджреНрдзрддрд╛ рдЪреБрдиреЗрдВ рдФрд░ рд╣рд▓ рдХрд░реЗрдВред рдЛрдгрд╛рддреНрдордХ, рджрд╢рдорд▓рд╡ рдФрд░ рднрд┐рдиреНрди рдЬреИрд╕реА рджрд┐рдЦрдиреЗ рд╡рд╛рд▓реА рджрд╢рдорд▓рд╡ рд╕рдВрдЦреНрдпрд╛рдПрдБ тАФ рд╕рднреА рд╕реНрд╡реАрдХрд╛рд░реНрдп рд╣реИрдВред рдпрджрд┐ A рдореЗрдВ рдкрдВрдХреНрддрд┐рдпреЛрдВ рдФрд░ рд╕реНрддрдВрднреЛрдВ рдХреА рд╕рдВрдЦреНрдпрд╛ n рдХреЗ рдмрд░рд╛рдмрд░ рд╣реИ рдФрд░ b рдХреА рд▓рдВрдмрд╛рдИ рднреА рдЙрддрдиреА рд╣реА рд╣реИ, рддреЛ рдЖрдкрдХреЛ рд╣рд▓ рдорд┐рд▓ рдЬрд╛рдПрдЧрд╛; рдЕрдиреНрдпрдерд╛ рдЯреВрд▓ рдЖрдпрд╛рдо рдХреЗ рдмреЗрдореЗрд▓ (dimension mismatch) рдХреА рд╕реВрдЪрдирд╛ рджреЗрдЧрд╛ред
рд╡рд┐рдзрд┐ рдХреА рд╡реНрдпрд╛рдЦреНрдпрд╛
рдпрд╣ рд╕реЙрд▓реНрд╡рд░ рдЖрдВрд╢рд┐рдХ рдкрд┐рд╡рдЯрд┐рдВрдЧ (partial pivoting) рдХреЗ рд╕рд╛рде рдЧреЙрд╕рд┐рдпрди рдЙрдиреНрдореВрд▓рди рдХрд░рддрд╛ рд╣реИ, рдЬреЛ рдЧрдгрд┐рддреАрдп рд░реВрдк рд╕реЗ LU рдЕрдкрдШрдЯрди PA = LU рдХреЗ рд╕рдорддреБрд▓реНрдп рд╣реИред рдкреНрд░рддреНрдпреЗрдХ рд╕реНрддрдВрдн рдХреЗ рд▓рд┐рдП рдпрд╣ рд╕рдмрд╕реЗ рдмрдбрд╝реЗ рдкрд░рд┐рдорд╛рдг рд╡рд╛рд▓реЗ рдЙрдкрд▓рдмреНрдз рдкрд┐рд╡рдЯ рдХреЛ рдЪреБрдирддрд╛ рд╣реИ рддрд╛рдХрд┐ рдЕрдВрдХрдЧрдгрд┐рдд рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рд░реВрдк рд╕реЗ рд╕реНрдерд┐рд░ рд░рд╣реЗ, рдкрд┐рд╡рдЯ рдХреЗ рдиреАрдЪреЗ рдХреА рдкреНрд░рд╡рд┐рд╖реНрдЯрд┐рдпреЛрдВ рдХреЛ рд╕рдорд╛рдкреНрдд рдХрд░рддрд╛ рд╣реИ, рдФрд░ рдлрд┐рд░ рдЕрдВрддрд┐рдо рдЕрдЬреНрдЮрд╛рдд рд╕реЗ рдКрдкрд░ рдХреА рдУрд░ рдкреНрд░рддрд┐рд╕реНрдерд╛рдкрди (back-substitution) рдХрд░рддрд╛ рд╣реИред рдкреНрд░рдХреНрд░рд┐рдпрд╛ рдХреЛ рд╕рдВрдХреНрд╖реЗрдк рдореЗрдВ рдЗрд╕ рдкреНрд░рдХрд╛рд░ рд▓рд┐рдЦрд╛ рдЬрд╛ рд╕рдХрддрд╛ рд╣реИ:
$$\begin{gathered} \mathbf{A}\,\mathbf{x} = \mathbf{b}, \qquad \mathbf{A} = \mathbf{L}\mathbf{U} \\[1.5em] \text{solve}\quad \left\{ \begin{aligned} \mathbf{L}\mathbf{y} &= \mathbf{b} \quad(\text{forward}) \\ \mathbf{U}\mathbf{x} &= \mathbf{y} \quad(\text{back substitution}) \end{aligned} \right. \end{gathered}$$рдпрджрд┐ рдХреЛрдИ рдкрд┐рд╡рдЯ рд╡реНрдпрд╛рд╡рд╣рд╛рд░рд┐рдХ рд░реВрдк рд╕реЗ рд╢реВрдиреНрдп рд╣реИ, рддреЛ рд╕рд╛рд░рдгрд┐рдХ рд╢реВрдиреНрдп рд╣реЛрддрд╛ рд╣реИ, рдореИрдЯреНрд░рд┐рдХреНрд╕ рд╡рд┐рдЪрд┐рддреНрд░ (singular) рдХрд╣рд▓рд╛рддрд╛ рд╣реИ, рдФрд░ рдирд┐рдХрд╛рдп рдХрд╛ рдХреЛрдИ рдЕрджреНрд╡рд┐рддреАрдп рд╣рд▓ рдирд╣реАрдВ рд╣реЛрддрд╛ тАФ рдРрд╕реЗ рдореЗрдВ рдЯреВрд▓ рд╢реВрдиреНрдп рд╕реЗ рднрд╛рдЧ рджреЗрдиреЗ рдХреЗ рдмрдЬрд╛рдп рдЗрд╕рдХреА рд╕реВрдЪрдирд╛ рджреЗ рджреЗрддрд╛ рд╣реИред рд╕рд╛рд░рдгрд┐рдХ рдХреА рдЧрдгрдирд╛ рдЗрд╕ рдкреНрд░рдХрд╛рд░ рд╣реЛрддреА рд╣реИ:
$$\det(\mathbf{A}) = (-1)^{s}\prod_{k=1}^{n} u_{kk}$$
рд╣рд▓ рдХрд┐рдпрд╛ рдЧрдпрд╛ рдЙрджрд╛рд╣рд░рдг
рдорд╛рди рд▓реАрдЬрд┐рдП \(2x + y - z = 8\), \(-3x - y + 2z = -11\), \(-2x + y + 2z = -3\)ред рддреЛ \(\mathbf{A} = [[2,1,-1],[-3,-1,2],[-2,1,2]]\) рдФрд░ \(\mathbf{b} = [8,-11,-3]\) рд╣реЛрдЧрд╛ред рдЙрдиреНрдореВрд▓рди рдХрд░рдиреЗ рдкрд░ \(x = 2\), \(y = 3\), \(z = -1\) рдорд┐рд▓рддрд╛ рд╣реИред рдкрд╣рд▓реЗ рд╕рдореАрдХрд░рдг рдХреА рдЬрд╛рдБрдЪ рдХрд░реЗрдВ:
$$2(2) + 3 - (-1) = 4 + 3 + 1 = 8$$рд╕рд╣реА рд╣реИред
рдЕрдХреНрд╕рд░ рдкреВрдЫреЗ рдЬрд╛рдиреЗ рд╡рд╛рд▓реЗ рдкреНрд░рд╢реНрди
рдЕрдЧрд░ рд╕рд╛рд░рдгрд┐рдХ рд╢реВрдиреНрдп рд╣реЛ рддреЛ рдХреНрдпрд╛ рд╣реЛрдЧрд╛? рдореИрдЯреНрд░рд┐рдХреНрд╕ рд╡рд┐рдЪрд┐рддреНрд░ (singular) рд╣реЛрддреА рд╣реИ, рдпрд╛рдиреА рд╕рдореАрдХрд░рдг рдЖрдкрд╕ рдореЗрдВ рдирд┐рд░реНрднрд░ рдпрд╛ рдЕрд╕рдВрдЧрдд рд╣реИрдВ; рдРрд╕реЗ рдореЗрдВ рдХреЛрдИ рдПрдХрд▓ рд╣рд▓ рдирд╣реАрдВ рд╣реЛрддрд╛, рдЗрд╕рд▓рд┐рдП рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рд╡рд┐рдЪрд┐рддреНрд░ рдореИрдЯреНрд░рд┐рдХреНрд╕ рдХреА рд╕реВрдЪрдирд╛ рджреЗрддрд╛ рд╣реИред
рдЖрдВрд╢рд┐рдХ рдкрд┐рд╡рдЯрд┐рдВрдЧ рдХреНрдпреЛрдВ? рд╕рдмрд╕реЗ рдмрдбрд╝реЗ рдкрд┐рд╡рдЯ рдХреЛ рдЪреБрдирдиреЗ рд╕реЗ рд░рд╛рдЙрдВрдбрд┐рдВрдЧ рддреНрд░реБрдЯрд┐ рдмрдврд╝рдиреЗ рд╕реЗ рдмрдЪрд╛рд╡ рд╣реЛрддрд╛ рд╣реИ, рдЬрд┐рд╕рд╕реЗ рдХрдард┐рди рдореИрдЯреНрд░рд┐рдХреНрд╕ рдХреЗ рд▓рд┐рдП рднреА рд╕рдЯреАрдХ рдкрд░рд┐рдгрд╛рдо рдорд┐рд▓рддреЗ рд╣реИрдВред
рдХреНрдпрд╛ рд╣рд▓ рдкреВрд░реНрдгрд╛рдВрдХ рдХреЗ рдЕрд▓рд╛рд╡рд╛ рднреА рд╣реЛ рд╕рдХрддрд╛ рд╣реИ? рд╣рд╛рдБред рд╣рд▓ рдлреНрд▓реЛрдЯрд┐рдВрдЧ-рдкреЙрдЗрдВрдЯ рдореЗрдВ рдЧрдгрдирд╛ рдХрд┐рдП рдЬрд╛рддреЗ рд╣реИрдВ рдФрд░ рджрд╢рдорд▓рд╡ рд╣реЛ рд╕рдХрддреЗ рд╣реИрдВ; рдкреНрд░рджрд░реНрд╢рди рдкрд░рд┐рд╢реБрджреНрдзрддрд╛ рдХреА рд╕реЗрдЯрд┐рдВрдЧ рдпрд╣ рдирд┐рдпрдВрддреНрд░рд┐рдд рдХрд░рддреА рд╣реИ рдХрд┐ рдХрд┐рддрдиреЗ рд╕рд╛рд░реНрдердХ рдЕрдВрдХ (significant figures) рджрд┐рдЦрд╛рдП рдЬрд╛рдПрдБред