Cyfrowe Przetwarzanie Sygnałów
2010.03.08
Analiza widmowa
05
dr inż. Jarosław Bułat
Ćwiczenie 1.
W poniższym programie wykonano analizę częstotliwościową różnymi metodami przebiegu
opisanego następującą zależnością:
x śąt źą=cosśą2 Ćą f tźąƒÄ…0.5cosśą2Ćą f tźą , gdzie : f =100Hz , f =102Hz
1 2 1 2
" Która metoda jest właściwa?
" Powtórz zadanie dla N = 1000 oraz zaprezentuj widmo sygnału dla obu przypadków na
jednym wykresie (oś X wykresu powinna być wyskalowana w częstotliwości unormowanej,
oÅ› Y w dB)
" Oblicz wartość bezwzględną współczynników transformacji odpowiadających częstotliwościom
f i f i porównaj z wartościami odczytanymi z wykresów.
1 2
" dla f =100Hz, f =102.5Hz i N=1000 wykonaj analizę poniższym programem, wytłumacz
1 2
dlaczego wartość maksymalna współczynników odpowiadających f zmniejszyła się chociaż
2
energia składowej sygnału f pozostała bez zmian (nieznacznie zmieniła się częstotliwość)
2
clear all; close all;
fs = 1000; Nx = 1000;
A1 = 1;
A2 = 0.5;
f1 = 100;
f2 = 102;
t = 1/fs*(0:Nx-1);
x = A1*cos( 2*pi*f1*t ) + A2*cos( 2*pi*f2*t );
figure, freqz( x ); % freqz
N = 500; % FFT
X = fft( x, N );
figure, plot( 20*log10(abs(X(1:end/2+1))) );
fvtool( x ); % fvtool
Ćwiczenie 2.
Dla sygnału z ćw. 1 o parametrach f =100Hz, f =102.5Hz wykonaj analizę częstotliwościową z przy
1 2
użyciu transformacji Fouriera interpolując widmo poprzez dodanie zer do sygnału x, w taki sposób
aby cały wektor miał długość 10x większą od oryginalnego. Następnie oblicz transformację Fouriera
DTFT typu ,,lupa'' (patrz ćwiczenie 3 z laboratorium 3) w zakresie f =90Hz, f =120Hz dla 1024
d g
współczynników. Oba wykresy wraz z widmem prostej transformacji FFT (ćw. 1) narysuj na jednym
wykresie (oÅ› X wyskalowana w hercach, oÅ› Y w dB).
Powtórz ćwiczenie stosując różne okna czasowe, tzn. nie analizuj interesujący nas sygnał tylko jego
iloczyn z wybraną funkcją okna (np.: Hanning, Hamming, Blackman, Kaiser, Bartlett, etc...). Zmień
amplitudę drugiej składowej sygnału na 0.0001. Dobierz funkcję okna oraz liczbę próbek sygnału,
zapewniające wymaganą w tym przypadku rozdzielczość częstotliwościową i amplitudową. Obie
składowe mają być w widmie widoczne.
Ćwiczenie 3.
Wyznacz widmo sygnału pseudolosowego n z rozdzielczością 1Hz dla przedziału 0-128Hz (f =1kHz).
s
clear all; close all;
Nn = 10000;
randn( 'state', 0 ); % inicjalizacja generatora liczb pseudolosowych
n = 0.3*randn( Nn, 1 );
Jakie powinno być widmo sygnału pseudolosowego o rozkładzie normalnym? Dlaczego widmo
takiego sygnału obliczone z jednego okna nie jest takie jakie być powinno? (Jest różne od
teoretycznego). Co należy zrobić aby przybliżyć go do widma teoretycznego? Zastosuj funkcję
spectrum() do wyznaczenia widma tego sygnału. Co ona liczy?
Dodaj ten szum do sygnału z ćw. 1: N=1000, f =100Hz, f =102Hz. Oblicz fft() i spectrum() dla
1 2
sygnału sumy. Jakie wnioski można wyciągnąć z tej analizy?
Ćwiczenie 4.
Wygeneruj sygnał spróbkowany częstotliwością f =1000Hz o czasie trwania t=1s opisany
s
następującą zależnością:
x śąt źą=cosśą2 Ćąśą f tƒÄ…0.5 df t2źąźąƒÄ…cosśą2 Ćąśąśą f ƒÄ… f źątƒÄ…0.5 df t2źąźą , f =100Hz , f =20Hz , df =100Hz
0 0
Następnie wykonaj analizę częstotliwościową jak w ćw. 1 oraz analizę czasowo-częstotliwościową za
pomocą funkcji spectrogram(...). Dobierz długość oraz rodzaj okna analizy w ten sposób aby oba
przebiegi o zmiennej częstotliwości były rozróżnialne. Kosztem czego się ją osiąga? Powtórz analizę
z wykorzystaniem funkcji spectrum(). Poniżej został przedstawiony fragment programu generujący
niewłaściwe wykresy - parametry okna zostały tak dobrane, że nie pozwalają na rozróżnienie obu
składowych sinusoidalnych.
[S,F,T,P] = spectrogram( x, boxcar(64), 60, [0:500], fs );
figure; imagesc( T, F, P ); colormap( jet ); colorbar;
xlabel( 't [s]' ); ylabel( 'f [Hz]' );
Ćwiczenie 5.
Dla podanego przebiegu, wykonaj analizę częstotliwościową i podaj parametry każdej ze składowych
sygnałów z których został złożony cały sygnał:
" czas rozpoczęcia i zakończenia każdej składowej
" częstotliwość początkową i końcową każdej składowej
" amplituda (w dziedzinie czasu) na początku i końcu każdej składowej
" poziom szumu w dB
Załóż, że parametry każdej składowej sygnału takiej jak: amplituda, częstotliwość zmieniają się
liniowo w czasie. Przyjmij częstotliwość próbkowania jako: f =1000Hz.
s