Chapter 21

Eigenvalues and Eigenvectors

[32/100%]
Show
LaTeX
Problem 21.1

Show that an n×nn \times n skew-symmetric matrix A\mathbf{A} has no nonzero eigenvalues.

?
Problem 21.2

Let A\mathbf{A} represent an n×nn \times n matrix, B\mathbf{B} a k×kk \times k matrix, and X\mathbf{X} an n×kn \times k matrix such that AX=XB\mathbf{A}\mathbf{X} = \mathbf{X}\mathbf{B}.

?
(a)

Show that C(X)C(\mathbf{X}) is an invariant subspace (of Rn×1\mathbb {R}^{n \times 1}) relative to A\mathbf{A}.

(b)

Show that if X\mathbf{X} is of full column rank, then every eigenvalue of B\mathbf{B} is an eigenvalue of A\mathbf{A}.

Problem 21.3

Let p(λ)p(\lambda ) represent the characteristic polynomial of an n×nn \times n matrix A\mathbf{A}, and let c0,c1,c2,…,cnc_{0}, c_{1}, c_{2}, \ldots , c_{n} represent the respective coefficients of the characteristic polynomial, so that

p(λ)=c0λ0+c1λ+c2λ2+⋯+cnλn=∑s=0ncsλs p(\lambda ) = c_{0}\lambda ^{0} + c_{1}\lambda + c_{2}\lambda ^{2} + \cdots + c_{n}\lambda ^{n} = \sum _{s=0}^{n} c_{s}\lambda ^{s}

(for λ∈R\lambda \in \mathbb {R}). Further, let P\mathbf{P} represent the n×nn \times n matrix obtained from p(λ)p(\lambda ) by formally replacing the scalar λ\lambda with the n×nn \times n matrix A\mathbf{A} (and by setting A0=In\mathbf{A}^{0} = \mathbf{I}_{n}). That is, let

P=c0I+c1A+c2A2+⋯+cnAn=∑s=0ncsAs. \mathbf{P} = c_{0}\mathbf{I} + c_{1}\mathbf{A} + c_{2}\mathbf{A}^{2} + \cdots + c_{n}\mathbf{A}^{n} = \sum _{s=0}^{n} c_{s}\mathbf{A}^{s}.

Show that P=0\mathbf{P} = \mathbf{0} (a result that is known as the Cayley-Hamilton theorem) by carrying out the following four steps.

?
(a)

Letting B(λ)=A−λI\mathbf{B}(\lambda ) = \mathbf{A} - \lambda \mathbf{I} and letting H(λ)\mathbf{H}(\lambda ) represent the adjoint matrix of B(λ)\mathbf{B}(\lambda ), show that (for λ∈R\lambda \in \mathbb {R})

H(λ)=K0+λK1+λ2K2+⋯+λn−1Kn−1, \mathbf{H}(\lambda ) = \mathbf{K}_{0} + \lambda \mathbf{K}_{1} + \lambda ^{2}\mathbf{K}_{2} + \cdots + \lambda ^{n-1}\mathbf{K}_{n-1},

where K0,K1,K2,…,Kn−1\mathbf{K}_{0}, \mathbf{K}_{1}, \mathbf{K}_{2}, \ldots , \mathbf{K}_{n-1} are n×nn \times n matrices (that do not vary with λ\lambda).

(b)

Letting T0=AK0\mathbf{T}_{0} = \mathbf{A}\mathbf{K}_{0}, Tn=−Kn−1\mathbf{T}_{n} = -\mathbf{K}_{n-1}, and (for s=1,…,n−1s = 1, \ldots , n-1) Ts=AKs−Ks−1\mathbf{T}_{s} = \mathbf{A}\mathbf{K}_{s} - \mathbf{K}_{s-1}, show that (for λ∈R\lambda \in \mathbb {R})

T0+λT1+λ2T2+⋯+λnTn=p(λ)In. \mathbf{T}_{0} + \lambda \mathbf{T}_{1} + \lambda ^{2}\mathbf{T}_{2} + \cdots + \lambda ^{n}\mathbf{T}_{n} = p(\lambda )\mathbf{I}_{n}.
(c)

