Section 7.2

Diagonalization by Similarity Transformations

[27/100%]
Show
LaTeX
Problem 7.2.1

Diagonalize A=(−8−61210)\mathbf{A} = \begin{pmatrix} -8 & -6 \\ 12 & 10 \end{pmatrix} with a similarity transformation, or else explain why A\mathbf{A} can't be diagonalized.

?
Problem 7.2.2
?
(a)

Verify that alg mult⁡A(λ)=geo mult⁡A(λ)\operatorname {alg\, mult}_{\mathbf{A}}(\lambda ) = \operatorname {geo\, mult}_{\mathbf{A}}(\lambda ) for each eigenvalue of

A=(−4−3−30−10665). \mathbf{A} = \begin{pmatrix} -4 & -3 & -3 \\ 0 & -1 & 0 \\ 6 & 6 & 5 \end{pmatrix}.
(b)

Find a nonsingular P\mathbf{P} such that P−1AP\mathbf{P}^{-1}\mathbf{A}\mathbf{P} is a diagonal matrix.

Problem 7.2.3

Show that similar matrices need not have the same eigenvectors by giving an example of two matrices that are similar but have different eigenspaces.

?
Problem 7.2.4

λ=2\lambda =2 is an eigenvalue for A=(321020−2−30)\mathbf{A} = \begin{pmatrix} 3 & 2 & 1 \\ 0 & 2 & 0 \\ -2 & -3 & 0 \end{pmatrix}. Find alg mult⁡A(λ)\operatorname {alg\, mult}_{\mathbf{A}}(\lambda ) as well as geo mult⁡A(λ)\operatorname {geo\, mult}_{\mathbf{A}}(\lambda ). Can you conclude anything about the diagonalizability of A\mathbf{A} from these results?

?
Problem 7.2.5

If B=P−1AP\mathbf{B} = \mathbf{P}^{-1}\mathbf{A}\mathbf{P}, explain why Bk=P−1AkP\mathbf{B}^k = \mathbf{P}^{-1}\mathbf{A}^k\mathbf{P}.

?
Problem 7.2.6

Compute lim⁡n→∞An\lim_{n\to \infty } \mathbf{A}^n for A=(7/51/5−11/2)\mathbf{A} = \begin{pmatrix} 7/5 & 1/5 \\ -1 & 1/2 \end{pmatrix}.

?
Problem 7.2.7

Let {x→1,x→2,…,x→t}\left\{ \overrightarrow {x}_1, \overrightarrow {x}_2, \ldots , \overrightarrow {x}_t\right\} be a set of linearly independent eigenvectors for An×n\mathbf{A}_{n\times n} associated with respective eigenvalues {λ1,λ2,…,λt}\left\{ \lambda_1, \lambda_2, \ldots , \lambda_t\right\}, and let X\mathbf{X} be any n×(n−t)n \times (n-t) matrix such that Pn×n=(x→1∣⋯∣x→t∣X)\mathbf{P}_{n\times n} = (\overrightarrow {x}_1 \mid \cdots \mid \overrightarrow {x}_t \mid \mathbf{X}) is nonsingular. Prove that if P−1=(y→1∗⋮y→t∗Y∗)\mathbf{P}^{-1} = \begin{pmatrix} \overrightarrow {y}_1^{*} \\ \vdots \\ \overrightarrow {y}_t^{*} \\ \mathbf{Y}^{*} \end{pmatrix}, where the y→i∗\overrightarrow {y}_i^{*}'s are rows and Y∗\mathbf{Y}^{*} is (n−t)×n(n-t) \times n, then {y→1∗,y→2∗,…,y→t∗}\left\{ \overrightarrow {y}_1^{*}, \overrightarrow {y}_2^{*}, \ldots , \overrightarrow {y}_t^{*}\right\} is a set of linearly independent left-hand eigenvectors associated with {λ1,λ2,…,λt}\left\{ \lambda_1, \lambda_2, \ldots , \lambda_t\right\}, respectively (i.e., y→i∗A=λiy→i∗\overrightarrow {y}_i^{*}\mathbf{A} = \lambda_i \overrightarrow {y}_i^{*}).

?
Problem 7.2.8

Let A\mathbf{A} be a diagonalizable matrix, and let ρ(⋆)\rho (\star ) denote the spectral radius. Prove that lim⁡k→∞Ak=0\lim_{k\to \infty } \mathbf{A}^k = \mathbf{0} if and only if ρ(A)<1\rho (\mathbf{A}) < 1.

