3500000 3000000 2500000 2000000 1500000 1000000 500000 О 0 25 50 75 100 125 150 175 200 Рис. 2.3. Спектр фрагмента одно и то же небольшое значение. Это означает, что применив фильтр высоких частот, мы подавим основную мощность сигнала. Найдем коэффициенты фильтра с помощью функции finvin. Длина фильтра М = 128. Определим параметр, задающий частоту отсечения. В данном примере частота стробирования Fr = 22050, а длина интервала, по которому найдены коэффициенты, равна N = 3*М = 381. Согласно (1.20). коэффициенту Фурье с номером 100 отвечает частота freq = Fr * 100/N = 5742. Воспользовавшись формулой (1.42). найдем параметр с = 2 * f req /Fr = 0.521, отвечающий этой частоте. Теперь для конструирования фильтра можно выбрать любое значение из интервала [0.521,1) Выберем значение равным 0.8 и найдем коэффициенты фильтра В= f i r w i n ( 1 2 7 , 0 . 8 , zero_pa s s=Fa l se ) Чтобы посмотреть на передаточную функцию построенного фильтра, используем функцию scipy.signal.freqz. В = sgn . f i rw in (M , 0 . 8 ,p a s s _ z e r o=F a l s e ) (W,H) = s g n . f r e q z ( B , 1) p i t . p l o t (F r / 2 *W / n p . p i , ab s (H ) , ’k ’ )
RkJQdWJsaXNoZXIy MTExODQxMg==