Technopedia Center
PMB University Brochure
Faculty of Engineering and Computer Science
S1 Informatics S1 Information Systems S1 Information Technology S1 Computer Engineering S1 Electrical Engineering S1 Civil Engineering

faculty of Economics and Business
S1 Management S1 Accountancy

Faculty of Letters and Educational Sciences
S1 English literature S1 English language education S1 Mathematics education S1 Sports Education
teknopedia

teknopedia

teknopedia

teknopedia

teknopedia
  • Registerasi
  • Brosur UTI
  • Kip Scholarship Information
  • Performance
  1. Weltenzyklopädie
  2. Inverse Iteration – Wikipedia
Inverse Iteration – Wikipedia
aus Wikipedia, der freien Enzyklopädie

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 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.

Beschreibung

[Bearbeiten | Quelltext bearbeiten]

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

Statt wie bei der Potenzmethode in jedem Schritt die Matrix mit einem Vektor zu multiplizieren, wird nun ein lineares Gleichungssystem gelöst, da ( A − θ I ) − 1 {\displaystyle (A-\theta I)^{-1}} {\displaystyle (A-\theta I)^{-1}} nicht explizit verfügbar ist. Diese Matrix ist schlechter konditioniert, je näher λ {\displaystyle \lambda } {\displaystyle \lambda } an θ {\displaystyle \theta } {\displaystyle \theta } liegt, allerdings hat der Fehler eine dominante Komponente in Richtung des gesuchten Eigenvektors, so dass das Verfahren praktisch nutzbar ist.

Algorithmus

[Bearbeiten | Quelltext bearbeiten]

Gegeben sei eine quadratische Matrix A ∈ C n × n {\displaystyle A\in \mathbb {C} ^{n\times n}} {\displaystyle A\in \mathbb {C} ^{n\times n}}, ein Startvektor x 0 ∈ R n {\displaystyle x_{0}\in \mathbb {R} ^{n}} {\displaystyle x_{0}\in \mathbb {R} ^{n}} und ein Shift θ ∈ R {\displaystyle \theta \in \mathbb {R} } {\displaystyle \theta \in \mathbb {R} } so dass ( A − θ I ) {\displaystyle (A-\theta I)} {\displaystyle (A-\theta I)} regulär ist. Der Startvektor kann bis auf eine Lebesgue-Nullmenge beliebig gewählt werden.

Für k = 1 , 2 , . . . {\displaystyle k=1,2,...} {\displaystyle k=1,2,...}

  1. q k = x k − 1 ‖ x k − 1 ‖ {\displaystyle q_{k}={\frac {x_{k-1}}{\|x_{k-1}\|}}} {\displaystyle q_{k}={\frac {x_{k-1}}{\|x_{k-1}\|}}}
  2. Löse ( A − θ I ) x k = q k {\displaystyle (A-\theta I)x_{k}=q_{k}} {\displaystyle (A-\theta I)x_{k}=q_{k}}

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

λ k = x k T A x k x k T x k {\displaystyle \lambda _{k}={\frac {x_{k}^{T}Ax_{k}}{x_{k}^{T}x_{k}}}} {\displaystyle \lambda _{k}={\frac {x_{k}^{T}Ax_{k}}{x_{k}^{T}x_{k}}}}

Erweiterungen

[Bearbeiten | Quelltext bearbeiten]

Wählt man in jedem Schritt über θ = λ k {\displaystyle \theta =\lambda _{k}} {\displaystyle \theta =\lambda _{k}} einen neuen Shift so erhält man die Rayleigh-Quotienten-Iteration.

Literatur

[Bearbeiten | Quelltext bearbeiten]
  • Gene H. Golub, Charles F. van Loan: Matrix Computations
  • James H. Wilkinson: The Algebraic Eigenvalue Problem
  • Hans-Rudolf Schwarz, Norbert Köckler: Numerische Mathematik. 5., überarbeitete Auflage. Teubner, Stuttgart u. a. 2004, ISBN 3-519-42960-8.
Abgerufen von „https://de.teknopedia.teknokrat.ac.id/w/index.php?title=Inverse_Iteration&oldid=178014568“
Kategorie:
  • Numerische lineare Algebra

  • indonesia
  • Polski
  • العربية
  • Deutsch
  • English
  • Español
  • Français
  • Italiano
  • مصرى
  • Nederlands
  • 日本語
  • Português
  • Sinugboanong Binisaya
  • Svenska
  • Українська
  • Tiếng Việt
  • Winaray
  • 中文
  • Русский
Sunting pranala
Pusat Layanan

UNIVERSITAS TEKNOKRAT INDONESIA | ASEAN's Best Private University
Jl. ZA. Pagar Alam No.9 -11, Labuhan Ratu, Kec. Kedaton, Kota Bandar Lampung, Lampung 35132
Phone: (0721) 702022
Email: pmb@teknokrat.ac.id