Section 7.1

Elementary Properties of Eigensystems

[24/100%]
Show
LaTeX
Problem 7.1.1

Determine the eigenvalues and eigenvectors for the following matrices.

A=(−10−71411),B=(21684148−8−32−18),C=(3−250140−15), \mathbf{A} = \begin{pmatrix} -10 & -7 \\ 14 & 11 \end{pmatrix}, \quad \mathbf{B} = \begin{pmatrix} 2 & 16 & 8 \\ 4 & 14 & 8 \\ -8 & -32 & -18 \end{pmatrix}, \quad \mathbf{C} = \begin{pmatrix} 3 & -2 & 5 \\ 0 & 1 & 4 \\ 0 & -1 & 5 \end{pmatrix}, D=(063−151−124),E=(300030003). \mathbf{D} = \begin{pmatrix} 0 & 6 & 3 \\ -1 & 5 & 1 \\ -1 & 2 & 4 \end{pmatrix}, \quad \mathbf{E} = \begin{pmatrix} 3 & 0 & 0 \\ 0 & 3 & 0 \\ 0 & 0 & 3 \end{pmatrix}.

Which, if any, are deficient in eigenvectors in the sense that there fails to exist a complete linearly independent set?

?
Problem 7.1.2

Without doing an eigenvalue--eigenvector computation, determine which of the following are eigenvectors for

A=(−9−6−2−4−8−6−3−12015853221712), \mathbf{A} = \begin{pmatrix} -9 & -6 & -2 & -4 \\ -8 & -6 & -3 & -1 \\ 20 & 15 & 8 & 5 \\ 32 & 21 & 7 & 12 \end{pmatrix},

and for those which are eigenvectors, identify the associated eigenvalue.

?
(a)

(−1,1,0,1)T(-1,1,0,1)^T.

(b)

(1,0,−1,0)T(1,0,-1,0)^T.

(c)

(−1,0,2,2)T(-1,0,2,2)^T.

(d)

(0,1,−3,0)T(0,1,-3,0)^T.

Problem 7.1.3

Explain why the eigenvalues of triangular and diagonal matrices

T=(t11t12⋯t1n0t22⋯t2n⋮⋮⋱⋮00⋯tnn)andD=(λ10⋯00λ2⋯0⋮⋮⋱⋮00⋯λn) \mathbf{T} = \begin{pmatrix} t_{11} & t_{12} & \cdots & t_{1n} \\ 0 & t_{22} & \cdots & t_{2n} \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & t_{nn} \end{pmatrix} \quad \text{and} \quad \mathbf{D} = \begin{pmatrix} \lambda _1 & 0 & \cdots & 0 \\ 0 & \lambda _2 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & \lambda _n \end{pmatrix}

are simply the diagonal entries---the tiit_{ii}'s and λi\lambda_i's.

?
Problem 7.1.4

For T=(AB0C)\mathbf{T} = \begin{pmatrix} \mathbf{A} & \mathbf{B} \\ \mathbf{0} & \mathbf{C} \end{pmatrix}, prove det⁡(T−λI)=det⁡(A−λI)det⁡(C−λI)\operatorname {det}\left(\mathbf{T}-\lambda \mathbf{I}\right) = \operatorname {det}\left(\mathbf{A}-\lambda \mathbf{I}\right)\operatorname {det}\left(\mathbf{C}-\lambda \mathbf{I}\right) to conclude that σ(AB0C)=σ(A)∪σ(C)\sigma \begin{pmatrix} \mathbf{A} & \mathbf{B} \\ \mathbf{0} & \mathbf{C} \end{pmatrix} = \sigma (\mathbf{A}) \cup \sigma (\mathbf{C}) for square A\mathbf{A} and C\mathbf{C}.

?
Problem 7.1.5

Determine the eigenvectors of D=diag⁡(λ1,λ2,…,λn)\mathbf{D} = \operatorname {diag}(\lambda_1, \lambda_2, \ldots , \lambda_n). In particular, what is the eigenspace associated with λi\lambda_i?

?
Problem 7.1.6

