Matrixpotenz

In der linearen Algebra bezeichnet die Matrixpotenz das Ergebnis einer wiederholten Matrixmultiplikation.

Die Potenz einer quadratischen Matrix A ∈ R n × n {\displaystyle A\in R^{n\times n}} über einem Halbring R {\displaystyle R} wird analog zur Potenz einer Zahl als wiederholte Multiplikation definiert. Ist A {\displaystyle A} eine quadratische Matrix, so ist

A 0 = E ; A 1 = A ; A 2 = A ⋅ A ; A 3 = A ⋅ A ⋅ A {\displaystyle A^{0}=E;\quad A^{1}=A;\quad A^{2}=A\cdot A;\quad A^{3}=A\cdot A\cdot A} usw.

Allgemein:

A n = A ⋅ A … A ⏟ n -mal {\displaystyle A^{n}=\underbrace {A\cdot A\dotsc A} _{n{\text{-mal}}}} .

Formal definiert man die Potenz rekursiv: Ist A {\displaystyle A} eine quadratische Matrix, so ist

• A 0 = E {\displaystyle A^{0}=E} die Einheitsmatrix und
• für alle k ∈ N = { 0 , 1 , 2 , … } {\displaystyle k\in \mathbb {N} =\{0,1,2,\dotsc \}} gilt A k + 1 = A k ⋅ A {\displaystyle A^{k+1}=A^{k}\cdot A} .

Es gelten die Potenzgesetze: Für alle n , m ∈ N {\displaystyle n,m\in \mathbb {N} } gilt

• A n + m = A n ⋅ A m {\displaystyle A^{n+m}=A^{n}\cdot A^{m}} ,
• A n ⋅ m = ( A n ) m {\displaystyle A^{n\cdot m}=\left(A^{n}\right)^{m}} .

Für invertierbare Matrizen sind auch Potenzen mit negativen ganzzahligen Exponenten definiert. Die Schreibweise A − 1 {\displaystyle A^{-1}} für die inverse Matrix kann auch als Matrixpotenz interpretiert werden. Für negative Exponenten − n {\displaystyle -n} , n ∈ N {\displaystyle n\in \mathbb {N} } , setzt man

A − n = ( A − 1 ) n {\displaystyle A^{-n}=\left(A^{-1}\right)^{n}} .

Matrixpotenzen mit nicht ganzzahligen Exponenten, beispielsweise die Quadratwurzel einer Matrix, können nur in Sonderfällen definiert werden.

In manchen Fällen kann die Matrixpotenz auf die Potenz von reellen Zahlen zurückgeführt werden. Lässt sich die Matrix A {\displaystyle A} diagonalisieren, existieren also eine reguläre Matrix T {\displaystyle T} und eine Diagonalmatrix D {\displaystyle D} mit A = T ⋅ D ⋅ T − 1 {\displaystyle A=T\cdot D\cdot T^{-1}} (d. h. A {\displaystyle A} ist ähnlich zu D {\displaystyle D} ), so gilt

A n = T ⋅ D n ⋅ T − 1   . {\displaystyle A^{n}=T\cdot D^{n}\cdot T^{-1}\ .}

Die Potenz einer Diagonalmatrix erhält man durch Potenzieren der Diagonalelemente. Sind die Diagonalelemente von D {\displaystyle D} (also die Eigenwerte von A {\displaystyle A} ) positiv, so bleiben obige Potenzgesetze auch für gebrochene Exponenten gültig.

Wenn sich eine Matrix nicht diagonalisieren lässt, so findet man eine sinnvolle Verallgemeinerung der Matrixpotenz über die binomische Reihe. Eine schnelle Berechnungsmethode für diese Verallgemeinerung erhält man über die Jordansche Normalform. Ist A = T ⋅ J ⋅ T − 1 {\displaystyle A=T\cdot J\cdot T^{-1}} eine Jordanzerlegung, dann gilt

