Householdertransformation schlägt fehl |
04.12.2010, 22:07 | Duedi | Auf diesen Beitrag antworten » | |||||
Householdertransformation schlägt fehl Programmiere gerade eine QR-Zerlegung mit Householderspiegelung und in einem Schritt scheint was nicht zu funktionieren (Werte gerundet, da in Scilab): Ich möchte die erste Zeile zum Einheitsvektor machen: Dann ist und Nun ist aber jetzt Die erste Spalte ist jetzt aber alles Andere als Einheitsvektor Weiß jemand Rat? Dank schonmal, Duedi |
|||||||
04.12.2010, 22:28 | Quadratzahl-Jan | Auf diesen Beitrag antworten » | |||||
Hi! So wie ich die HHT kenne, müßte v = a_1 - ||a_1|| * (1,0,0) sein. Hast du vergessen dieses Vielfache des Einheitsvektors abzuziehen ? Oder ist das so eure Methode? a_1 soll hier der erste Spaltenvektor der Matrix A sein, also a_1 = (0,1,2) |
|||||||
04.12.2010, 22:32 | Duedi | Auf diesen Beitrag antworten » | |||||
Genau, das müsste ich eigentlich abziehen, aber (bei uns) mit dem Signum des (1,1)-Eintrags von A gewichtet, was in dem Fall 0 ist. |
|||||||
04.12.2010, 23:06 | Quadratzahl-Jan | Auf diesen Beitrag antworten » | |||||
OK. Zur sicherheit: kann es sein, dass man Sign(0)=1 gesetzt hat? Meine mich an sowas erinnern zu können. Gucke aber noch mal weiter, gehe davon aus, dass es nicht das Problem löst. (Du hast ja schon geschrieben, dass sign(0)=0. aber ein Versuch ist es wert ) |
|||||||
04.12.2010, 23:48 | Duedi | Auf diesen Beitrag antworten » | |||||
Tatsächlich, das löst einen Großteil der Probleme, vielen Dank! (Sind nur noch ein paar Implementierungsfehler, die ich ausmerzen muss) |
|||||||
05.12.2010, 14:10 | tigerbine | Auf diesen Beitrag antworten » | |||||
|
|||||||
Anzeige | |||||||
|
|||||||
05.12.2010, 23:36 | Duedi | Auf diesen Beitrag antworten » | |||||
Hm ich schaue mir morgen mal an, warum das bei dir klappt und bei mir nicht. (Eigentlich habe ich es mit zum Laufen gebracht, jedenfalls sind die Ergebnisse von Maple gleich). Danke auf jeden Fall! |
|
Verwandte Themen
Die Beliebtesten » |
Die Größten » |
|
Die Neuesten » |
|