Rekonstruktion von "Wellenfunktions"-Phasen aus |ψ(x)||ψ(x)||\psi(x)| und |ψ~(p)||ψ~(p)||\tilde \psi(p)|

Betrachten Sie eine "Wellenfunktion" ψ ( x ) , die eine Fourier-Transformation aufweist ψ ~ ( p )

Angenommen, wir wüssten es für jeden x , | ψ ( x ) | 2 , und das wissen wir, für jeden p , | ψ ~ ( p ) | 2 .

Haben wir genug Informationen, um die "Wellenfunktion" zu rekonstruieren? ψ ( x ) , dh erhalten Sie die Phase der "Wellenfunktion" für alle x (Bis auf eine globale Phase interessieren uns nur die relativen Phasen zwischen den x )?

Lesern könnte auch dieses verwandte Papier gefallen .
Dies hängt mit dem Phasenwiederherstellungsproblem bei der Beugung zusammen. Dort haben Sie die Intensität des gebeugten Lichts - analog zu Ihrer | ψ ~ ( p ) | 2 , aber anstatt den Modul von zu kennen ψ Überall kennen Sie die Phase des Analogen – die Dichte der Streuer. Insbesondere wissen Sie, dass die Phase Null ist, da die Dichte real ist. Ihr Problem ist also im Grunde wie das Phasenwiederherstellungsproblem, außer dass Sie den Modul von kennen ψ statt Phase. Da die Informationen vergleichbar sind, würde ich vermuten, dass dies lösbar ist.
Verwandte Frage zu MO.SE: mathoverflow.net/q/188931/13917
Diese Frage wurde übrigens ursprünglich von Pauli gestellt
Hier ist die Referenz: W. Pauli, in Quantentheorie, herausgegeben von H. Geiger und K. Scheel, Handbuch der Physik Bd. 24 (Springer, Berlin, 1933), Pt. 1, p. 98

Antworten (4)

Nein. Betrachten Sie irgendeinen Zustand mit einer Impulswellenfunktion, die symmetrisch um Null ist. Seine normquadratischen Wahrscheinlichkeitsverteilungen im Ort-Raum und im Impuls-Raum werden durch die Zeitumkehr nicht verändert, obwohl dies bei der Wellenfunktion eindeutig der Fall ist.

Hier ist ein explizites Beispiel. Nehmen Sie das vier Gaußsche Wellenpaket der mittleren Positionen x 0 oder x 0 , mittlere Impulse p 0 oder p 0 , und räumliche Ausbreitung σ :

ψ ( ± x , ± p ) ( x ) e ( x x x 0 ) ) 2 / 4 σ 2 ich x ( ± p p 0 ) .

(Hier, ± x und ± p sind zwei binäre Variablen, die die Werte Plus oder Minus annehmen. x und p sind ihre Gegensätze.) Betrachten Sie nun die beiden Superpositionen

ϕ a w a j = ψ ( + , + ) + ψ ( , ) ,

ϕ t Ö w a r d = ψ ( + , ) + ψ ( , + ) .

ϕ a w a j ist die Überlagerung zweier Wellenpakete, die durch getrennt sind 2 x 0 sich mit relativer Geschwindigkeit voneinander entfernen 2 p 0 / m . ϕ t Ö w a r d ist dasselbe mit den Paketen, die sich aufeinander zu bewegen. Das kann man überprüfen | ϕ t Ö w a r d ( x ) | 2 = | ϕ a w a j ( x ) | 2 und | ϕ ~ t Ö w a r d ( p ) | 2 = | ϕ ~ a w a j ( p ) | 2 .

Ich weiß nicht, ob es andere Beispiele als die Zeitumkehr gibt.

Ein intuitiver dimensionaler Grund, warum es nicht funktionieren konnte: ein Zustandsvektor in C N + 1 wird durch 2N reelle Koordinaten beschrieben (eine komplexe Dimension ist irrelevant), ebenso wie seine Fourier-Transformation. Wenn wir nur die normalisierten quadratischen Moduln der Komponenten betrachten, haben wir auch 2N reelle Zahlen, also wenn diese tatsächlich unabhängig wären, sollten wir in der Lage sein, den ursprünglichen Vektor zu rekonstruieren.

Sie sind jedoch nicht unabhängig: Die bekannteste Abhängigkeit zwischen dem quadratischen Modul einer Funktion und dem ihrer Fourier-Transformation ist die Heisenberg-Unschärferelation (von der Analoga in einer endlichdimensionalen Umgebung existieren, wenn Sie in einer Umgebung bleiben wollen, in der Dimension Zählen ist einfach).

Ein weiteres wird durch das Paley-Wiener-Theorem bereitgestellt, das impliziert, dass eine kompakt unterstützte Funktion eine Fourier-Transformation hat, die auf keiner offenen Menge identisch Null ist.

Die kleinsten Gegenbeispiele treten in einem 2-dimensionalen Zustandsraum auf: die Fourier-Transformation von ( a b ) ist ( a + b a b ) (bis auf eine multiplikative Konstante), also zB ( 1 ich ) und ( 1 ich ) beschreiben unterschiedliche Zustände, aber sie haben gleiche Moduli, ebenso wie ihre Fourier-Transformationen.