Note: This result also holds for nondiagonalizable matrices, as demonstrated elsewhere in the text.

?
Problem 7.2.9

Apply the technique used to prove Schur's triangularization theorem to construct an orthogonal matrix P\mathbf{P} such that PTAP\mathbf{P}^T\mathbf{A}\mathbf{P} is upper triangular for A=(13−916−11)\mathbf{A} = \begin{pmatrix} 13 & -9 \\ 16 & -11 \end{pmatrix}.

?
Problem 7.2.10

Verify the Cayley--Hamilton theorem for A=(1−4−48−11−8−885)\mathbf{A} = \begin{pmatrix} 1 & -4 & -4 \\ 8 & -11 & -8 \\ -8 & 8 & 5 \end{pmatrix}.

Hint: This is the matrix from Example 7.2.1 on p. 507.

?
Problem 7.2.11

Since each row sum in the following symmetric matrix A\mathbf{A} is 44, it's clear that x→=(1,1,1,1)T\overrightarrow {x} = (1,1,1,1)^T is both a right-hand and left-hand eigenvector associated with λ=4∈σ(A)\lambda =4 \in \sigma (\mathbf{A}). Use the deflation technique (rescale x→\overrightarrow {x} to unit length y→\overrightarrow {y}, build an orthogonal matrix P=[y→∣X]\mathbf{P} = [\overrightarrow {y} \mid \mathbf{X}] having y→\overrightarrow {y} as its first column, and then the remaining eigenvalues of A\mathbf{A} are exactly the eigenvalues of the smaller matrix XTAX\mathbf{X}^T\mathbf{A}\mathbf{X}) to determine the remaining eigenvalues of

A=(1021021121101102). \mathbf{A} = \begin{pmatrix} 1 & 0 & 2 & 1 \\ 0 & 2 & 1 & 1 \\ 2 & 1 & 1 & 0 \\ 1 & 1 & 0 & 2 \end{pmatrix}.
?
Problem 7.2.12

Explain why AGi=GiA=λiGi\mathbf{A}\mathbf{G}_i = \mathbf{G}_i\mathbf{A} = \lambda_i \mathbf{G}_i for the spectral projector Gi\mathbf{G}_i associated with the eigenvalue λi\lambda_i of a diagonalizable matrix A\mathbf{A}.

?
Problem 7.2.13

Prove that A=c→n×1d→1×nT\mathbf{A} = \overrightarrow {c}_{n\times 1}\overrightarrow {d}_{1\times n}^T is diagonalizable if and only if d→Tc→≠0\overrightarrow {d}^T\overrightarrow {c} \ne 0.

?
Problem 7.2.14

Prove that A=(W00Z)\mathbf{A} = \begin{pmatrix} \mathbf{W} & \mathbf{0} \\ \mathbf{0} & \mathbf{Z} \end{pmatrix} is diagonalizable if and only if Ws×s\mathbf{W}_{s\times s} and Zt×t\mathbf{Z}_{t\times t} are each diagonalizable.

?
Problem 7.2.15

Prove that if AB=BA\mathbf{A}\mathbf{B} = \mathbf{B}\mathbf{A}, then A\mathbf{A} and B\mathbf{B} can be simultaneously triangularized by a unitary similarity transformation---i.e., U∗AU=T1\mathbf{U}^{*}\mathbf{A}\mathbf{U} = \mathbf{T}_1 and U∗BU=T2\mathbf{U}^{*}\mathbf{B}\mathbf{U} = \mathbf{T}_2 for some unitary matrix U\mathbf{U}.

Hint: Recall Exercise 7.1.20 along with the development of Schur's triangularization theorem.

?
Problem 7.2.16

For diagonalizable matrices, prove that AB=BA\mathbf{A}\mathbf{B} = \mathbf{B}\mathbf{A} if and only if A\mathbf{A} and B\mathbf{B} can be simultaneously diagonalized---i.e., P−1AP=D1\mathbf{P}^{-1}\mathbf{A}\mathbf{P} = \mathbf{D}_1 and P−1BP=D2\mathbf{P}^{-1}\mathbf{B}\mathbf{P} = \mathbf{D}_2 for some P\mathbf{P}.

