Spectral Decomposition

deejayosamu·2026년 1월 16일

다변량통계

목록 보기
4/11

Let ARp×pA \in \mathbb{R}^{p \times p} be a symmetric matrix with eigenvalues λ1λp\lambda_1 \geq \cdots \geq \lambda_p and eigenvectors ξ1ξp(Rp)\xi_1 \geq \cdots \geq \xi_p (\in \mathbb{R}^p)
Then for P=(ξ1,,ξp)Rp×pP=(\xi_1,\cdots,\xi_p) \in \mathbb{R}^{p \times p} and Λ=diag(λj)Rp×p\Lambda=diag(\lambda_j) \in \mathbb{R}^{p \times p}
A=i=1pλiξiξiT=PΛPTA=\sum_{i=1}^{p} \lambda_i \xi_i \xi_i^T=P \Lambda P^T
pf)
Note that Aξi=λiξi (i=1,,p)A \xi_i=\lambda_i \xi_i \space (i=1,\cdots,p)
Since PP is orthogonal, A=PΛPTA=P \Lambda P^T

  • Remarks

    Let ARp×pA \in \mathbb{R}^{p \times p} be symmetric with eigenvalues λ1,,λp\lambda_1, \cdots, \lambda_p

    A1A^{-1} has eigenvalues 1/λ1,,1/λp1/\lambda_1, \cdots, 1/\lambda_p (eigenvector는 동일)
    Ax=λx=>x=λA1x=>λ1x=A1x\because Ax=\lambda x => x=\lambda A^{-1}x => \lambda^{-1}x=A^{-1}x

    tr(A)=i=1pλitr(A)=\sum_{i=1}^{p} \lambda_i and A=i=1pλi|A|=\prod_{i=1}^p \lambda_i
    \because Using spectral decomposition, tr(A)=tr(PΛPT)=tr(Λ)tr(A)=tr(P \Lambda P^T)=tr(\Lambda) and A=PΛPT=PΛP1=Λ=i=1pλi|A|=|P \Lambda P^T|=|P| |\Lambda| |P^{-1}|=|\Lambda|=\prod_{i=1}^p \lambda_i

    rank(A)rank(A) is the number of nonzero eigenvalues of AA
    \because Let kk be the number of nonzero eigenvalues of AA
    rank
    ++ Let AA be positive definite

    ④ Eigenvalues of AA are positive
    \because By spectral decomposition, we have xTAx=xTPΛPTx=j=1pλj(ξjTx)2x^T A x=x^T P \Lambda P^T x=\sum_{j=1}^{p} \lambda_j (\xi_j^T x)^2
    If λj0\lambda_{j^*} \leq 0 for some jj^*, x=ξjx=\xi_{j^*} gives xTAx=λj0x^T Ax=\lambda_{j^*} \leq 0 => contradiction

    A1/2A^{1/2} is given by A1/2=PΛ1/2PTA^{1/2}=P \Lambda^{1/2} P^T where Λ1/2=diag(λ1,,λp)\Lambda^{1/2}=diag(\sqrt\lambda_1, \cdots, \sqrt \lambda_p)

Theorem)
Suppose ARp×pA \in \mathbb{R}^{p \times p} is nonnegative definite with eigenvalues λ1λp0\lambda_1 \geq \cdots \geq \lambda_p \geq 0
Then, for x0x \neq 0, λpxTAxxTxλ1\lambda_p \leq \frac{x^T Ax}{x^T x} \leq \lambda_1
pf)
theorem_pf

0개의 댓글