Logarithm of 2x2 Matrix: A General Guide

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
2 replies · 9K views
jhendren
Messages
34
Reaction score
0
How do you take the log of a 2x2 matrix in general where all entries are non-zero
 
Physics news on Phys.org
jhendren said:
How do you take the log of a 2x2 matrix in general where all entries are non-zero

Hey jhendren and welcome to the forums.

Have you tried either (a) diagonalizing the matrix or (b) using an operator expansion based on the Taylor series of the logarithm?

The first one is based on the eigen-decomposition and the second one is based on the operator algebra results for functions of a linear operator.
 
If y= ln(x) then [itex]x= e^y[/itex]. So to define the logarithm is to define the exponential and vice-versa. And it is easier to work with the exponential. Its Taylor series is [itex]\sum_{n=0}^\infty x^n/n![/itex].

One can show that if matrix x is "diagonalizable", that is, if there exist a matrix P such that [itex]x= PDP^{-1}[/itex] where D is a diagonal matrix, then that is [itex]\sum_{n= 0}^\infty (PDP^{-1})^n/n!= P\left(\sum_{n=0}^\infty D^n\right)P^{-1}[/itex]. And [itex]D^n[/itex] is just the diagonal matrix with the nth powers of the diagonal elements of D on it diagonal. That reduces to [itex]e^x= Pe^DP^{-1}[/itex] where, now, [itex]e^D[/itex] is the diagonal matrix having the exponentials of the diagonal elements of D on its diagonal.

If x is not diagonalizable, it can still be written in "Jordan Normal Form" but the exponential of that is trickier.

If, for example,
[tex]D= \begin{bmatrix}a & 0 \\ 0 & b\end{bmatrix}[/tex]
then
[tex]e^D= \begin{bmatrix}e^a & 0 \\ 0 & e^b\end{bmatrix}[/tex].

If A is the "Jordan Normal Form", written as
[tex]A= \begin{bmatrix}a & 1 \\ 0 & a\end{bmatrix}[/tex]
then it is easy to show that
[tex]A^n= \begin{bmatrix}a^n & na^{n-1} & 0 & a^n\end{bmatrix}[/tex]
so that
[tex]e^A= \begin{bmatrix}\sum a^n/n! & sum na^{n-1}/n!\\ 0 \sum a^n/n!\end{bmatrix}= \begin{bmatrix}\sum a^n/n! & sum a^{n-1}{n!} \\ 0 & a^n/n! \end{bmatrix}[/tex]
[tex]e^A= \begin{bmatrix}e^a & e^a \\ 0 & e^a\end{bmatrix}[/tex]
 
Last edited by a moderator: