Skrypt APSC - MARS

ROZDZIAŁ 7. FILTRY CYFROWE 139Ts=1/ f s ; % obliczamy okres próbkowaniawga = 2/Ts∗ tan ((2∗ pi ∗ f g / f s ) / 2 ) ; % p u l s a c j a graniczna f . analogowegof ga = wga /(2∗ pi ) ; % c z ę s t o t l i w o ś ć graniczna f . analog .wa = 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% p r o j e k t f i l t r u analogowego ( Czebyszewa , 2−go rzędu )[ ba , aa]=cheby1 (2 ,R, wga , ’ s ’ ) ;% charakter ystyka c z ę s t o t l i w o ś c i o w a f . analogowegof i g u r e (1)[ Ha ,w]= f r e q s ( ba , aa , wa ) ;hold onp l o t (wa/(2∗ pi ) ,20∗ log10 ( abs (Ha ) ) ) ;p r i n t s y s ( ba , aa , ’ s ’ ) % postać t r a n s m i t a n c j i ( zmienna s )a=ba ( 3 ) ;b=aa ( 2 ) ;c=aa ( 3 ) ;% podstawienie wg przykładuTS2 = (2/ Ts ) ^ 2 ; % podstawienie upraszczając e% o b l i c z e n i e współczynników f i l t r u cyfrowegob0 = a /(TS2+b∗2/Ts+c ) ;b1 = 2∗a /(TS2+b∗2/Ts+c ) ;b2 = a /(TS2+b∗2/Ts+c ) ;a0 = 1 ;a1 = (2∗ c−2∗TS2 )/( TS2+b∗2/Ts+c ) ;a2 = (TS2−b∗2/Ts+c )/( TS2+b∗2/Ts+c ) ;% wektor współczynników ( l i c z n i k t r a n s m i t a n c j i z )bd = [ b0 b1 b2 ] ;% wektor współczynników ( mianownik t r a n s m i t a n c j i z )ad = [ a0 a1 a2 ] ;% postać t r a n s m i t a n c j i cyfrowej ( zmienna z )p r i n t s y s (bd , ad , ’ z ’ )% charakter ystyka c z ę s t o t l i w o ś c i o w a f . 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 on