SectionPEEProperties of Eigenvalues and Eigenvectors

The previous section introduced eigenvalues and eigenvectors, and concentrated on their existence and determination. This section will be more about theorems, and the various properties eigenvalues and eigenvectors enjoy. Like a good 4×100 meter relay, we will lead-off with one of our better theorems and save the very best for the anchor leg.

SubsectionBPEBasic Properties of Eigenvalues

If $$p=1\text{,}$$ then the set $$S=\set{\vect{x}_1}$$ is linearly independent since eigenvectors are nonzero (Definition EEM), so assume for the remainder that $$p\geq 2\text{.}$$

We will prove this result by contradiction (Proof Technique CD). Suppose to the contrary that $$S$$ is a linearly dependent set. Define $$S_i=\set{\vectorlist{x}{i}}$$ and let $$k$$ be an integer such that $$S_{k-1}=\set{\vectorlist{x}{k-1}}$$ is linearly independent and $$S_k=\set{\vectorlist{x}{k}}$$ is linearly dependent. We have to ask if there is even such an integer $$k\text{?}$$ First, since eigenvectors are nonzero, the set $$\set{\vect{x}_1}$$ is linearly independent. Since we are assuming that $$S=S_p$$ is linearly dependent, there must be an integer $$k\text{,}$$ $$2\leq k\leq p\text{,}$$ where the sets $$S_i$$ transition from linear independence to linear dependence (and stay that way). In other words, $$\vect{x}_k$$ is the vector with the smallest index that is a linear combination of just vectors with smaller indices.

Since $$\set{\vectorlist{x}{k}}$$ is a linearly dependent set there must be scalars, $$\scalarlist{a}{k}\text{,}$$ not all zero (Definition LI), so that

\begin{gather*} \zerovector=\lincombo{a}{x}{k}\text{.} \end{gather*}

Then,

\begin{align*} \zerovector &=\left(A-\lambda_kI_n\right)\zerovector&& \knowl{./knowl/theorem-ZVSM.html}{\text{Theorem ZVSM}}\\ &=\left(A-\lambda_kI_n\right)\left(\lincombo{a}{x}{k}\right)&& \knowl{./knowl/definition-RLD.html}{\text{Definition RLD}}\\ &=\left(A-\lambda_kI_n\right)a_1\vect{x}_1+ \cdots+ \left(A-\lambda_kI_n\right)a_k\vect{x}_k&& \knowl{./knowl/theorem-MMDAA.html}{\text{Theorem MMDAA}}\\ &=a_1\left(A-\lambda_kI_n\right)\vect{x}_1+ \cdots+ a_k\left(A-\lambda_kI_n\right)\vect{x}_k&& \knowl{./knowl/theorem-MMSMM.html}{\text{Theorem MMSMM}}\\ &=a_1\left(A\vect{x}_1-\lambda_kI_n\vect{x}_1\right)+ \cdots+ a_k\left(A\vect{x}_k-\lambda_kI_n\vect{x}_k\right)&& \knowl{./knowl/theorem-MMDAA.html}{\text{Theorem MMDAA}}\\ &=a_1\left(A\vect{x}_1-\lambda_k\vect{x}_1\right)+ \cdots+ a_k\left(A\vect{x}_k-\lambda_k\vect{x}_k\right)&& \knowl{./knowl/theorem-MMIM.html}{\text{Theorem MMIM}}\\ &=a_1\left(\lambda_1\vect{x}_1-\lambda_k\vect{x}_1\right)+ \cdots+ a_k\left(\lambda_k\vect{x}_k-\lambda_k\vect{x}_k\right)&& \knowl{./knowl/definition-EEM.html}{\text{Definition EEM}}\\ &=a_1\left(\lambda_1-\lambda_k\right)\vect{x}_1+ \cdots+ a_k\left(\lambda_k-\lambda_k\right)\vect{x}_k&& \knowl{./knowl/theorem-MMDAA.html}{\text{Theorem MMDAA}}\\ &=a_1\left(\lambda_1-\lambda_k\right)\vect{x}_1+ \cdots+ a_{k-1}\left(\lambda_{k-1}-\lambda_k\right)\vect{x}_{k-1}+ a_k\left(0\right)\vect{x}_k&& \knowl{./knowl/property-AICN.html}{\text{Property AICN}}\\ &=a_1\left(\lambda_1-\lambda_k\right)\vect{x}_1+ \cdots+ a_{k-1}\left(\lambda_{k-1}-\lambda_k\right)\vect{x}_{k-1}+ \zerovector&& \knowl{./knowl/theorem-ZSSM.html}{\text{Theorem ZSSM}}\\ &=a_1\left(\lambda_1-\lambda_k\right)\vect{x}_1+ \cdots+ a_{k-1}\left(\lambda_{k-1}-\lambda_k\right)\vect{x}_{k-1}&& \knowl{./knowl/property-Z.html}{\text{Property Z}}\text{.} \end{align*}

This equation is a relation of linear dependence on the linearly independent set $$\set{\vectorlist{x}{k-1}}\text{,}$$ so the scalars must all be zero. That is, $$a_i\left(\lambda_i-\lambda_k\right)=0$$ for $$1\leq i\leq k-1\text{.}$$ However, we have the hypothesis that the eigenvalues are distinct, so $$\lambda_i\neq\lambda_k$$ for $$1\leq i\leq k-1$$ and so $$\lambda_i-\lambda_k\neq 0\text{.}$$ So Theorem ZPZT implies $$a_i=0$$ for $$1\leq i\leq k-1\text{.}$$

This reduces the original relation of linear dependence on $$\set{\vectorlist{x}{k}}$$ to the simpler equation $$a_k\vect{x}_k=\zerovector\text{.}$$ By Theorem SMEZV we conclude that $$a_k=0$$ or $$\vect{x}_k=\zerovector\text{.}$$ Eigenvectors are never the zero vector (Definition EEM), so $$a_k=0\text{.}$$ So now all of the scalars $$a_i\text{,}$$ $$1\leq i\leq k$$ are zero, contradicting their introduction as the scalars creating a nontrivial relation of linear dependence on the set $$\set{\vectorlist{x}{k}}\text{.}$$ With a contradiction in hand, we conclude that $$S$$ must be linearly independent.

