Ćw2 pl, AGH Mechatronika, mgr sem.1, identyfikacja układów mechatronicznych
[ Pobierz całość w formacie PDF ]
//-->Identyfikacja układów mechatronicznychProjekt #23. PróbkowanieW tym zadaniu będziesz przerzedzał próbkowany sygnał np. obniżając próbkowanie oryginalnego sygnału wsygnał niedopróbkowany. Przerzedzanie wykonaj poprzez wyciągnięcie do drugiej wartości z oryginalnegosygnału. Gdy sygnały są dyskretne czas pomiędzy dwoma kolejnymi pomiarami jest nazwany „okresempróbkowania” a jego odwrotność „częstotliwością próbkowania”.a)Stwórz sygnał:x = 2*cos(2*pi*0.1*k) + 2*cos(2*pi*0.4*k);dlak=1:Ns; Ns = 5000;Zawiera on dwa harmoniczne komponenty z różnymi częstotliwościami.Oblicz transformatę Fouriera dla tego sygnałuh = fft(x);Wykreśl wartość bezwzględnąhw skali powiązanej z częstotliwością próbkowaniaFs = 1kHztak by na osiczęstotliwości zawrzeć wartości między 0 i częstotliwość Nyquista.figure(1), plot(f,abs(h(1:Ns/2)))Przerzedź sygnał poprzez stworzenie wektoraxdzawierającego co drugi element wektora x.xd = x(1:2:end);Określ widmo sygnałuxdużywając poleceniaffti zaprezentuj wyniki dla połowicznej częstotliwościpróbkowaniaFs = 500Hznafigure(2)w taki sam sposób jak wcześniej (zauważ iż częstotliwośćNyquista oraz liczba próbek jest zredukowana więc powinieneś stworzyć nowy wektor częstotliwościfd).Porównaj widma z obu figur i wyciągnij wnioski.b)Powtórzcałe zadanie a) dla sygnału:x = 2*cos(2*pi*0.1*k) - 2*cos(2*pi*0.4*k);dla k=1:Ns; Ns = 5000;Zaobserwuj zmiany w widmie na figurze 2.Sprawozdanie powinno zawierać:Figury 1 i 2 z widmami dla przypadków a i b.Skomentuj wyniki: porównaj różnice w widmach. Które częstotliwości zostały pominięte poprzerzedzeniu sygnału i dlaczego?Dlaczego widmo z przerzedzonym sygnałem wygląda tak jak wygląda? (np. które częstotliwościsą stłumione i do jakiego poziomu? Co stało się z amplitudami?)
[ Pobierz całość w formacie PDF ]