Чтобы сравнить эффективность двух интерполяторов, один CIC-фильтр и другой КИХ-фильтр, используют gain функция, чтобы настроить CIC-фильтр выходная амплитуда, чтобы совпадать с КИХ фильтрует выходную амплитуду.

Запустите путем создания набора входных данных, синусоидального сигнала x.

fs = 1000;          % Input sampling frequency.
t = 0:1/fs:1.5;     % Signal length = 1501 samples.
x = sin(2*pi*10*t); % Amplitude = 1 sinusoid.
x = x';

Спроектируйте каскад двух dsp.FIRInterpolator объекты с полным коэффициентом интерполяции 4.

l = 4; % Interpolation factor for FIR filter.
d = fdesign.interpolator(l);
firInterp = design(d,'multistage','SystemObject',true);

Запустите данные через интерполятор.

yfir = firInterp(x);

Спроектируйте dsp.CICInterpolator объект с коэффициентом интерполяции установил 4, дифференциальный набор задержки к 1, и количество набора разделов к 4.

r = 4; % Interpolation factor for the CIC filter.
d = fdesign.interpolator(r,'cic');
cicInterp = design(d,'multisection','SystemObject',true);

Запустите те же данные через фильтр.

ycic = cicInterp(x);

Используйте gain функция, чтобы настроить CIC-фильтр выходная амплитуда, чтобы совпадать с КИХ фильтрует выходную амплитуду.

gaincic = gain(cicInterp);
subplot(211);
plot([yfir; double(ycic)]);
subplot(212)
plot([yfir; double(ycic)/gain(cicInterp)]);

После исправления для усиления, вызванного интерполятором CIC, во втором подграфике, вы видите, что КИХ-фильтр и CIC-фильтр обеспечивают почти идентичную интерполяцию.

Это усиление равняется усилению последнего раздела CIC-фильтра. Чтобы подтвердить, откорректируйте КИХ-амплитуду фильтра с помощью gain($\mathrm{cicInterp}$,$2\mathit{N}$). Если $\mathit{N}$ количество интегратора и разделы расчески CIC-фильтра, затем $2\mathit{N}$ последний раздел CIC-фильтра. $\mathit{N}$ дан cicInterp.NumSections.

Второй подграфик показывает, что КИХ-фильтр и CIC-фильтр обеспечивают почти идентичную интерполяцию, когда усиление коррекции равняется усилению последнего раздела CIC-фильтра.

subplot(212);
plot([yfir; double(ycic)/gain(cicInterp,2*cicInterp.NumSections)]);