There is a simple connection between the eigenvalues of a matrix and whether or not the matrix is nonsingular.

We have the following equivalences:

\begin{align*} A\text{ is singular}&\iff\text{there exists }\vect{x}\neq\zerovector,\,A\vect{x}=\zerovector&& \knowl{./knowl/definition-NM.html}{\text{Definition NM}}\\ &\iff\text{there exists }\vect{x}\neq\zerovector,\,A\vect{x}=0\vect{x}&& \knowl{./knowl/theorem-ZSSM.html}{\text{Theorem ZSSM}}\\ &\iff\lambda=0\text{ is an eigenvalue of A}&& \knowl{./knowl/definition-EEM.html}{\text{Definition EEM}}\text{.} \end{align*}

With an equivalence about singular matrices we can update our list of equivalences about nonsingular matrices.

The equivalence of the first and last statements is Theorem SMZE, reformulated by negating each statement in the equivalence. So we are able to improve on Theorem NME7 with this addition.

SageNME8.Nonsingular Matrix Equivalences, Round 8.

Zero eigenvalues are another marker of singular matrices. We illustrate with two matrices, the first nonsingular, the second not.

A = matrix(QQ, [[ 1,  0,  2,  1,  1, -3,  4, -6],
[ 0,  1,  0,  5,  3, -8,  0,  6],
[ 1,  1,  4,  4,  2, -8,  4, -2],
[ 0,  1,  6,  0, -2, -1,  1,  0],
[-1,  1,  0,  1,  0, -2, -1,  5],
[ 0, -1, -6, -3,  1,  4,  1, -5],
[-1, -1, -2, -3, -2,  7, -4,  2],
[ 0,  2,  0,  0, -2, -2,  0,  6]])
A.is_singular()

0 in A.eigenvalues()

B = matrix(QQ, [[ 2, -1,  0,  4, -6, -2,  8, -2],
[-1,  1, -1, -4, -1,  7, -7, -2],
[-1,  1,  1,  0,  7, -5, -6,  7],
[ 2, -1, -1,  3, -5, -2,  5,  2],
[ 1, -2,  0,  2, -1, -3,  8,  0],
[-1,  1,  0, -2,  6, -1, -8,  5],
[ 2, -1, -1,  2, -7,  1,  7, -3],
[-1,  0,  1,  0,  6, -4, -2,  4]])
B.is_singular()

0 in B.eigenvalues()


Certain changes to a matrix change its eigenvalues in a predictable way.

Let $$\vect{x}\neq\zerovector$$ be one eigenvector of $$A$$ for $$\lambda\text{.}$$ Then

\begin{align*} \left(\alpha A\right)\vect{x}&=\alpha\left(A\vect{x}\right)&& \knowl{./knowl/theorem-MMSMM.html}{\text{Theorem MMSMM}}\\ &=\alpha\left(\lambda\vect{x}\right)&& \knowl{./knowl/definition-EEM.html}{\text{Definition EEM}}\\ &=\left(\alpha\lambda\right)\vect{x}&& \knowl{./knowl/property-SMAC.html}{\text{Property SMAC}}\text{.} \end{align*}

So $$\vect{x}\neq\zerovector$$ is an eigenvector of $$\alpha A$$ for the eigenvalue $$\alpha\lambda\text{.}$$

Unfortunately, there are not parallel theorems about the sum or product of arbitrary matrices. But we can prove a similar result for powers of a matrix.

Let $$\vect{x}\neq\zerovector$$ be one eigenvector of $$A$$ for $$\lambda\text{.}$$ Then we proceed by induction on $$s$$ (Proof Technique I). First, for $$s=0\text{,}$$ employing Theorem MMIM and Property OC to establish the base case,

\begin{gather*} A^s\vect{x}=A^0\vect{x}=I_n\vect{x}=\vect{x}=1\vect{x}=\lambda^0\vect{x}=\lambda^s\vect{x}\text{.} \end{gather*}

So $$\lambda^s$$ is an eigenvalue of $$A^s$$ in this special case. For the induction step, we assume the theorem is true for $$s\text{,}$$ and find

\begin{align*} A^{s+1}\vect{x}&=A^sA\vect{x}\\ &=A^s\left(\lambda\vect{x}\right)&& \knowl{./knowl/definition-EEM.html}{\text{Definition EEM}}\\ &=\lambda\left(A^s\vect{x}\right)&& \knowl{./knowl/theorem-MMSMM.html}{\text{Theorem MMSMM}}\\ &=\lambda\left(\lambda^s\vect{x}\right)&&\text{Induction Hypothesis}\\ &=\left(\lambda\lambda^s\right)\vect{x}&& \knowl{./knowl/property-SMAC.html}{\text{Property SMAC}}\\ &=\lambda^{s+1}\vect{x}\text{.} \end{align*}

So $$\vect{x}\neq\zerovector$$ is an eigenvector of $$A^{s+1}$$ for $$\lambda^{s+1}\text{,}$$ and induction tells us the theorem is true for all $$s\geq 0\text{.}$$

While we cannot prove that the sum of two arbitrary matrices behaves in any reasonable way with regard to eigenvalues, we can work with the sum of dissimilar powers of the same matrix. We have already seen two connections between eigenvalues and polynomials, in the proof of Theorem EMHE and the characteristic polynomial (Definition CP). Our next theorem strengthens this connection.

Let $$\vect{x}\neq\zerovector$$ be one eigenvector of $$A$$ for $$\lambda\text{,}$$ and write $$q(x)=a_0+a_1x+a_2x^2+\cdots+a_mx^m\text{.}$$ Then