Hint: If A\mathbf{A} and B\mathbf{B} commute, then so do P−1AP=(λ1I00D)\mathbf{P}^{-1}\mathbf{A}\mathbf{P} = \begin{pmatrix} \lambda_1\mathbf{I} & \mathbf{0} \\ \mathbf{0} & \mathbf{D} \end{pmatrix} and P−1BP=(WXYZ)\mathbf{P}^{-1}\mathbf{B}\mathbf{P} = \begin{pmatrix} \mathbf{W} & \mathbf{X} \\ \mathbf{Y} & \mathbf{Z} \end{pmatrix}.

?
Problem 7.2.17

Explain why the following ``proof'' of the Cayley--Hamilton theorem is not valid: p(λ)=det⁡(A−λI)  ⟹  p(A)=det⁡(A−AI)=det⁡(0)=0p(\lambda ) = \operatorname {det}\left(\mathbf{A}-\lambda \mathbf{I}\right) \implies p(\mathbf{A}) = \operatorname {det}\left(\mathbf{A}-\mathbf{A}\mathbf{I}\right) = \operatorname {det}\left(\mathbf{0}\right) = 0.

?
Problem 7.2.18

Show that the eigenvalues of the finite difference matrix

A=(2−1−12−1⋱⋱⋱−12−1−12)n×n \mathbf{A} = \begin{pmatrix} 2 & -1 & & & \\ -1 & 2 & -1 & & \\ & \ddots & \ddots & \ddots & \\ & & -1 & 2 & -1 \\ & & & -1 & 2 \end{pmatrix}_{n\times n}

are λj=4sin⁡2jπ2(n+1)\lambda_j = 4\sin^2\dfrac {j\pi }{2(n+1)}, 1≤j≤n1 \le j \le n.

Example 7.2.5 establishes that the eigenvalues of this finite difference matrix are λj=2−2cos⁡jπn+1\lambda_j = 2 - 2\cos \dfrac {j\pi }{n+1}, 1≤j≤n1 \le j \le n.

?
Problem 7.2.19

Let N=(01⋱⋱⋱10)n×n\mathbf{N} = \begin{pmatrix} 0 & 1 & & \\ & \ddots & \ddots & \\ & & \ddots & 1 \\ & & & 0 \end{pmatrix}_{n\times n}.

?
(a)

Show that λ∈σ(N+NT)\lambda \in \sigma (\mathbf{N}+\mathbf{N}^T) if and only if iλ∈σ(N−NT)\mathrm{i}\lambda \in \sigma (\mathbf{N}-\mathbf{N}^T).

Example 7.2.5 shows that the eigenvalues of N+NT\mathbf{N}+\mathbf{N}^T are λj=2cos⁡(jπ/(n+1))\lambda_j = 2\cos (j\pi /(n+1)), 1≤j≤n1 \le j \le n.

(b)

Explain why N+NT\mathbf{N}+\mathbf{N}^T is nonsingular if and only if nn is even.

(c)

Evaluate det⁡(N−NT)/det⁡(N+NT)\operatorname {det}\left(\mathbf{N}-\mathbf{N}^T\right)/\operatorname {det}\left(\mathbf{N}+\mathbf{N}^T\right) when nn is even.

Problem 7.2.20

A Toeplitz matrix having the form

C=(c0cn−1cn−2⋯c1c1c0cn−1⋯c2c2c1c0⋯c3⋮⋮⋮⋱⋮cn−1cn−2cn−3⋯c0)n×n \mathbf{C} = \begin{pmatrix} c_0 & c_{n-1} & c_{n-2} & \cdots & c_1 \\ c_1 & c_0 & c_{n-1} & \cdots & c_2 \\ c_2 & c_1 & c_0 & \cdots & c_3 \\ \vdots & \vdots & \vdots & \ddots & \vdots \\ c_{n-1} & c_{n-2} & c_{n-3} & \cdots & c_0 \end{pmatrix}_{n\times n}

is called a circulant matrix. If p(x)=c0+c1x+⋯+cn−1xn−1p(x) = c_0+c_1 x+\cdots +c_{n-1}x^{n-1}, and if {1,ξ,ξ2,…,ξn−1}\left\{ 1,\xi ,\xi^2,\ldots ,\xi^{n-1}\right\} are the nnth roots of unity, then the results of Exercise 5.8.12 (p. 379) insure that

FnCFn−1=(p(1)0⋯00p(ξ)⋯0⋮⋮⋱⋮00⋯p(ξn−1)) \mathbf{F}_n\mathbf{C}\mathbf{F}_n^{-1} = \begin{pmatrix} p(1) & 0 & \cdots & 0 \\ 0 & p(\xi ) & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & p(\xi ^{n-1}) \end{pmatrix}

