Eigenvector and Eigenvalues

Theorem. Given a symmetric matrix A∈Rn×nA \in \mathbb{R}^{n \times n} that is positive semi-definitive, there exists an orthogonal matrix V∈Rn×nV \in \mathbb{R}^{n \times n}, and a diagonal matrix Σ∈Rn×n\Sigma \in \mathbb{R}^{n \times n} with all entries non-negative, such that

A=VΣVT. A = V \Sigma V^T.

We need the following fact for the proof.

Fact. If a function f:S→Rf: S \rightarrow \mathbb{R} is continuous on a closed set SS, then it achieves maximum value at some point of SS.

Proof. Consider the maximization problem:

max⁡12xTAxs.t.xTx=1 \begin{aligned} &\max & \frac{1}{2} x^T A x \\ &s.t. & x^T x = 1 \end{aligned}

Question to ponder: Explain the geometric meaning of the term xTAxx^T A x.

The corresponding Lagrange function is given by

L(x,λ)=12xTAx−λ(xTx−1) L(x, \lambda) = \frac{1}{2} x^T A x - \lambda (x^T x - 1)

Taking derivative with respect to xx gives

∂L∂x=Ax−λx. \frac{\partial L}{\partial x} = Ax - \lambda x.

Existence of maximum solution: Note that {x:xTx=1}\{ x : x^T x = 1 \} is defined the boundary of the unit cycle, which is a closed set. Combined with the fact that f≐(1/2)⋅xTAxf \doteq (1 / 2) \cdot x^T A x is a continuous function (with respect to xx), it takes its maximum value at some point in {x:xTx=1}\{ x : x^T x = 1 \}. Denote this point v1∈{x:xTx=1}v_1 \in \{ x : x^T x = 1 \}.

Gradient at the maximum point: At this point, it holds that

∂L∂x∣x=v1=Av1−λv1=0, \frac{\partial L}{\partial x} \mid_{x = v_1} = Av_1 - \lambda v_1 = 0,

which implies that v1v_1 is an eigenvector of AA. Denote its corresponding eigenvalue λ1\lambda_1.

Question to ponder: Prove that λ1\lambda_1 is the largest eigenvalue.

After finding v1v_1, consider a new optimization problem:

max⁡12xTAxs.t.xTx=1v1Tx=0 \begin{aligned} &\max & \frac{1}{2} x^T A x \\ &s.t. & x^T x = 1 \\ & & v_1^T x = 0 \end{aligned}

That is, we want to find a new unit vector that is orthogonal to v1v_1 while maximizing 12xTAx\frac{1}{2} x^T A x.

Existence of maximum solution. By the lemma and the facts:

  1. {x:xTx=1}∩{x:v1Tx=0}\{x : x^T x = 1\} \cap \{ x : v_1^T x = 0 \} is closed.
  2. 12xTAx\frac{1}{2} x^T A x is continuous.

Denote this maximum point v2v_2. Now consider the new Lagrange function for the optimization

L(x,λ)=12xTAx−λ(xTx−1)−β(v1Tx−0). L(x, \lambda) = \frac{1}{2} x^T A x - \lambda (x^T x - 1) - \beta(v_1^T x - 0).

Gradient at the maximum point: At this point, it holds that

∂L∂x∣x=v2=Av2−λv2−βv1=0. \frac{\partial L}{\partial x} \mid_{x = v_2} = Av_2 - \lambda v_2 - \beta v_1 = 0.

Left multiplying v1Tv_1^T gives

(v1TA)v2−λv1Tv2−βv1Tv1=(Av1)Tv2−β=λ1v1Tv2−β=−β=0, (v_1^T A)v_2 - \lambda v_1^T v_2 - \beta v_1^T v_1 = (A v_1)^T v_2 - \beta = \lambda_1 v_1^T v_2 - \beta = -\beta = 0,

where v1TA=(Av1)Tv_1^T A = (A v_1)^T follows from the symmetry of AA, i.e., A=ATA = A^T.

Therefore,

Av2−λv2=0. Av_2 - \lambda v_2 = 0.

It follows that v2v_2 is an eigenvector of AA. Denote its associated eigenvalue λ2\lambda_2.

In a similar manner, we can define

