[Jools 0]

Привет всем.

 

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

 

H(f) = (abs(sin(pi*M*f) / sin(pi*f/R)))^2N, где

R - децимирующий фактор

N - порядок

M - дифференциальная задержка

 

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

 

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

[hobgoblin 0]

H(f) = (abs(sin(pi*M*f) / sin(pi*f/R)))^2N, где

R - децимирующий фактор

N - порядок

M - дифференциальная задержка

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

f is the frequency relative to the low sampling rate fs/R

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

[petrov 6]

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

[Jools 0]

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

 

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 6]

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

 

Да конечно.