Show that, for s=0,1,…,ns = 0, 1, \ldots , n, Ts=csI\mathbf{T}_{s} = c_{s}\mathbf{I}.

(d)

Show that

P=T0+AT1+A2T2+⋯+AnTn=0. \mathbf{P} = \mathbf{T}_{0} + \mathbf{A}\mathbf{T}_{1} + \mathbf{A}^{2}\mathbf{T}_{2} + \cdots + \mathbf{A}^{n}\mathbf{T}_{n} = \mathbf{0}.
Remark.
?

It follows from a fundamental result on adjoint matrices that (for λ∈R\lambda \in \mathbb {R}) B(λ)H(λ)=∣B(λ)∣In=p(λ)In\mathbf{B}(\lambda )\mathbf{H}(\lambda ) = \left|\mathbf{B}(\lambda )\right|\mathbf{I}_{n} = p(\lambda )\mathbf{I}_{n}.

Problem 21.4

Let c0,c1,…,cn−1,cnc_{0}, c_{1}, \ldots , c_{n-1}, c_{n} represent the respective coefficients of the characteristic polynomial p(λ)p(\lambda ) of an n×nn \times n matrix A\mathbf{A} [so that p(λ)=c0+c1λ+⋯+cn−1λn−1+cnλnp(\lambda ) = c_{0} + c_{1}\lambda + \cdots + c_{n-1}\lambda^{n-1} + c_{n}\lambda^{n} (for λ∈R\lambda \in \mathbb {R})]. Using the result of Exercise 3 (the Cayley-Hamilton theorem), show that if A\mathbf{A} is nonsingular, then c0≠0c_{0} \neq 0, and

A−1=−(1/c0)(c1I+c2A+⋯+cnAn−1). \mathbf{A}^{-1} = -(1/c_{0})(c_{1}\mathbf{I} + c_{2}\mathbf{A} + \cdots + c_{n}\mathbf{A}^{n-1}).
?
Problem 21.5

Show that if an n×nn \times n matrix B\mathbf{B} is similar to an n×nn \times n matrix A\mathbf{A}, then (1) Bk\mathbf{B}^{k} is similar to Ak\mathbf{A}^{k} (k=2,3,…k = 2, 3, \ldots) and (2) B′\mathbf{B}' is similar to A′\mathbf{A}'.

?
Problem 21.6

Show that if an n×nn \times n matrix B\mathbf{B} is similar to an (n×nn \times n) idempotent matrix, then B\mathbf{B} is idempotent.

?
Problem 21.7

Let A=(1001)\mathbf{A} = \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix} and B=(1101)\mathbf{B} = \begin{pmatrix} 1 & 1 \\ 0 & 1 \end{pmatrix}. Show that B\mathbf{B} has the same rank, determinant, trace, and characteristic polynomial as A\mathbf{A}, but that, nevertheless, B\mathbf{B} is not similar to A\mathbf{A}.

?
Problem 21.8

Expand on the result of Exercise 7 by showing (for an arbitrary positive integer nn) that for an n×nn \times n matrix B\mathbf{B} to be similar to an n×nn \times n matrix A\mathbf{A} it is not sufficient for B\mathbf{B} to have the same rank, determinant, trace, and characteristic polynomial as A\mathbf{A}.

?
Problem 21.9

Let A\mathbf{A} represent an n×nn \times n matrix, B\mathbf{B} a k×kk \times k matrix, and X\mathbf{X} an n×kn \times k matrix such that AX=XB\mathbf{A}\mathbf{X} = \mathbf{X}\mathbf{B}. Show that if X\mathbf{X} is of full column rank, then there exists an orthogonal matrix Q\mathbf{Q} such that Q′AQ=(T11T120T22)\mathbf{Q}'\mathbf{A}\mathbf{Q} = \begin{pmatrix} \mathbf{T}_{11} & \mathbf{T}_{12} \\ \mathbf{0} & \mathbf{T}_{22} \end{pmatrix}, where T11\mathbf{T}_{11} is a k×kk \times k matrix that is similar to B\mathbf{B}.

?
Problem 21.10

