Masteries · Jason GlowneyArithmoí Book 1

Masteries

1.1.13 Matrix Powers and Logarithms

The ex­po­nen­tial of a matrix MM, is de­fined as the fol­low­ing in­fi­nite series

eM=∑n=0∞Mnn!=N,\begin{array}{r} e^{M} = \sum_{n = 0}^{\infty}\frac{M^{n}}{n!} = N, \end{array}
(1.1.64)

which con­verges to the re­sult­ing matrix NN. We can define the matrix M, as the matrix log­a­rithm of N, since eM=N.e^{M} = N.In terms of the ex­is­tence of the matrix log­a­rithm, com­plex ma­tri­ces are straight­for­ward, where the log­a­rithm does exist if the matrix is in­vert­ible. Log­a­rithms may not nec­es­sar­i­ly be unique, but if there are no eigen­val­ues ∈ℝ−\in \mathbb{R}^{-}, a unique matrix log­a­rithm, termed the prin­ci­pal log­a­rithm, lies in region z∈ℂ|−π<z<π|z\in \mathbb{C}| -\pi < z < \pi|. For real ma­tri­ces, a matrix log­a­rithm exists if that matrix is in­vert­ible, and for all neg­a­tive eigen­val­ues −λ- \lambda, the Jordan blocks are in even num­bers.

In Lie theory, we have the ex­po­nen­tial map­ping from the Lie al­ge­bra 𝔤\mathfrak{g}, to its Lie group e𝔤→Ge^{\mathfrak{g}} \rightarrow G, and the cor­re­spond­ing log­a­rithm map­ping from the Lie group GG, to the Lie al­ge­bra ln(G)→𝔤.{\ ln}(G)\rightarrow \mathfrak{g}.In this book, the matrix groups that we con­sid­er are di­ag­o­nal­iz­able, and we can take ad­van­tage of this to com­pute matrix ex­po­nen­tials and matrix log­a­rithms. We can take an ini­tial di­ag­o­nal­iz­able Lie al­ge­bra matrix 𝔤\mathfrak{g}, and ex­po­nen­ti­ate it

𝔤=V−1𝔡V,\begin{array}{r} \mathfrak{g} =V^{- 1}\mathfrak{d}V, \end{array}
(1.1.65)

for VV, a matrix formed by eigen­vec­tor columns of 𝔤\mathfrak{g}, and 𝔡\mathfrak{d}, the di­ag­o­nal­ized matrix of its eigen­val­ues. We can now com­pute both the ex­po­nen­tial and the re­sult­ing log­a­rithm, as

G=Ve𝔡V−1\begin{array}{r} G = Ve^{\mathfrak{d}}V^{- 1} \end{array}
(1.1.66)

and

𝔤=ln⁡(G)=V∙ln(e𝔡)∙V−1.\begin{array}{r} \mathfrak{g} =\ln(G) = V \bullet ln\left( e^{\mathfrak{d}} \right) \bullet V^{- 1}. \end{array}
(1.1.67)

With the matrix ex­po­nen­tial and the matrix log­a­rithm in mind, we now have the tools to con­sid­er matrix powers. For a com­plex matrix M, raised to a com­plex power, we can com­pute com­plex powers of MM, as

Ma+ib=eln⁡(M)(a+ib) .\begin{array}{r} M^{a + ib} = e^{\ln(M)(a + ib)\ }. \end{array}
(1.1.68)

Taking the di­ag­o­nal­iz­able matrix MM

M=(0ii0),\begin{array}{r} M = \begin{pmatrix} 0 & i \\ i & 0 \end{pmatrix}, \end{array}
(1.1.69)

we find that

