Logarytm macierzy

Logarytm macierzy – inna macierz taka, że wykładnicza macierz drugiej macierzy jest równa oryginalnej macierzy. A zatem jest to uogólnienie logarytmu skalarnego, a w pewnym sensie także odwrotna funkcja macierzy wykładniczej. Nie wszystkie macierze mają logarytm i te macierze, które posiadają logarytm, mogą mieć więcej niż jeden logarytm. Badanie logarytmów macierzy prowadzi do teorii Liego. Od kiedy macierz ma logarytm, jest w grupie Liego i logarytm jest odpowiednim elementem algebry Liego.

Definicja

Wykładnicza forma macierzy $A$ jest określona:

$e^{𝐀}\equiv \sum _{n=0}^{\mathrm{\infty }}\frac{{𝐀}^n}{n!}.$

Biorąc pod uwagę macierz $B,$ macierz $A$ jest logarytmem macierzy $B,$ jeżeli $e^{𝐀}=B.$ Często logarytmy matrycy nie są unikatowe, tak jak logarytmy liczb zespolonych, co wyjaśniono poniżej.

Przykład: Logarytm obrotów w płaszczyźnie

Obroty w płaszczyźnie to prosty przykład. Obrót kąta $\alpha$ wokół macierzy jest reprezentowany przez macierz $2\times 2$

$A=\left(\begin{array}{c}\cos (\alpha )-\sin (\alpha )\\ \sin (\alpha )\cos (\alpha )\end{array}\right).$

Dla każdej liczby naturalnej $n,$ macierz

$B_n=(\alpha +2\pi n)\left(\begin{array}{cc}0& -1\\ 1& 0\end{array}\right),$

to logarytm macierzy $A.$ Zatem macierz $A$ ma nieskończenie wiele logarytmów. To wynika z tego, że kąt obrotu jest ustalony tylko do wielokrotności $2\pi .$

W języku teorii Liego obroty macierzy $A$ są elementami grupy Liego SO(2). Odpowiadające im logarytmy B są elementami algebry Liego SO(2), która składa się z wszystkich macierzy skośno-symetrycznych.

Macierz

$\left(\begin{array}{cc}0& 1\\ -1& 0\end{array}\right)$

jest generatorem algebry Liego SO(2).

Czy istnieje?

Pytanie, czy macierz ma logarytm ma najprostszą odpowiedź, gdy rozpatrywana jest w złożonej konfiguracji. Macierz ma logarytm wtedy i tylko wtedy, gdy jest odwracalna. Logarytm nie jest wyjątkowy, ale jeśli matryca nie ma negatywnych rzeczywistych wartości własnych, to ma wyjątkowy logarytm, którego wartości własne zawierają się w klamrze: $\{z\in C.|-\pi <\mathrm{Im}z<\pi \}.$ Ten logarytm jest znany jako logarytm główny.

Szablon:Fakt.

Własności

Jeśli $A$ i $B$ są dodatnio określone oraz $A$ i $B$ są przemienne, czyli $AB=BA,$ to

$AB=e^{\ln (A)+\ln (B)}.$

Dla każdej odwracalnej macierzy, $A^{-1}=e^{-\ln (A)}.$

Obliczanie logarytmu macierzy diagonalnej

Sposób znalezienia $\ln (A)$ diagonalnej macierzy $A$ jest następujący:

• Znajdź macierz $V$ wektorów własnych $A$ (każda kolumna $V$ jest wektorem własnym $A$).
• Znajdź odwrotność $V^{-1}$ z $V.$
• Niech $A^{\prime }=V^{-1}AV.$

• Wtedy $A^{\prime }$ będzie macierzą diagonalną, której elementy przekątnej są wartościami własnymi macierzy $A.$
• Zamień każdy element przekątnej $A$ na jego logarytmu naturalny, w celu uzyskania $\ln A^{\prime }.$
• Wtedy $\ln A=V(\ln A^{\prime })V^{-1}.$

To, że logarytm A może być macierzą zespoloną, nawet jeśli A jest rzeczywiste, wynika z faktu, że macierz z pozytywnymi i rzeczywistymi wartościami może mieć negatywne lub nawet złożone wartości własne (dotyczy to przykładowo macierzy rotacji). Brak unikatowości logarytmu macierzy wynika z braku jednoznaczności logarytmu liczby zespolonej.

Obliczanie logarytmu niediagonalnej macierzy

