Анатолий, здесь все очень просто.
АЧХ децимируюещй системы (DDC в наших случаях) мы смотрим на частоте дискретизации АЦП.
С первым фильтром fir1 все понятно, он работает до децимации на 2.
Чтобы привести частоте дискретизации АЦП АЧХ второго фильтра fir2, который стоит после децимации на 2, мы должны его импульсную характристику h2 интероплировать на 2.
После этого берем свертку от h1 и интерполированной h2, и затем АЧХ от результата свертки.

Общее правило такое, если перед каким-то фильтром есть децимация на N, то для приведения его к частоте дискретизации АЦП, надо его импульсную характеристику интерполировать на N.
И так с каждым каскадом DDC. Потом свертка всех интероплированных результатов.


h1 = [-1 4 10 4 -1] ./ 16;
h2 = [0.00142192172987358 -0.00222268508122243 -0.0137433922997260 -0.0322604508637608 -0.0458124362883413 -0.0346645776223707 0.0163109682969986 0.102976693990851 0.197188326792324 0.258848499145519 0.258848499145519 0.197188326792324 0.102976693990851 0.0163109682969986 -0.0346645776223707 -0.0458124362883413 -0.0322604508637608 -0.0137433922997260 -0.00222268508122243 0.00142192172987358];

h_stage1 = conv(h1, h1); % Каскад 1 (без децимации

h_stage2 = h2;
h_stage2 = upsample(h_stage2, 2); % С учетом децимации на 2.

% Результирующая импульсная характеристика.
h_ddc = conv(h_stage1, h_stage2);
fvtool(h_ddc);

Большое спасибо за квалифицированный ответ! А я пытался наоборот - децимировать h1 и естественно ничего не получалось. Поехал дальше.

p.s.
А вот эта строчка зачем?

Зачем сворачивать самого с собой?

h=conv(h1,8*upsample(h2,8));

дальше freqz

Вот это ответ! Проще не бывает. Спасибо!

В моем примере стояло два фильтра h1 перед децимацией на 2, забыл это упомянуть.


Децимировать h надо если перед фильтром есть интеропляция, например когда делаете DUC.

Спасибо! Я понял свою ошибку.

Вот такая сквозная АЧХ моего DDC получилась. Спасибо всем участникам обсуждения!

Скажите,пожалуйста, можно ли это сделать в частотной области?

Можно.

Вот что-то читаю и никак не могу понять...


Полосу сигнала в 50МГц при частоте дискретизации 75МSPS??? Где-то ошибка?


Тут имелось ввиду "задержка" а не порядок? Порядок фильтра обеспечивает подавление лепестков, насколько я это понимаю (http://www.dsplib.ru/content/cic/cic.html). Т.е. кол-во нулей при разных порядках фильтра остается тем же, отличается только амплитуды "лепестков" в АЧХ. Разве не так?

Первый вопрос снят... Наверняка имелось ввиду полоса с центром на 0Гц.