\begin{align*} q(A)\vect{x}&=\left(a_0A^0+a_1A^1+a_2A^2+\cdots+a_mA^m\right)\vect{x}\\ &=(a_0A^0)\vect{x}+(a_1A^1)\vect{x}+(a_2A^2)\vect{x}+\cdots+(a_mA^m)\vect{x}&& \knowl{./knowl/theorem-MMDAA.html}{\text{Theorem MMDAA}}\\ &=a_0(A^0\vect{x})+a_1(A^1\vect{x})+a_2(A^2\vect{x})+\cdots+a_m(A^m\vect{x})&& \knowl{./knowl/theorem-MMSMM.html}{\text{Theorem MMSMM}}\\ &=a_0(\lambda^0\vect{x})+a_1(\lambda^1\vect{x})+a_2(\lambda^2\vect{x})+\cdots+a_m(\lambda^m\vect{x})&& \knowl{./knowl/theorem-EOMP.html}{\text{Theorem EOMP}}\\ &=(a_0\lambda^0)\vect{x}+(a_1\lambda^1)\vect{x}+(a_2\lambda^2)\vect{x}+\cdots+(a_m\lambda^m)\vect{x}&& \knowl{./knowl/property-SMAC.html}{\text{Property SMAC}}\\ &=\left(a_0\lambda^0+a_1\lambda^1+a_2\lambda^2+\cdots+a_m\lambda^m\right)\vect{x}&& \knowl{./knowl/property-DSAC.html}{\text{Property DSAC}}\\ &=q(\lambda)\vect{x}\text{.} \end{align*}

So $$\vect{x}\neq 0$$ is an eigenvector of $$q(A)$$ for the eigenvalue $$q(\lambda)\text{.}$$

In Example ESMS4 the $$4\times 4$$ symmetric matrix

\begin{equation*} C= \begin{bmatrix} 1 & 0 & 1 & 1\\ 0 & 1 & 1 & 1\\ 1 & 1 & 1 & 0\\ 1 & 1 & 0 & 1 \end{bmatrix} \end{equation*}

is shown to have the three eigenvalues $$\lambda=3,\,1,\,-1\text{.}$$ Suppose we wanted a $$4\times 4$$ matrix that has the three eigenvalues $$\lambda=4,\,0,\,-2\text{.}$$ We can employ Theorem EPM by finding a polynomial that converts $$3$$ to $$4\text{,}$$ $$1$$ to $$0\text{,}$$ and $$-1$$ to $$-2\text{.}$$ Such a polynomial is called an interpolating polynomial, and in this example we can use

\begin{equation*} r(x)=\frac{1}{4}\left(x^2+4x-5\right)=\frac{1}{4}x^2+x-\frac{5}{4}\text{.} \end{equation*}

We will not discuss how to concoct this polynomial, but a text on numerical analysis should provide the details. For now, simply verify that $$r(3)=4\text{,}$$ $$r(1)=0$$ and $$r(-1)=-2\text{.}$$

Now compute

\begin{align*} r(C)&=\frac{1}{4}C^2+C-\frac{5}{4}I_4\\ &= \frac{1}{4} \begin{bmatrix} 3 & 2 & 2 & 2\\ 2 & 3 & 2 & 2\\ 2 & 2 & 3 & 2\\ 2 & 2 & 2 & 3 \end{bmatrix} + \begin{bmatrix} 1 & 0 & 1 & 1\\ 0 & 1 & 1 & 1\\ 1 & 1 & 1 & 0\\ 1 & 1 & 0 & 1 \end{bmatrix} -\frac{5}{4} \begin{bmatrix} 1 & 0 & 0 & 0\\ 0 & 1 & 0 & 0\\ 0 & 0 & 1 & 0\\ 0 & 0 & 0 & 1 \end{bmatrix} = \frac{1}{2} \begin{bmatrix} 1 & 1 & 3 & 3\\ 1 & 1 & 3 & 3\\ 3 & 3 & 1 & 1\\ 3 & 3 & 1 & 1 \end{bmatrix}\text{.} \end{align*}

Theorem EPM tells us that if $$r(x)$$ transforms the eigenvalues in the desired manner, then $$r(C)$$ will have the desired eigenvalues. You can check this by computing the eigenvalues of $$r(C)$$ directly. Furthermore, notice that the multiplicities are the same, and the eigenspaces of $$C$$ and $$r(C)$$ are identical.

Inverses and transposes also behave predictably with regard to their eigenvalues.

Notice that since $$A$$ is assumed nonsingular, $$\inverse{A}$$ exists by Theorem NI, but more importantly, $$\lambda^{-1}=1/\lambda$$ does not involve division by zero since Theorem SMZE prohibits this possibility.

Let $$\vect{x}\neq\zerovector$$ be one eigenvector of $$A$$ for $$\lambda\text{.}$$ Suppose $$A$$ has size $$n\text{.}$$ Then

\begin{align*} \inverse{A}\vect{x}&=\inverse{A}(1\vect{x})&& \knowl{./knowl/property-OC.html}{\text{Property OC}}\\ &=\inverse{A}(\frac{1}{\lambda}\lambda\vect{x})&& \knowl{./knowl/property-MICN.html}{\text{Property MICN}}\\ &=\frac{1}{\lambda}\inverse{A}(\lambda\vect{x})&& \knowl{./knowl/theorem-MMSMM.html}{\text{Theorem MMSMM}}\\ &=\frac{1}{\lambda}\inverse{A}(A\vect{x})&& \knowl{./knowl/definition-EEM.html}{\text{Definition EEM}}\\ &=\frac{1}{\lambda}(\inverse{A}A)\vect{x}&& \knowl{./knowl/theorem-MMA.html}{\text{Theorem MMA}}\\ &=\frac{1}{\lambda}I_n\vect{x}&& \knowl{./knowl/definition-MI.html}{\text{Definition MI}}\\ &=\frac{1}{\lambda}\vect{x}&& \knowl{./knowl/theorem-MMIM.html}{\text{Theorem MMIM}}\text{.} \end{align*}