in which Fn\mathbf{F}_n is the Fourier matrix of order nn. Verify these facts for the circulant below by computing its eigenvalues and eigenvectors directly:

C=(1010010110100101). \mathbf{C} = \begin{pmatrix} 1 & 0 & 1 & 0 \\ 0 & 1 & 0 & 1 \\ 1 & 0 & 1 & 0 \\ 0 & 1 & 0 & 1 \end{pmatrix}.
?
Problem 7.2.21

Suppose that (λ,x→)(\lambda , \overrightarrow {x}) and (μ,y→∗)(\mu , \overrightarrow {y}^{*}) are right-hand and left-hand eigenpairs for A∈Rn×n\mathbf{A} \in \mathbb {R}^{n\times n}---i.e., Ax→=λx→\mathbf{A}\overrightarrow {x}=\lambda \overrightarrow {x} and y→∗A=μy→∗\overrightarrow {y}^{*}\mathbf{A}=\mu \overrightarrow {y}^{*}. Explain why y→∗x→=0\overrightarrow {y}^{*}\overrightarrow {x}=0 whenever λ≠μ\lambda \ne \mu.

?
Problem 7.2.22

Consider A∈Rn×n\mathbf{A} \in \mathbb {R}^{n\times n}.

?
(a)

Show that if A\mathbf{A} is diagonalizable, then there are right-hand and left-hand eigenvectors x→\overrightarrow {x} and y→∗\overrightarrow {y}^{*} associated with λ∈σ(A)\lambda \in \sigma (\mathbf{A}) such that y→∗x→≠0\overrightarrow {y}^{*}\overrightarrow {x} \ne 0 so that we can make y→∗x→=1\overrightarrow {y}^{*}\overrightarrow {x}=1.

(b)

Show that not every right-hand and left-hand eigenvector x→\overrightarrow {x} and y→∗\overrightarrow {y}^{*} associated with λ∈σ(A)\lambda \in \sigma (\mathbf{A}) must satisfy y→∗x→≠0\overrightarrow {y}^{*}\overrightarrow {x} \ne 0.

(c)

Show that (a) need not be true when A\mathbf{A} is not diagonalizable.

Problem 7.2.23

Consider A∈Rn×n\mathbf{A} \in \mathbb {R}^{n\times n} with λ∈σ(A)\lambda \in \sigma (\mathbf{A}).

?
(a)

Prove that if λ\lambda is simple, then y→∗x→≠0\overrightarrow {y}^{*}\overrightarrow {x} \ne 0 for every pair of respective right-hand and left-hand eigenvectors x→\overrightarrow {x} and y→∗\overrightarrow {y}^{*} associated with λ\lambda regardless of whether or not A\mathbf{A} is diagonalizable.

Hint: Use the core-nilpotent decomposition (p. 397).

(b)

Show that y→∗x→=0\overrightarrow {y}^{*}\overrightarrow {x}=0 is possible when λ\lambda is not simple.

Problem 7.2.24

For A∈Rn×n\mathbf{A} \in \mathbb {R}^{n\times n} with σ(A)={λ1,λ2,…,λk}\sigma (\mathbf{A}) = \left\{ \lambda_1,\lambda_2,\ldots ,\lambda_k\right\}, show A\mathbf{A} is diagonalizable if and only if Rn=N(A−λ1I)⊕N(A−λ2I)⊕⋯⊕N(A−λkI)\mathbb {R}^{n} = N(\mathbf{A}-\lambda_1\mathbf{I}) \oplus N(\mathbf{A}-\lambda_2\mathbf{I}) \oplus \cdots \oplus N(\mathbf{A}-\lambda_k\mathbf{I}).

Hint: Recall Exercise 5.9.14.

?
Problem 7.2.25

Schur's triangularization theorem insures that every square matrix A\mathbf{A} is unitarily similar to an upper-triangular matrix---say, U∗AU=T\mathbf{U}^{*}\mathbf{A}\mathbf{U}=\mathbf{T}. But even when A\mathbf{A} is real, U\mathbf{U} and T\mathbf{T} may have to be complex if A\mathbf{A} has some complex eigenvalues. However, the matrices (and the arithmetic) can be constrained to be real by settling for a block-triangular result with 2×22\times 2 or scalar entries on the diagonal. Prove that for each A∈Rn×n\mathbf{A} \in \mathbb {R}^{n\times n} there exists an orthogonal matrix P∈Rn×n\mathbf{P} \in \mathbb {R}^{n\times n} and real matrices Bij\mathbf{B}_{ij} such that