Show that if 00 is an eigenvalue of an n×nn \times n matrix A\mathbf{A}, then its algebraic multiplicity is greater than or equal to n−rk⁡(A)n - \operatorname {rk}\left(\mathbf{A}\right).

?
Problem 21.11

Let A\mathbf{A} represent an n×nn \times n matrix. Show that if a scalar λ\lambda is an eigenvalue of A\mathbf{A} of algebraic multiplicity γ\gamma, then rk⁡(A−λI)≥n−γ\operatorname {rk}\left(\mathbf{A} - \lambda \mathbf{I}\right) \geq n - \gamma.

?
Problem 21.12

Let γ1\gamma_{1} represent the algebraic multiplicity and ν1\nu_{1} the geometric multiplicity of 00 when 00 is regarded as an eigenvalue of an n×nn \times n (singular) matrix A\mathbf{A}. And let γ2\gamma_{2} represent the algebraic multiplicity and ν2\nu_{2} the geometric multiplicity of 00 when 00 is regarded as an eigenvalue of A2\mathbf{A}^{2}. Show that if ν1=γ1\nu_{1} = \gamma_{1}, then ν2=γ2=ν1\nu_{2} = \gamma_{2} = \nu_{1}.

?
Problem 21.13

Let x1\mathbf{x}_{1} and x2\mathbf{x}_{2} represent eigenvectors of an n×nn \times n matrix A\mathbf{A}, and let c1c_{1} and c2c_{2} represent nonzero scalars. Under what circumstances is the vector x=c1x1+c2x2\mathbf{x} = c_{1}\mathbf{x}_{1} + c_{2}\mathbf{x}_{2} an eigenvector of A\mathbf{A}?

?
Problem 21.14

Let A\mathbf{A} represent an n×nn \times n matrix, and suppose that there exists an n×nn \times n nonsingular matrix Q\mathbf{Q} such that Q−1AQ=D\mathbf{Q}^{-1}\mathbf{A}\mathbf{Q} = \mathbf{D} for some diagonal matrix D={di}\mathbf{D} = \{ d_{i}\}. Further, for i=1,…,ni = 1, \ldots , n, let ri′\mathbf{r}_{i}' represent the iith row of Q−1\mathbf{Q}^{-1}. Show

?
(a)

that A′\mathbf{A}' is diagonalized by (Q−1)′(\mathbf{Q}^{-1})',

(b)

that the diagonal elements of D\mathbf{D} are the (not necessarily distinct) eigenvalues of A′\mathbf{A}', and

(c)

that r1,…,rn\mathbf{r}_{1}, \ldots , \mathbf{r}_{n} are eigenvectors of A′\mathbf{A}' (with ri\mathbf{r}_{i} corresponding to the eigenvalue did_{i}).

Problem 21.15

Show that if an n×nn \times n nonsingular matrix A\mathbf{A} is diagonalized by an n×nn \times n nonsingular matrix Q\mathbf{Q}, then A−1\mathbf{A}^{-1} is also diagonalized by Q\mathbf{Q}.

?
Problem 21.16

Let A\mathbf{A} represent an n×nn \times n matrix whose spectrum comprises kk eigenvalues λ1,…,λk\lambda_{1}, \ldots , \lambda_{k} with algebraic multiplicities γ1,…,γk\gamma_{1}, \ldots , \gamma_{k}, respectively, that sum to nn. Show that A\mathbf{A} is diagonalizable if and only if, for i=1,…,ki = 1, \ldots , k, rk⁡(A−λiI)=n−γi\operatorname {rk}\left(\mathbf{A} - \lambda_{i}\mathbf{I}\right) = n - \gamma_{i}.

?
Problem 21.17

