Tridiagonalmatrix

In der linearen Algebra ist eine Tridiagonalmatrix (auch Dreibandmatrix) eine quadratische Matrix, die nur in der Hauptdiagonalen und in den beiden ersten Nebendiagonalen Einträge ungleich null enthält. Tridiagonalmatrizen treten in der Numerik recht häufig auf, zum Beispiel bei der Berechnung von kubischen Splines, bei der Diskretisierung der zweiten Ableitung auf eindimensionalen Gebieten (insbesondere bei Sturm-Liouville-Problemen), bei der Berechnung von orthogonalen Polynomen und Funktionensystemen (etwa bei der Berechnung von Besselfunktionen) und bei Krylow-Unterraum-Verfahren basierend auf Dreitermrekursionen.

Eine Matrix T ∈ C n × n {\displaystyle T\in \mathbb {C} ^{n\times n}} heißt tridiagonal, wenn sie die folgende Form hat:

T = ( t 1 , 1 t 1 , 2 0 … 0 t 2 , 1 t 2 , 2 t 2 , 3 ⋱ ⋮ 0 t 3 , 2 ⋱ ⋱ 0 ⋮ ⋱ ⋱ ⋱ t n − 1 , n 0 … 0 t n , n − 1 t n , n ) {\displaystyle T={\begin{pmatrix}t_{1,1}&t_{1,2}&0&\dots &0\\t_{2,1}&t_{2,2}&t_{2,3}&\ddots &\vdots \\0&t_{3,2}&\ddots &\ddots &0\\\vdots &\ddots &\ddots &\ddots &t_{n-1,n}\\0&\dots &0&t_{n,n-1}&t_{n,n}\end{pmatrix}}}

Es gilt also t i j = 0 {\displaystyle t_{ij}=0} für alle | i − j | > 1 {\displaystyle |i-j|>1} . Eine Tridiagonalmatrix heißt unreduziert oder irreduzibel, wenn die Elemente in den Nebendiagonalen alle ungleich Null sind, das heißt t i j ≠ 0 {\displaystyle t_{ij}\not =0} für alle | i − j | = 1 {\displaystyle |i-j|=1} gilt. Sind die Haupt- und Nebendiagonaleinträge konstant, gilt also t 1 , 1 = … = t n , n {\displaystyle t_{1,1}=\ldots =t_{n,n}} , t 1 , 2 = … = t n − 1 , n {\displaystyle t_{1,2}=\ldots =t_{n-1,n}} und t 2 , 1 = … = t n , n − 1 {\displaystyle t_{2,1}=\ldots =t_{n,n-1}} , so spricht man von einer Tridiagonal-Toeplitz-Matrix.

Eine Tridiagonalmatrix ist sowohl ein Spezialfall einer Bandmatrix als auch einer Hessenbergmatrix. Eine diagonaldominante Tridiagonalmatrix ist immer regulär.

Lineare Gleichungssysteme mit einer Tridiagonalmatrix lassen sich mit einem Aufwand von O(n) effizient lösen. Entweder mit dem sehr schnellen Thomas-Algorithmus oder bei Stabilitätsproblemen mit Hilfe des Gauß-Verfahrens mit Pivotisierung. Gleichungssysteme mit Tridiagonalmatrizen können also selbst bei vergleichsweise großer Dimension mittels eines direkten Lösers berechnet werden.

• Gerhard Opfer: Numerische Mathematik für Anfänger. Eine Einführung für Mathematiker, Ingenieure und Informatiker. 4., durchgesehene Auflage. Vieweg, Braunschweig u. a. 2002, ISBN 3-528-37265-6.

• Dreiecksmatrix
• Bidiagonalmatrix
• Blocktridiagonalmatrix