Hallo,
ich habe in Matlab einen iterativen Algorithmus, den s.g. Gerchberg-Saxton-Alg., implementiert. Nun bin ich in einer Arbeit darauf gestossen, dass man die Aufgabenstellung auch anders formulieren kann, nämlich als Eigenwertproblem: S.52 unten/S.53 oben
Nun stehe ich vor zweierlei Problemen:1. ich kenne b…
Hi,
also in matlab nimmst Du die eigs-Methode (auf keinen Fall eig, das ist ein schlechter Zufallszahlengenerator).
Da ist alles drin, was Du erwähnt hast, Householder, um das Problem auf Hessenbergform zu bringen, danach implizite Multishift-Iteration. Halt das, was in Lapack implementiert ist und noch etwas Tuning oben drauf.
Gruß, Lutz