Спроектируйте фильтр по умолчанию

Создайте dsp.FIRRateConverter объект. Объект использует объединенную антиобработку изображений и сглаживание КИХ-фильтра. По умолчанию этот фильтр создан с использованием designMultirateFIR функция. Функция проектирует фильтр на основе коэффициентов преобразования уровня, которые вы задаете, и хранит коэффициенты в Numerator свойство. Для коэффициента интерполяции 3 и фактора децимации 2, объект проектирует коэффициенты с помощью designMultirateFIR(3,2).

firrc = dsp.FIRRateConverter(3,2);

Передискретизируйте на коэффициент 3/2

Передискретизируйте сигнал на коэффициент 3/2.

y = firrc(x);

Постройте исходные и передискретизируемые сигналы. Для того, чтобы построить два сигнала на том же графике, получите параметры аффинного преобразования времени. Используйте эти параметры, чтобы вычислить выходные шаги расчета. ${\mathit{n}}^{\mathrm{th}}$ выходной шаг расчета равняется $\mathrm{scale}\times \mathit{n}+\mathrm{delay}$. Значения ввода и вывода совпадают каждые 3 выходных выборки и любая входная выборка, вследствие скорости преобразования 3/2.

scale = firrc.InterpolationFactor/firrc.DecimationFactor;
delay = length(firrc.Numerator)/(2*firrc.DecimationFactor);

% Calculate output times for vector y in input units
ty = (((0:length(y)-1)-delay)/scale)/sine.SampleRate;

stem(tx, x,'filled',MarkerSize=4); hold on;
stem(ty, y);
hold off;
xlim([0.0 0.0145])
ylim([-1.5 1.5])
legend('Original input','Resampled');

Передискретизируйте на коэффициент 5/3 в Автоматическом Режиме Создания фильтра

Теперь измените коэффициент интерполяции в 5 и фактор децимации к 3. Для создания фильтра, которое будет обновлено автоматически на основе новых коэффициентов преобразования уровня, устанавливает NumeratorSource свойство к 'Auto'. Альтернативно, можно передать ключевое слово 'auto' как вход при создании объекта. Объект затем действует в автоматическом режиме создания фильтра. Каждый раз существует изменение в коэффициентах преобразования уровня, обновления объекта создание фильтра соответственно.

release(firrc)
firrc.NumeratorSource = 'Auto';
firrc.InterpolationFactor = 5;
firrc.DecimationFactor = 3

firrc = 
  dsp.FIRRateConverter with properties:

    InterpolationFactor: 5
       DecimationFactor: 3
        NumeratorSource: 'Auto'
           DesignMethod: 'Kaiser'

  Show all properties

Чтобы получить доступ к коэффициентам фильтра в автоматическом режиме создания фильтра, введите firrc.Numerator в командной строке MATLAB.

Передискретизируйте сигнал с обновленными значениями преобразования уровня.

yAuto = firrc(x);

Постройте исходные и передискретизируемые сигналы. Повторно вычислите параметры аффинного преобразования времени, поскольку коэффициенты преобразования уровня изменились. Обратите внимание, что ввод и вывод совпадает каждые 3 входных выборки и каждые 5 выходных выборок, вследствие 5/3 коэффициента преобразования.

scale = firrc.InterpolationFactor/firrc.DecimationFactor;
delay = length(firrc.Numerator)/(2*firrc.DecimationFactor);

% Calculate output times for vector yAuto in input units
tyAuto = (((0:length(yAuto)-1)-delay)/scale)/sine.SampleRate;

stem(tx, x,'filled',MarkerSize=4); hold on;
stem(tyAuto, yAuto,'r');
hold off;
xlim([0.0 0.015])
ylim([-1.5 1.5])
legend('Original input','Resampled');

Задайте модель интерполяции сигнала

В режиме автоматического проектирования можно также задать базовую модель интерполяции сигнала D/A через DesignMethod свойство.

Установите DesignMethod к 'linear' и измените коэффициент интерполяции в 11.

release(firrc)
firrc.DesignMethod = 'linear';
firrc.InterpolationFactor = 11;

Передискретизируйте сигнал с помощью модели линейной интерполяции.

yLinear = firrc(x);

Постройте исходные и передискретизируемые сигналы. Выходные выборки лежат на кусочной линейной кривой. Обратите внимание, что ввод и вывод совпадает каждые три входных выборки и каждые 11 выходных выборок как ожидалось от отношения 11/3.

scale = firrc.InterpolationFactor/firrc.DecimationFactor;
delay = (length(firrc.Numerator)-1)/(2*firrc.DecimationFactor);

% Calculate output times for vector yLinear in input units
tyLinear = (((0:length(yLinear)-1)-delay)/scale)/sine.SampleRate;

stem(tx, x, 'filled',MarkerSize=4); hold on;
stem(tyLinear, yLinear);
plot(tyLinear, yLinear,Color=[1 0 0 0.3]);
hold off;
xlim([0.0 0.009])
ylim([-1.5 1.5])
legend('Original input','Resampled');