Question #81954

Solve the following
0.6x + 0.8y + 0.1z = 1
1.1 x + 0.4y + 0.3z= 0.2
x + y + 2z =0.5
by LU decomposition method and find the inverse of the coefficient matrix

Expert's answer

Answer on Question #81954 — Math — Linear Algebra

Question

Solve the following


0.6x+0.8y+0.1z=10.6x + 0.8y + 0.1z = 11.1x+0.4y+0.3z=0.21.1x + 0.4y + 0.3z = 0.2x+y+2z=0.5x + y + 2z = 0.5


by LU decomposition method and find the inverse of the coefficient matrix

Solution


A=(0.60.80.11.10.40.3112);b=(10.20.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}


Swap some rows in matrices.


A=(1121.10.40.30.60.80.1);b=(0.50.21)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}


LU = A


L=(100101);U=(101001);L = \begin{pmatrix} 1 & 0 & 0 \\ * & 1 & 0 \\ * & * & 1 \end{pmatrix}; \quad U = \begin{pmatrix} 1 & * & * \\ 0 & 1 & * \\ 0 & 0 & 1 \end{pmatrix};(1121.10.40.30.60.80.1)(0.50.21)R2=R2+(1.1)R1;R3=R3+(0.6)R1;\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(11200.71.900.21.1)(0.50.350.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};U=(11200.71.9001.643);b=(0.50.350.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}L=(1001.1100.61)L = \begin{pmatrix} 1 & 0 & 0 \\ 1.1 & 1 & 0 \\ 0.6 & * & 1 \end{pmatrix}(11200.71.900.21.1)(0.50.350.7)R3=R3+(0.20.7)R2(11200.71.9001.643)(0.50.350.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}L=(1001.1100.60.28571)L = \begin{pmatrix} 1 & 0 & 0 \\ 1.1 & 1 & 0 \\ 0.6 & -0.2857 & 1 \end{pmatrix}U=(11200.71.9001.643);b=(0.50.350.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}{x+y+2z=0.50.7y1.9z=0.351.643z=0.6\begin{cases} x + y + 2z = 0.5 \\ -0.7y - 1.9z = -0.35 \\ -1.643z = 0.6 \end{cases}


So, we have: Z=0.6/(1.643)=0.365Z = 0.6 / (-1.643) = -0.365

Y=0.35+1.9(0.365)0.7=1.4907Y = \frac {- 0.35 + 1.9 * (- 0.365)}{- 0.7} = 1.4907X=0.52(0.365)1.4907=0.2607X = 0.5 - 2 * (-0.365) - 1.4907 = -0.2607


Find the inverse of the coefficient matrix


(0.60.80.11001.10.40.3010112001);\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);


1. R1=R1/0.6R_{1}^{*} = R_{1} / 0.6

(11.330.1671.67001.10.40.3010112001);\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);


2. R2=R21.1R1;R3=R3R1R_{2}^{*} = R_{2} - 1.1 * R_{1}; R_{3}^{*} = R_{3} - R_{1}

(11.330.1671.670001.0670.11671.831000.331.831.6701);\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);


3. R2=R2/(1.067)R_{2}^{*} = R_{2} / (-1.067)

(11.330.1671.6700010.1091.7190.9375000.331.831.6701);\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);


4. R1=R11.33R2;R3=R3+0.33R2R_{1}^{*} = R_{1} - 1.33 * R_{2}; R_{3}^{*} = R_{3} + 0.33 * R_{2}

(100.31250.6251.250010.1091.7190.93750001.7971.0940.31251);\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);


5. R3=R31.797R_{3}^{*} = \frac{R_{3}}{1.797}

(100.31250.6251.250010.1091.7190.937500010.6090.1740.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);


6. R1=R10.3125R3;R2=R2+0.109R3R_{1}^{*} = R_{1} - 0.3125 * R_{3}; R_{2}^{*} = R_{2} + 0.109 * R_{3}

(1000.4351.3040.1740101.650.9570.0610010.6090.1740.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);


Answer: Z=0.365Z = -0.365; Y=1.4907Y = 1.4907; X=0.2607X = -0.2607;


(0.4351.3040.1741.650.9570.0610.6090.1740.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).


Answer provided by https://www.AssignmentExpert.com

Need a fast expert's response?

Submit order

and get a quick answer at the best price

for any assignment or question with DETAILED EXPLANATIONS!

LATEST TUTORIALS
APPROVED BY CLIENTS