STATLAB/Zadanie zaliczeniowe3

1. Przygotowanie danych

Nagraj przy użyciu mikrofonu kilka następujących po sobie dźwięków o różnej wysokości (np. cztery) zaśpiewanych lub zagranych na dowolnym instrumencie muzycznym. Dźwięk zapisz w formacie „raw”; zanotuj częstotliwość próbkowania.

2. Analiza czasowa widma

Dla wybranego rodzaju okna czasowego i ustalonej szerokości okna [math]N_w \ll N[/math] ([math]N[/math] = długość sygnału), oblicz krótkoczasową transformatę Fouriera (STFT), czyli taką dwu-wymiarową tablicę (macierz numpy)

[math] Z[k][n_0] = \frac{\left\vert X_{n_0}[k] \right\vert^2}{N_w} \,, [/math]

że [math]X_{n_0}[/math] jest dyskretną transformatą Fouriera okienkowanego wycinka sygnału od [math]n_0-\frac{N_w}{2}[/math] o szerokości [math]N_w[/math] (czyli wyśrodkowanego w [math]n_0[/math]).

Proszę pamiętać o okienkowaniu (periodogram) wycinka sygnału przed policzeniem DTF.

3. Prezentacja

Macierz wykreśl przy pomocy funkcji pylab.imshow, podając zakresy poszczególnych osi przy użyciu parametru

[math] \textrm{extent} = (t_{min}, t_{max}, f_{min}, f_{max}) \,. [/math]

4. Interpretacja wyniku

Powtórz punkty 2 i 3 dla kilku różnych szerokości okien [math]N_w[/math] i zinterpretuj widoczne różnice.