Построение частотной хар-ки CIC Опции

[Jools]

Привет всем.

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

H(f) = (abs(sin(pi*M*f) / sin(pi*f/R)))^2N, где
R - децимирующий фактор
N - порядок
M - дифференциальная задержка

Но получить по нему АЧХ не получается. Может направите меня в правильное русло? Только, чур, Matlab типа hcic = mfilt.cicdecim(R,M,N) не предлагать.

Заранее спасибо.

[hobgoblin]

H(f) = (abs(sin(pi*M*f) / sin(pi*f/R)))^2N, где
R - децимирующий фактор
N - порядок
M - дифференциальная задержка

Это Вы из статьи Хогенауэра взяли? Там вроде бы надо брать f не от -0.5 до +0.5, а от -R*0.5 до +R*0.5. Поскольку сказано, что

f is the frequency relative to the low sampling rate fs/R

Тогда все должно получиться