Algorytm przedstawiony powyżej nie działa na niediagonalnych macierzach, takich jak

$\left[\begin{array}{cc}1& 1\\ 0& 1\end{array}\right].$

Dla takich matryc trzeba znaleźć ich rozkład Jordana i, zamiast obliczania logarytmów elementów przekątnej jak wyżej, trzeba obliczyć logarytm bloków Jordana.

Ten ostatni jest tworzony, dzięki zauważeniu, że można napisać blok Jordan jako

$B=\left(\begin{array}{ccccc}\lambda 1& 0& 0& \cdots & 0\\ 0& \lambda 1& 0& \cdots & 0\\ 0& 0& \lambda 1\cdots & 0\\ ⋮& ⋮& ⋮& \ddots & \ddots & ⋮\\ 0& 0& 0& 0& \lambda 1\\ 0& 0& 0& 0& 0& \lambda \end{array}\right)=\lambda \left(\begin{array}{cccccc}1& {\lambda }^{-1}& 0& 0& \cdots & 0\\ 0& 1& {\lambda }^{-1}& 0& \cdots & 0\\ 0& 0& 1& {\lambda }^{-1}\cdots & 0\\ ⋮& ⋮& ⋮& \ddots & \ddots & ⋮\\ 0& 0& 0& 0& 1& {\lambda }^{-1}\\ 0& 0& 0& 0& 0& 1\end{array}\right)=\lambda (I+K)$

gdzie $K$ jest macierzą z zerami na i poniżej głównej przekątnej (liczba $\lambda$ jest niezerowa z powodu założenia, że macierz którego logarytm próbuje się obliczyć jest odwracalna).

Następnie, dzięki rozwinięciu Mercatora

$\ln (1+x)=x-\frac{x^2}{2}+\frac{x^3}{3}-\frac{x^4}{4}+\cdot$

staje się

$\ln B=\ln (\lambda (I+K))=\ln (\lambda I)+\ln (I+K)=\ln (\lambda )I+K-K^2/2+K^3/3-K^4/4+\dots$

To rozwinięcie na ogół nie jest zbieżne dla każdej macierzy $K,$ ponieważ nie jest zbieżne dla każdej liczby rzeczywistej z wartością bezwzględną większą od jedności, to jednak $K$ jest macierzą nilpotenta, tak więc cykl rzeczywiście ma skończoną liczbę wyrazów ($K^m$ jest równe zero, jeśli jest wymiarem $k$).

Stosując to odkrycie, otrzymujemy

$\ln \left[\begin{array}{cc}1& 1\\ 0& 1\end{array}\right]=\left[\begin{array}{cc}0& 1\\ 0& 0\end{array}\right].$

Funkcjonalna analiza perspektywiczna

Szablon:Fakt.

Korzystanie z narzędzi holomorficznego rachunku funkcyjnego, biorąc pod uwagę funkcję holomorficzną $f(z),$ określoną na zbiorze otwartym w płaszczyźnie zespolonej i ograniczony operator liniowy $T$ można obliczyć $f(t)$ tak długo, jak $f(z)$ jest określona na Widmie $T.$

Funkcja $F(z)=\ln (z)$ może być zdefiniowana w dowolnym prostym połączonym zbiorze otwartym w płaszczyźnie zespolonej nie zawierającej początku i jest holomorficzny w takiej domenie. Oznacza to, że można określić $\ln (T),$ pod warunkiem, że widmo $T$ nie zawiera początku i istnieje ścieżka począwszy od środka do nieskończoności nie przekraczająca spektrum $T$ (na przykład, w przypadku widma z $T,$ którym jest koło z początkiem wewnątrz koła, jest możliwe określenie $\ln (T)$).

Szablon:Fakt.

Bibliografia

• Gantmacher, Felix R.(1959) Teoria matryc 1, Nowy York, Chelsea, s. 239–241.
• Culver, Walter J. (1966), O istnieniu i unikatowości rzeczywistych logarytmów macierzy, „Proceedings of the American Mathematical Society” 17 (5): s. 1146–1151, doi:10.1090/S0002-9939-1966-0202740-6, ISSN 0002-9939.
• Higham, Nicholas (2008), Functions of Matrices. Theory and Computation, SIAM, Szablon:ISBN.
• Engo, Kenth (June 2001), On the BCH-formula in so(3), „BIT Numerical Mathematics” 41 (3): 629-632, doi: 10.1023/A:1021979515229, ISSN 0006-3835.