Skrypt APSC - MARS

ROZDZIAŁ 7. FILTRY CYFROWE 146nych o częstotliwościach 10Hz i 131Hz i o różnych amplitudach. Tym razemzaprojektowaliśmy filtr IIR o charakterystyce Butterwortha, piątego rzędu, oczęstotliwości granicznej 50Hz. Filtracja odbywa się za pomocą funkcji filter.Na koniec, za pomocą transformaty DFT, sporządzany jest wykres widma sygnałuwejściowego i przefiltrowanego, na którym widać, że składowa sygnału oczęstotliwości 131Hz została wytłumiona.f s =1000; % c z ę s t o t l i w o ś ć próbkowaniaL=250; % l i c z b a próbekA1=5; % amplituda 1fx1 =10; % c z ę s t o t l i w o ś ć 1A2=3; % amplituda 2fx2 =131; % c z ę s t o t l i w o ś ć 2N=5; % rząd f i l t r uf g =50; % c z ę s t o t l i w o ś ć graniczna f i l t r uwg = 2∗ f g / f s ; % p r z e l i c z e n i e na c z ę s t o t l i w o ś ć znormalizowanądt=1/ f s ; % okres próbkowaniat=dt ∗ ( 0 : L−1); % wektor czasu% złożony sygnał testowy :x=A1∗ s i n (2∗ pi ∗ fx1 ∗ t)+A2∗ s i n (2∗ pi ∗ fx2 ∗ t ) ;f i g u r e (1)p l o t ( t , x ) ; g r i d ; y l a b e l ( ’ x ( t ) ’ ) ; x l a b e l ( ’ t [ s ] ’ ) ;t i t l e ( ’ Sygnał wejściowy x ( t ) ’ )[ b , a]= butter (N, wg ) ; % p r o j e k t f i l t r uy=f i l t e r (b , a , x ) ; % f i l t r a c j a% sygnał po f i l t r a c j i :f i g u r e (2)p l o t ( t , y ) ; g r i d ; y l a b e l ( ’ y ( t ) ’ ) ; x l a b e l ( ’ t [ s ] ’ ) ;t i t l e ( ’ Sygnał p r z e f i l t r o w a n y y ( t ) ’ )% DFT powyższych sygnałów :NFFT = 2^nextpow2 (L ) ; % Długość 2^xX = f f t ( x ,NFFT)/L ;Y = f f t ( y ,NFFT)/L ;f = f s /2∗ l i n s p a c e (0 , 1 ,NFFT/2+1);% Rysowanie spektrum c z ę s t o t l i w o ś c i :f i g u r e (3)p l o t ( f , 2∗ abs (X( 1 :NFFT/2+1)) , ’ r −−’, f , 2∗ abs (Y( 1 :NFFT/2+1)) , ’b− ’)t i t l e ( ’ Spektrum c z ę s t o t l i w o ś c i o w e sygnałów x ( t ) i y ( t ) ’ )x l a b e l ( ’ f [ Hz ] ’ )y l a b e l ( ’ Amplituda ’ )legend ( ’ |X( f ) | ’ , ’ |Y( f ) | ’ )g r i d on