Skip to content

Generalized Inverse


Definition

Let AA be a real m×nm\times n matrix. A generalized inverse of AA is any n×mn\times m matrix AIA^{I} satisfying

AAIA=A.AA^{I}A=A.

The matrix products are defined because AIA^{I} has the shape of A⊤A^{\top}. No condition on mm, nn or the rank of AA is imposed.


Interpretation

An inverse matrix exists only for square nonsingular matrices. The condition AAIA=AAA^{I}A=A keeps only what is needed to solve equations: whenever b=Ax\boldsymbol{b}=A\boldsymbol{x} is reachable, AIbA^{I}\boldsymbol{b} is a vector that AA maps back to b\boldsymbol{b}, because A(AIAx)=AxA(A^{I}A\boldsymbol{x})=A\boldsymbol{x}.

Because one equation is a weak requirement, a generalized inverse is in general far from unique, and the notation AIA^{I} means "some generalized inverse". Adding the three further Penrose conditions singles out one of them, the Moore-Penrose inverse A+A^+. A common mistake is to assume properties of A+A^+, such as symmetry of AAIAA^{I}, for an arbitrary AIA^{I}.

One-sided inverses are special cases: if LA=InLA=I_n (left inverse, full column rank) then ALA=AALA=A; likewise for a right inverse.


Key Results

  • Existence: every matrix has a generalized inverse. If A=P(IrOOO)QA=P\begin{pmatrix}I_r&O\\O&O\end{pmatrix}Q with P,QP,Q invertible, then the generalized inverses are exactly Q−1(IrXYZ)P−1Q^{-1}\begin{pmatrix}I_r&X\\Y&Z\end{pmatrix}P^{-1} with X,Y,ZX,Y,Z arbitrary.
  • Uniqueness: AIA^{I} is unique if and only if AA is square and invertible, and then AI=A−1A^{I}=A^{-1}.
  • AAIAA^{I} and AIAA^{I}A are idempotent, rank⁡(AAI)=rank⁡(AIA)=rank⁡A\operatorname{rank}(AA^{I})=\operatorname{rank}(A^{I}A)=\operatorname{rank}A, and rank⁡AI≥rank⁡A\operatorname{rank}A^{I}\ge\operatorname{rank}A.
  • Consistency: Ax=bA\boldsymbol{x}=\boldsymbol{b} has a solution if and only if AAIb=bAA^{I}\boldsymbol{b}=\boldsymbol{b}.
  • General solution: if the system is consistent, its solutions are exactly x=AIb+(In−AIA)z\boldsymbol{x}=A^{I}\boldsymbol{b}+(I_n-A^{I}A)\boldsymbol{z}, z∈Rn\boldsymbol{z}\in\mathbb{R}^n, an affine subspace.

Examples

Let A=(1111)A=\begin{pmatrix}1&1\\1&1\end{pmatrix} (rank 11, singular).

Two different generalized inverses. For G=(1000)G=\begin{pmatrix}1&0\\0&0\end{pmatrix}: AG=(1010)AG=\begin{pmatrix}1&0\\1&0\end{pmatrix} and AGA=(1111)=AAGA=\begin{pmatrix}1&1\\1&1\end{pmatrix}=A. The matrix G′=(0001)G'=\begin{pmatrix}0&0\\0&1\end{pmatrix} works too, and so does 14A\tfrac{1}{4}A, since A3=4AA^3=4A. Note that AGAG is idempotent but not symmetric.

Solving a system. For b=(2,2)⊤\boldsymbol{b}=(2,2)^{\top}: AGb=(2,2)⊤=bAG\boldsymbol{b}=(2,2)^{\top}=\boldsymbol{b}, so the system is consistent. With I−GA=(0−101)I-GA=\begin{pmatrix}0&-1\\0&1\end{pmatrix},

x=Gb+(I−GA)z=(2,0)⊤+z2(−1,1)⊤,\boldsymbol{x}=G\boldsymbol{b}+(I-GA)\boldsymbol{z} =(2,0)^{\top}+z_2(-1,1)^{\top},

which is the line x1+x2=2x_1+x_2=2. For b=(1,0)⊤\boldsymbol{b}=(1,0)^{\top}: AGb=(1,1)⊤≠bAG\boldsymbol{b}=(1,1)^{\top}\neq\boldsymbol{b}, so there is no solution.

The solution line x1 + x2 = 2 of the example system, containing the two particular solutions Gb = (2,0) and G'b = (0,2) and the direction (I - GA)z along the line.

The figure draws the consistent example with b=(2,2)⊤\boldsymbol{b}=(2,2)^{\top}. The two generalized inverses GG and G′G' give different particular solutions, Gb=(2,0)⊤G\boldsymbol{b}=(2,0)^{\top} and G′b=(0,2)⊤G'\boldsymbol{b}=(0,2)^{\top}, but both lie on the same line x1+x2=2x_1+x_2=2. Adding (I−GA)z=z2(−1,1)⊤(I-GA)\boldsymbol{z}=z_2(-1,1)^{\top} moves along that line, so the choice of AIA^{I} changes the starting point, never the solution set.


Prerequisites



Sources

Aliases: g-inverse, generalised inverse, {1}-inverse, inner inverse, 1-inverse Keywords: A^I, A^-, A A^I A = A, consistency condition, x = A^I b + (I - A^I A) z, conditional inverse Symbols: element of, real numbers

Last updated: 22 October 2026 19:04