Чтобы сравнить эффективность двух интерполяторов, одного 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)]);