So $$\vect{x}\neq\zerovector$$ is an eigenvector of $$\inverse{A}$$ for the eigenvalue $$\frac{1}{\lambda}\text{.}$$

The proofs of the theorems above have a similar style to them. They all begin by grabbing an eigenvalue-eigenvector pair and adjusting it in some way to reach the desired conclusion. You should add this to your toolkit as a general approach to proving theorems about eigenvalues.

So far we have been able to reserve the characteristic polynomial for strictly computational purposes. However, sometimes a theorem about eigenvalues can be proved easily by employing the characteristic polynomial (rather than using an eigenvalue-eigenvector pair). The next theorem is an example of this.

Suppose $$A$$ has size $$n\text{.}$$ Then

\begin{align*} \charpoly{A}{x} &=\detname{A-xI_n}&& \knowl{./knowl/definition-CP.html}{\text{Definition CP}}\\ &=\detname{\transpose{\left(A-xI_n\right)}}&& \knowl{./knowl/theorem-DT.html}{\text{Theorem DT}}\\ &=\detname{\transpose{A}-\transpose{\left(xI_n\right)}}&& \knowl{./knowl/theorem-TMA.html}{\text{Theorem TMA}}\\ &=\detname{\transpose{A}-x\transpose{I_n}}&& \knowl{./knowl/theorem-TMSM.html}{\text{Theorem TMSM}}\\ &=\detname{\transpose{A}-xI_n}&& \knowl{./knowl/definition-IM.html}{\text{Definition IM}}\\ &=\charpoly{\transpose{A}}{x}&& \knowl{./knowl/definition-CP.html}{\text{Definition CP}}\text{.} \end{align*}

So $$A$$ and $$\transpose{A}$$ have the same characteristic polynomial, and by Theorem EMRCP, their eigenvalues are identical and have equal algebraic multiplicities. Notice that what we have proved here is a bit stronger than the stated conclusion in the theorem.

If a matrix has only real entries, then the computation of the characteristic polynomial (Definition CP) will result in a polynomial with coefficients that are real numbers. Complex numbers could result as roots of this polynomial, but they are roots of quadratic factors with real coefficients, and as such, come in conjugate pairs. The next theorem proves this, and a bit more, without mentioning the characteristic polynomial.

We have

\begin{align*} A\conjugate{\vect{x}}&=\conjugate{A}\conjugate{\vect{x}}&& A\text{ has real entries}\\ &=\conjugate{A\vect{x}}&& \knowl{./knowl/theorem-MMCC.html}{\text{Theorem MMCC}}\\ &=\conjugate{\lambda\vect{x}}&& \knowl{./knowl/definition-EEM.html}{\text{Definition EEM}}\\ &=\conjugate{\lambda}\conjugate{\vect{x}}&& \knowl{./knowl/theorem-CRSM.html}{\text{Theorem CRSM}}\text{.} \end{align*}

So $$\conjugate{\vect{x}}$$ is an eigenvector of $$A$$ for the eigenvalue $$\conjugate{\lambda}\text{.}$$

This phenomenon is amply illustrated in Example CEMS6, where the four complex eigenvalues come in two pairs, and the two basis vectors of the eigenspaces are complex conjugates of each other. Theorem ERMCP can be a time-saver for computing eigenvalues and eigenvectors of real matrices with complex eigenvalues, since the conjugate eigenvalue and eigenspace can be inferred from the theorem rather than computed.

SubsectionMEMultiplicities of Eigenvalues

A polynomial of degree $$n$$ will have exactly $$n$$ roots. From this fact about polynomial equations we can say more about the algebraic multiplicities of eigenvalues.

We will prove a more general result by induction (Proof Technique I). Then the theorem will be true as a special case. We will carefully state this result as a proposition indexed by $$m\text{,}$$ $$m\geq 1\text{.}$$

$$P(m)\text{:}$$ Suppose that $$A$$ is an $$m\times m$$ matrix whose entries are complex numbers or linear polynomials in the variable $$x$$ of the form $$c-x\text{,}$$ where $$c$$ is a complex number. Suppose further that there are exactly $$k$$ entries that contain $$x$$ and that no row or column contains more than one such entry. Then, when $$k=m\text{,}$$ $$\detname{A}$$ is a polynomial in $$x$$ of degree $$m\text{,}$$ with leading coefficient $$\pm 1\text{,}$$ and when $$k\lt m\text{,}$$ $$\detname{A}$$ is a polynomial in $$x$$ of degree $$k$$ or less.

Base Case: Suppose $$A$$ is a $$1\times 1$$ matrix. Then its determinant is equal to the lone entry (Definition DM). When $$k=m=1\text{,}$$ the entry is of the form $$c-x\text{,}$$ a polynomial in $$x$$ of degree $$m=1$$ with leading coefficient $$-1\text{.}$$ When $$k\lt m\text{,}$$ then $$k=0$$ and the entry is simply a complex number, a polynomial of degree $$0\leq k\text{.}$$ So $$P(1)$$ is true.

Induction Step: Assume $$P(m)$$ is true, and that $$A$$ is an $$(m+1)\times(m+1)$$ matrix with $$k$$ entries of the form $$c-x\text{.}$$ There are two cases to consider.

Suppose $$k=m+1\text{.}$$ Then every row and every column will contain an entry of the form $$c-x\text{.}$$ Suppose that for the first row, this entry is in column $$t\text{.}$$ Compute the determinant of $$A$$ by an expansion about this first row (Definition DM). The term associated with entry $$t$$ of this row will be of the form

\begin{equation*} (c-x)(-1)^{1+t}\detname{\submatrix{A}{1}{t}}\text{.} \end{equation*}