Prove that 0∈σ(A)0 \in \sigma (\mathbf{A}) if and only if A\mathbf{A} is a singular matrix.

?
Problem 7.1.7

Explain why it's apparent that

An×n=(n11⋯11n1⋯111n⋯1⋮⋮⋮⋱⋮111⋯n) \mathbf{A}_{n\times n} = \begin{pmatrix} n & 1 & 1 & \cdots & 1 \\ 1 & n & 1 & \cdots & 1 \\ 1 & 1 & n & \cdots & 1 \\ \vdots & \vdots & \vdots & \ddots & \vdots \\ 1 & 1 & 1 & \cdots & n \end{pmatrix}

doesn't have a zero eigenvalue, and hence why A\mathbf{A} is nonsingular.

?
Problem 7.1.8

Explain why the eigenvalues of A∗A\mathbf{A}^{*}\mathbf{A} and AA∗\mathbf{A}\mathbf{A}^{*} are real and nonnegative for every A∈Cm×n\mathbf{A} \in \mathcal{C}^{m\times n}.

Hint: Consider ∥Ax→∥22/∥x→∥22\left\| \mathbf{A}\overrightarrow {x}\right\|_2^2 / \left\| \overrightarrow {x}\right\|_2^2. When are the eigenvalues of A∗A\mathbf{A}^{*}\mathbf{A} and AA∗\mathbf{A}\mathbf{A}^{*} strictly positive?

?
Problem 7.1.9
?
(a)

If A\mathbf{A} is nonsingular, and if (λ,x→)(\lambda , \overrightarrow {x}) is an eigenpair for A\mathbf{A}, show that (λ−1,x→)(\lambda^{-1}, \overrightarrow {x}) is an eigenpair for A−1\mathbf{A}^{-1}.

(b)

For all α∉σ(A)\alpha \notin \sigma (\mathbf{A}), prove that x→\overrightarrow {x} is an eigenvector of A\mathbf{A} if and only if x→\overrightarrow {x} is an eigenvector of (A−αI)−1(\mathbf{A}-\alpha \mathbf{I})^{-1}.

Problem 7.1.10
?
(a)

Show that if (λ,x→)(\lambda , \overrightarrow {x}) is an eigenpair for A\mathbf{A}, then (λk,x→)(\lambda^k, \overrightarrow {x}) is an eigenpair for Ak\mathbf{A}^k for each positive integer kk.

(b)

If p(x)=α0+α1x+α2x2+⋯+αkxkp(x) = \alpha_0 + \alpha_1 x + \alpha_2 x^2 + \cdots + \alpha_k x^k is any polynomial, then we define p(A)p(\mathbf{A}) to be the matrix

p(A)=α0I+α1A+α2A2+⋯+αkAk. p(\mathbf{A}) = \alpha _0 \mathbf{I} + \alpha _1 \mathbf{A} + \alpha _2 \mathbf{A}^2 + \cdots + \alpha _k \mathbf{A}^k.

Show that if (λ,x→)(\lambda , \overrightarrow {x}) is an eigenpair for A\mathbf{A}, then (p(λ),x→)(p(\lambda ), \overrightarrow {x}) is an eigenpair for p(A)p(\mathbf{A}).

Problem 7.1.11

Explain why Gerschgorin's theorem implies that

A=(10−200120−410−100500) \mathbf{A} = \begin{pmatrix} 1 & 0 & -2 & 0 \\ 0 & 12 & 0 & -4 \\ 1 & 0 & -1 & 0 \\ 0 & 5 & 0 & 0 \end{pmatrix}

must have at least two real eigenvalues. Corroborate this fact by computing the eigenvalues of A\mathbf{A}.

Gerschgorin's theorem: if the nn Gerschgorin circles Gi={z∈C∣∣z−aii∣≤∑j≠i∣aij∣}G_i = \left\{ z \in \mathcal{C} \mid \left|z - a_{ii}\right| \le \sum_{j \ne i} \left|a_{ij}\right|\right\} associated with the rows of an n×nn \times n matrix A\mathbf{A} are such that kk of them form a connected region that is disjoint from the remaining n−kn-k circles, then exactly kk of the eigenvalues of A\mathbf{A} lie in the connected region formed by the kk circles, and the remaining n−kn-k eigenvalues lie in the region formed by the union of the other n−kn-k circles.

