Question #52680

For the following matrices, check whether the reexists an invertible matrix P such that P−1 A Pis diagonal. When such a P exists,find P.
i)A=( 0 1 -3)
( 2 -1 6)
( 1 -1 4)
(ii) B is equal to. ( 0 0 -2)
( 1 2 1)
( 1 0 3)

b)Find the inverse of the matrix B in part (a) by using Cayley-Hamilton theorem.(3)c)Using the fact that det(AB)=det(A)det(B)for any two matrices A and B,prove the identity (a2+b2) (c2+d2)=(ac−bd)2+(ad+bc)2

Expert's answer

Answer on Question #52680 – Math – Linear Algebra

a) For the following matrices, check whether there exists an invertible matrix PP such that P−1APP^{-1}AP is diagonal. When such a PP exists, find PP.

(i) A=(01−32−161−14)A = \begin{pmatrix} 0 & 1 & -3 \\ 2 & -1 & 6 \\ 1 & -1 & 4 \end{pmatrix}

(ii) B=(00−2121103)B = \begin{pmatrix} 0 & 0 & -2 \\ 1 & 2 & 1 \\ 1 & 0 & 3 \end{pmatrix}

b) Find the inverse of the matrix BB in part (a) by using Cayley-Hamilton theorem.

c) Using the fact that det⁡(AB)=det⁡(A)det⁡(B)\det(AB) = \det(A)\det(B) for any two matrices AA and BB, prove the identity (a2+b2)(c2+d2)=(ac−bd)2+(ad+bc)2(a^2 + b^2)(c^2 + d^2) = (ac - bd)^2 + (ad + bc)^2

Solution

a) Such matrix PP exists if there are no zero vectors among eigenvectors of given matrix (AA or BB)

(i) The columns of matrix PP are given by the eigenvectors of the matrix AA. Let's first find its eigenvalues. The characteristic polynomial of AA:


p(λ)=det⁡(λI3−A)=det⁡(λ−13−2λ+1−6−11λ−4)=λ3−3λ2+7λ−5p(\lambda) = \det(\lambda I_3 - A) = \det \begin{pmatrix} \lambda & -1 & 3 \\ -2 & \lambda + 1 & -6 \\ -1 & 1 & \lambda - 4 \end{pmatrix} = \lambda^3 - 3\lambda^2 + 7\lambda - 5


Since eigenvalues satisfy the equation p(λ)=0p(\lambda) = 0, we obtain the following eigenvalues


λ1,2,3=1\lambda_{1,2,3} = 1


The eigenvectors are

λ=1\lambda = 1:


e1=(110),e2=(30−1),e3=(000)e_1 = \begin{pmatrix} 1 \\ 1 \\ 0 \end{pmatrix}, \qquad e_2 = \begin{pmatrix} 3 \\ 0 \\ -1 \end{pmatrix}, \qquad e_3 = \begin{pmatrix} 0 \\ 0 \\ 0 \end{pmatrix}


Therefore there is no such matrix for AA.

(ii) The columns of matrix PP are given by the eigenvectors of the matrix BB. Let's first find its eigenvalues. The characteristic polynomial of BB:


p(λ)=det⁡(λI3−B)=det⁡(λ02−1λ−2−1−10λ−3)=λ3−3λ2+7λ−5p(\lambda) = \det(\lambda I_3 - B) = \det \begin{pmatrix} \lambda & 0 & 2 \\ -1 & \lambda - 2 & -1 \\ -1 & 0 & \lambda - 3 \end{pmatrix} = \lambda^3 - 3\lambda^2 + 7\lambda - 5


Since eigenvalues satisfy the equation p(λ)=0p(\lambda) = 0, we obtain the following eigenvalues


λ1,2,3=1,2,2\lambda_{1,2,3} = 1, 2, 2


The eigenvectors are


λ=1:\lambda = 1:e1=(−211)e_1 = \begin{pmatrix} -2 \\ 1 \\ 1 \end{pmatrix}λ=2:\lambda = 2:e2=(−101),e3=(010)e_2 = \begin{pmatrix} -1 \\ 0 \\ 1 \end{pmatrix}, \qquad e_3 = \begin{pmatrix} 0 \\ 1 \\ 0 \end{pmatrix}


Therefore


P=(−2−10101110)P = \begin{pmatrix} -2 & -1 & 0 \\ 1 & 0 & 1 \\ 1 & 1 & 0 \end{pmatrix}


b) The characteristic polynomial of BB is given by


p(λ)=det⁡(λI3−B)p(\lambda) = \det(\lambda I_3 - B)


The Cayley-Hamilton theorem states that


p(B)=0p(B) = 0


Since