The submatrix $$\submatrix{A}{1}{t}$$ is an $$m\times m$$ matrix with $$k=m$$ terms of the form $$c-x\text{,}$$ no more than one per row or column. By the induction hypothesis, $$\detname{\submatrix{A}{1}{t}}$$ will be a polynomial in $$x$$ of degree $$m$$ with coefficient $$\pm 1\text{.}$$ So this entire term is then a polynomial of degree $$m+1$$ with leading coefficient $$\pm 1\text{.}$$

The remaining terms (which constitute the sum that is the determinant of $$A$$) are products of complex numbers from the first row with cofactors built from submatrices that lack the first row of $$A$$ and lack some column of $$A\text{,}$$ other than column $$t\text{.}$$ As such, these submatrices are $$m\times m$$ matrices with $$k=m-1\lt m$$ entries of the form $$c-x\text{,}$$ no more than one per row or column. Applying the induction hypothesis, we see that these terms are polynomials in $$x$$ of degree $$m-1$$ or less. Adding the single term from the entry in column $$t$$ with all these others, we see that $$\detname{A}$$ is a polynomial in $$x$$ of degree $$m+1$$ and leading coefficient $$\pm 1\text{.}$$

The second case occurs when $$k\lt m+1\text{.}$$ Now there is a row of $$A$$ that does not contain an entry of the form $$c-x\text{.}$$ We consider the determinant of $$A$$ by expanding about this row (Theorem DER), whose entries are all complex numbers. The cofactors employed are built from submatrices that are $$m\times m$$ matrices with either $$k$$ or $$k-1$$ entries of the form $$c-x\text{,}$$ no more than one per row or column. In either case, $$k\leq m\text{,}$$ and we can apply the induction hypothesis to see that the determinants computed for the cofactors are all polynomials of degree $$k$$ or less. Summing these contributions to the determinant of $$A$$ yields a polynomial in $$x$$ of degree $$k$$ or less, as desired.

Definition CP tells us that the characteristic polynomial of an $$n\times n$$ matrix is the determinant of a matrix having exactly $$n$$ entries of the form $$c-x\text{,}$$ no more than one per row or column. As such we can apply $$P(n)$$ to see that the characteristic polynomial has degree $$n\text{.}$$

By the definition of the algebraic multiplicity (Definition AME), we can factor the characteristic polynomial as

\begin{equation*} \charpoly{A}{x}=c(x-\lambda_1)^{\algmult{A}{\lambda_1}}(x-\lambda_2)^{\algmult{A}{\lambda_2}}(x-\lambda_3)^{\algmult{A}{\lambda_3}}\cdots(x-\lambda_k)^{\algmult{A}{\lambda_k}} \end{equation*}

where $$c$$ is a nonzero constant. (We could prove that $$c=(-1)^{n}\text{,}$$ but we do not need that specificity right now. See Exercise PEE.T30) The left-hand side is a polynomial of degree $$n$$ by Theorem DCP and the right-hand side is a polynomial of degree $$\sum_{i=1}^{k}\algmult{A}{\lambda_i}\text{.}$$ So the equality of the polynomials' degrees gives the equality $$\sum_{i=1}^{k}\algmult{A}{\lambda_i}=n\text{.}$$

Since $$\lambda$$ is an eigenvalue of $$A\text{,}$$ there is an eigenvector of $$A$$ for $$\lambda\text{,}$$ $$\vect{x}\text{.}$$ Then $$\vect{x}\in\eigenspace{A}{\lambda}\text{,}$$ so $$\geomult{A}{\lambda}\geq 1\text{,}$$ since we can extend $$\set{\vect{x}}$$ into a basis of $$\eigenspace{A}{\lambda}$$ (Theorem ELIS).

To show $$\geomult{A}{\lambda}\leq\algmult{A}{\lambda}$$ is the most involved portion of this proof. To this end, let $$g=\geomult{A}{\lambda}$$ and let $$\vectorlist{x}{g}$$ be a basis for the eigenspace of $$\lambda\text{,}$$ $$\eigenspace{A}{\lambda}\text{.}$$ Construct another $$n-g$$ vectors, $$\vectorlist{y}{n-g}\text{,}$$ so that

\begin{equation*} \set{\vectorlist{x}{g},\,\vectorlist{y}{n-g}} \end{equation*}

is a basis of $$\complex{n}\text{.}$$ This can be done by repeated applications of Theorem ELIS.

Finally, define a matrix $$S$$ by

\begin{equation*} S=[\vect{x}_1|\vect{x}_2|\vect{x}_3|\ldots|\vect{x}_g|\vect{y}_1|\vect{y}_2|\vect{y}_3|\ldots|\vect{y}_{n-g}]=[\vect{x}_1|\vect{x}_2|\vect{x}_3|\ldots|\vect{x}_g|R] \end{equation*}

where $$R$$ is an $$n\times(n-g)$$ matrix whose columns are $$\vectorlist{y}{n-g}\text{.}$$ The columns of $$S$$ are linearly independent by design, so $$S$$ is nonsingular (Theorem NMLIC) and therefore invertible (Theorem NI).

Then

\begin{align*} \matrixcolumns{e}{n}&=I_n\\ &=\inverse{S}S\\ &=\inverse{S}[\vect{x}_1|\vect{x}_2|\vect{x}_3|\ldots|\vect{x}_g|R]\\ &=[\inverse{S}\vect{x}_1|\inverse{S}\vect{x}_2|\inverse{S}\vect{x}_3|\ldots|\inverse{S}\vect{x}_g|\inverse{S}R]\text{.} \end{align*}

So

\begin{equation*} \inverse{S}\vect{x}_i=\vect{e}_i\quad 1\leq i\leq g\tag{*}\text{.} \end{equation*}

Preparations in place, we compute the characteristic polynomial of $$A\text{,}$$

