Wahrscheinlichkeitsverteilung für eine verrauschte Sinuswelle

12

Ich möchte eine Wahrscheinlichkeitsverteilung von Abtastpunkten aus einer oszillierenden Funktion analytisch berechnen, wenn ein Messfehler vorliegt. Ich habe bereits die Wahrscheinlichkeitsverteilung für den Teil "ohne Rauschen" berechnet (ich werde dies am Ende setzen), aber ich kann nicht herausfinden, wie "Rauschen" eingeschlossen werden soll.

Numerische Schätzung

Stellen Sie sich zur Verdeutlichung vor, es gibt eine Funktion der Sie während eines einzelnen Zyklus zufällig Punkte auswählen. Wenn Sie die Punkte in einem Histogramm ablegen, erhalten Sie etwas, das mit der Verteilung zusammenhängt.y(x)=sin(x)

Ohne Lärm

Zum Beispiel ist hier die und das entsprechende Histogrammsin(x)

Geben Sie hier die Bildbeschreibung ein

Mit Lärm

Wenn nun ein Messfehler auftritt, ändert sich die Form des Histogramms (und daher denke ich, dass die zugrunde liegende Verteilung). Beispielsweise

Geben Sie hier die Bildbeschreibung ein

Analytische Berechnung

Hoffentlich habe ich Sie davon überzeugt, dass es einen Unterschied zwischen den beiden gibt. Jetzt werde ich aufschreiben, wie ich den Fall "ohne Rauschen" berechnet habe:

Ohne Lärm

y(x)=sin(x)

Wenn dann die Zeiten, zu denen wir abtasten, gleichmäßig verteilt sind, muss die Wahrscheinlichkeitsverteilung für erfüllen:y

P(y)dy=dx2π

dann seit

dxdy=ddy(arcsin(y))=11y2

und so

P(y)=12π1y2

Dies passt bei entsprechender Normalisierung zu dem im Fall "kein Rauschen" erzeugten Histogramm.

Mit Lärm

Meine Frage lautet also: Wie kann ich Rauschen analytisch in die Verteilung einbeziehen? Ich denke, es ist so etwas wie das clevere Kombinieren der Verteilungen oder das Einbeziehen von Rauschen in die Definition von , aber ich habe keine Ideen und Möglichkeiten, mich vorwärts zu bewegen, sodass Hinweise / Tipps oder sogar empfohlene Lektüre viel sind geschätzt.y(x)

Greg
quelle

Antworten:

10

Dies hängt davon ab, wie der Rauschprozess strukturiert ist.

Y

XixYi|Xi=xisin(xi)Yg

ϵiN(0,σ2)f(ϵ)=12πσexp(ϵ22σ2)fg

fY+ϵ(z)=(fg)(z)=fY(y)fϵ(zy)dy=fY(zw)fϵ(w)dw

Geben Sie hier die Bildbeschreibung ein

(Diese Faltung wurde numerisch durchgeführt. Ich weiß nicht, wie nachvollziehbar dieses Integral in diesem Beispiel ist, weil ich es nicht versucht habe.)

Glen_b - Monica neu starten
quelle
Wunderbares Zeug, mir fehlte die Idee der "Faltung", wie Ihre Zahlen zeigen, dass dies genau richtig ist. Ich muss nur die Integration versuchen. Danke
Greg
2
Sie finden es vielleicht unlösbar, aber es ist normalerweise nicht schwierig, das Ergebnis numerisch zu approximieren.
Glen_b -State Monica
0

Ich denke, der abgeleitete Ausdruck für P (x) ist um den Faktor zwei versetzt. Eine gleichmäßig verteilte Abtastzeit entspricht einer gleichmäßig verteilten Phase über das Intervall -pi, pi. Die trigonometrische Funktion verteilt die Wahrscheinlichkeit über das y-Intervall {-1,1}. Die Integration von P (y) über dieses Intervall muss = 1 sein, nicht 2, wie mit Ihrem obigen Integranden erhalten. Ich denke P (y) = 1 / (pi Sqrt (1-y ^ 2))

Jeff Patterson
quelle
Es könnte gut sein, deshalb habe ich "mit angemessener Normalisierung" gesagt, da ich zu faul war, um darüber nachzudenken. Vielen Dank.
Greg