Просто реализуйте функцию фильтра, как описано здесь . Это похоже на простой фильтр линии задержки, где b
и a
- коэффициенты, а x
- входные данные:
a[1]y[n] = b[1]x[n]+b[2]x[n-1]+...+b[N]x[n-B+1]-a[2]y[n-1]-...-a[N]y[n-N+1]
Вам нужно только позаботиться о еще не существующих индексах, которые обозначены в квадратных скобках. a
и b
- векторы длины N
.