logo

Exponential and Logarithm of Matrices 📂Matrix Algebra

Exponential and Logarithm of Matrices

Definition

Aiming to generalize the exponential function exp\exp and the logarithmic function log\log to matrices.

Matrix Exponential

The generalization of the exponential function to matrices exp:Cn×nCn×n\exp : \mathbb{C}^{n \times n} \to \mathbb{C}^{n \times n} is defined for a matrix ACn×nA \in \mathbb{C}^{n \times n} as follows: expA:=k=0Akk! \exp A := \sum_{k=0}^{\infty} \frac{A^{k}}{k!} The expression expA\exp A can also be simply written as eAe^{A}, and it is called the matrix exponential of AA.

Matrix Logarithm

The generalization of the logarithmic function to matrices log:Cn×nCn×n\log : \mathbb{C}^{n \times n} \to \mathbb{C}^{n \times n} is defined such that if there exists a matrix AA that satisfies expA=B\exp A = B with BB, then we denote AA as logB\log B, and it is called the matrix logarithm of BB.

Theorem 1

If two matrices A,BCn×nA, B \in \mathbb{C}^{n \times n} satisfy AB=BAAB = BA, then the following holds: logAB=logA+logB2πiU(logA+logB) \log AB = \log A + \log B - 2 \pi i \mathcal{U} \left( \log A + \log B \right) In particular, for the eigenvalues μk\mu_{k} of AA and the eigenvalues νk\nu_{k} of BB, we obtain the following corollary: logAB=logA+logB    k=1,,n:argμk+argνk(π,π] \log AB = \log A + \log B \iff \forall k = 1 , \cdots , n : \arg \mu_{k} + \arg \nu_{k} \in ( - \pi , \pi ] Here arg\arg is the argument of a complex number, and U\mathcal{U} is the matrix unwinding function defined as follows: U(A):=AlogeA2πi,ACn×n \mathcal{U} (A) := {\frac{ A - \log e^{A} }{ 2 \pi i }} \qquad , A \in \mathbb{C}^{n \times n}

Explanation

The matrix exponential and logarithm formally extend the exponential and logarithm defined in complex numbers. As seen, the matrix exponential is defined using the power series of the matrix.

In the Space of Hermitian and Positive Definite Matrices

Space of Hermitian and Positive Definite Matrices:

  • The Hermitian matrix space: The set of all Hermitian matrices of size n×nn \times n is represented as follows. For a scalar field R\mathbb{R}, Hn\mathbb{H}_{n} is a vector space. Hn:={ACn×n:A=A} \mathbb{H}_{n} := \left\{ A \in \mathbb{C}^{n \times n} : A = A^{\ast} \right\}
  • The set of positive definite matrices: The set of all positive definite matrices of size n×nn \times n is denoted as Pn\mathbb{P}_{n}. PnHn\mathbb{P}_{n} \subset \mathbb{H}_{n} is a convex cone of Hn\mathbb{H}_{n}.

Instead of the space of all matrices, considering the Hermitian matrix space Hn\mathbb{H}_{n} as the domain and the convex cone of positive definite matrices Pn\mathbb{P}_{n} as the codomain, exp:HnPn\exp : \mathbb{H}_{n} \to \mathbb{P}_{n} is a bijection, with its inverse being log:PnHn\log : \mathbb{P}_{n} \to \mathbb{H}_{n}. Particularly, these are diffeomorphisms.

HP.png

See Also


  1. Aprahamian, M., & Higham, N. J. (2014). The matrix unwinding function, with an application to computing the matrix exponential. SIAM Journal on Matrix Analysis and Applications, 35(1), 88-109. https://doi.org/10.1137/130920137 Lemma 3.12 ↩︎