Let A\mathbf{A} represent an n×nn \times n symmetric matrix with not-necessarily-distinct eigenvalues d1,…,dnd_{1}, \ldots , d_{n} that have been ordered so that d1≤d2≤⋯≤dnd_{1} \leq d_{2} \leq \cdots \leq d_{n}. And, let Q\mathbf{Q} represent an n×nn \times n orthogonal matrix such that Q′AQ=\diag(d1,…,dn)\mathbf{Q}'\mathbf{A}\mathbf{Q} = \diag (d_{1}, \ldots , d_{n}) --- the existence of which is guaranteed. Further, for m=2,…,n−1m = 2, \ldots , n-1, define Sm={x∈Rn×1:x≠0,Qm′x=0}S_{m} = \{ \mathbf{x} \in \mathbb {R}^{n \times 1} : \mathbf{x} \neq \mathbf{0}, \mathbf{Q}_{m}'\mathbf{x} = \mathbf{0}\} and Tm={x∈Rn×1:x≠0,Pm′x=0}T_{m} = \{ \mathbf{x} \in \mathbb {R}^{n \times 1} : \mathbf{x} \neq \mathbf{0}, \mathbf{P}_{m}'\mathbf{x} = \mathbf{0}\}, where Qm=(q1,…,qm−1)\mathbf{Q}_{m} = (\mathbf{q}_{1}, \ldots , \mathbf{q}_{m-1}) and Pm=(qm+1,…,qn)\mathbf{P}_{m} = (\mathbf{q}_{m+1}, \ldots , \mathbf{q}_{n}). Show that, for m=2,…,n−1m = 2, \ldots , n-1,

dm=min⁡x∈Smx′Axx′x=max⁡x∈Tmx′Axx′x. d_{m} = \min _{\mathbf{x} \in S_{m}} \frac{\mathbf{x}'\mathbf{A}\mathbf{x}}{\mathbf{x}'\mathbf{x}} = \max _{\mathbf{x} \in T_{m}} \frac{\mathbf{x}'\mathbf{A}\mathbf{x}}{\mathbf{x}'\mathbf{x}}.
?
Problem 21.18

Let A\mathbf{A} represent an n×nn \times n symmetric matrix, and adopt the following notation: d1,…,dnd_{1}, \ldots , d_{n} are the not-necessarily-distinct eigenvalues of A\mathbf{A}, q1,…,qn\mathbf{q}_{1}, \ldots , \mathbf{q}_{n} are orthonormal eigenvectors corresponding to d1,…,dnd_{1}, \ldots , d_{n}, respectively, Q=(q1,…,qn)\mathbf{Q} = (\mathbf{q}_{1}, \ldots , \mathbf{q}_{n}), {λ1,…,λk}\{ \lambda_{1}, \ldots , \lambda_{k}\} is the spectrum of A\mathbf{A}; and, for j=1,…,kj = 1, \ldots , k, Sj={i:di=λj}S_{j} = \{ i : d_{i} = \lambda_{j}\}, Ej=∑i∈Sjqiqi′\mathbf{E}_{j} = \sum_{i \in S_{j}} \mathbf{q}_{i}\mathbf{q}_{i}', and Qj=(qi1,…,qiνj)\mathbf{Q}_{j} = (\mathbf{q}_{i_{1}}, \ldots , \mathbf{q}_{i_{\nu_{j}}}), where i1,…,iνji_{1}, \ldots , i_{\nu_{j}} denote the elements of SjS_{j}.

?
(a)

Show that the matrices E1,…,Ek\mathbf{E}_{1}, \ldots , \mathbf{E}_{k}, which appear in the spectral decomposition A=∑j=1kλjEj\mathbf{A} = \sum_{j=1}^{k} \lambda_{j}\mathbf{E}_{j}, have the following properties:

  1. E1+⋯+Ek=I\mathbf{E}_{1} + \cdots + \mathbf{E}_{k} = \mathbf{I};

  2. E1,…,Ek\mathbf{E}_{1}, \ldots , \mathbf{E}_{k} are nonnull, symmetric, and idempotent;

  3. for t≠j=1,…,kt \neq j = 1, \ldots , k, EtEj=0\mathbf{E}_{t}\mathbf{E}_{j} = \mathbf{0}; and

  4. rk⁡(E1)+⋯+rk⁡(Ek)=n\operatorname {rk}\left(\mathbf{E}_{1}\right) + \cdots + \operatorname {rk}\left(\mathbf{E}_{k}\right) = n.

(b)

Take F1,…,Fr\mathbf{F}_{1}, \ldots , \mathbf{F}_{r} to be n×nn \times n nonnull idempotent matrices such that F1+⋯+Fr=I\mathbf{F}_{1} + \cdots + \mathbf{F}_{r} = \mathbf{I}. And, suppose that, for some distinct scalars τ1,…,τr\tau_{1}, \ldots , \tau_{r},

A=τ1F1+⋯+τrFr. \mathbf{A} = \tau _{1}\mathbf{F}_{1} + \cdots + \tau _{r}\mathbf{F}_{r}.

Show that r=kr = k and that there exists a permutation t1,…,trt_{1}, \ldots , t_{r} of the first rr positive integers such that (for j=1,…,rj = 1, \ldots , r) τj=λtj\tau_{j} = \lambda_{t_{j}} and Fj=Etj\mathbf{F}_{j} = \mathbf{E}_{t_{j}}.

Problem 21.19

Let A\mathbf{A} represent an n×nn \times n symmetric matrix, and let d1,…,dnd_{1}, \ldots , d_{n} represent the (not-necessarily-distinct) eigenvalues of A\mathbf{A}. Show that lim⁡k→∞Ak=0\lim_{k \rightarrow \infty } \mathbf{A}^{k} = \mathbf{0} if and only if, for i=1,…,ni = 1, \ldots , n, ∣di∣<1\left|d_{i}\right| < 1.

?
Problem 21.20
?
(a)

Show that if 00 is an eigenvalue of an n×nn \times n not-necessarily-symmetric matrix A\mathbf{A}, then it is also an eigenvalue of A+\mathbf{A}^{+} and that the geometric multiplicity of 00 is the same when it is regarded as an eigenvalue of A+\mathbf{A}^{+} as when it is regarded as an eigenvalue of A\mathbf{A}.

(b)

Show (via an example) that the reciprocals of the nonzero eigenvalues of a square nonsymmetric matrix A\mathbf{A} are not necessarily eigenvalues of A+\mathbf{A}^{+}.

Problem 21.21

Show that, for any positive integer nn that is divisible by 2, there exists an n×nn \times n orthogonal matrix that has no eigenvalues.

?
Remark.
?

Find a 2×22 \times 2 orthogonal matrix Q\mathbf{Q} that has no eigenvalues, and then consider the block-diagonal matrix \diag(Q,Q,…,Q)\diag (\mathbf{Q}, \mathbf{Q}, \ldots , \mathbf{Q}).

Problem 21.22

Let Q\mathbf{Q} represent an n×nn \times n orthogonal matrix, and let p(λ)p(\lambda ) represent the characteristic polynomial of Q\mathbf{Q}. Show that (for λ≠0\lambda \neq 0)

p(λ)=±λnp(1/λ). p(\lambda ) = \pm \lambda ^{n} p(1/\lambda ).
?
Problem 21.23

Let A\mathbf{A} represent an n×nn \times n matrix, and suppose that the scalar 1 is an eigenvalue of A\mathbf{A} of geometric multiplicity ν\nu. Show that ν≤rk⁡(A)\nu \leq \operatorname {rk}\left(\mathbf{A}\right) and that if ν=rk⁡(A)\nu = \operatorname {rk}\left(\mathbf{A}\right), then A\mathbf{A} is idempotent.

?
Problem 21.24

Let A\mathbf{A} represent an n×nn \times n nonsingular matrix. And, let λ\lambda represent an eigenvalue of A\mathbf{A}, and x\mathbf{x} represent an eigenvector of A\mathbf{A} corresponding to λ\lambda. Show that ∣A∣/λ\left|\mathbf{A}\right|/\lambda is an eigenvalue of adj⁡(A)\operatorname {adj}(\mathbf{A}) and that x\mathbf{x} is an eigenvector of adj⁡(A)\operatorname {adj}(\mathbf{A}) corresponding to ∣A∣/λ\left|\mathbf{A}\right|/\lambda.

?
Problem 21.25

Let A\mathbf{A} represent an n×nn \times n matrix, and let p(λ)p(\lambda ) represent the characteristic polynomial of A\mathbf{A}. And, let λ1,…,λk\lambda_{1}, \ldots , \lambda_{k} represent the distinct eigenvalues of A\mathbf{A}, and γ1,…,γk\gamma_{1}, \ldots , \gamma_{k} represent their respective algebraic multiplicities, so that (for all λ\lambda)

p(λ)=(−1)nq(λ)∏j=1k(λ−λj)γj p(\lambda ) = (-1)^{n} q(\lambda ) \prod _{j=1}^{k} (\lambda - \lambda _{j})^{\gamma _{j}}

for some polynomial q(λ)q(\lambda ) (of degree n−∑j=1kγjn - \sum_{j=1}^{k} \gamma_{j}) that has no real roots. Further, define B=A−λ1UV′\mathbf{B} = \mathbf{A} - \lambda_{1}\mathbf{U}\mathbf{V}', where U=(u1,…,uγ1)\mathbf{U} = (\mathbf{u}_{1}, \ldots , \mathbf{u}_{\gamma_{1}}) is an n×γ1n \times \gamma_{1} matrix whose columns u1,…,uγ1\mathbf{u}_{1}, \ldots , \mathbf{u}_{\gamma_{1}} are (not necessarily linearly independent) eigenvectors of A\mathbf{A} corresponding to λ1\lambda_{1} and where V=(v1,…,vγ1)\mathbf{V} = (\mathbf{v}_{1}, \ldots , \mathbf{v}_{\gamma_{1}}) is an n×γ1n \times \gamma_{1} matrix such that V′U\mathbf{V}'\mathbf{U} is diagonal.

?
(a)

Show that the characteristic polynomial, say r(λ)r(\lambda ), of B\mathbf{B} is such that (for all λ\lambda)

r(λ)=(−1)nq(λ)∏i=1γ1[λ−(1−vi′ui)λ1]∏j=2k(λ−λj)γj.(E.1) r(\lambda ) = (-1)^{n} q(\lambda ) \prod _{i=1}^{\gamma _{1}} [\lambda - (1 - \mathbf{v}_{i}'\mathbf{u}_{i})\lambda _{1}] \prod _{j=2}^{k} (\lambda - \lambda _{j})^{\gamma _{j}}. \tag {E.1}
(b)

Show that in the special case where U′V=cI\mathbf{U}'\mathbf{V} = c\mathbf{I} for some nonzero scalar cc, the distinct eigenvalues of B\mathbf{B} are either λ2,…,λs−1,λs,λs+1,…,λk\lambda_{2}, \ldots , \lambda_{s-1}, \lambda_{s}, \lambda_{s+1}, \ldots , \lambda_{k} with algebraic multiplicities γ2,…,γs−1,γs+γ1,γs+1,…,γk\gamma_{2}, \ldots , \gamma_{s-1}, \gamma_{s} + \gamma_{1}, \gamma_{s+1}, \ldots , \gamma_{k}, respectively, or (1−c)λ1,λ2,…,λk(1-c)\lambda_{1}, \lambda_{2}, \ldots , \lambda_{k} with algebraic multiplicities γ1,γ2,…,γk\gamma_{1}, \gamma_{2}, \ldots , \gamma_{k}, respectively [depending on whether or not (1−c)λ1=λs(1-c)\lambda_{1} = \lambda_{s} for some ss (2≤s≤k2 \leq s \leq k)].

(c)

Show that, in the special case where γ1=1\gamma_{1} = 1, (1) u1\mathbf{u}_{1} is an eigenvector of B\mathbf{B} corresponding to the eigenvalue (1−v1′u1)λ1(1 - \mathbf{v}_{1}'\mathbf{u}_{1})\lambda_{1} and (2) for any eigenvector x\mathbf{x} of B\mathbf{B} corresponding to an eigenvalue λ\lambda [other than (1−v1′u1)λ1(1 - \mathbf{v}_{1}'\mathbf{u}_{1})\lambda_{1}], the vector

x−λ1(λ1−λ)−1(v1′x)u1 \mathbf{x} - \lambda _{1}(\lambda _{1} - \lambda )^{-1}(\mathbf{v}_{1}'\mathbf{x})\mathbf{u}_{1}

is an eigenvector of A\mathbf{A} corresponding to λ\lambda.

Remark.
?

Since the left and right sides of equality (E.1) are polynomials (in λ\lambda), it suffices to show that they are equal for all λ\lambda other than λ1,…,λk\lambda_{1}, \ldots , \lambda_{k}.

Problem 21.26

Let A\mathbf{A} represent an m×nm \times n matrix of rank rr. And, take P\mathbf{P} to be any m×mm \times m orthogonal matrix and D1\mathbf{D}_{1} to be any r×rr \times r nonsingular diagonal matrix such that

P′AA′P=(D12000). \mathbf{P}'\mathbf{A}\mathbf{A}'\mathbf{P} = \begin{pmatrix} \mathbf{D}_{1}^{2} & \mathbf{0} \\ \mathbf{0} & \mathbf{0} \end{pmatrix}.

Further, partition P\mathbf{P} as P=(P1,P2)\mathbf{P} = (\mathbf{P}_{1}, \mathbf{P}_{2}), where P1\mathbf{P}_{1} has rr columns, and let Q=(Q1,Q2)\mathbf{Q} = (\mathbf{Q}_{1}, \mathbf{Q}_{2}), where Q1=A′P1D1−1\mathbf{Q}_{1} = \mathbf{A}'\mathbf{P}_{1}\mathbf{D}_{1}^{-1} and where Q2\mathbf{Q}_{2} is any n×(n−r)n \times (n-r) matrix such that Q1′Q2=0\mathbf{Q}_{1}'\mathbf{Q}_{2} = \mathbf{0}. Show that

P′AQ=(D1000). \mathbf{P}'\mathbf{A}\mathbf{Q} = \begin{pmatrix} \mathbf{D}_{1} & \mathbf{0} \\ \mathbf{0} & \mathbf{0} \end{pmatrix}.
?
Problem 21.27

Let A\mathbf{A} represent an m×nm \times n matrix. And, let P\mathbf{P} represent an m×mm \times m orthogonal matrix, Q\mathbf{Q} an n×nn \times n orthogonal matrix, and D1={si}\mathbf{D}_{1} = \{ s_{i}\} an r×rr \times r diagonal matrix with (strictly) positive diagonal elements such that P′AQ=(D1000)\mathbf{P}'\mathbf{A}\mathbf{Q} = \begin{pmatrix} \mathbf{D}_{1} & \mathbf{0} \\ \mathbf{0} & \mathbf{0} \end{pmatrix}; denote by p1,…,pr\mathbf{p}_{1}, \ldots , \mathbf{p}_{r} the first through rrth columns of P\mathbf{P} and by q1,…,qr\mathbf{q}_{1}, \ldots , \mathbf{q}_{r} the first through rrth columns of Q\mathbf{Q}; and define α1,…,αk\alpha_{1}, \ldots , \alpha_{k} to be the distinct values represented among s1,…,srs_{1}, \ldots , s_{r}. Then, the singular value decomposition of A\mathbf{A} is A=∑j=1kαjUj\mathbf{A} = \sum_{j=1}^{k} \alpha_{j}\mathbf{U}_{j}, where (for j=1,…,kj = 1, \ldots , k) Uj=∑i∈Ljpiqi′\mathbf{U}_{j} = \sum_{i \in L_{j}} \mathbf{p}_{i}\mathbf{q}_{i}' with Lj={i:si=αj}L_{j} = \{ i : s_{i} = \alpha_{j}\}. Show that the matrices U1,…,Uk\mathbf{U}_{1}, \ldots , \mathbf{U}_{k}, which appear in the singular value decomposition, are such that UjUj′Uj=Uj\mathbf{U}_{j}\mathbf{U}_{j}'\mathbf{U}_{j} = \mathbf{U}_{j} (for j=1,…,kj = 1, \ldots , k) and Ut′Uj=0\mathbf{U}_{t}'\mathbf{U}_{j} = \mathbf{0} and UtUj′=0\mathbf{U}_{t}\mathbf{U}_{j}' = \mathbf{0} (for t≠j=1,…,kt \neq j = 1, \ldots , k).

?
Problem 21.28

Let A\mathbf{A} represent an m×nm \times n matrix. And, as in Exercise 27, take P\mathbf{P} to be an m×mm \times m orthogonal matrix, Q\mathbf{Q} an n×nn \times n orthogonal matrix, and D1\mathbf{D}_{1} an r×rr \times r nonsingular diagonal matrix such that P′AQ=(D1000)\mathbf{P}'\mathbf{A}\mathbf{Q} = \begin{pmatrix} \mathbf{D}_{1} & \mathbf{0} \\ \mathbf{0} & \mathbf{0} \end{pmatrix}. Further, partition P\mathbf{P} and Q\mathbf{Q} as P=(P1,P2)\mathbf{P} = (\mathbf{P}_{1}, \mathbf{P}_{2}) and Q=(Q1,Q2)\mathbf{Q} = (\mathbf{Q}_{1}, \mathbf{Q}_{2}), where each of the matrices P1\mathbf{P}_{1} and Q1\mathbf{Q}_{1} has rr columns. Show that C(A)=C(P1)C(\mathbf{A}) = C(\mathbf{P}_{1}) and that N(A)=C(Q2)N(\mathbf{A}) = C(\mathbf{Q}_{2}).

?
Problem 21.29

Let A1,…,Ak\mathbf{A}_{1}, \ldots , \mathbf{A}_{k} represent n×nn \times n not-necessarily-symmetric matrices, each of which is diagonalizable. Show that if A1,…,Ak\mathbf{A}_{1}, \ldots , \mathbf{A}_{k} commute in pairs, then A1,…,Ak\mathbf{A}_{1}, \ldots , \mathbf{A}_{k} are simultaneously diagonalizable.

?
Problem 21.30

Let V\mathbf{V} represent an n×nn \times n symmetric nonnegative definite matrix, X\mathbf{X} an n×pn \times p matrix of rank rr, and d\mathbf{d} a p×1p \times 1 vector. Using the results of Exercise 19.11 (or otherwise), show that each of the following three conditions is necessary and sufficient for the vector X(X′X)−d\mathbf{X}(\mathbf{X}'\mathbf{X})^{-}\mathbf{d} to be a solution, for every d∈C(X′)\mathbf{d} \in C(\mathbf{X}'), to the problem of minimizing the quadratic form a′Va\mathbf{a}'\mathbf{V}\mathbf{a} (in a\mathbf{a}) subject to X′a=d\mathbf{X}'\mathbf{a} = \mathbf{d}:

?
(a)

there exists an orthogonal matrix that simultaneously diagonalizes V\mathbf{V} and PX\mathbf{P}_{\mathbf{X}};

(b)

there exists a subset of rr orthonormal eigenvectors of V\mathbf{V} that is a basis for C(X)C(\mathbf{X});

(c)

there exists a subset of rr eigenvectors of V\mathbf{V} that is a basis for C(X)C(\mathbf{X}).

Problem 21.31

Let A\mathbf{A} represent an n×nn \times n symmetric matrix, and let B\mathbf{B} represent an n×nn \times n symmetric positive definite matrix. And, let λmax⁡\lambda_{\max } and λmin⁡\lambda_{\min } represent, respectively, the largest and smallest roots of ∣A−λB∣\left|\mathbf{A} - \lambda \mathbf{B}\right|. Show that

λmin⁡≤x′Axx′Bx≤λmax⁡ \lambda _{\min } \leq \frac{\mathbf{x}'\mathbf{A}\mathbf{x}}{\mathbf{x}'\mathbf{B}\mathbf{x}} \leq \lambda _{\max }

for every nonnull vector x\mathbf{x} in Rn\mathbb {R}^{n}.

?
Problem 21.32

Let A\mathbf{A} represent an n×nn \times n symmetric matrix, and let B\mathbf{B} represent an n×nn \times n symmetric positive definite matrix. Show that A−B\mathbf{A} - \mathbf{B} is nonnegative definite if and only if all nn (not necessarily distinct) roots of ∣A−λB∣\left|\mathbf{A} - \lambda \mathbf{B}\right| are greater than or equal to 1 and is positive definite if and only if all nn roots are (strictly) greater than 1.

?