Реализация нерекурсивного фильтра CIC Опции

[soldat_shveyk]

Скажите, а как правильно посмотреть суммарную АЧХ такой цепочки : fir1-децимация-fir2?

Анатолий, здесь все очень просто.
АЧХ децимируюещй системы (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);

[_Anatoliy]

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

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

Код

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

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

[thermit]

Потихоньку двигаюсь дальше. Скажите, а как правильно посмотреть суммарную АЧХ такой цепочки : fir1-децимация-fir2? Допустим в качестве fir1 выступает фильтр с h1=conv(conv(conv(conv(p,p),p),p),p);
Сигнал с его выхода децимирую на 8 и подаю на второй фильтр с ИХ h2.

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

дальше freqz

[_Anatoliy]

h=conv(h1,8*upsample(h2,8));
дальше freqz

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

[soldat_shveyk]

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

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

А я пытался наоборот - децимировать h1 и естественно ничего не получалось.

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

[_Anatoliy]

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

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

[_Anatoliy]

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

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

 

[Кнкн]

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

дальше freqz

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

[thermit]

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

Можно.

[honeycomb0]

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

Мне нужно продецимировать сигнал полосой 50МГц в 16 раз из 1200 MSPS в 75 MSPS.

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

порядок определяется количеством нулей передаточной функции.

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

[honeycomb0]

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

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