Wielandt-Iteration

Wielandt-Iteration

Die inverse Iteration ist ein numerisches Verfahren zur Berechnung von Eigenwerten von Matrizen. Sie ist eine Variante der von-Mises-Iteration, mit deren Hilfe allerdings beliebige Eigenwerte berechnet werden können. Das Verfahren wurde 1944 von Helmut Wielandt bei der Stabilitätsanalyse von Strukturen, die kleine Störungen bekannter Systeme sind, eingeführt. In diesem Fall sind gute Approximationen für die relevanten Eigenwerte bekannt, und man erhält rasche Konvergenz.

Inhaltsverzeichnis

Beschreibung

Ist λ ein Eigenwert von A und x der zugehörige Eigenvektor, so ist λ − θ ein Eigenwert von (A − θI) zum Eigenvektor x, wobei I die Einheitsmatrix ist. Des Weiteren ist dann \frac{1}{\lambda-\theta} ein Eigenwert von (A − θI) − 1 zum Eigenvektor x. Ist λ nun der Eigenwert von A, der θ am nächsten liegt, so ist \frac{1}{\lambda-\theta} der betragsmäßig größte Eigenwert von (A − θI) − 1. Wendet man nun auf (A − θI) − 1 die Potenzmethode an, so konvergiert xk gegen den Eigenvektor zum Eigenwert λ von A, der θ am nächsten liegt.

Algorithmus

Gegeben sei eine quadratische Matrix A\in\mathbb{R}^{n\times n}, ein Startvektor x_0\in\mathbb{R}^n und ein Shift \theta\in\mathbb{R} so dass (A − θI) regulär ist.

Für k = 1,2,....

  1. q_k=\frac{x_{k-1}}{\|x_{k-1}\|}
  2. Löse (A − θI)xk = qk

Über den Rayleigh-Quotienten erhält man eine Näherung für den zugehörigen Eigenwert.

\lambda_k=\frac{x_k^TAx_k}{x_k^Tx_k}

Erweiterungen

Wählt man in jedem Schritt über θ = λk einen neuen Shift so erhält man die Rayleigh-Quotienten-Iteration.

Literatur

  • Gene H. Golub / Charles F. van Loan Matrix Computations
  • J.H. Wilkinson The Algebraic Eigenvalue Problem

Wikimedia Foundation.

Игры ⚽ Поможем сделать НИР

Schlagen Sie auch in anderen Wörterbüchern nach:

  • Helmut Wielandt — (* 19. Dezember 1910 in Niedereggenen; † 14. Februar 2001 in Schliersee) war ein deutscher Mathematiker. Sein Hauptarbeitsgebiet war die Gruppentheorie, speziell die Theorie der Permutationsgruppen. Wielandts Beweis der für di …   Deutsch Wikipedia

  • Inverse Iteration — Die inverse Iteration ist ein numerisches Verfahren zur Berechnung von Eigenwerten und Eigenvektoren von Matrizen. Sie ist eine Variante der von Mises Iteration, mit deren Hilfe allerdings beliebige Eigenwerte berechnet werden können. Das… …   Deutsch Wikipedia

  • Orr-Sommerfeld-Gleichung — Dieser Artikel befasst sich mit dem Anwachsen kleiner Störungen in einem Strömungsfeld. Für weitere Bedeutungen des Begriffs Lineare Stabilitätstheorie siehe Stabilitätstheorie. Die Lineare Stabilitätstheorie (kurz LST) beschreibt in einer… …   Deutsch Wikipedia

  • Squire-Gleichung — Dieser Artikel befasst sich mit dem Anwachsen kleiner Störungen in einem Strömungsfeld. Für weitere Bedeutungen des Begriffs Lineare Stabilitätstheorie siehe Stabilitätstheorie. Die Lineare Stabilitätstheorie (kurz LST) beschreibt in einer… …   Deutsch Wikipedia

  • Lineare Stabilitätstheorie — Die Lineare Stabilitätstheorie (kurz LST) beschreibt physikalisch in einer Strömung das Anwachsen wellenförmiger Störungen mit kleiner Amplitude. Durch Vorgabe eines stationären Strömungsfeldes lassen sich Anfachungsraten und Form der linearen… …   Deutsch Wikipedia

  • Liste des articles de mathematiques — Projet:Mathématiques/Liste des articles de mathématiques Cette page recense les articles relatifs aux mathématiques, qui sont liés aux portails de mathématiques, géométrie ou probabilités et statistiques via l un des trois bandeaux suivants  …   Wikipédia en Français

  • Hellmuth Kneser — (* 16. April 1898 in Dorpat; † 23. August 1973 in Tübingen) war ein deutscher Mathematiker. Hellmuth Kneser, ca. 1930. Inhaltsverzeichnis …   Deutsch Wikipedia

  • Liste mathematischer Sätze — Inhaltsverzeichnis A B C D E F G H I J K L M N O P Q R S T U V W X Y Z A Satz von Abel Ruffini: eine allgemeine Polynomgleichung vom …   Deutsch Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”