p(λ)=∣λ02−1λ−2−1−10λ−3∣=λ3−5λ2+8λ−4,p(\lambda) = \begin{vmatrix} \lambda & 0 & 2 \\ -1 & \lambda - 2 & -1 \\ -1 & 0 & \lambda - 3 \end{vmatrix} = \lambda^3 - 5\lambda^2 + 8\lambda - 4,


we obtain


p(B)=B3−5B2+8B−4I3=0p(B) = B^3 - 5B^2 + 8B - 4I_3 = 0


Multiplying by B−1B^{-1} we obtain the following


B2−5B+8I3−4B−1=0B^2 - 5B + 8I_3 - 4B^{-1} = 0


Therefore


B−1=14(B2−5B+8I3)B^{-1} = \frac{1}{4}(B^2 - 5B + 8I_3)


Since


B2=(00−2121103)(00−2121103)=(−20−6343307),B ^ {2} = \left( \begin{array}{c c c} 0 & 0 & - 2 \\ 1 & 2 & 1 \\ 1 & 0 & 3 \end{array} \right) \left( \begin{array}{c c c} 0 & 0 & - 2 \\ 1 & 2 & 1 \\ 1 & 0 & 3 \end{array} \right) = \left( \begin{array}{c c c} - 2 & 0 & - 6 \\ 3 & 4 & 3 \\ 3 & 0 & 7 \end{array} \right),


we obtain


B−1=14((−20−6343307)−5(00−2121103)+8(100010001))==14(604−22−2−200)=(3/201−1/21/2−1/2−1/200)\begin{array}{l} B ^ {- 1} = \frac {1}{4} \left(\left( \begin{array}{c c c} - 2 & 0 & - 6 \\ 3 & 4 & 3 \\ 3 & 0 & 7 \end{array} \right) - 5 \left( \begin{array}{c c c} 0 & 0 & - 2 \\ 1 & 2 & 1 \\ 1 & 0 & 3 \end{array} \right) + 8 \left( \begin{array}{c c c} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{array} \right)\right) = \\ = \frac {1}{4} \left( \begin{array}{c c c} 6 & 0 & 4 \\ - 2 & 2 & - 2 \\ - 2 & 0 & 0 \end{array} \right) = \left( \begin{array}{c c c} 3 / 2 & 0 & 1 \\ - 1 / 2 & 1 / 2 & - 1 / 2 \\ - 1 / 2 & 0 & 0 \end{array} \right) \\ \end{array}


c) Let's consider two matrices:


A=(abdc),B=(adbc)A = \left( \begin{array}{c c} a & b \\ d & c \end{array} \right), \qquad B = \left( \begin{array}{c c} a & d \\ b & c \end{array} \right)


Then


AB=(abdc)(adbc)=(a2+b2ad+bcad+bcc2+d2)A B = \left( \begin{array}{c c} a & b \\ d & c \end{array} \right) \left( \begin{array}{c c} a & d \\ b & c \end{array} \right) = \left( \begin{array}{c c} a ^ {2} + b ^ {2} & a d + b c \\ a d + b c & c ^ {2} + d ^ {2} \end{array} \right)


Since


det⁡(AB)=det⁡(A)det⁡(B),\det (A B) = \det (A) \det (B),


we obtain


det⁡(a2+b2ad+bcad+bcc2+d2)=det⁡(abdc)det⁡(adbc)\det \left( \begin{array}{c c} a ^ {2} + b ^ {2} & a d + b c \\ a d + b c & c ^ {2} + d ^ {2} \end{array} \right) = \det \left( \begin{array}{c c} a & b \\ d & c \end{array} \right) \det \left( \begin{array}{c c} a & d \\ b & c \end{array} \right)(a2+b2)(c2+d2)−(ad+bc)2=(ac−bd)(ac−bd)(a ^ {2} + b ^ {2}) (c ^ {2} + d ^ {2}) - (a d + b c) ^ {2} = (a c - b d) (a c - b d)(a2+b2)(c2+d2)=(ac−bd)2+(ad+bc)2(a ^ {2} + b ^ {2}) (c ^ {2} + d ^ {2}) = (a c - b d) ^ {2} + (a d + b c) ^ {2}


**Answer:**

a)

(i) Such a matrix does not exist

(ii) (−2−10101110)\left( \begin{array}{ccc} - 2 & -1 & 0\\ 1 & 0 & 1\\ 1 & 1 & 0 \end{array} \right)

b) (3/201−1/21/2−1/2−1/200)\left( \begin{array}{ccc}3 / 2 & 0 & 1\\ -1 / 2 & 1 / 2 & -1 / 2\\ -1 / 2 & 0 & 0 \end{array} \right)

LATEST TUTORIALS
APPROVED BY CLIENTS