Jaka jest faza i wielkość odpowiedzi białego szumu?

16

Chciałbym stworzyć biały szum w dziedzinie częstotliwości, a następnie przekształcić go w domenę czasu za pomocą Pythona. Aby zrozumieć problem, po prostu wygenerowałem biały szum w dziedzinie czasu i przekształciłem go w domenę częstotliwości:

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()

Nie patrzę wcale tak, jak się spodziewałem: Wróży wykres białego szumu Pytania:

  • Czy biały szum nie powinien mieć płaskiej odpowiedzi? (równe kwoty dla wszystkich częstotliwości)
  • Jaki jest związek między odchyleniem standardowym (1 w moim przykładzie) a wielkością i fazą?

Z góry dziękuję!

Uffe
źródło

Odpowiedzi:

21

Czy biały szum nie powinien mieć płaskiej odpowiedzi? (równe kwoty dla wszystkich częstotliwości)

Oczekiwać odpowiedzi wielkość białego szumu jest płaska (to co JasonR nazywa widmowa gęstość mocy). Żaden konkretny przypadek sekwencji białego szumu nie będzie miał dokładnie płaskiej odpowiedzi (to, co JasonR nazywa widmem mocy).

W rzeczywistości transformacja Fouriera białego szumu jest ... białym szumem!

Jaki jest związek między odchyleniem standardowym (1 w moim przykładzie) a wielkością i fazą?

Nie będzie związku między odchyleniem standardowym a fazą. Jeśli chodzi o wielkość, załóżmy, że n(t) jest stacjonarnym białym szumem ze średnią zerową i odchyleniem standardowymσ . Następnie autokorelacja (kowariancja) to:

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

Zatem gęstość widmowa mocy wynosi tylko σ2) (chociaż dla czasu dyskretnego nastąpi skalowanie na podstawie czasu trwania sygnału).


Pytania z komentarza:

  1. Kiedy mówisz, że transformata Fouriera jest również białym szumem, jak mogę zmierzyć std-dev, gdy transformacja jest złożona? Prawdziwa, wyimaginowana część czy kombinacja?

Załóżmy, że nasz szum jest czasem dyskretnym i wynosi n[m] (średnia zero, szum Gaussa, biały szum z wariancją σ2) ). Następnie transformacja jest:

N.[k]=m=0M.-1n[m]mi-jot2)πmk/M.=m=0M.-1n[m]sałata(2)πmk/M.)+jotn[m]grzech(2)πmk/M.)

a oczekiwana wartość to:

mi[N.[k]]=mi[m=0M.-1n[m]mi-jot2)πmk/M.]=m=0M.-1mi[n[m]]mi-jot2)πmk/M.=0

Wariancja części rzeczywistej jest dana przez:

mi[(N.[k])2)]=mi[m=0M.-1n[m]sałata(2)πmk/M.)p=0M.-1n[p]sałata(2)πpk/M.)]=mi[m=0M.-1p=0M.-1n[m]n[p]δ[n-p]sałata(2)πmk/M.)sałata(2)πpk/M.)]=m=0M.-1mi[n[m]2)]sałata2)(2)πmk/M.)=σ2)m=0M.-1sałata2)(2)πmk/M.)=σ2)(M.2)+sałata(M.+1)2)πk/M.grzech(2)πM.k/M.)2)grzech(2)πk/M.)   )=σ2)M.2)

Wierzę, że wyobrażona część będzie zachowywać się w ten sam sposób.

  1. Czy mógłbyś mi wyjaśnić, w jaki sposób czas trwania sygnału odnosi się do gęstości widmowej mocy (dla dyskretnych sytuacji czasowych)

Uważam, że (w oparciu o powyższe wyprowadzenie) gęstość widmowa mocy (oczekiwana wartość kwadratu DFT) będzie skalowana liniowo w zależności od czasu trwania.

  1. Jeśli std-dev nie wpływa na fazę, co decyduje o amplitudzie 3 stopni i rodzaju rozkładu (wydaje się raczej jednolita niż normalna)

Sprawdź tabelę na stronie 2 tego pliku PDF . mówi, że argument (faza) współczynników będzie równomiernie rozłożony, jak sam twierdzisz. Zrzut ekranu z tabeli poniżej.

wprowadź opis zdjęcia tutaj

Peter K.
źródło
3
W szczególności dwie koncepcje, które OP wprowadza w błąd, to gęstość widmowa mocy białego szumu i widmo mocy jednej konkretnej realizacji losowego procesu szumu białego.
Jason R
Dzięki! Mam kilka pytań uzupełniających. 1: Kiedy powiesz, że transformata Fouriera jest również białym szumem, jak mogę zmierzyć std-dev, gdy transformacja jest złożona? Prawdziwa, wyimaginowana część czy kombinacja? 2: Czy mógłbyś mi wyjaśnić, w jaki sposób czas trwania sygnału odnosi się do gęstości widmowej mocy (dla dyskretnych sytuacji czasowych) 3: Jeśli std-dev nie wpływa na fazę, co określa amplitudę 3 stopni i rodzaj dystrybucja (wydaje się raczej jednolita niż normalna)
Uffe
πσ2)M.2)
To jest bieżący link do dokumentu PDF wymienionego powyżej ( radarsp.weebly.com/uploads/2/1/4/7/21471216/dft_of_noise.pdf ), który jest uszkodzony.
Gość
@ Gość Dzięki! W przyszłości po prostu spróbuj edytować odpowiedź za pomocą nowego linku. Nie przejdzie bezpośrednio, ponieważ będzie musiał zostać sprawdzony przez użytkownika o większej liczbie powtórzeń, ale dostanie się tam (i otrzyma +2 powtórzenia).
Peter K.