Suppose A (size n) has n linearly independent eigenvectors v1,…,vn with eigenvalues λ1,…,λn. Stack the eigenvectors as columns of V and eigenvalues on the diagonal of Λ:
V=[v1v2⋯vn],Λ=diag(λ1,…,λn)
Look at what A does to all columns at once:
AV=[Av1⋯Avn]=[λ1v1⋯λnvn]=VΛ
So AV=VΛ. If V is invertible (needs independent eigenvectors), right-multiply... actually multiply by V−1:
WHY huge: raising a diagonal matrix to a power is just powering the diagonal entries — O(n) instead of repeated O(n3) multiplies. This is how Markov chains find steady states and how we analyze long-term behavior.
Why do we require det(A−λI)=0? → so (A−λI)v=0 has a nonzero solution (matrix singular).
What must V contain and satisfy? → columns = eigenvectors; must be linearly independent (invertible).
Formula for Ak? → VΛkV−1.
When is V orthogonal? → when A is symmetric (spectral theorem).
Recall Explain to a 12-year-old (Feynman)
Imagine a big rubber sheet. When you pull on it, most drawn arrows on it swing to point somewhere new. But a few magic arrows keep pointing the same way — they just get longer or shorter. Those magic arrows are eigenvectors, and "how much longer" is the eigenvalue. If you know all the magic arrows, you understand everything the pull does: just stretch along each magic arrow. Eigendecomposition is writing down that list of magic arrows and their stretch amounts.
Dekho, ek matrix A basically ek machine hai jo vectors ko transform karti hai — kisi ko ghumati hai, kisi ko stretch karti hai. Lekin kuch special vectors hote hain jinko A sirf lamba ya chhota karti hai, direction bilkul same rehti hai. Yehi hote hain eigenvectors, aur kitna stretch hua woh number hai eigenvalueλ. Equation simple hai: Av=λv.
Ab eigenvalue nikalne ke liye hum Av=λv ko rearrange karte hain: (A−λI)v=0. Yahan I isliye ghusaya kyunki λ ek scalar hai aur usko matrix banana zaroori tha subtract karne ke liye. Iss system ka nonzero solution tabhi milega jab matrix singular ho, matlab det(A−λI)=0. Iss equation ke roots hi eigenvalues hain. Phir har λ ko wapas daal ke null space solve karo — woh eigenvector hai.
Jab tumhare paas n independent eigenvectors ho, unhe columns me rakh do V me, aur eigenvalues ko diagonal Λ me. Toh A=VΛV−1. Isko yaad rakho "Enter, Stretch, Exit" — pehle eigen-coordinates me jao (V−1), phir har axis ko stretch karo (Λ), phir wapas aao (V). Sabse bada faayda: Ak=VΛkV−1, matlab power sirf diagonal pe lagao, bahut fast!
Yeh cheez ML me har jagah hai — PCA me covariance matrix (jo symmetric hoti hai) ka eigendecomposition hi principal directions deta hai. Ek warning: har matrix decompose nahi hoti — agar independent eigenvectors kam pad jaayein (defective matrix), toh SVD use karna padta hai. Bas yeh dhyaan rakho aur tum set ho.