\begin{align*} \charpoly{A}{x}&=\detname{A-xI_n}&& \knowl{./knowl/definition-CP.html}{\text{Definition CP}}\\ &=1\detname{A-xI_n}&& \knowl{./knowl/property-OCN.html}{\text{Property OCN}}\\ &=\detname{I_n}\detname{A-xI_n}&& \knowl{./knowl/definition-DM.html}{\text{Definition DM}}\\ &=\detname{\inverse{S}S}\detname{A-xI_n}&& \knowl{./knowl/definition-MI.html}{\text{Definition MI}}\\ &=\detname{\inverse{S}}\detname{S}\detname{A-xI_n}&& \knowl{./knowl/theorem-DRMM.html}{\text{Theorem DRMM}}\\ &=\detname{\inverse{S}}\detname{A-xI_n}\detname{S}&& \knowl{./knowl/property-CMCN.html}{\text{Property CMCN}}\\ &=\detname{\inverse{S}\left(A-xI_n\right)S}&& \knowl{./knowl/theorem-DRMM.html}{\text{Theorem DRMM}}\\ &=\detname{\inverse{S}AS-\inverse{S}xI_nS}&& \knowl{./knowl/theorem-MMDAA.html}{\text{Theorem MMDAA}}\\ &=\detname{\inverse{S}AS-x\inverse{S}I_nS}&& \knowl{./knowl/theorem-MMSMM.html}{\text{Theorem MMSMM}}\\ &=\detname{\inverse{S}AS-x\inverse{S}S}&& \knowl{./knowl/theorem-MMIM.html}{\text{Theorem MMIM}}\\ &=\detname{\inverse{S}AS-xI_n}&& \knowl{./knowl/definition-MI.html}{\text{Definition MI}}\\ &=\charpoly{\inverse{S}AS}{x}&& \knowl{./knowl/definition-CP.html}{\text{Definition CP}}\text{.} \end{align*}

What can we learn then about the matrix $$\inverse{S}AS\text{?}$$

\begin{align*} \inverse{S}AS&=\inverse{S}A[\vect{x}_1|\vect{x}_2|\vect{x}_3|\ldots|\vect{x}_g|R]\\ &=\inverse{S}[A\vect{x}_1|A\vect{x}_2|A\vect{x}_3|\ldots|A\vect{x}_g|AR]&& \knowl{./knowl/definition-MM.html}{\text{Definition MM}}\\ &=\inverse{S}[\lambda\vect{x}_1|\lambda\vect{x}_2|\lambda\vect{x}_3|\ldots|\lambda\vect{x}_g|AR]&& \knowl{./knowl/definition-EEM.html}{\text{Definition EEM}}\\ &=[\inverse{S}\lambda\vect{x}_1|\inverse{S}\lambda\vect{x}_2|\inverse{S}\lambda\vect{x}_3|\ldots|\inverse{S}\lambda\vect{x}_g|\inverse{S}AR]&& \knowl{./knowl/definition-MM.html}{\text{Definition MM}}\\ &=[\lambda\inverse{S}\vect{x}_1|\lambda\inverse{S}\vect{x}_2|\lambda\inverse{S}\vect{x}_3|\ldots|\lambda\inverse{S}\vect{x}_g|\inverse{S}AR]&& \knowl{./knowl/theorem-MMSMM.html}{\text{Theorem MMSMM}}\\ &=[\lambda\vect{e}_1|\lambda\vect{e}_2|\lambda\vect{e}_3|\ldots|\lambda\vect{e}_g|\inverse{S}AR]&& \inverse{S}S=I_n,\text{ (*) above} \end{align*}

Now imagine computing the characteristic polynomial of $$A$$ by computing the characteristic polynomial of $$\inverse{S}AS$$ using the form just obtained. The first $$g$$ columns of $$\inverse{S}AS$$ are all zero, save for a $$\lambda$$ on the diagonal. So if we compute the determinant by expanding about the first column, successively, we will get successive factors of $$(\lambda-x)\text{.}$$ More precisely, let $$T$$ be the square matrix of size $$n-g$$ that is formed from the last $$n-g$$ rows and last $$n-g$$ columns of $$\inverse{S}AR\text{.}$$ Then

\begin{equation*} \charpoly{A}{x}=\charpoly{\inverse{S}AS}{x}=(\lambda-x)^g\charpoly{T}{x}\text{.} \end{equation*}

This says that $$(x-\lambda)$$ is a factor of the characteristic polynomial at least $$g$$ times, so the algebraic multiplicity of $$\lambda$$ as an eigenvalue of $$A$$ is greater than or equal to $$g$$ (Definition AME). In other words,

\begin{equation*} \geomult{A}{\lambda}=g\leq\algmult{A}{\lambda} \end{equation*}

as desired.

Theorem NEM says that the sum of the algebraic multiplicities for all the eigenvalues of $$A$$ is equal to $$n\text{.}$$ Since the algebraic multiplicity is a positive quantity, no single algebraic multiplicity can exceed $$n$$ without the sum of all of the algebraic multiplicities doing the same.

Suppose that $$A$$ has $$k$$ distinct eigenvalues, $$\scalarlist{\lambda}{k}\text{.}$$ Then

\begin{align*} k&=\sum_{i=1}^{k}1\\ &\leq\sum_{i=1}^{k}\algmult{A}{\lambda_i}&& \knowl{./knowl/theorem-ME.html}{\text{Theorem ME}}\\ &=n&& \knowl{./knowl/theorem-NEM.html}{\text{Theorem NEM}}\text{.} \end{align*}

SubsectionEHMEigenvalues of Hermitian Matrices

Recall that a matrix is Hermitian (or self-adjoint) if $$A=\adjoint{A}$$ (Definition HM). In the case where $$A$$ is a matrix whose entries are all real numbers, being Hermitian is identical to being symmetric (Definition SYM). Keep this in mind as you read the next two theorems. Their hypotheses could be changed to “suppose $$A$$ is a real symmetric matrix.”

Let $$\vect{x}\neq\zerovector$$ be one eigenvector of $$A$$ for the eigenvalue $$\lambda\text{.}$$ Then

