Wie ist die Phasen- und Amplitudenantwort von weißem Rauschen?

16

Ich möchte weißes Rauschen im Frequenzbereich erzeugen und es dann mit Python in den Zeitbereich umwandeln. Um das Problem zu verstehen, habe ich einfach weißes Rauschen im Zeitbereich erzeugt und es in den Freq-Bereich umgewandelt:

import scipy.signal as sg
import numpy as np
import matplotlib.pyplot as plt

e = np.random.normal(0,1,1e3)
E = sg.fft(e)

plt.figure("Bode plot")
plt.subplot(211)
plt.title("Magitude")
plt.plot(abs(E))
plt.subplot(212)
plt.title("Phase")
plt.plot(np.angle(E))
plt.show()

Ich sehe überhaupt nicht so aus, wie ich es erwartet hatte: Bode-Plot von weißem Rauschen Fragen:

  • Soll weißes Rauschen nicht eine flache Amplitudenantwort haben? (gleiche Beträge für alle Frequenzen)
  • Wie ist die Beziehung zwischen der Standardabweichung (1 in meinem Beispiel) und der Größe und Phase?

Danke im Voraus!

Uffe
quelle

Antworten:

21

Soll weißes Rauschen nicht eine flache Amplitudenantwort haben? (gleiche Beträge für alle Frequenzen)

Die erwartete Amplitudenantwort von weißem Rauschen ist flach (JasonR nennt dies die spektrale Leistungsdichte). Eine bestimmte Instanz einer Sequenz mit weißem Rauschen hat keine präzise flache Reaktion (dies ist der Kommentar von JasonR, der als Leistungsspektrum bezeichnet wird).

Tatsächlich ist die Fourier-Transformation von weißem Rauschen ... weißes Rauschen!

Wie ist die Beziehung zwischen der Standardabweichung (1 in meinem Beispiel) und der Größe und Phase?

Es gibt keine Beziehung zwischen der Standardabweichung und der Phase. Angenommen, ist stationäres weißes Rauschen mit dem Mittelwert Null und der Standardabweichung . Dann ist die Autokorrelation (Kovarianz):n(t)σ

Rnn(τ)=E[n(t)n(t+τ)]=σ2δ(τ)

Die spektrale Leistungsdichte beträgt also nur (obwohl es für die diskrete Zeit eine Skalierung auf der Grundlage der Signaldauer gibt).σ2


Fragen von Kommentar:

  1. Wenn Sie sagen, dass die Fourier-Transformation auch weißes Rauschen ist, wie kann ich dann den Standardwert messen, wenn die Transformation komplex ist? Real-, Imaginärteil oder irgendeine Kombination?

Angenommen, unser Rauschen ist zeitdiskret und ist (Mittelwert Null, Gaußsches, weißes Rauschen mit Varianz ). Dann ist die Transformation:n[m]σ2

N[k]=m=0M1n[m]ej2πmk/M=m=0M1n[m]cos(2πmk/M)+jn[m]sin(2πmk/M)

und der erwartete Wert ist:

E[N[k]]=E[m=0M1n[m]ej2πmk/M]=m=0M1E[n[m]]ej2πmk/M=0

Die Varianz des Realteils ergibt sich aus:

E[(N[k])2]=E[m=0M1n[m]cos(2πmk/M)p=0M1n[p]cos(2πpk/M)]=E[m=0M1p=0M1n[m]n[p]δ[np]cos(2πmk/M)cos(2πpk/M)]=m=0M1E[n[m]2]cos2(2πmk/M)=σ2m=0M1cos2(2πmk/M)=σ2(M2+cos(M+1)2πk/Msin(2πMk/M)2sin(2πk/M)   )=σ2M2

Ich glaube, der Imaginärteil wird sich genauso verhalten.

  1. Könnten Sie mir bitte erklären, wie die Dauer des Signals mit der spektralen Leistungsdichte zusammenhängt (für zeitdiskrete Situationen)?

Ich glaube, dass (basierend auf der obigen Ableitung) die spektrale Leistungsdichte (der erwartete Wert des Quadrats der DFT) linear als die Dauer skaliert.

  1. Wenn die Phase nicht vom std-dev beeinflusst wird, was bestimmt die 3-Grad-Amplitude und die Art der Verteilung (scheint eher einheitlich als normal zu sein)

Schauen Sie sich die Tabelle auf Seite 2 dieser PDF-Datei an . Sie besagt, dass das Argument (die Phase) der Koeffizienten gleichmäßig verteilt wird, wie Sie angeben. Screenshot der folgenden Tabelle.

Bildbeschreibung hier eingeben

Peter K.
quelle
3
Insbesondere sind die beiden Konzepte, die das OP verwirrt, die Leistungsspektraldichte von weißem Rauschen und das Leistungsspektrum einer bestimmten Realisierung eines Zufallsprozesses mit weißem Rauschen.
Jason R
Vielen Dank! Ich habe einige Anschlussfragen. 1: Wenn Sie sagen, dass die Fouriertransformation auch weißes Rauschen ist, wie kann ich dann den std-dev messen, wenn die Transformation komplex ist? Real-, Imaginärteil oder irgendeine Kombination? 2: Könnten Sie mir bitte erklären, wie die Dauer des Signals mit der spektralen Leistungsdichte zusammenhängt (für diskrete Zeitsituationen)? Verteilung (scheint eher einheitlich als normal zu sein)
Uffe
Perfekt! Dies erklärt alle Funktionen des Frequenzgangs. Mir ist jetzt klar, dass die Amplitude der Phase nicht 3, sondern war und dass die Standardabweichung des Real- und Imaginärteils einfach von der Anzahl der Elemente im Vektor abhängt. Ich kann auch experimentell bestätigen, dass der Imaginärteil ebenfalls eine Varianz von . σ 2 Mπσ2M2
Uffe
Dies ist der aktuelle Link zum oben genannten PDF-Dokument ( radarsp.weebly.com/uploads/2/1/4/7/21471216/dft_of_noise.pdf ), der fehlerhaft ist.
Gast
@ Gast Danke! Versuchen Sie in Zukunft einfach, die Antwort mit dem neuen Link zu bearbeiten. Es wird nicht direkt ausgeführt, da es von einem Benutzer mit höheren Wiederholungen überprüft werden muss, aber es wird dort ankommen (und Sie erhalten dabei +2 Wiederholungen).
Peter K.