?
Problem 7.1.12

If A\mathbf{A} is nilpotent (Ak=0\mathbf{A}^k = \mathbf{0} for some kk), explain why trace⁡(A)=0\operatorname {trace}(\mathbf{A}) = 0.

Hint: What is σ(A)\sigma (\mathbf{A})?

?
Problem 7.1.13

If x→1,x→2,…,x→k\overrightarrow {x}_1, \overrightarrow {x}_2, \ldots , \overrightarrow {x}_k are eigenvectors of A\mathbf{A} associated with the same eigenvalue λ\lambda, explain why every nonzero linear combination

v→=α1x→1+α2x→2+⋯+αnx→n \overrightarrow {v} = \alpha _1 \overrightarrow {x}_1 + \alpha _2 \overrightarrow {x}_2 + \cdots + \alpha _n \overrightarrow {x}_n

is also an eigenvector for A\mathbf{A} associated with the eigenvalue λ\lambda.

?
Problem 7.1.14

Explain why an eigenvector for a square matrix A\mathbf{A} cannot be associated with two distinct eigenvalues for A\mathbf{A}.

?
Problem 7.1.15

Suppose σ(An×n)=σ(Bn×n)\sigma (\mathbf{A}_{n\times n}) = \sigma (\mathbf{B}_{n\times n}). Does this guarantee that A\mathbf{A} and B\mathbf{B} have the same characteristic polynomial?

?
Problem 7.1.16

Construct 2×22 \times 2 examples to prove the following statements.

?
(a)

λ∈σ(A)\lambda \in \sigma (\mathbf{A}) and μ∈σ(B)\centernot  ⟹  λ+μ∈σ(A+B)\mu \in \sigma (\mathbf{B}) \centernot \implies \lambda +\mu \in \sigma (\mathbf{A}+\mathbf{B}).

(b)

λ∈σ(A)\lambda \in \sigma (\mathbf{A}) and μ∈σ(B)\centernot  ⟹  λμ∈σ(AB)\mu \in \sigma (\mathbf{B}) \centernot \implies \lambda \mu \in \sigma (\mathbf{A}\mathbf{B}).

Problem 7.1.17

Suppose that {λ1,λ2,…,λn}\left\{ \lambda_1, \lambda_2, \ldots , \lambda_n\right\} are the eigenvalues for An×n\mathbf{A}_{n\times n}, and let (λk,c→)(\lambda_k, \overrightarrow {c}) be a particular eigenpair.

?
(a)

For λ∉σ(A)\lambda \notin \sigma (\mathbf{A}), explain why (A−λI)−1c→=c→/(λk−λ)(\mathbf{A}-\lambda \mathbf{I})^{-1}\overrightarrow {c} = \overrightarrow {c}/(\lambda_k-\lambda ).

(b)

For an arbitrary vector d→n×1\overrightarrow {d}_{n\times 1}, prove that the eigenvalues of A+c→d→T\mathbf{A}+\overrightarrow {c}\overrightarrow {d}^T agree with those of A\mathbf{A} except that λk\lambda_k is replaced by λk+d→Tc→\lambda_k + \overrightarrow {d}^T\overrightarrow {c}.

(c)

How can d→\overrightarrow {d} be selected to guarantee that the eigenvalues of A+c→d→T\mathbf{A}+\overrightarrow {c}\overrightarrow {d}^T and A\mathbf{A} agree except that λk\lambda_k is replaced by a specified number μ\mu?

Problem 7.1.18

Suppose that A\mathbf{A} is a square matrix.

?
(a)

Explain why A\mathbf{A} and AT\mathbf{A}^T have the same eigenvalues.

(b)

Explain why λ∈σ(A)  ⟺  λ‾∈σ(A∗)\lambda \in \sigma (\mathbf{A}) \iff \overline{\lambda } \in \sigma (\mathbf{A}^{*}).

