Answer on Question #81954 — Math — Linear Algebra
Question
Solve the following
0.6 x + 0.8 y + 0.1 z = 1 0.6x + 0.8y + 0.1z = 1 0.6 x + 0.8 y + 0.1 z = 1 1.1 x + 0.4 y + 0.3 z = 0.2 1.1x + 0.4y + 0.3z = 0.2 1.1 x + 0.4 y + 0.3 z = 0.2 x + y + 2 z = 0.5 x + y + 2z = 0.5 x + y + 2 z = 0.5
by LU decomposition method and find the inverse of the coefficient matrix
Solution
A = ( 0.6 0.8 0.1 1.1 0.4 0.3 1 1 2 ) ; b = ( 1 0.2 0.5 ) A = \begin{pmatrix} 0.6 & 0.8 & 0.1 \\ 1.1 & 0.4 & 0.3 \\ 1 & 1 & 2 \end{pmatrix}; \quad b = \begin{pmatrix} 1 \\ 0.2 \\ 0.5 \end{pmatrix} A = ⎝ ⎛ 0.6 1.1 1 0.8 0.4 1 0.1 0.3 2 ⎠ ⎞ ; b = ⎝ ⎛ 1 0.2 0.5 ⎠ ⎞
Swap some rows in matrices.
A = ( 1 1 2 1.1 0.4 0.3 0.6 0.8 0.1 ) ; b = ( 0.5 0.2 1 ) A = \begin{pmatrix} 1 & 1 & 2 \\ 1.1 & 0.4 & 0.3 \\ 0.6 & 0.8 & 0.1 \end{pmatrix}; \quad b = \begin{pmatrix} 0.5 \\ 0.2 \\ 1 \end{pmatrix} A = ⎝ ⎛ 1 1.1 0.6 1 0.4 0.8 2 0.3 0.1 ⎠ ⎞ ; b = ⎝ ⎛ 0.5 0.2 1 ⎠ ⎞
LU = A
L = ( 1 0 0 ∗ 1 0 ∗ ∗ 1 ) ; U = ( 1 ∗ ∗ 0 1 ∗ 0 0 1 ) ; L = \begin{pmatrix} 1 & 0 & 0 \\ * & 1 & 0 \\ * & * & 1 \end{pmatrix}; \quad U = \begin{pmatrix} 1 & * & * \\ 0 & 1 & * \\ 0 & 0 & 1 \end{pmatrix}; L = ⎝ ⎛ 1 ∗ ∗ 0 1 ∗ 0 0 1 ⎠ ⎞ ; U = ⎝ ⎛ 1 0 0 ∗ 1 0 ∗ ∗ 1 ⎠ ⎞ ; ( 1 1 2 1.1 0.4 0.3 0.6 0.8 0.1 ) ( 0.5 0.2 1 ) → R 2 ∗ = R 2 + ( − 1.1 ) ∗ R 1 ; R 3 ∗ = R 3 + ( − 0.6 ) ∗ R 1 ; → \begin{pmatrix} 1 & 1 & 2 \\ 1.1 & 0.4 & 0.3 \\ 0.6 & 0.8 & 0.1 \end{pmatrix} \begin{pmatrix} 0.5 \\ 0.2 \\ 1 \end{pmatrix} \to R_2^* = R_2 + (-1.1) * R_1; \quad R_3^* = R_3 + (-0.6) * R_1; \quad \to ⎝ ⎛ 1 1.1 0.6 1 0.4 0.8 2 0.3 0.1 ⎠ ⎞ ⎝ ⎛ 0.5 0.2 1 ⎠ ⎞ → R 2 ∗ = R 2 + ( − 1.1 ) ∗ R 1 ; R 3 ∗ = R 3 + ( − 0.6 ) ∗ R 1 ; → ( 1 1 2 0 − 0.7 − 1.9 0 0.2 − 1.1 ) ( 0.5 − 0.35 0.7 ) ; \begin{pmatrix} 1 & 1 & 2 \\ 0 & -0.7 & -1.9 \\ 0 & 0.2 & -1.1 \end{pmatrix} \begin{pmatrix} 0.5 \\ -0.35 \\ 0.7 \end{pmatrix}; ⎝ ⎛ 1 0 0 1 − 0.7 0.2 2 − 1.9 − 1.1 ⎠ ⎞ ⎝ ⎛ 0.5 − 0.35 0.7 ⎠ ⎞ ; U = ( 1 1 2 0 − 0.7 − 1.9 0 0 − 1.643 ) ; b = ( 0.5 − 0.35 0.6 ) U = \begin{pmatrix} 1 & 1 & 2 \\ 0 & -0.7 & -1.9 \\ 0 & 0 & -1.643 \end{pmatrix}; \quad b = \begin{pmatrix} 0.5 \\ -0.35 \\ 0.6 \end{pmatrix} U = ⎝ ⎛ 1 0 0 1 − 0.7 0 2 − 1.9 − 1.643 ⎠ ⎞ ; b = ⎝ ⎛ 0.5 − 0.35 0.6 ⎠ ⎞ L = ( 1 0 0 1.1 1 0 0.6 ∗ 1 ) L = \begin{pmatrix} 1 & 0 & 0 \\ 1.1 & 1 & 0 \\ 0.6 & * & 1 \end{pmatrix} L = ⎝ ⎛ 1 1.1 0.6 0 1 ∗ 0 0 1 ⎠ ⎞ ( 1 1 2 0 − 0.7 − 1.9 0 0.2 − 1.1 ) ( 0.5 − 0.35 0.7 ) → R 3 ∗ = R 3 + ( 0.2 0.7 ) ∗ R 2 → ( 1 1 2 0 − 0.7 − 1.9 0 0 − 1.643 ) ( 0.5 − 0.35 0.6 ) \begin{pmatrix} 1 & 1 & 2 \\ 0 & -0.7 & -1.9 \\ 0 & 0.2 & -1.1 \end{pmatrix} \begin{pmatrix} 0.5 \\ -0.35 \\ 0.7 \end{pmatrix} \to R_3^* = R_3 + \begin{pmatrix} 0.2 \\ 0.7 \end{pmatrix} * R_2 \to \begin{pmatrix} 1 & 1 & 2 \\ 0 & -0.7 & -1.9 \\ 0 & 0 & -1.643 \end{pmatrix} \begin{pmatrix} 0.5 \\ -0.35 \\ 0.6 \end{pmatrix} ⎝ ⎛ 1 0 0 1 − 0.7 0.2 2 − 1.9 − 1.1 ⎠ ⎞ ⎝ ⎛ 0.5 − 0.35 0.7 ⎠ ⎞ → R 3 ∗ = R 3 + ( 0.2 0.7 ) ∗ R 2 → ⎝ ⎛ 1 0 0 1 − 0.7 0 2 − 1.9 − 1.643 ⎠ ⎞ ⎝ ⎛ 0.5 − 0.35 0.6 ⎠ ⎞ L = ( 1 0 0 1.1 1 0 0.6 − 0.2857 1 ) L = \begin{pmatrix} 1 & 0 & 0 \\ 1.1 & 1 & 0 \\ 0.6 & -0.2857 & 1 \end{pmatrix} L = ⎝ ⎛ 1 1.1 0.6 0 1 − 0.2857 0 0 1 ⎠ ⎞ U = ( 1 1 2 0 − 0.7 − 1.9 0 0 − 1.643 ) ; b = ( 0.5 − 0.35 0.6 ) U = \begin{pmatrix} 1 & 1 & 2 \\ 0 & -0.7 & -1.9 \\ 0 & 0 & -1.643 \end{pmatrix}; \quad b = \begin{pmatrix} 0.5 \\ -0.35 \\ 0.6 \end{pmatrix} U = ⎝ ⎛ 1 0 0 1 − 0.7 0 2 − 1.9 − 1.643 ⎠ ⎞ ; b = ⎝ ⎛ 0.5 − 0.35 0.6 ⎠ ⎞ { x + y + 2 z = 0.5 − 0.7 y − 1.9 z = − 0.35 − 1.643 z = 0.6 \begin{cases}
x + y + 2z = 0.5 \\
-0.7y - 1.9z = -0.35 \\
-1.643z = 0.6
\end{cases} ⎩ ⎨ ⎧ x + y + 2 z = 0.5 − 0.7 y − 1.9 z = − 0.35 − 1.643 z = 0.6
So, we have: Z = 0.6 / ( − 1.643 ) = − 0.365 Z = 0.6 / (-1.643) = -0.365 Z = 0.6/ ( − 1.643 ) = − 0.365
Y = − 0.35 + 1.9 ∗ ( − 0.365 ) − 0.7 = 1.4907 Y = \frac {- 0.35 + 1.9 * (- 0.365)}{- 0.7} = 1.4907 Y = − 0.7 − 0.35 + 1.9 ∗ ( − 0.365 ) = 1.4907 X = 0.5 − 2 ∗ ( − 0.365 ) − 1.4907 = − 0.2607 X = 0.5 - 2 * (-0.365) - 1.4907 = -0.2607 X = 0.5 − 2 ∗ ( − 0.365 ) − 1.4907 = − 0.2607
Find the inverse of the coefficient matrix
( 0.6 0.8 0.1 1 0 0 1.1 0.4 0.3 0 1 0 1 1 2 0 0 1 ) ; \left( \begin{array}{cccc} 0.6 & 0.8 & 0.1 & 1 & 0 & 0 \\ 1.1 & 0.4 & 0.3 & 0 & 1 & 0 \\ 1 & 1 & 2 & 0 & 0 & 1 \end{array} \right); ⎝ ⎛ 0.6 1.1 1 0.8 0.4 1 0.1 0.3 2 1 0 0 0 1 0 0 0 1 ⎠ ⎞ ;
1. R 1 ∗ = R 1 / 0.6 R_{1}^{*} = R_{1} / 0.6 R 1 ∗ = R 1 /0.6
( 1 1.33 0.167 1.67 0 0 1.1 0.4 0.3 0 1 0 1 1 2 0 0 1 ) ; \left( \begin{array}{cccc} 1 & 1.33 & 0.167 & 1.67 & 0 & 0 \\ 1.1 & 0.4 & 0.3 & 0 & 1 & 0 \\ 1 & 1 & 2 & 0 & 0 & 1 \end{array} \right); ⎝ ⎛ 1 1.1 1 1.33 0.4 1 0.167 0.3 2 1.67 0 0 0 1 0 0 0 1 ⎠ ⎞ ;
2. R 2 ∗ = R 2 − 1.1 ∗ R 1 ; R 3 ∗ = R 3 − R 1 R_{2}^{*} = R_{2} - 1.1 * R_{1}; R_{3}^{*} = R_{3} - R_{1} R 2 ∗ = R 2 − 1.1 ∗ R 1 ; R 3 ∗ = R 3 − R 1
( 1 1.33 0.167 1.67 0 0 0 − 1.067 0.1167 − 1.83 1 0 0 − 0.33 1.83 − 1.67 0 1 ) ; \left( \begin{array}{cccc} 1 & 1.33 & 0.167 & 1.67 & 0 & 0 \\ 0 & -1.067 & 0.1167 & -1.83 & 1 & 0 \\ 0 & -0.33 & 1.83 & -1.67 & 0 & 1 \end{array} \right); ⎝ ⎛ 1 0 0 1.33 − 1.067 − 0.33 0.167 0.1167 1.83 1.67 − 1.83 − 1.67 0 1 0 0 0 1 ⎠ ⎞ ;
3. R 2 ∗ = R 2 / ( − 1.067 ) R_{2}^{*} = R_{2} / (-1.067) R 2 ∗ = R 2 / ( − 1.067 )
( 1 1.33 0.167 1.67 0 0 0 1 − 0.109 1.719 − 0.9375 0 0 − 0.33 1.83 − 1.67 0 1 ) ; \left( \begin{array}{cccc} 1 & 1.33 & 0.167 & 1.67 & 0 & 0 \\ 0 & 1 & -0.109 & 1.719 & -0.9375 & 0 \\ 0 & -0.33 & 1.83 & -1.67 & 0 & 1 \end{array} \right); ⎝ ⎛ 1 0 0 1.33 1 − 0.33 0.167 − 0.109 1.83 1.67 1.719 − 1.67 0 − 0.9375 0 0 0 1 ⎠ ⎞ ;
4. R 1 ∗ = R 1 − 1.33 ∗ R 2 ; R 3 ∗ = R 3 + 0.33 ∗ R 2 R_{1}^{*} = R_{1} - 1.33 * R_{2}; R_{3}^{*} = R_{3} + 0.33 * R_{2} R 1 ∗ = R 1 − 1.33 ∗ R 2 ; R 3 ∗ = R 3 + 0.33 ∗ R 2
( 1 0 0.3125 − 0.625 1.25 0 0 1 − 0.109 1.719 − 0.9375 0 0 0 1.797 − 1.094 − 0.3125 1 ) ; \left( \begin{array}{cccc} 1 & 0 & 0.3125 & -0.625 & 1.25 & 0 \\ 0 & 1 & -0.109 & 1.719 & -0.9375 & 0 \\ 0 & 0 & 1.797 & -1.094 & -0.3125 & 1 \end{array} \right); ⎝ ⎛ 1 0 0 0 1 0 0.3125 − 0.109 1.797 − 0.625 1.719 − 1.094 1.25 − 0.9375 − 0.3125 0 0 1 ⎠ ⎞ ;
5. R 3 ∗ = R 3 1.797 R_{3}^{*} = \frac{R_{3}}{1.797} R 3 ∗ = 1.797 R 3
( 1 0 0.3125 − 0.625 1.25 0 0 1 − 0.109 1.719 − 0.9375 0 0 0 1 − 0.609 − 0.174 0.556 ) ; \left( \begin{array}{cccc} 1 & 0 & 0.3125 & -0.625 & 1.25 & 0 \\ 0 & 1 & -0.109 & 1.719 & -0.9375 & 0 \\ 0 & 0 & 1 & -0.609 & -0.174 & 0.556 \end{array} \right); ⎝ ⎛ 1 0 0 0 1 0 0.3125 − 0.109 1 − 0.625 1.719 − 0.609 1.25 − 0.9375 − 0.174 0 0 0.556 ⎠ ⎞ ;
6. R 1 ∗ = R 1 − 0.3125 ∗ R 3 ; R 2 ∗ = R 2 + 0.109 ∗ R 3 R_{1}^{*} = R_{1} - 0.3125 * R_{3}; R_{2}^{*} = R_{2} + 0.109 * R_{3} R 1 ∗ = R 1 − 0.3125 ∗ R 3 ; R 2 ∗ = R 2 + 0.109 ∗ R 3
( 1 0 0 − 0.435 1.304 − 0.174 0 1 0 1.65 − 0.957 0.061 0 0 1 − 0.609 − 0.174 0.556 ) ; \left( \begin{array}{cccc} 1 & 0 & 0 & -0.435 & 1.304 & -0.174 \\ 0 & 1 & 0 & 1.65 & -0.957 & 0.061 \\ 0 & 0 & 1 & -0.609 & -0.174 & 0.556 \end{array} \right); ⎝ ⎛ 1 0 0 0 1 0 0 0 1 − 0.435 1.65 − 0.609 1.304 − 0.957 − 0.174 − 0.174 0.061 0.556 ⎠ ⎞ ;
Answer: Z = − 0.365 Z = -0.365 Z = − 0.365 ; Y = 1.4907 Y = 1.4907 Y = 1.4907 ; X = − 0.2607 X = -0.2607 X = − 0.2607 ;
( − 0.435 1.304 − 0.174 1.65 − 0.957 0.061 − 0.609 − 0.174 0.556 ) . \left( \begin{array}{cccc} -0.435 & 1.304 & -0.174 \\ 1.65 & -0.957 & 0.061 \\ -0.609 & -0.174 & 0.556 \end{array} \right). ⎝ ⎛ − 0.435 1.65 − 0.609 1.304 − 0.957 − 0.174 − 0.174 0.061 0.556 ⎠ ⎞ .
Answer provided by https://www.AssignmentExpert.com