STATLAB/Zadanie zaliczeniowe3
Spis treści
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.
