полосовой фильтр не уважая отсечку

Question

Я использую этот фильтр в Python:

def bandpass_firwin(ntaps, lowcut, highcut, fs, window='hamming'): 
    nyq = 0.5 * fs 
    taps = firwin(ntaps, [lowcut, highcut], nyq=nyq, pass_zero=False, 
        window=window, scale=False) 

где мои ntaps = 128; lowcut = 0.7; highcut = 4; fs = 61

Я отфильтровываю свой сигнал, который имеет 610 образцов, отобранных на 61 Гц (так что это 10 секунд).

Когда я смотрю на спектр сигнала, который был отфильтрован с помощью полосового фильтра этого, я вижу это:

РЕЕК в этом спектре при 0,61 Гц. Коэффициент не находится в диапазоне от 0,7 до 4 Гц.

Как это возможно? & Как я могу предотвратить это?

Answer 1

Ваш фильтр не является магии - существуют внутренние ограничения на пропускную способность. Попробуйте использовать больше кранов, если вам действительно нужно плотное обрезание.

Чем больше кранов вы используете, тем больше вам нужно подумать о краевых эффектах и ​​о том, как вы их обрабатываете (поскольку предположения края все больше и больше проникают в блок данных). Возможно, вы хотите, чтобы плавный спуск по краю? Или зеркало и повторение данных? Возможно, вы можете полностью игнорировать его ...

Другой метод свертки заключается в фильтрации непосредственно в частотной области путем простого умножения на требуемый спектр. Это накладывает предположение о том, что ваш сигнал повторяется, хотя вы можете изменить это, расширив свой сигнал по своему усмотрению. Если вы хотите знать поддержку эквивалентного FIR-фильтра, возьмите IFFT окна, и вы увидите, насколько начало временного блока будет мазать в начале.