logo

Pseudo Inverse Matrix 📂Matrix Algebra

Pseudo Inverse Matrix

Overview

The Pseudoinverse Matrix generalizes the concept of an inverse matrix to matrices ARm×nA \in \mathbb{R}^{m \times n} that are not square matrices because their rows and columns are not of the same size. It refers to a matrix that acts as an ’effective’ inverse matrix for them. If a matrix transformation TA:N(A)C(A)T_{A} : \mathcal{N} (A) \to \mathcal{C} (A) satisfies
TAx=Ax T_{A} \mathbf{x} = A \mathbf{x} for all xN(A)\mathbf{x} \in \mathcal{N} (A)^{\perp}, then TAT_{A} becomes bijective. This can be seen as narrowing down the codomain of TAT_{A} to forcefully make it surjective, and an inverse transformation TA1:C(A)N(A)T_{A}^{-1} : \mathcal{C} (A) \to \mathcal{N} (A) exists through which the pseudoinverse can be defined.

Definition1

Given a matrix ARm×nA \in \mathbb{R}^{m \times n} and matrix transformation TA:N(A)C(A)T_{A} : \mathcal{N} (A) \to \mathcal{C} (A), for every vector y=y1+y2Rm\mathbf{y} = \mathbf{y}_{1} + \mathbf{y}_{2} \in \mathbb{R}^{m} if the matrix AA^{\dagger} satisfies the following, AA^{\dagger} is called the Moore-Penrose Pseudoinverse of AA. Ay=TA1y1 A^{\dagger} \mathbf{y} = T_{A}^{-1} \mathbf{y}_{1}


  • Where y1C(A)\mathbf{y}_{1} \in \mathcal{C}(A) and y2C(A)\mathbf{y}_{2} \in \mathcal{C}(A)^{\perp}.
  • \dagger is represented by a symbol that looks like a dagger pointed downwards, and indeed, it is read as Dagger. In mathematical physics, it might also represent a conjugate transpose matrix, but the context should help prevent confusion.

Theorem

The pseudoinverse of a matrix ARm×nA \in \mathbb{R}^{m \times n} is calculated as follows. A=limδ0(ATA+δ2I)1AT=limδ0AT(ATA+δ2I)1 \begin{align*} A^{\dagger} =& \lim_{\delta \to 0} \left( A^{T} A + \delta^{2} I \right)^{-1} A^{T} \\ =& \lim_{\delta \to 0} A^{T} \left( A^{T} A + \delta^{2} I \right)^{-1} \end{align*}


  1. 김상동, 김필수, 신병춘, 이용훈. (2012). 수치행렬해석: p78. ↩︎