Как вычисляется общая АЧХ CIC и КИХ в связке, MATLAB/Octave Опции

[honeycomb0]

Вот у меня есть CIC фильтр после которого идет КИХ. Понятно, что для общего АЧХ этой связки надо перемножить АЧХ фильтров. С АЧХ КИХ фильтра вопросов нет. А вот с АЧХ CIC у меня сплошная путаница. Сам СIC у меня очень простой (т.е. первого порядка, четырехкратная децимация) и описан он очень примитивно:

CODE

R = 4; % четырехкратная децимация
D = 2; % задержка
delayBuffer = zeros(1,D);
Acum = 0;
xn = zeros(1,128); xn(1) = 1; % входной сигнал, Дельта функция
output = zeros(1,(length(xn) / R));

decim = 0;
n = 1;

% Сначала интегрируем
for ii=1:length(xn)
Acum = Acum + xn(ii);
decim = decim + 1;

% ... потом децимируем
if(decim == R)
% ... тут гребеньчатый фильтр
Comb = Acum - delayBuffer(end);
delayBuffer(2:end) = delayBuffer(1:end-1);
delayBuffer(1) = Acum;
output(n) = Comb;
n = n + 1;
decim = 0;
end
end

Для АЧХ этого фильтра я делаю freqz(output) и получаю (см. первый график). И тут у меня начались "проблемы в понимании". Неужели АЧХ не должно выглядить как-то вот так (см. второй график)? Я догадываюсь что АЧХ на выходе "свернулось" в один лепесток из-за прореживания дециматора - это так? Как тогда понимать, что при расчете CIC фильтра, учитывают уровни подавления последующих лепестков (за счет увеличения порядка фильтра), но при этом на выходе фильтра имеем АХЧ с одним лепестком...? А так же, уровень усиления в обоих случаях разнится Очевидно, я что-то не так понял в CIC...

Сообщение отредактировал honeycomb0 - Sep 8 2016, 10:22

Эскизы прикрепленных изображений