Ja, Ihr Gegenbeispiel ist eine diskrete Version von mir. (Ich kombiniere unsere Antworten gerne, obwohl ich nicht weiß, wie das gemacht wird.) Eine natürliche Folgefrage ist, ob das Paar (|f(x)|^2, |f(p)|^2) genügt, um einen Zustand bis zur Zeitumkehr eindeutig zu bestimmen . Basierend auf Ihrem Dimensionsargument vermute ich, dass dies der Fall ist (möglicherweise mit zusätzlichen unbestimmten diskreten Symmetrien).
sehr klar - fantastische Einfachheit! @JessRiedel Ich stimme zu, dass eine eindeutige Bestimmung bis zur Zeitumkehr wahrscheinlich erscheint. Werde noch mal drüber nachdenken - hört sich so an als wärst du es auch.

Kommentare zur Frage (v1):

I) Rekonstruktion von Phasen aus Modul 1 | f ( x ) | eines Signals f ( x ) und Modul | f ~ ( k ) | seines Fourier-transformierten (FT) Signals

(1) f ~ ( k )   :=   1 2 π R d x   e ich k x f ( x )

ist ein interessantes und wahrscheinlich gut untersuchtes technisches Problem, entweder für kontinuierliche oder diskrete Fourier-Transformation.

II) Beispiel: Ein Gauß-Signal

(2) f ( x )   =   EIN e a 2 x 2 + b x , EIN , a , b C , R e ( a )   >   0 ,

mit Fouriertransformation

(3) f ~ ( k )   =   EIN a exp [ ( k + ich b ) 2 2 a ]   =   EIN a exp [ a ¯ 2 | a | 2 ( k + ich b ) 2 ]   =   EIN a exp [ a ¯ 2 | a | 2 ( k 2 + 2 ich b k b 2 ) ] ;

mit Modul

(4) | f ( x ) |   =   | EIN | exp [ R e ( a ) 2 x 2 + R e ( b ) x ] ,

und

(5) | f ~ ( k ) |   =   | EIN | | a | exp [ R e ( a ) k 2 2 ich m ( a ¯ b ) k R e ( a ¯ b 2 ) 2 | a | 2 ] ,

beziehungsweise. Interessant ist, dass wenn man zusätzlich weiß, dass das Signal von Gaußscher Form ist (2), dann lässt sich aus (4) und (5) die Konstante rekonstruieren a bis auf evtl. eine Zeichenmehrdeutigkeit von ich m ( a ) ; die Konstante EIN bis zu einer Phase; und die Konstante b ist einzigartig für die gegebene Wahl von a .

III) Das obige Gaußsche Beispiel weckt die Hoffnung, dass der Modul eines Signals und der Modul seiner FT ausreichend komplementäre Informationen sind, so dass eine Rekonstruktion bis zu möglicherweise einer endlichen Anzahl von selbstkonsistenten Lösungen und modulo einer gesamten globalen Phase möglich ist.

IV) Wir spekulieren, dass es in der Praxis möglich sein könnte, ein selbstkonsistentes Signal aus seinem Modulus und dem Modulus seiner FT über einen iterativen Festkommaalgorithmus zu rekonstruieren 2 : Erste Fourier-Transformation des bloßen Moduls des Signals; als nächstes die Phasen des Ergebnisses mit dem anfänglich gegebenen FT-Modul multiplizieren; dann inverse FT; als nächstes die Phasen des Ergebnisses mit dem anfänglich gegebenen Modul multiplizieren; dann FT; und so weiter, bis eine selbstkonsistente Festpunktkonfiguration auf jeder Seite der FT erreicht ist.

--

1 In der Signalanalyse wird der Modul auch als Amplitude, Magnitude oder Absolutwert bezeichnet.

2 Update: Es stellt sich heraus, dass dieser Algorithmus existiert und als Gerchberg-Saxton-Algorithmus bekannt ist (Huttipp: WetSavannaAnimal alias Rod Vance).

Was die Rekonstruktion betrifft, liegen Sie genau richtig - was Sie "spekuliert" haben, ist real und als Gerchberg-Saxton-Algorithmus bekannt . Ich kenne die Theorie der Konvergenz solcher Dinge nicht - wir haben keine Kontraktionsabbildung, sondern eine Isometrie für die FT - aber ich weiß, dass so etwas in der Praxis sehr gut funktioniert, und ich bin mir sicher, dass die Konvergenz hätte studiert worden. Ich wäre überrascht, wenn es nicht ziemlich universell wäre - ich war verblüfft, wie einfach es ist, dies selbst mit sehr lauten Signalen zum Laufen zu bringen (Isometrie hilft hier eindeutig).

Antwort ist nein. Betrachten Sie zwei Wellen: ψ 1 ( x , t ) = ψ ( x , t ) e ich a t   und   ψ 2 ( x , t ) = ψ ( x , t ) e ich β t . Sicherlich sind diese Wellen anders, aber | ψ 1 ( x , t ) | = | ψ 2 ( x , t ) | = | ψ ( x , t ) | . Aus physikalischen Messungen kann man also keine Welle im Detail bestimmen.

Globale Phasen sind irrelevant, und jedes Verfahren zum Abrufen von Phasen befasst sich mit relativen Phasen, wie es in dieser Frage angegeben ist.