Wurzelzieher

Inhalt

Householdertransformation
Definition und Eigenschaften

Konstruktion einer spezifischen Spiegelung

&a-2\,(a+\lambda e)\,\frac{/ &a-2\,(a+\lambda e)\,/ &-\lambda e\\[0.5em]

Anwendung: QR-Zerlegung

  

Pseudocode

Literatur

 

 

Householdertransformation

Definition und Eigenschaften

Illustration der Householder-Transformation in 2D

Die Spiegel-Hyperebene kann durch einen Normalenvektor v, also einen Vektor, der orthogonal zur Hyperebene ist, definiert werden. Ist v als Spaltenvektor gegeben und I die Einheitsmatrix, dann wird die oben beschriebene lineare Abbildung durch die folgende Matrix dargestellt:

Dabei bezeichnet vT die Transponierte des Spaltenvektors v, also einen Zeilenvektor. Der Nenner vT v ist das Skalarprodukt von v mit sich selbst, vvT das dyadische Produkt, eine Matrix. Die Matrix beschreibt die Orthogonalprojektion auf die durch v gegebene Richtung.Ist v auf die Länge eins normiert, also vT v = 1, so vereinfacht sich die Formel zu

H = I - 2vvT

Die Spiegelungseigenschaft ersieht man daraus, dass


,

wobei das Skalarprodukt bezeichnet. Der Term entspricht dabei dem Abstand des Punktes x zur Hyperebene . Der Vektor x wird also in zwei orthogonale Anteile zerlegt, wobei der erste Anteil in der Hyperebene liegt und der zweite ein Vielfaches des Vektors v ist. Unter der Spiegelung wird der Anteil in der Ebene invariant gelassen, der Anteil in Richtung v, also senkrecht zur Ebene, wird "umgeklappt", also nun abgezogen statt addiert.

Die Householder-Matrix hat folgende Eigenschaften:

  • Sie ist symmetrisch: H = HT
  • Sie ist orthogonal: H-1 = HT
  • Sie ist involutorisch: H2 = I (Dies folgt aus der Symmetrie und der Orthogonalität.)
  • Sie hat die Eigenwerte -1 (mit Vielfachheit 1) und 1 (mit Vielfachheit n-1).
  • Matrix-Vektor-Multplikationen mit H sind schnell berechenbar.

 

 

 

 

Copyright- und Lizenzinformationen: Diese Seite basiert auf dem Artikel Householdertransformation aus der freien Enzyklοpädιe Wιkιpedιa und steht unter der Lizenz Creative Commons CC-BY-SA 3.0 Unported (Kurzfassung). Liste der Autoren

Anbieterkennzeichnung

 



Load: 10; Render: 0; Total: 10