4.5.42 · Maths › Linear Algebra (Full)
Ek normal inverse A − 1 tab hi exist karta hai jab A square aur invertible ho. Lekin real life mein kaafi matrices rectangular hoti hain (m × n ), ya square-lekin-singular. Hum phir bhi A ko "undo" karna chahte hain — jitna best ho sake. ==Pseudoinverse A + == ek inverse ka best possible substitute hai: yeh least-squares, minimum-norm answer deta hai A x = b ke liye, chahe exact inverse exist kare ya na kare.
YEH KYUN IMPORTANT HAI: har overdetermined system (equations zyada, unknowns kam — data fitting) aur har underdetermined system (unknowns zyada, equations kam — many solutions) ko A + solve karta hai.
Definition Moore–Penrose pseudoinverse
Kisi bhi real matrix A ∈ R m × n ke liye, ==pseudoinverse A + ∈ R n × m == woh unique matrix hai jo ye chaar Penrose conditions satisfy karti hai:
A A + A = A
A + A A + = A +
( A A + ) ⊤ = A A + (matlab A A + symmetric hai)
( A + A ) ⊤ = A + A (matlab A + A symmetric hai)
Jab A invertible hoti hai, to A + = A − 1 . Toh yeh inverse ko generalize karta hai.
Ek key geometric fact: A A + ek orthogonal projection hai A ke column space par, aur A + A ek orthogonal projection hai A ke row space par.
Hum A x = b solve karna chahte hain. Agar b column space mein nahi hai, to exact solution exist nahi hoti — toh hum error ∥ A x − b ∥ 2 minimize karte hain (least squares).
Ab A x = b ke infinitely many solutions hain. Hum woh solution choose karte hain jiska norm ∥ x ∥ sabse chhota ho — Kyun? Yeh unique "no wasted energy" solution hai, jo nullspace ke perpendicular hai.
Worked example Example 1 — tall matrix, least squares
A = 1 1 1 0 1 2 , b = 6 0 0 . t = 0 , 1 , 2 par ek line y = c + d t fit karo.
Step 1: A ⊤ A = ( 3 3 3 5 ) . Kyun? Full column rank ⇒ Case 1 use karo.
Step 2: A ⊤ b = ( 6 0 ) .
Step 3: Solve karo ( 3 3 3 5 ) x = ( 6 0 ) . Determinant = 6 .
x = 6 1 ( 5 − 3 − 3 3 ) ( 6 0 ) = ( 5 − 3 )
Best-fit line y = 5 − 3 t . Yeh answer kyun, exact kyun nahi? Koi bhi line ( 0 , 6 ) , ( 1 , 0 ) , ( 2 , 0 ) se nahi guzarti; A + least-squares line deta hai.
Worked example Example 2 — fat matrix, minimum norm
A = ( 1 1 ) , b = ( 2 ) . Solve karo x 1 + x 2 = 2 sabse chhote ∥ x ∥ ke saath.
Step 1: Full row rank ⇒ Case 2. A A ⊤ = 2 .
Step 2: A + = A ⊤ ( A A ⊤ ) − 1 = ( 1 1 ) ⋅ 2 1 = ( 1/2 1/2 ) .
Step 3: x = A + b = ( 1 1 ) . Kyun? Sabhi ( x 1 , x 2 ) jinka sum 2 hai, unme se ( 1 , 1 ) ki length minimum hai — yeh nullspace direction ( 1 , − 1 ) ke perpendicular hai.
Worked example Example 3 — rank-deficient via SVD
A = ( 2 0 0 0 ) . Singular values σ 1 = 2 , σ 2 = 0 .
Σ + = ( 1/2 0 0 0 ) (2 ko invert karo, 0 ko chhodho). U = V = I ke saath:
A + = ( 1/2 0 0 0 )
Zero kyun waise hi rehta hai? Doosri direction nullspace mein hai — invert karne ke liye koi information nahi.
Common mistake Common errors ko steel-man karna
Mistake 1: "A + = ( A ⊤ A ) − 1 A ⊤ hamesha."
Kyun sahi lagta hai: yeh formula har least-squares textbook mein hai. Catch: iske liye full column rank chahiye, warna A ⊤ A singular hai aur invert nahi ho sakta. Fix: rank-deficient matrices ke liye SVD formula use karo.
Mistake 2: "A + A = I hamesha."
Kyun sahi lagta hai: ek inverse cheezein undo karta hai. Sach: A + A = I sirf full column rank ke liye hota hai. Generally A + A ek projection hai, identity nahi (iska eigenvalue 0 bhi ho sakta hai).
Mistake 3: "Zero singular values ko bhi invert karo (ek tiny number use karo)."
Kyun tempting lagta hai: awkward zero se bachna ho. Danger: 1/ ϵ bahut bada ho jaata hai, noise enormously amplify ho jaati hai. Fix: truly zero (ya negligible) σ i ko 0 par hi map karo.
Recall Feynman: 12-saal ke bacche ko samjhao
Ek vending machine socho. Ek normal inverse ek aisi machine hai jo, snack dene par, exactly wahi coins wapas deti hai jo tune daale the. Lekin kuch machines tuti hoti hain — woh information squish kar deti hain toh coins perfectly wapas nahi mil sakte, ya bahut saari coin-combinations ek hi snack deti hain. Pseudoinverse ek sabse fair refund machine hai: agar perfect refund possible nahi, to woh answer deti hai jo sabse close ho (least error); agar kaafi refunds kaam karte hain, to woh sabse chhota, simplest wala deti hai (least wasted coins). Yeh kabhi crash nahi karti, chahe machine kaisi bhi ho.
Mnemonic Teeno formulas yaad rakho
"Tall Tops, Fat Bottoms, SVD always."
Tall (full column rank): ( A ⊤ A ) − 1 top-left par → A + = ( A ⊤ A ) − 1 A ⊤ .
Fat (full row rank): A ⊤ upar, inverse neeche → A + = A ⊤ ( A A ⊤ ) − 1 .
SVD (V Σ + U ⊤ ) kisi bhi case ko handle karta hai — U , V flip karo, nonzero σ invert karo.
Chaar Penrose conditions kya hain? A A + A = A ; A + A A + = A + ; ( A A + ) ⊤ = A A + ; ( A + A ) ⊤ = A + A .
Full column rank ke liye Pseudoinverse? A + = ( A ⊤ A ) − 1 A ⊤ (ek left inverse, A + A = I ).
Full row rank ke liye Pseudoinverse? A + = A ⊤ ( A A ⊤ ) − 1 (ek right inverse, A A + = I ).
A + ke liye universal SVD formula?A + = V Σ + U ⊤ , har nonzero σ i ko invert karo, zeros ko zeros rehne do.
A A + geometrically kya represent karta hai?A ke column space par orthogonal projection.
Least squares se A ⊤ A x = A ⊤ b kyun milta hai? ∇∥ A x − b ∥ 2 = 0 set karne par; residual column space ke perpendicular hota hai.
A + = A − 1 kab hota hai?Jab A square aur invertible ho.
Fat-matrix pseudoinverse kaun sa problem solve karta hai? A x = b ke infinitely many solutions mein se minimum-norm solution.
Zero singular values ko invert kyun nahi karte? 1/0 undefined hai aur tiny values noise ko blow up kar deti hain; un directions mein koi information nahi hoti.
Inverse A^-1 needs square invertible