PUB - Publikationen an der Universität Bielefeld

Der Lanczosalgorithmus ist ein iteratives Verfahren zur Berechnung (bzw. Approximation) einiger Eigenwerte und Eigenvektoren großer, dünn besetzter symmetrischer Matrizen. Aufbauend auf Ergebnissen von Christopher Conway Paige wird das numerische Verhalten dieses Verfahrens insbesondere hinsichtlich der Bildung von Cluster von Ritzwerten analysiert: Es stellt sich heraus, dass sich Cluster von Ritzwerten nur um einen Eigenwert bilden können. Darüber hinaus berechnet der Lanczosalgorithmus (bzw. das cg-Verfahren) implizit eine Folge von Gaußquadraturapproximationen an ein durch den Startvektor des Verfahrens und den Eigenwerten der Matrix, für die die Eigenwertapproximationen bestimmt werden, definiertes Riemann-Stieltjes Integral. Es wird die so genannte Stabilisierung der Gewichte der Gaußquadratur sowohl beim Rechnen in exakter als auch in endlicher Arithmetik nachgewiesen. Ein Lanczos-ähnlicher Algorithmus für unitäre Matrizen zeigt numerisch ein analoges Verhalten wie der symmetrische Lanczosalgorithmus. Es wird der analoge Zusammenhang zwischen Orthogonalitätsverlust unter den Lanczosvektoren und der Konvergenz von Ritzwerten aufgezeigt.