PTAP=(B11B12⋯B1k0B22⋯B2k⋮⋮⋱⋮00⋯Bkk),where Bjj is 1×1 or 2×2. \mathbf{P}^T\mathbf{A}\mathbf{P} = \begin{pmatrix} \mathbf{B}_{11} & \mathbf{B}_{12} & \cdots & \mathbf{B}_{1k} \\ \mathbf{0} & \mathbf{B}_{22} & \cdots & \mathbf{B}_{2k} \\ \vdots & \vdots & \ddots & \vdots \\ \mathbf{0} & \mathbf{0} & \cdots & \mathbf{B}_{kk} \end{pmatrix}, \quad \text{where } \mathbf{B}_{jj} \text{ is } 1\times 1 \text{ or } 2\times 2.

If Bjj=[λj]\mathbf{B}_{jj} = [\lambda_j] is 1×11\times 1, then λj∈σ(A)\lambda_j \in \sigma (\mathbf{A}), and if Bjj\mathbf{B}_{jj} is 2×22\times 2, then σ(Bjj)={λj,λj‾}⊆σ(A)\sigma (\mathbf{B}_{jj}) = \left\{ \lambda_j, \overline{\lambda_j}\right\} \subseteq \sigma (\mathbf{A}).

?
Problem 7.2.26

When A∈Rn×n\mathbf{A} \in \mathbb {R}^{n\times n} is diagonalizable by a similarity transformation S\mathbf{S}, then S\mathbf{S} may have to be complex if A\mathbf{A} has some complex eigenvalues. Analogous to Exercise 7.2.25, we can stay in the realm of real numbers by settling for a block-diagonal result with 1×11\times 1 or 2×22\times 2 entries on the diagonal. Prove that if A∈Rn×n\mathbf{A} \in \mathbb {R}^{n\times n} is diagonalizable with real eigenvalues {ρ1,…,ρr}\left\{ \rho_1,\ldots ,\rho_r\right\} and complex eigenvalues {λ1,λ1‾,λ2,λ2‾,…,λt,λt‾}\left\{ \lambda_1,\overline{\lambda_1},\lambda_2,\overline{\lambda_2},\ldots ,\lambda_t,\overline{\lambda_t}\right\} with 2t+r=n2t+r=n, then there exists a nonsingular P∈Rn×n\mathbf{P} \in \mathbb {R}^{n\times n} and Bj\mathbf{B}_j's ∈R2×2\in \mathbb {R}^{2\times 2} such that

P−1AP=(D0⋯00B1⋯0⋮⋮⋱⋮00⋯Bt),where D=(ρ10⋯00ρ2⋯0⋮⋮⋱⋮00⋯ρr), \mathbf{P}^{-1}\mathbf{A}\mathbf{P} = \begin{pmatrix} \mathbf{D} & \mathbf{0} & \cdots & \mathbf{0} \\ \mathbf{0} & \mathbf{B}_1 & \cdots & \mathbf{0} \\ \vdots & \vdots & \ddots & \vdots \\ \mathbf{0} & \mathbf{0} & \cdots & \mathbf{B}_t \end{pmatrix}, \quad \text{where } \mathbf{D} = \begin{pmatrix} \rho _1 & 0 & \cdots & 0 \\ 0 & \rho _2 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & \rho _r \end{pmatrix},

and where Bj\mathbf{B}_j has eigenvalues λj\lambda_j and λj‾\overline{\lambda_j}.

?
Problem 7.2.27

For A∈Cn×n\mathbf{A} \in \mathcal{C}^{n\times n}, prove that x→∗Ax→=0\overrightarrow {x}^{*}\mathbf{A}\overrightarrow {x}=0 for all x→∈Cn×1  ⟹  A=0\overrightarrow {x} \in \mathcal{C}^{n\times 1} \implies \mathbf{A}=\mathbf{0}. Show that x→TAx→=0\overrightarrow {x}^T\mathbf{A}\overrightarrow {x}=0 for all x→∈Rn×1\centernot  ⟹  A=0\overrightarrow {x} \in \mathbb {R}^{n\times 1} \centernot \implies \mathbf{A}=\mathbf{0}, even if A\mathbf{A} is real.

?