from m a t p l o t l i b import pyp lo t as p i t M = 2000 LimPrint = 40 Arg = range(LimPrint) ArgF = np.arange(M/5) * 2 . * n p .p i * 5 0 . /N Sn = np . s in(ArgF) В = s g n . f i r w i n ( l l , 0 . 9 ) FSig = s g n . I f i l t e r ( В , 1 . ,Sn) p i t . p l o t (A r g ,F s i g [ : L imPr in t ] , ’k : ’ , Arg ,Sn[ : L imP r in t ] , ’k ’ ) Соответствующий график представлен на рис. 3.2. Заметим, что длина FIR фильтра равна 2 * 5 + 1, поэтому фазовый сдвиг произошел на 5 отсчетов. Горизонтальный участок в начале отфильтрованного сигнала возникает из-за нулевых начальных значений сигнала, используемых функцией lfilter. Рис. 3.2. Фазовый сдвиг после фильтрации сигнала. Сплошная линия — исходный сигнал, двоеточия — сигнал после фильтрации
RkJQdWJsaXNoZXIy MTExODQxMg==