max⁡12xTAxs.t.xTx=1v1Tx=0v2Tx=0 \begin{aligned} &\max & \frac{1}{2} x^T A x \\ &s.t. & x^T x = 1 \\ & & v_1^T x = 0 \\ & & v_2^T x = 0 \\ \end{aligned}

There exists an optimal point v3v_3. The derivative of Lagrange function at this point is

∂L∂x∣x=v3=Av3−λv3−β1v1−β2v2=0, \frac{\partial L}{\partial x} \mid_{x = v_3} = Av_3 - \lambda v_3 - \beta_1 v_1 - \beta_2 v_2 = 0,

for some parameters λ,β1,β2\lambda, \beta_1, \beta_2. Left multiplying v1Tv_1^T or v2Tv_2^T gives β1=0\beta_1 = 0 and β2=0\beta_2 = 0 so that v3v_3 is an eigenvector of AA.

Following the procedure, we can find nn eigenvectors v1,v2,…,vnv_1, v_2, \ldots, v_n of AA. Define

Σ=[λ1λ2…λn],V=[v1,v2,…,vn]. \Sigma = \begin{bmatrix} \lambda_1 & & \\ &\lambda_2 & \\ & & \ldots \\ & & & \lambda_n \end{bmatrix}, \qquad V = \begin{bmatrix} v_1, v_2, \ldots, v_n \end{bmatrix}.

Then

AV=VΣ. A V = V \Sigma.

Hence,

A=VΣV−1=VΣVT. A = V \Sigma V^{-1} = V \Sigma V^T.

□\square

Singular Vector Decomposition

In general, a matrix A∈Rm×nA \in \mathbb{R}^{m \times n} is considered a linear mapping from Rn\mathbb{R}^n to Rm\mathbb{R}^m. It takes a vector x∈Rnx \in \mathbb{R}^n to a new one Ax∈RmAx \in \mathbb{R}^m. We are interested in the following question.

Is there a vector xx, whose length ratio,

max⁡x∈Rn∥Ax∥∥x∥, \max_{x \in \mathbb{R}^n } \frac{ \Vert Ax \Vert }{ \Vert x \Vert },

is maximized by the mapping?

As ∥x∥\Vert x \Vert is a number,

∥Ax∥∥x∥=∥Ax∥x∥∥. \frac{ \Vert Ax \Vert }{ \Vert x \Vert } = \left\Vert A \frac{x}{\Vert x \Vert} \right\Vert .

Hence, it suffices to consider unit vectors, and the following optimization problem:

max⁡12xTATAxs.t.xTx=1. \begin{aligned} &\max & \frac{1}{2} x^T A^T Ax \\ &s.t. & x^T x = 1 \end{aligned}.

In this case, the objective function is 0.5⋅∥Ax∥2=0.5⋅(Ax)T(Ax)=0.5⋅xTATAx0.5 \cdot \Vert Ax \Vert^2 = 0.5 \cdot (Ax)^T (Ax) = 0.5 \cdot x^T A^T A x.

Following similar process of previous section, we can find a set of eigenvectors v1,v2,…,vnv_1, v_2, \ldots, v_n of ATAA^T A. Define U∈Rm×nU \in \mathbb{R}^{m \times n}

U=AVΣ−1=AV[1/λ11/λ2…1/λn]. U = A V { \sqrt{\Sigma^{-1} } } = AV \begin{bmatrix} \sqrt { 1 / \lambda_1} & & \\ &\sqrt { 1 / \lambda_2} & \\ & & \ldots \\ & & &\sqrt {1/\lambda_n} \end{bmatrix}.

Then UΣVT=AU { \sqrt \Sigma } V^T = A, which is called the singular vector decomposition of AA.

Finally, we can verify that the columns of UU are perpendicular: ∀i≠j,i,j∈[n]\forall i \neq j, i, j \in [n],

uiTuj=(1λiAVi)T(1λjAVj)=1λi1λjviTATAvj=0, u_i^T u_j = \left(\sqrt{\frac{1}{\lambda_i}} A V_i \right)^T \left( \sqrt{\frac{1}{\lambda_j}} A V_j \right) = \sqrt {\frac{1}{\lambda_i}} \sqrt{\frac{1}{\lambda_j}} v_i^T A^T A v_j = 0,

as viv_i and vjv_j are orthogonal eigenvectors of ATAA^T A.