Нужно рассчитать БИХ фильтр на определенную частоту. Делаю так:
Нули помещаю в точки 1,i0; -1,i0; Полюса в точки R*exp(i*Q) и R*exp(-i*Q).

где Q = 2pi*(f/fs) f- нужная частота, fs-частота выборок.
R = 1-fp/fs*pi fp - полоса пропускания

Получаю передаточную функцию (1-(Z^-2))/(1-(Z^-1)*2*R*cos(Q)+(Z^-2)*R^2)

и с помощью маткада коэффициент усиления на резонансной частоте
K = 2/(1-R)*SQR(1-cos(Q)^2)

Вопрос - что неправильно?