V=(121212−12),   V−1=1det⁡(V)adj(V)=(121212−12),\begin{array}{r} V = \begin{pmatrix} \frac{1}{\sqrt{2}} & \frac{1}{\sqrt{2}} \\ \frac{1}{\sqrt{2}} & - \frac{1}{\sqrt{2}} \end{pmatrix},\ \ \ V^{- 1} = \frac{1}{\det(V)}adj(V) = \begin{pmatrix} \frac{1}{\sqrt{2}} & \frac{1}{\sqrt{2}} \\ \frac{1}{\sqrt{2}} & - \frac{1}{\sqrt{2}} \end{pmatrix}, \end{array}
(1.1.70)

for adj(V)adj(V), rep­re­sent­ing the ad­ju­gate of the VV-matrix, where Vadj(V)Vadj(V) can be seen to equal the prod­uct det⁡(V)I\det(V)I. We then have

D=V−1MV=(i00−i),\begin{array}{r} D = V^{- 1}MV = \begin{pmatrix} i & 0 \\ 0 & - i \end{pmatrix}, \end{array}
(1.1.71)

for DD, the di­ag­o­nal matrix of eigen­val­ues from MM, with which, we can find ln⁡(M)\ln(M)

ln⁡(M) =V(ln⁡D)V−1=(0iπ2iπ20).\begin{array}{r} \ln(M)\ = V(\ln{D)}V^{- 1} = \begin{pmatrix} 0 & \frac{i\pi}{2} \\ \frac{i\pi}{2} & 0 \end{pmatrix}. \end{array}
(1.1.71)

With the result for ln⁡(M)\ln(M), we have what is re­quired to cal­cu­late the matrix MM, raised to a com­plex power Ma+ibM^{a + ib}, by the fol­low­ing

Ma+ib=eln⁡(M)∙(a+ib)→ea∙ln⁡(M)+ib∙ln⁡(M),\begin{array}{r} M^{a + ib} = e^{\ln(M) \bullet (a + ib)} \rightarrow e^{{a \bullet \ln}{(M) + ib \bullet \ln(M)}}, \end{array}
(1.1.72)

find­ing

Ma+ib=(0ii0)a+ib=ea∙(0iπ2iπ20)+ib∙(0iπ2iπ20)=ea(0iπ2iπ20)e(0−bπ2−bπ20).\begin{array}{r} M^{a + ib} = \begin{pmatrix} 0 & i \\ i & 0 \end{pmatrix}^{a + ib}{= e}^{a \bullet \begin{pmatrix} 0 & \frac{i\pi}{2} \\ \frac{i\pi}{2} & 0 \end{pmatrix} + ib \bullet \begin{pmatrix} 0 & \frac{i\pi}{2} \\ \frac{i\pi}{2} & 0 \end{pmatrix}} = e^{a\begin{pmatrix} 0 & \frac{i\pi}{2} \\ \frac{i\pi}{2} & 0 \end{pmatrix}}e^{\begin{pmatrix} 0 & - \frac{b\pi}{2} \\ - \frac{b\pi}{2} & 0 \end{pmatrix}}. \end{array}
(1.1.73)

A more straight­for­ward ex­am­ple would be to cal­cu­late

eMa+ib=eln⁡(eM)∙(a+ib)=eM∙(a+ib)=e(0−b+ia−b+ia0)=(cosh⁡(−b+ai)sinh⁡(−b+ai)sinh⁡(−b+ai)cosh⁡(−b+ai)).\begin{array}{r} {e^{M}}^{a + ib} = e^{\ln\left( e^{M} \right) \bullet (a + ib)} = e^{M \bullet (a + ib)} = e^{\begin{pmatrix} 0 & - b + ia \\ - b + ia & 0 \end{pmatrix}} = \begin{pmatrix} \cosh( - b + ai) & \sinh( - b + ai) \\ \sinh( - b + ai) & \cosh( - b + ai) \end{pmatrix}. \end{array}
(1.1.74)

The math­e­mat­ics in­volved in rais­ing a number or matrix to imag­i­nary/com­plex powers is an im­por­tant con­cept to know. We use it when con­sid­er­ing the angle spaces later in the book, which are not bound by the con­straints of geo­met­ric spaces as they are nor­mal­ly treat­ed.