Построение частотной хар-ки 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) не предлагать.

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

[petrov]

От реализации CIC можно абстрагироваться. Построить АЧХ обычного FIR фильтра можете? АЧХ FIR фильтра с прямоугольной импульсной характеристикой? АЧХ последовательного соединения таких фильтров? Затем понимаете как при децимации наложение спектров происходит?

[Jools]

Вот вроде разобрался. Если кому интересно:

clear all;

Fs = 40000;

M = 1; % Differencial delay
N = 2; % Number of section
R = 8; % Decimation factor
f = 0:1/Fs:0.5; % Normalize frequency axis

H = (sin(pi*f*R*M)./sin(pi*f))/R;
Hk = H.^N;

plot(f,20*log10(abs(Hk)));

АЧХ последовательного соединения таких фильтров?

Если на одной частоте дискретизации, то перемножением АЧХ. Я прав?

[petrov]

Если на одной частоте дискретизации, то перемножением АЧХ. Я прав?

Да конечно.