\begin{align*} \left(\lambda - \conjugate{\lambda}\right)\innerproduct{\vect{x}}{\vect{x}} &=\lambda\innerproduct{\vect{x}}{\vect{x}} - \conjugate{\lambda}\innerproduct{\vect{x}}{\vect{x}}&& \knowl{./knowl/property-DCN.html}{\text{Property DCN}}\\ &=\innerproduct{\vect{x}}{\lambda\vect{x}} - \innerproduct{\lambda\vect{x}}{\vect{x}}&& \knowl{./knowl/theorem-IPSM.html}{\text{Theorem IPSM}}\\ &=\innerproduct{\vect{x}}{A\vect{x}} - \innerproduct{A\vect{x}}{\vect{x}}&& \knowl{./knowl/definition-EEM.html}{\text{Definition EEM}}\\ &=\innerproduct{A\vect{x}}{\vect{x}} - \innerproduct{A\vect{x}}{\vect{x}}&& \knowl{./knowl/theorem-HMIP.html}{\text{Theorem HMIP}}\\ &=0\text{.} \end{align*}

Since $$\vect{x}\neq\zerovector\text{,}$$ by Theorem PIP we know $$\innerproduct{\vect{x}}{\vect{x}}\neq 0\text{.}$$ Then by Theorem ZPZT, $$\lambda - \conjugate{\lambda}=0\text{,}$$ and so $$\lambda = \conjugate{\lambda}\text{.}$$ If a complex number is equal to its conjugate, then it has a complex part equal to zero, and therefore is a real number.

Notice the key step of this proof is the ability to pitch a Hermitian matrix from one side of the inner product to the other.

Look back and compare Example ESMS4 and Example CEMS6. In Example CEMS6 the matrix has only real entries, yet the characteristic polynomial has roots that are complex numbers, and so the matrix has complex eigenvalues. However, in Example ESMS4, the matrix has only real entries, but is also symmetric, and hence Hermitian. So by Theorem HMRE, we were guaranteed eigenvalues that are real numbers.

In many physical problems, a matrix of interest will be real and symmetric, or Hermitian. Then if the eigenvalues are to represent physical quantities of interest, Theorem HMRE guarantees that these values will not be complex numbers.

The eigenvectors of a Hermitian matrix also enjoy a pleasing property that we will exploit later.

Let $$\vect{x}$$ be an eigenvector of $$A$$ for $$\lambda$$ and let $$\vect{y}$$ be an eigenvector of $$A$$ for a different eigenvalue $$\rho\text{.}$$ So we have $$\lambda-\rho\neq 0\text{.}$$ Then

\begin{align*} \left(\lambda-\rho\right)\innerproduct{\vect{x}}{\vect{y}} &=\lambda\innerproduct{\vect{x}}{\vect{y}}-\rho\innerproduct{\vect{x}}{\vect{y}}&& \knowl{./knowl/property-DCN.html}{\text{Property DCN}}\\ &=\innerproduct{\conjugate{\lambda}\vect{x}}{\vect{y}}-\innerproduct{\vect{x}}{\rho\vect{y}}&& \knowl{./knowl/theorem-IPSM.html}{\text{Theorem IPSM}}\\ &=\innerproduct{\lambda\vect{x}}{\vect{y}}-\innerproduct{\vect{x}}{\rho\vect{y}}&& \knowl{./knowl/theorem-HMRE.html}{\text{Theorem HMRE}}\\ &=\innerproduct{A\vect{x}}{\vect{y}}-\innerproduct{\vect{x}}{A\vect{y}}&& \knowl{./knowl/definition-EEM.html}{\text{Definition EEM}}\\ &=\innerproduct{A\vect{x}}{\vect{y}}-\innerproduct{A\vect{x}}{\vect{y}}&& \knowl{./knowl/theorem-HMIP.html}{\text{Theorem HMIP}}\\ &=0\text{.} \end{align*}

Because $$\lambda$$ and $$\rho$$ are presumed to be different, $$\lambda-\rho\neq 0\text{,}$$ and Theorem ZPZT implies that $$\innerproduct{\vect{x}}{\vect{y}}=0\text{.}$$ In other words, $$\vect{x}$$ and $$\vect{y}$$ are orthogonal vectors according to Definition OV.

Notice again how the key step in this proof is the fundamental property of a Hermitian matrix (Theorem HMIP) — the ability to swap $$A$$ across the two arguments of the inner product. Notice too, that we can always apply the Gram-Schmidt procedure (Theorem GSP) to any basis of any eigenspace, along with scaling the resulting the orthogonal vectors by their norm to arrive at an orthonormal basis of the eigenspace. For a Hermitian matrix, pairs of eigenvectors from different eigenspaces are also orthogonal. If we dumped all these basis vectors into one big set it would be an orthonormal set. We will build on these results and continue to see some more interesting properties in Section OD.

1.

How can you identify a nonsingular matrix just by looking at its eigenvalues?

2.

How many different eigenvalues may a square matrix of size $$n$$ have?

3.

What is amazing about the eigenvalues of a Hermitian matrix and why is it amazing?

ExercisesPEEExercises

T10.

Suppose that $$A$$ is a square matrix. Prove that the constant term of the characteristic polynomial of $$A$$ is equal to the determinant of $$A\text{.}$$

Solution.

Suppose that the characteristic polynomial of $$A$$ is

\begin{equation*} \charpoly{A}{x}=a_0+a_1x+a_2x^2+\dots+a_nx^n\text{.} \end{equation*}

Then

\begin{align*} a_0 &=a_0+a_1(0)+a_2(0)^2+\dots+a_n(0)^n\\ &=\charpoly{A}{0}\\ &=\detname{A-0I_n}&& \knowl{./knowl/definition-CP.html}{\text{Definition CP}}\\ &=\detname{A}\text{.} \end{align*}

T20.

Suppose that $$A$$ is a square matrix. Prove that a single vector may not be an eigenvector of $$A$$ for two different eigenvalues.

Solution.

