Drazin inverse explained
In mathematics, the Drazin inverse, named after Michael P. Drazin, is a kind of generalized inverse of a matrix.
Let A be a square matrix. The index of A is the least nonnegative integer k such that rank(Ak+1) = rank(Ak). The Drazin inverse of A is the unique matrix AD that satisfies
Ak+1AD=Ak, ADAAD=AD, AAD=ADA.
It's not a generalized inverse in the classical sense, since
in general.
, then
.
- If A is a block diagonal matrix
A=\begin{bmatrix}B&0\\
0&N
\end{bmatrix}
where
is invertible with inverse
and
is a
nilpotent matrix, then
AD=\begin{bmatrix}B-1&0\\
0&0
\end{bmatrix}
- Drazin inversion is invariant under conjugation. If
is the Drazin inverse of
, then
is the Drazin inverse of
.
- The Drazin inverse of a matrix of index 0 or 1 is called the group inverse or -inverse and denoted A#. The group inverse can be defined, equivalently, by the properties AA#A = A, A#AA# = A#, and AA# = A#A.
- A projection matrix P, defined as a matrix such that P2 = P, has index 1 (or 0) and has Drazin inverse PD = P.
- If A is a nilpotent matrix (for example a shift matrix), then
The hyper-power sequence is
Ai+1:=Ai+Ai\left(I-AAi\right);
for convergence notice that
For
or any regular
with
chosen such that
\left\|A0-A0AA0\right\|<\left\|A0\right\|
the sequence tends to its Drazin inverse,
Drazin inverses in categories
A study of Drazin inverses via category-theoretic techniques, and a notion of Drazin inverse for a morphism of a category, has been recently initiated by Cockett, Pacaud Lemay and Srinivasan. This notion is a generalization of the linear algebraic one, as there is a suitably defined category
having morphisms matrices
with complex entries; a Drazin inverse for the matrix M amounts to a Drazin inverse for the corresponding morphism in
.
Jordan normal form and Jordan-Chevalley decomposition
As the definition of the Drazin inverse is invariant under matrix conjugations, writing
, where J is in Jordan normal form, implies that
. The Drazin inverse is then the operation that maps invertible Jordan blocks to their inverses, and nilpotent Jordan blocks to zero.
where
is semisimple and
is nilpotent and both operators commute. The two terms can be block diagonalized with blocks corresponding to the kernel and cokernel of
. The Drazin inverse in the same basis is then defined to be zero on the kernel of
, and equal to the inverse of
on the cokernel of
.
See also
References
- Drazin. M. P.. Pseudo-inverses in associative rings and semigroups. The American Mathematical Monthly. 65. 1958. 506–514. 2308576. 10.2307/2308576. 7.
- 10.1016/S0096-3003(03)00786-0 . Generalized inverse A(2)T,S and a rank equation. 2004. Zheng. Bing. Bapat. R.B. Applied Mathematics and Computation. 155. 2. 407.
- Drazin Inverses in Categories . 2024 . Cockett . Robin . Pacaud Lemay . Jean-Simon . Srinivasan . Priyaa Varshinee . 2402.18226 . math.CT.
External links