Skrypt APSC - MARS

ROZDZIAŁ 7. FILTRY CYFROWE 138wa = 2∗ pi ∗ [ 0 : f s /1000: f s / 2 ] ; % wektor p u l s a c j i analogowejwd = ( 0 : 1 / f s : 1 ) ; % wektor p u l s a c j i cyfrowej[ ba , aa]=cheby1 (N,R, wg , ’ s ’ ) ; % p r o j e k t f i l t r uf i g u r e (1)c l fhold on% odpowiedź impulsowa f i l t r u analogowegoimpulse ( ba , aa ,TMAX) ;% p r o j e k t f i l t r u cyfrowego[ bd , ad ] = impinvar ( ba , aa , f s ) ;g = f s ; % skalowanie amplitudy% odpowiedź impulsowa f i l t r u cyfrowegoimpz ( g∗bd , ad , [ ] , f s ) ;legend ( [ ’ h( t ) ’ ; ’ h(n ) ’ ] ) ;hold o f ff i g u r e (2)c l fhold on% charakter ystyka amplitudowo−c z ę s t o t l i w o ś c i o w a% f i l t r u analogowego[ Ha ,w]= f r e q s ( ba , aa , wa ) ;p l o t (wa/(2∗ pi ) ,20∗ log10 ( abs (Ha ) ) ) ;% i f i l t r u cyfrowego[ Hd, wd]= f r e q z ( bd , ad , 5 1 2 ) ;p l o t (wd/(2∗ pi )∗ f s , 20∗ log10 ( abs (Hd ) ) , ’ : ’ ) ;x l a b e l ( ’ f [ Hz ] ’ ) ;y l a b e l ( ’ |H| [ dB ] ’ ) ;legend ( [ ’H( s ) ’ ; ’H( z ) ’ ] ) ;a x i s ( [ 0 500 −40 5 ] ) ;g r i d on7.7.2 Metoda transformacji biliniowejMatlab daje kilka możliwości zaprojektowania filtrów tą metodą. Oto skrypt pokazującymożliwie dużo etapów obliczeniowych, czyli obliczający współczynnikifiltru IIR „na piechotę”:f s = 1000; % c z ę s t o t l i w o ś ć próbkowaniaf g = 100; % c z ę s t o t l i w o ś ć graniczna f i l t r uR = 1 ;% zafalowania , dB