A n = T ⋅ J n ⋅ T − 1 {\displaystyle A^{n}=T\cdot J^{n}\cdot T^{-1}}

Ist der Exponent eine ganze Zahl, so lässt sich die Matrixpotenz effizient mit binärer Exponentiation berechnen. Die Einschränkungen an den Zahlenbereich der Matrixelemente sind gering:

• Ist der Exponent nicht-negativ, so müssen die Matrixelemente in einem Ring liegen.
• Ist der Exponent negativ, so müssen die Matrixelemente in einem Körper liegen.

Ist der Zahlenbereich der Matrixelemente algebraisch abgeschlossen, kann man also darin beliebige algebraische Gleichungen lösen, so kann der Exponent auch rational sein und die Matrixpotenz kann über die Jordansche Normalform von A {\displaystyle A} auf Potenzen von skalaren Werten zurückgeführt werden, siehe oben.

Mittels der Matrixpotenz lassen sich Polynome auch für Matrizen definieren. Ein Beispiel dafür ist z. B. das Minimalpolynom. Genauso kann man auch Potenzreihen für Matrizen definieren, die wichtigsten Reihen sind dabei der Matrixlogarithmus, das Matrixexponential sowie die Neumann-Reihe.

Durch geeignete Wahl des zugrunde liegenden Halbrings R {\displaystyle R} lässt sich das Finden der kürzesten Pfade in einem Graphen auf die Berechnung einer Potenz der Adjazenzmatrix des Graphen zurückführen. Die Min-Plus-Matrixmultiplikation erhält man, indem man als Trägermenge von R {\displaystyle R} die erweiterten reellen Zahlen R ∗ = R + ∪ { ∞ } {\displaystyle \mathbb {R} ^{\ast }=\mathbb {R} ^{+}\cup \{\infty \}} wählt. Die Addition in R {\displaystyle R} entspricht dann der Minimumbildung in R ∗ {\displaystyle \mathbb {R} ^{\ast }} und die Multiplikation in R {\displaystyle R} der Addition in R ∗ {\displaystyle \mathbb {R} ^{\ast }} , wobei man x + ∞ = ∞ + x = ∞ {\displaystyle x+\infty =\infty +x=\infty } setzt. Die absorbierende Null in R {\displaystyle R} ist dann ∞ {\displaystyle \infty } , während das Einselement in R {\displaystyle R} durch die Zahl 0 {\displaystyle 0} dargestellt wird. Ist nun K ∈ R n × n {\displaystyle K\in R^{n\times n}} die Kostenmatrix eines Graphen mit n {\displaystyle n} Knoten, dann ist D = ∑ k = 0 n K k {\displaystyle \textstyle D=\sum _{k=0}^{n}K^{k}} die zugehörige Entfernungsmatrix mit den Längen der kürzesten Pfade zwischen allen Knoten des Graphen. Da die Addition in R {\displaystyle R} idempotent ist, ist D = ( 1 + K ) n {\displaystyle D=(1+K)^{n}} .

• In der theoretischen Ökonomie bzw. Biologie werden Matrixpotenzen zur Analyse langfristiger Populationsentwicklungen eingesetzt, beispielsweise unter Nutzung einer Leslie-Matrix.^[1]
• Des Weiteren gibt es Anwendungen bei der Stereobasisverbreiterung.

• Peter Knabner, Wolf Barth: Lineare Algebra. Grundlagen und Anwendungen (= Springer-Lehrbuch). Springer Spektrum, Berlin u. a. 2013, ISBN 978-3-642-32185-6. 
• Gilbert Strang: Lineare Algebra. Eine Einführung für Studienanfänger (= Springer-Lehrbuch). Springer, Berlin 2003, ISBN 3-540-43949-8 (englisch: Introduction to linear algebra. Übersetzt von Michael Dellnitz). 

1. ↑ Populationsentwicklung. (PDF; 72 kB) Archiviert vom Original am 4. März 2018; abgerufen am 27. Februar 2022. , Archivlink