Hint: Recall Exercise 6.1.8.

(c)

Do these results imply that λ∈σ(A)  ⟺  λ‾∈σ(A)\lambda \in \sigma (\mathbf{A}) \iff \overline{\lambda } \in \sigma (\mathbf{A}) when A\mathbf{A} is a square matrix of real numbers?

(d)

A nonzero row vector y→∗\overrightarrow {y}^{*} is called a left-hand eigenvector for A\mathbf{A} whenever there is a scalar μ∈C\mu \in \mathcal{C} such that y→∗(A−μI)=0→\overrightarrow {y}^{*}(\mathbf{A}-\mu \mathbf{I}) = \overrightarrow {0}. Explain why μ\mu must be an eigenvalue for A\mathbf{A} in the ``right-hand'' sense of the term when A\mathbf{A} is a square matrix of real numbers.

Problem 7.1.19

Consider matrices Am×n\mathbf{A}_{m\times n} and Bn×m\mathbf{B}_{n\times m}.

?
(a)

Explain why AB\mathbf{A}\mathbf{B} and BA\mathbf{B}\mathbf{A} have the same characteristic polynomial if m=nm=n.

Hint: Recall Exercise 6.2.16.

(b)

Explain why the characteristic polynomials for AB\mathbf{A}\mathbf{B} and BA\mathbf{B}\mathbf{A} can't be the same when m≠nm \ne n, and then explain why σ(AB)\sigma (\mathbf{A}\mathbf{B}) and σ(BA)\sigma (\mathbf{B}\mathbf{A}) agree, with the possible exception of a zero eigenvalue.

Problem 7.1.20

If AB=BA\mathbf{A}\mathbf{B} = \mathbf{B}\mathbf{A}, prove that A\mathbf{A} and B\mathbf{B} have a common eigenvector.

Hint: For λ∈σ(A)\lambda \in \sigma (\mathbf{A}), let the columns of X\mathbf{X} be a basis for N(A−λI)N(\mathbf{A}-\lambda \mathbf{I}) so that (A−λI)BX=0(\mathbf{A}-\lambda \mathbf{I})\mathbf{B}\mathbf{X} = \mathbf{0}. Explain why there exists a matrix P\mathbf{P} such that BX=XP\mathbf{B}\mathbf{X} = \mathbf{X}\mathbf{P}, and then consider any eigenpair for P\mathbf{P}.

?
Problem 7.1.21

For fixed matrices Pm×m\mathbf{P}_{m\times m} and Qn×n\mathbf{Q}_{n\times n}, let T\mathbf{T} be the linear operator on Cm×n\mathcal{C}^{m\times n} defined by T(A)=PAQ\mathbf{T}(\mathbf{A}) = \mathbf{P}\mathbf{A}\mathbf{Q}.

?
(a)

Show that if x→\overrightarrow {x} is a right-hand eigenvector for P\mathbf{P} and y→∗\overrightarrow {y}^{*} is a left-hand eigenvector for Q\mathbf{Q}, then x→y→∗\overrightarrow {x}\overrightarrow {y}^{*} is an eigenvector for T\mathbf{T}.

(b)

Explain why trace⁡(T)=trace⁡(P)trace⁡(Q)\operatorname {trace}(\mathbf{T}) = \operatorname {trace}(\mathbf{P})\operatorname {trace}(\mathbf{Q}).

Problem 7.1.22

Let D=diag⁡(λ1,λ2,…,λn)\mathbf{D} = \operatorname {diag}(\lambda_1, \lambda_2, \ldots , \lambda_n) be a diagonal real matrix such that λ1<λ2<⋯<λn\lambda_1 < \lambda_2 < \cdots < \lambda_n, and let v→n×1\overrightarrow {v}_{n\times 1} be a column of real nonzero numbers.

?
(a)

Prove that if α\alpha is real and nonzero, then λi\lambda_i is not an eigenvalue for D+αv→v→T\mathbf{D}+\alpha \overrightarrow {v}\overrightarrow {v}^T. Show that the eigenvalues of D+αv→v→T\mathbf{D}+\alpha \overrightarrow {v}\overrightarrow {v}^T are in fact given by the solutions of the secular equation f(ξ)=0f(\xi )=0 defined by

