Привет всем.

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

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

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

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

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

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

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

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

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)));




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

Да конечно.