Laboratorium TSS cw4


ÿþPolitechnika BiaBostocka WydziaB Mechaniczny Instrukcja do zaj laboratoryjnych Temat wiczenia: Analiza czstotliwo[ciowa z wykorzystaniem transformaty Fouriera Numer wiczenia: 4 Laboratorium z przedmiotu: Teoria sygnaBów i systemów Kod: BiaBystok 2005 1. Wprowadzenie 2 3 Gsto[ widmowa mocy Na podstawie widma sygnaBu X(f) mo|emy wyznaczy gsto[ widmow mocy Pxx(f) okre[lon nastpujco: - gsto[ widmowa mocy wBasnej: 2 * X( f ) X( f )X ( f ) Pxx( f ) = lim = lim (5) T ’!" T ’!" T T lub inaczej przy N fragmentach sygnaBu o czasie trwania To: 2 2 N X ( f ) X ( f ) 1 1 Pxx( f ) = lim = (6) " i N’!" N To N To i=1 - gsto[ widmowa mocy wzajemnej: N 1 X( f )Y*( f ) 1 Xi( f )Yi( f ) Pxy( f ) = lim = (7) " To N’!" N To N i=1 W praktyce czsto korzysta si z realizacji gsto[ci widmowej mocy obserwujc pro- ces N-krotnie w czasie jego trwania To. Ilo[ fragmentów i ich transformat Xi(f), i = 1, ..., N nale|y dobra tak, aby ocena widma byBa stacjonarna ze wzgldu na N. Stosowanie okna pomiarowego Pomiar ka|dego sygnaBu odbywa si w skoDczonym czasie a zatem ma pocztek i ko- niec (równie| sygnaB wygenerowany programowo przez system komputerowy). Zadane tak warunki nie zapewniaj peBnej reprezentacji sygnaBu w czasie (sygnaB istnieje równie| poza wyznaczonymi momentem jego pocztku i koDca). Zatem mówi si o wycinku sygnaBu a nie o jego caBkowitej reprezentacji a co z tego wynika o istnieniu okna pomiarowego. Podczas 4 analizy takiego sygnaBu okno pomiarowe wprowadza znieksztaBcenia, ujawniajce si gBównie przy analizie czstotliwo[ciowej a przyczyn takiego stanu rzeczy jest fakt, |e niezale|nie od woli osoby analizujcej sygnaB, transformacji do przestrzeni czstotliwo[ci podlega iloczyn sygnaBu i okna, w którym dokonano pomiaru. Na rys. 1 pokazano sygnaB rzeczywisty Æ(t) oraz funkcj okna prostoktnego w(t) (w(t) = 1 dla  T/2 d" t d" T/2 , w(t) = 0 dla t <  T/2 i t > T/2). Rys. 1. Reprezentacja graficzna sygnaBu rzeczywistego w prostoktnym oknie pomiarowym Transformata Fouriera okna prostoktnego opisuje posta (8): sin(ÀfT ) W ( f ) = T (8) ÀfT lub te| sin(À ( f - fi )T ) W ( f ) = T (9) À ( f - fi )T gdzie: T  dBugo[ okna, f  czstotliwo[ bie|ca, fi  dyskretna warto[ czstotliwo[ci [rodkowej transformaty okna (dowolne warto[ci poBo|enia pr|ków w badanym widmie). Ró|nica (f  fi) wyra|a odstp czstotliwo[ci bie|cej od [rodkowej. WpByw okna po- miarowego ujawnia si gBównie w miejscach wystpowania pr|ków czstotliwo[ciowych w widmie badanego sygnaBu. Reprezentacj graficzn funkcji W(f) pokazano na rys. 2. 5 Rys. 2. Reprezentacja graficzna transformaty Fouriera prostoktnego okna pomiarowego W wyniku uwzgldnienia okna pomiarowego w(t), na podstawie pierwotnego sygnaBu Æ(t), otrzymuje si sygnaB f(t): f (t) = Õ(t) Å" w(t) . (10) Wówczas transformata F(f) okre[lona jest zale|no[ci: " F(t) = (11) +"Õ(t) Å" w(t) Å" e-2Àftdt . -" Zapisujc mno|enie splotowe w postaci analitycznej i podstawiajc (9) otrzymuje si: " T sin(À ( f - fi )T ) F( f ) = ¦ ( f ) df (12) +" À ( f - fi )T -" a zatem: F( f ) `"¦ ( f ) (13) PrzykBadem powstania znieksztaBcenia widma wywoBanym istnieniem okna jest fakt, |e przy zmniejszaniu si odlegBo[ci pomidzy czstotliwo[ciami analizowanego sygnaBu f1 i f2 wykresy funkcji W(f)  rozpite na tych czstotliwo[ciach zaczn si sumowa. Wystpuj wówczas znieksztaBcenia jako[ciowe (jedna skBadowa zamiast dwóch) i ilo[ciowe (wiksza amplituda). Oba bBdy wynikaj z utraty rozró|nialno[ci skBadowych widmowych. W celu polepszenia rozró|nialno[ci i zmniejszenia zakBóceD widma wynikajcych z obcicia sygnaBu przez okno pomiarowe stosuje si ró|nego rodzaju okna korekcyjne. Nieza- le|nie od typu zastosowanego okna korekcyjnego, osignita w ten sposób modyfikacja pole- ga na zBagodzeniu pionowych obci sygnaBu na pocztku i na koDcu pomiaru. Wyró|nia si trygonometryczne, potgowe i wykBadnicze okna korekcyjne. Najcz[ciej wykorzystywane s okna trygonometryczne a w[ród nich najbardziej znane s okna, które przyjBy nazwy od na- zwisk ich twórców: - okno von Hanna (okno Hanninga): 6 2Àt T wHn (t) = 0,5 + 0,5 Å" cos( ) dla t d" T 2 i T wHn (t) = 0 dla t > 2 - okno Hamminga: 2Àt wHm (t) = a + (1- a) Å" cos( ) dla t d" T T i , gdzie: a "< 0,5; 1> wHm (t) = 0 dla t > T Okno Hanninga jest najbardziej rozpowszechnionym oknem korekcyjnym i w syste- mach obróbki sygnaBów bywa obligatoryjnie umieszczane w oprogramowaniu oraz czsto w cz[ci sprztowej. Na rys. 3 przedstawiono porównanie transformaty prostoktnego okna pomiarowego W(f) i transformaty korekcyjnego okna Hanninga WHn(f). Oczekuje si tutaj, |e przy zastoso- waniu okna korekcyjnego nastpi znaczne zmniejszenie znieksztaBcenia widma, wynikajce z istnienia listków bocznych. Rys. 3. Porównanie transformaty okna prostoktnego i okna Hanninga Funkcje programu MATLAB do wykorzystania w wiczeniu: Funkcja Przeznaczenie fft Dyskretna transformata Fouriera abs Obliczenie moduBu (warto[ci bezwzgldnej) psd Estymacja gsto[ci widmowej mocy metod Welcha hann Okno Hanna (Hanning) hamming Okno Hamminga (Hamming) rectwin Okno prostoktne triang Okno trójktne blackman Okno typu Blackman 7 2. Wykonanie wiczenia a) wygenerowa sygnaBy x1(t), x2(t) i x3(t) opisane nastpujco: x1( t ) = A11 Å" cos( 2 Å"À Å" f11 Å" t ), x2( t ) = A21 Å" cos( 2 Å"À Å" f21 Å" t ) + A22 Å" cos( 2 Å"À Å" f22 Å" t ), x3( t ) = A31 Å" cos( 2 Å"À Å" f31 Å" t ) + A32 Å" cos( 2 Å"À Å" f32 Å" t ) + A33 Å" cos( 2 Å"À Å" f33 Å" t ); gdzie: f11, f21, f22, f31, f32, f33  czstotliwo[ci kolejnych skBadowych, A11, A21, A22, A31, A32, A33  amplitudy kolejnych skBadowych, b) wyznaczy widma X1(f), X2(f) i X3(f) (funkcja fft(..)), c) wykre[li widma amplitudowe | X1(f)|, | X2)| i X3f)|, d) dokona skalowania osi czstotliwo[ci zgodnie z zale|no[ci: i fi = Å" f , gdzie: i = 0 .. N-1 (N  ilo[ danych wektora sygnaBu), p N fp  czstotliwo[ próbkowania. e) dokona skalowania osi amplitudy zgodnie z zale|no[ci (4), e) wykre[li widma amplitudowe wyskalowane w warto[ciach amplitudy i czstotliwo[ci, f) wyznaczy estymacj gsto[ci widmowej mocy sygnaBu x1(t), x2(t) i x3(t) metod Welcha (funkcja psd(..)), g) wyznaczy gsto[ widmow mocy sygnaBu x1(t), x2(t) i x3(t) przy wykorzystaniu widma X1(f), X2(f) i X3(f), h) wyznaczy widma sygnaBu x1(t), x2(t) i x3(t) z zastosowaniem okna Hanninga, Hamminga, okna typu Blackman, okna prostoktnego i trójktnego, i) wykre[li widma amplitudowe sygnaBów z zastosowaniem okien, j) porówna widmo amplitudowe sygnaBu x1(t) z oknem prostoktnym o czasie trwania bd- cym caBkowit wielokrotno[ci okresu T1=1/f11 z widmem amplitudowym tego samego sygnaBu z oknem prostoktnym o czasie trwania nie bdcym caBkowit wielokrotno[ci okresu T1, k) porówna widmo amplitudowe sygnaBu x1(t) z oknem prostoktnym o czasie trwania nie bdcym caBkowit wielokrotno[ci okresu T1 z widmem amplitudowym z zastosowaniem okna Hanninga i Hamminga o takiej samej dBugo[ci jak przy oknie prostoktnym. 8 Literatura 1. Szabatin J.: Podstawy teorii sygnaBów. WKA, Warszawa 2003. 2. Pasko M., Walczak J.: Teoria sygnaBów. Wydaw. Politechniki Zlskiej, Gliwice 1999. 3. Izydorczyk J., PBonka G.: teoria sygnaBów. Wydaw. HELION, Gliwice 1999. 4. Kurowski W.: Podstawy teoretyczne komputerowego miernictwa systemów mecha- nicznych. Wydaw. Politechniki BiaBostockiej, BiaBystok 1994. 5. Smyczek J.: Teoria sygnaBów i informacji. Cz. I., Wydaw. Politechniki Aódzkiej, Aódz 1991. 6. Wojnar A.: Teoria sygnaBów. WNT, Warszawa 1988. 9

Wyszukiwarka

Podobne podstrony:
Laboratorium TSS cw5
Laboratorium TSS cw1
Rola laboratoriów w świetle wymagań systemów zarządzania jakoscią
Laboratorium 3
O cw4
Ćwiczenie laboratoryjne nr 6 materiały
Windows 2 Laboratorium 4b
Chemia żywnosciCwiczenie laboratoryjne nr 1 wyodrebnianie i badanie własciwosci fizykochemicznych b
Laboratorium 3
LABORATORIUM CHEMIA I WYTRZYMALOSC MATERIALOW sprawko 1
Ćw4 Instrukcje CASE i IF
Ustawa o medycznej diagnostyce laboratoryjnej

więcej podobnych podstron