f(ξ)=1+α∑i=1nvi2λi−ξ. f(\xi ) = 1 + \alpha \sum _{i=1}^n \frac{v_i^2}{\lambda _i - \xi }.

For n=4n=4 and α>0\alpha >0, verify that the graph of f(ξ)f(\xi ) has nn vertical asymptotes at the λi\lambda_i's, with ff increasing from 00 to +∞+\infty on (λi,λi+1)(\lambda_i, \lambda_{i+1}) and a further branch beyond λn\lambda_n approaching the horizontal asymptote at 11, so that each interval (λi,λi+1)(\lambda_i, \lambda_{i+1}) and the ray (λn,∞)(\lambda_n, \infty ) contains exactly one root ξi\xi_i of ff, and thereby conclude that the eigenvalues of D+αv→v→T\mathbf{D}+\alpha \overrightarrow {v}\overrightarrow {v}^T interlace with those of D\mathbf{D}.

Figure 7.1.4Figure 7.1.4

(b)

Verify that (D−ξiI)−1v→(\mathbf{D}-\xi_i\mathbf{I})^{-1}\overrightarrow {v} is an eigenvector for D+αv→v→T\mathbf{D}+\alpha \overrightarrow {v}\overrightarrow {v}^T that is associated with the eigenvalue ξi\xi_i.

Problem 7.1.23

Let λ1,…,λn\lambda_1, \ldots , \lambda_n be the roots of the polynomial p(λ)=λn+c1λn−1+c2λn−2+⋯+cnp(\lambda ) = \lambda^n + c_1\lambda^{n-1} + c_2\lambda^{n-2} + \cdots + c_n, and let τk=λ1k+λ2k+⋯+λnk\tau_k = \lambda_1^k + \lambda_2^k + \cdots + \lambda_n^k. Newton's identities say ck=−(τ1ck−1+τ2ck−2+⋯+τk−1c1+τk)/kc_k = -(\tau_1 c_{k-1} + \tau_2 c_{k-2} + \cdots + \tau_{k-1}c_1 + \tau_k)/k. Derive these identities by executing the following steps:

?
(a)

Show p′(λ)=p(λ)∑i=1n(λ−λi)−1p'(\lambda ) = p(\lambda )\sum_{i=1}^n (\lambda -\lambda_i)^{-1} (logarithmic differentiation).

(b)

Use the geometric series expansion for (λ−λi)−1(\lambda -\lambda_i)^{-1} to show that for ∣λ∣>max⁡i∣λi∣\left|\lambda \right| > \max_i \left|\lambda_i\right|,

∑i=1n1λ−λi=nλ+τ1λ2+τ2λ3+⋯ . \sum _{i=1}^n \frac{1}{\lambda -\lambda _i} = \frac{n}{\lambda } + \frac{\tau _1}{\lambda ^2} + \frac{\tau _2}{\lambda ^3} + \cdots .
(c)

Combine these two results, and equate like powers of λ\lambda.

Problem 7.1.24

Let the characteristic equation for A\mathbf{A} be given by λn+c1λn−1+c2λn−2+⋯+cn=0\lambda^n + c_1\lambda^{n-1} + c_2\lambda^{n-2} + \cdots + c_n = 0, and define a sequence by taking B0=I\mathbf{B}_0 = \mathbf{I} and

Bk=−trace⁡(ABk−1)kI+ABk−1for k=1,2,…,n. \mathbf{B}_k = -\frac{\operatorname {trace}(\mathbf{A}\mathbf{B}_{k-1})}{k}\mathbf{I} + \mathbf{A}\mathbf{B}_{k-1} \quad \text{for } k = 1, 2, \ldots , n.

Prove that for each kk,

ck=−trace⁡(ABk−1)k. c_k = -\frac{\operatorname {trace}(\mathbf{A}\mathbf{B}_{k-1})}{k}.

Hint: Use Newton's identities (Exercise 7.1.23), and recall Exercise 7.1.10(a).

?