Suppose that the vector $$\vect{x}\neq\zerovector$$ is an eigenvector of $$A$$ for the two eigenvalues $$\lambda$$ and $$\rho\text{,}$$ where $$\lambda\neq\rho\text{.}$$ Then $$\lambda-\rho\neq 0\text{,}$$ and we also have

\begin{align*} \zerovector &=A\vect{x}-A\vect{x}&& \knowl{./knowl/property-AIC.html}{\text{Property AIC}}\\ &=\lambda\vect{x}-\rho\vect{x}&& \knowl{./knowl/definition-EEM.html}{\text{Definition EEM}}\\ &=(\lambda-\rho)\vect{x}&& \knowl{./knowl/property-DSAC.html}{\text{Property DSAC}}\text{.} \end{align*}

By Theorem SMEZV, either $$\lambda-\rho=0$$ or $$\vect{x}=\zerovector\text{,}$$ which are both contradictions.

T21.

Suppose that $$A$$ is a square matrix of size $$n$$ with an eigenvalue $$\lambda\text{.}$$ For $$\lambda\in\complexes\text{,}$$ prove that $$\lambda+\alpha$$ is an eigenvalue of the matrix $$A+\alpha I_n\text{.}$$

1. Construct a proof that begins by employing a vector $$\vect{x}$$ that is an eigenvector of $$A$$ for $$\lambda\text{.}$$

2. Construct a proof that establishes that $$\lambda+\alpha$$ is a root of the characteristic polynomial of $$A+\alpha I_n\text{.}$$

T22.

Suppose that $$U$$ is a unitary matrix with eigenvalue $$\lambda\text{.}$$ Prove that $$\lambda$$ has modulus 1, i.e. $$\modulus{\lambda}=1\text{.}$$ This says that all of the eigenvalues of a unitary matrix lie on the unit circle of the complex plane.

T30.

Theorem DCP tells us that the characteristic polynomial of a square matrix of size $$n$$ has degree $$n\text{.}$$ By suitably augmenting the proof of Theorem DCP prove that the coefficient of $$x^n$$ in the characteristic polynomial is $$(-1)^n\text{.}$$

T50.

Theorem EIM says that if $$\lambda$$ is an eigenvalue of the nonsingular matrix $$A\text{,}$$ then $$\frac{1}{\lambda}$$ is an eigenvalue of $$\inverse{A}\text{.}$$ Write an alternate proof of this theorem using the characteristic polynomial and without making reference to an eigenvector of $$A$$ for $$\lambda\text{.}$$

Solution.

Since $$\lambda$$ is an eigenvalue of a nonsingular matrix, $$\lambda\neq 0$$ (Theorem SMZE). $$A$$ is invertible (Theorem NI), and so $$-\lambda A$$ is invertible (Theorem MISM). Thus $$-\lambda A$$ is nonsingular (Theorem NI) and $$\detname{-\lambda A}\neq 0$$ (Theorem SMZD). We have

\begin{align*} \charpoly{\inverse{A}}{\frac{1}{\lambda}} &=\detname{\inverse{A}-\frac{1}{\lambda}I_n}&& \knowl{./knowl/definition-CP.html}{\text{Definition CP}}\\ &=1\detname{\inverse{A}-\frac{1}{\lambda}I_n}&& \knowl{./knowl/property-OCN.html}{\text{Property OCN}}\\ &=\frac{1}{\detname{-\lambda A}} \detname{-\lambda A}\detname{\inverse{A}-\frac{1}{\lambda}I_n}&& \knowl{./knowl/property-MICN.html}{\text{Property MICN}}\\ &=\frac{1}{\detname{-\lambda A}} \detname{\left(-\lambda A\right)\left(\inverse{A}-\frac{1}{\lambda}I_n\right)}&& \knowl{./knowl/theorem-DRMM.html}{\text{Theorem DRMM}}\\ &=\frac{1}{\detname{-\lambda A}} \detname{-\lambda A\inverse{A}-\left(-\lambda A\right)\frac{1}{\lambda}I_n}&& \knowl{./knowl/theorem-MMDAA.html}{\text{Theorem MMDAA}}\\ &=\frac{1}{\detname{-\lambda A}} \detname{-\lambda I_n-\left(-\lambda A\right)\frac{1}{\lambda}I_n}&& \knowl{./knowl/definition-MI.html}{\text{Definition MI}}\\ &=\frac{1}{\detname{-\lambda A}} \detname{-\lambda I_n+\lambda\frac{1}{\lambda}AI_n}&& \knowl{./knowl/theorem-MMSMM.html}{\text{Theorem MMSMM}}\\ &=\frac{1}{\detname{-\lambda A}} \detname{-\lambda I_n+1AI_n}&& \knowl{./knowl/property-MICN.html}{\text{Property MICN}}\\ &=\frac{1}{\detname{-\lambda A}} \detname{-\lambda I_n+AI_n}&& \knowl{./knowl/property-OCN.html}{\text{Property OCN}}\\ &=\frac{1}{\detname{-\lambda A}} \detname{-\lambda I_n+A}&& \knowl{./knowl/theorem-MMIM.html}{\text{Theorem MMIM}}\\ &=\frac{1}{\detname{-\lambda A}} \detname{A-\lambda I_n}&& \knowl{./knowl/property-ACM.html}{\text{Property ACM}}\\ &=\frac{1}{\detname{-\lambda A}} \charpoly{A}{\lambda}&& \knowl{./knowl/definition-CP.html}{\text{Definition CP}}\\ &=\frac{1}{\detname{-\lambda A}}\,0&& \knowl{./knowl/theorem-EMRCP.html}{\text{Theorem EMRCP}}\\ &=0&& \knowl{./knowl/property-ZCN.html}{\text{Property ZCN}}\text{.} \end{align*}

So $$\frac{1}{\lambda}$$ is a root of the characteristic polynomial of $$\inverse{A}$$ and so is an eigenvalue of $$\inverse{A}\text{.}$$ This proof is due to Sara Bucht.