Проект дробных КИХ-фильтров задержки

Дробные фильтры задержки нацелены на сдвиг цифровой последовательности значением нецелого числа посредством интерполяции и передискретизации, объединенной в один фильтр свертки. Этот пример демонстрирует разработку и реализацию дробных КИХ-инструментов использования фильтров задержки, доступных в DSP System Toolbox™.

Задержитесь как система свертки

Целочисленные задержки задержек

Рассмотрите задержку цифрового сигнала, y[n]=x[n-D] где D является целым числом. Эта операция может быть представлена как фильтр свертки y=h*x, с конечной импульсной характеристикой h[n]=δ[n-D]. Соответствующая передаточная функция H(z)=z-D, и frequeny ответ H(ω)=e-iωD. Программно, можно реализовать такой целочисленный фильтр задержки с помощью следующего кода MATLAB®.

% Create the FIR
D = 3; % Delay value
h = [zeros(1,D) 1]
h = 1×4

     0     0     0     1

Переключите последовательность путем пропущения его через КИХ h. Отметьте начальные нули в начале выхода, те показывают начальное условие, которое свойственно к таким фильтрам.

x = (1:10)';
dfir = dsp.FIRFilter(h);
y = dfir(x)'
y = 1×10

     0     0     0     1     2     3     4     5     6     7

Задержки нецелого числа через Интерполяцию D/A

Задержка последовательности x[n-D] не задан каждый раз, когда D не является целым числом. Чтобы сделать такие дробные задержки разумными, нужно добавить промежуточный этап интерполяции D/A, чтобы произвести выход на contniuum. Таким образом, y[n]=xnˆ(n-D) где xnˆ обозначает некоторую интерполяцию D/A входа sequnce x. D/A, интерполирующий функциюxnˆ мог зависеть от с n и мог считаться представлением базовой модели аналогового сигнала от который последовательность x был произведен. Эта стратегия используется в других проблемах передискретизации, таких как преобразование уровня.

Этот пример покажет дробные фильтры задержки с помощью двух моделей интерполяции, обе из которых предлагаются как часть Системного Тулбокса DSP.

  1. Находящаяся в sinc модель интерполяции, которая использует bandlimited реконструкцию для xnˆ.

  2. Основанная на Лагранже модель интерполяции, которая использует полиномиальную реконструкцию для xnˆ.

Bandlimited дробные фильтры задержки

Интерполяционная формула Шенноновского Уиттекера xˆ(t)=kx[k]sinc(t-k) модели bandlimited сигналы. Таким образом, промежуточное преобразование D/A xˆ bandlimited реконструкция входной последовательности. Для значения задержки D, дробной задержки y[n]=xˆ(n-D), который использует то же самое xˆ для каждого n, может быть представлен как фильтр свертки. Этот фильтр называется идеалом bandlimited дробным фильтром задержки, и его импульсная характеристика

hD[k]=sinc(k-D).

Соответствующей частотной характеристикой (то есть, DTFT) дают HD(ω)=e-iωD.

Причинное КИХ-приближение идеального фильтра сдвига Bandlimited

Идеал sinc фильтр сдвига, описанный в предыдущем разделе, является фильтром все-передачи (i.e. |Hd(ω)|=1), но это имеет бесконечную и непричинную импульсную характеристику hD. В MATLAB это не может быть представлено как вектор, а скорее в зависимости от индекса k.

% Ideal Filter sequence
D = 0.4;
hIdeal = @(k) sinc(k-D); 

В практических и вычислительных целях идеальный фильтр может быть усеченным на конечном окне индекса по стоимости некоторой потери пропускной способности. Для целевого значения задержки D и желаемая длина N, окно индексов kудовлетворение |D-k|N2 симметрично о D, и получает основной лепесток идеального фильтра. Для D=i0+FD где 0FD1 и целое число i0, явные индексы окна{i0-N-12,,i0+N2}. Целое число i0 упоминается как целочисленная задержка и может быть выбран произвольно. Чтобы сделать КИХ причинным, установить i0=N-12, таким образом, окно индекса {0,,N-1}. Код ниже изображает объяснение позади причинного КИХ-приближения.

% FIR approximation with causal shift
N = 6;
idxWindow = (-floor((N-1)/2):floor(N/2))';
i0 = -idxWindow(1); % Causal latency
hApprox = hIdeal(idxWindow);
plot_causal_fir("sinc",D,N,i0,hApprox,hIdeal);

Figure contains 2 axes. Axes 1 with title sinc-based Fractional Delay Non-Causal FIR, Truncated on [-2:3] contains 6 objects of type stem, patch, scatter, line, constantline. These objects represent Truncated FIR, Anti-Causal Indices, Causal Indices, Ideal Filter. Axes 2 with title Causal Shift of sinc-based Fractional Delay FIR, Truncated on [0:5] contains 3 objects of type stem, patch, constantline. These objects represent Truncated Causal FIR, All Causal Indices.

Усечение фильтра sinc вызывает пульсацию в частотной характеристике, которая может быть обращена путем применения весов {wk} (такие как Кайзер или Хэмминг) к КИХ-коэффициентам.

Наконец, получившаяся КИХ-модель приближения идеала bandlimited дробный фильтр задержки приведена ниже.

h[k]=wkhd[k]={wksinc(k-FD-i0)0kN-10otherwise 

Можно спроектировать такой фильтр с помощью designFracDelayFIR функционируйте и dsp.VariableFractionalDelay Системный объект в 'FIR' режим, оба из которых используют веса окна Кайзера.

Основанные на Лагранже дробные фильтры задержки

Основанные на Лагранже дробные фильтры задержки используют polynmial, соответствующий на движущемся окне входных выборок. Таким образом, xnˆ(t) полином некоторой фиксированной степени K. Как находящиеся в sinc фильтры задержки, основанные на Лагранже фильтры задержки могут быть сформулированы как причинная КИХ-свертка (i.e. y=h*x) из длины N=K+1, и поддерживаемый на окне индекса{-N-12,N2}. Так же к находящейся в sinc модели, примените причинную задержку i0=N-12. Учитывая дробную задержку FD, КИХ-коэффициенты h[0],,h[K] из (причинный переключенный) Лагранжев фильтр задержки может быть получен путем решения системы линейного equaions, как записано ниже. Те уравнения определяют стандартный Лагранжев многочлен подходящая проблема.

k=0Ktnkh[k]=(FD)n,n=0,,K

Здесь, t0,,tK перечислимые индексы демонстрационного окна. Реализация является прямой.

% Filter parameters
FD = 0.4;
K = 7;    % Polynomial degree
N = K+1;  % FIR Length
idxWindow = (-floor((N-1)/2):floor(N/2))';

% Define and solve Lagrange interpolation equations
V = idxWindow.^(0:K); % Vandermonde structure
C = FD.^(0:K);

hLagrange = C/V;  % Solve for the coefficients
i0 = -idxWindow(1); % Causal latency
plot_causal_fir("Lagrange",FD,N,i0,hLagrange);

Figure contains 2 axes. Axes 1 with title Lagrange-based Fractional Delay Non-Causal FIR, Supported on [-3:4] contains 4 objects of type stem, patch, constantline. These objects represent FIR, Anti-Causal Indices, Causal Indices. Axes 2 with title Causal Shift of Lagrange-based Fractional Delay FIR, Supported on [0:7] contains 3 objects of type stem, patch, constantline. These objects represent Causal FIR, All Causal Indices.

Эта модель может быть реализована как КИХ-фильтр прямой формы, если значение задержки, FD фиксируется, или использование структуры Фэрроу, если значение задержки варьируется. Существует раздел ниже специализированного implementaiton Лагранжевой интерполяции с помощью dsp.VariableFractionalDelay в 'Farrow' режим.

Разработка и реализация находящиеся в sinc Дробные КИХ-Фильтры Задержки

Следующий раздел является особым вниманием при разработке и реализовании находящихся в sinc дробных фильтров задержки.

Функция designFracDelayFIR в режиме проектирования длины

Функциональный designFracDelayFIR обеспечивает простой интерфейс, чтобы спроектировать дробный КИХ-фильтр задержки значения задержки FD и длины N.

FD = 0.32381;
N = 10;
h = designFracDelayFIR(FD,N)
h = 1×10

    0.0046   -0.0221    0.0635   -0.1664    0.8198    0.3926   -0.1314    0.0552   -0.0200    0.0042

Сама реализация может быть сделана как стандартный КИХ-фильтр, такой как dsp.FIRFilter Системный объект.

% Create an FIR filter object
fdfir = dsp.FIRFilter(h); 

Задержите сигнал путем пропущения его через спроектированный фильтр.

% Generate some input
n = (1:100)';
x = gen_input_signal(n);

% Filter the input signal
y = fdfir(x);
plot_sequences(n,x, n,y);
legend('Filter Output','Original Sequence')
title('Raw Filter Output v.s. Input Sequence')

Figure contains an axes. The axes with title Raw Filter Output v.s. Input Sequence contains 2 objects of type scatter. These objects represent Filter Output, Original Sequence.

release(fdfir);

Заметьте, что фактическая задержка фильтра не FD, а скорее FD+i0 из-за причинной целочисленной задержки i0. Та задержка возвращена в designFracDelayFIR функция также. Переключите график входной последовательности FD+i0 выровнять фильтр выход с ожидаемым результатом.

[h,i0] = designFracDelayFIR(FD,N);
y = fdfir(x);
plot_sequences(n+i0+FD,x, n,y);
legend('Filter Output','Input Sequence (shifted by FD+i0)')
title('Filter Output v.s. Time Adjusted Input Sequence')

Figure contains an axes. The axes with title Filter Output v.s. Time Adjusted Input Sequence contains 2 objects of type scatter. These objects represent Filter Output, Input Sequence (shifted by FD+i0).

Обратите внимание на то, что переключенные входные маркеры, расположенные в (n+FD+i0,x[n]) обычно не совпадайте с выходными демонстрационными маркерами (n,y[n]), потому что n+FD+i0 падения на значениях нецелого числа на оси X, тогда как n является целым числом. Скорее переключенные входные выборки падают приблизительно на линию, соединяющую каждые два последовательных выхода выборки.

plot_sequences(n+i0+FD,x, n,y,'with line');
legend('Filter Output','Input Sequence (shifted by FD+i0)')
title('Output Samples v.s. Shifted Input Samples ')
xlim([20,30])

Figure contains an axes. The axes with title Output Samples v.s. Shifted Input Samples contains 3 objects of type scatter, line. These objects represent Filter Output, Input Sequence (shifted by FD+i0).

dsp.VariableFractionalDelay Системный объект in 'FIR' Режим

Так же к designFracDelayFIR, dsp.VariableFractionalDelay объект может также спроектировать находящиеся в sinc фильтры задержки, когда используется с 'FIR' режим интерполяции. Начните путем создания экземпляра Системного объекта. КИХ-длина всегда даже и задана как поясной параметр.

vfd_fir = dsp.VariableFractionalDelay('InterpolationMethod','FIR','FilterHalfLength',N/2);
i0_vfd_fir = vfd_fir.FilterHalfLength;    % Interger latency

Передайте желаемую дробную задержку как второй входной параметр к вызову объектов. Убедитесь, что значение задержки, которое вы задаете, включает целочисленную задержку.

y = vfd_fir(x,i0+FD);
release(vfd_fir)
plot_sequences(n+i0+FD,x, n,y);
legend('Filter Output','Input Sequence (shifted by FD+i0)')
title('dsp.VariableFractionalDelay in FIR Mode')

Figure contains an axes. The axes with title dsp.VariableFractionalDelay in FIR Mode contains 2 objects of type scatter. These objects represent Filter Output, Input Sequence (shifted by FD+i0).

Сравнение designFracDelayFIR и dsp.VariableFractionalDelay в 'КИХ-' Режиме

Оба designFracDelayFIR и dsp.VariableFractionalDelay в 'FIR' режим обеспечивает находящиеся в sinc дробные фильтры задержки, но их реализации отличаются.

  • dsp.VariableFractionalDelay аппроксимирует значение задержки рациональным числом FDkL до некоторого допуска, и затем производит дробную задержку как k-ую фазу (длинного) фильтра интерполяции длины L. Это требует увеличенного использования памяти и дает к менее точной задержке.

  • В отличие от этого, designFracDelayFIR генерирует КИХ-коэффициенты непосредственно, вместо того, чтобы произвести их от более длинного КИХ. Это дает точное дробное значение задержки и стоит меньшей памяти.

  • designFracDelayFIR имеет простой функциональный интерфейс, возвращающий КИХ-коэффициенты, оставляя реализацию фильтра пользователь. dsp.VariableFractionalDelay Системный объект, предназначенный, чтобы инкапсулировать создание фильтра и реализацию полностью.

Использование designFractionalDelayFIR предпочтен по dsp.VariableFractionalDelay в 'FIR' режим для его простоты, лучшей эффективности и КПД. На рисунке ниже, фильтр спроектирован с dsp.VariableFractionalDelay имеет более короткую пропускную способность, и ее групповая задержка выключена ~0.02 от номинальной стоимости.

% Obtain the FIR coefficients from the dsp.VariableFractionalDelay object
h_vfd_fir = vfd_fir([1;zeros(31,1)],i0_vfd_fir+FD);
release(vfd_fir);
plot_freq_and_gd(h,i0,[],"designFracDelayFIR", h_vfd_fir,i0_vfd_fir,[],"dsp.VariableFractionalDelay FIR mode");
hold on;
yline(FD,'DisplayName','Target Fractional Delay');
ylim([-0.1,0.4])

Figure contains 2 axes. Axes 1 with title Gain response contains 2 objects of type line. These objects represent Gain ResponsedesignFracDelayFIR, Gain Response dsp.VariableFractionalDelay FIR mode. Axes 2 with title Group Delay (i0 adjusted) contains 3 objects of type line, constantline. These objects represent Group Delay designFracDelayFIR, Group Delay dsp.VariableFractionalDelay FIR mode, Target Fractional Delay.

Разработка и реализация основанные на Лагранже фильтры задержки

Основанный на Лагранже дробный фильтр задержки является в вычислительном отношении дешевым и может быть реализован эффективно с помощью структуры Фэрроу. Фильтр Фэрроу является специальным типом КИХ, который реализован с помощью только элементарные алгебраические операции, такие как скалярные сложения и умножение. В отличие от находящихся в sinc проектов, фильтры Фэрроу не требуют, чтобы специализированные функции (такие как sinc или функция Бесселя) вычислили КИХ-коэффициенты задержки. Это делает Фэрроу, которую дробная задержка фильтрует особенно простой реализовать на базовом оборудовании.

На оборотной стороне основанные на Лагранже фильтры задержки ограничиваются младшими разрядами, из-за очень нестабильного характера знатных полиномиальных приближений. Это обычно заканчивается с более низкой пропускной способностью, при сравнении с находящимся в sinc фильтром.

Системный объект dsp.VariableFractionalDelay в 'Farrow' mode

Используйте системный объект dsp.VariableFractionalDelay в 'Farrow' режим, чтобы создать и реализовать Фэрроу задерживает фильтры. Начните путем создания экземпляра системного объекта:

vfd = dsp.VariableFractionalDelay('InterpolationMethod','Farrow','FilterLength',8);
i0var = floor(vfd.FilterLength/2)  % Interger latency of the filter
i0var = 4

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

y = vfd(x,i0var+FD);
plot_sequences(n+i0var+FD,x, n,y);
legend('Farrow Fractional Delay Output','Input Sequence (shifted by FD+i0)')
title('dsp.VariableFractionalDelay in Farrow Mode')

Figure contains an axes. The axes with title dsp.VariableFractionalDelay in Farrow Mode contains 2 objects of type scatter. These objects represent Farrow Fractional Delay Output, Input Sequence (shifted by FD+i0).

Можно также варьироваться дробные значения задержки. Код ниже работает с системами координат 20 выборок при увеличении значения задержки с каждой системой координат. Отметьте увеличение задержки выходного графика, соответствуя изменениям в значениях задержки.

release(vfd)
FDs = i0var+5*(0:0.2:0.8); % Fractional delays vector
xsource = dsp.SignalSource(x,20);
ysink = dsp.AsyncBuffer;
for FD=FDs
    xk = xsource();
    yk = vfd(xk, FD);
    write(ysink,yk);
end
y = read(ysink);

plot_sequences(n+i0var,x, n,y);
legend('Variable Fractional Delay Output','Original Sequence (shifted by i0)')
title('dsp.VariableFractionalDelay in Farrow Mode, Varying Delay')

Figure contains an axes. The axes with title dsp.VariableFractionalDelay in Farrow Mode, Varying Delay contains 2 objects of type scatter. These objects represent Variable Fractional Delay Output, Original Sequence (shifted by i0).

Пропускная способность КИХ дробные фильтры задержки: анализ и проектирование

Более длинные фильтры дают лучшее приближение идеального фильтра задержки. Действительно, в терминах необработанных квадратичных норм это имеет место. Однако нам нужна метрика, которая более практически значима, такова как пропускная способность. Функциональный designFracDelayFIR меры объединили пропускную способность, которая задана как частотный диапазон, в котором и усиление и групповая задержка в 1% их номинальной стоимости. Измеренная объединенная пропускная способность может быть получена как возвращаемое значение designFracDelayFIR функция. Сравните фильтр длины 16 (синих) с фильтром длины 256 (красных) на рисунке ниже. Как ожидалось, более длинный фильтр имеют значительно более высокую объединенную пропускную способность.

FD = 0.3;
N1 = 16;
N2 = 256;
[h1,i1,bw1] = designFracDelayFIR(FD, N1);
[h2,i2,bw2] = designFracDelayFIR(FD, N2);

plot_freq_and_gd(h1,i1,bw1,"N="+num2str(N1), h2,i2,bw2,"N="+num2str(N2));
ylim([-0.2,0.6])

Figure contains 2 axes. Axes 1 with title Gain response contains 4 objects of type line, constantline. These objects represent Gain ResponseN=16, Gain Response N=256, Combined Bandwidth N=16, Combined Bandwidth N=256. Axes 2 with title Group Delay (i0 adjusted) contains 4 objects of type line, constantline. These objects represent Group Delay N=16, Group Delay N=256, Combined Bandwidth N=16, Combined Bandwidth N=256.

Функция designFracDelayFIR в режиме проектирования пропускной способности

Режим проектирования пропускной способности designFracDelayFIR может определить необходимую длину для данной пропускной способности. Задайте значение задержки и желаемую целевую пропускную способность как входные параметры к функции, и функция найдет соответствующую длину.

FD = 0.3;
bwLower = 0.9; % Target bandwidth lower limit
[h,i0fixed,bw] = designFracDelayFIR(FD,bwLower);
fdfir = dsp.FIRFilter(h);
info(fdfir)
ans = 6x35 char array
    'Discrete-Time FIR Filter (real)    '
    '-------------------------------    '
    'Filter Structure  : Direct-Form FIR'
    'Filter Length     : 52             '
    'Stable            : Yes            '
    'Linear Phase      : No             '

Обратите внимание на то, что bwLower просто нижняя граница для объединенной пропускной способности. Функция возвращает фильтр, объединенная пропускная способность которого является, по крайней мере, значением, заданным в bwLow.

Искажение в высоких сигналах пропускной способности

В этом разделе мы сравниваем эффективность двух точек проекта (длинный sinc v.s. короткий Лагранж) с высоким входом пропускной способности. dsp.VariableFractionalDelay в предыдущем разделе структура Фэрроу с 8 степенями, эффективно КИХ длины 9. Фильтр получен designFracDelayFIR(FD,0.9) имеет продолжительность 52 выборок. Соединение двух КИХ-частотных характеристик на том же графике демонстрирует различие в пропускной способности между двумя.

release(vfd);
hvar = vfd([1;zeros(31,1)],i0var+FD);
plot_freq_and_gd(h,i0fixed,bw,"Sinc-based", hvar,i0var,[],"Farrow");
ylim([-0.2,0.6])

Figure contains 2 axes. Axes 1 with title Gain response contains 2 objects of type line. These objects represent Gain ResponseSinc-based, Gain Response Farrow. Axes 2 with title Group Delay (i0 adjusted) contains 2 objects of type line. These objects represent Group Delay Sinc-based, Group Delay Farrow.

Примените два фильтра на высокий сигнал пропускной способности, как сравнено на рисунке ниже. Sinc на левом столбце, Фэрроу справа. Временной интервал на верхней строке, частота на нижней части. Результаты, как ожидалось:

  • Дольше sinc фильтр имеет более высокую пропускную способность. Более короткий фильтр Фэрроу имеет более низкую пропускную способность.

  • Искажение сигнала является фактически несуществующим использованием дольше sinc фильтр, но легко примечательный в более коротком фильтре Фэрроу.

  • Более высокая точность прибывает за счет более длинной задержки: приблизительно 25 выборок v.s. только 4 в более коротком фильтре.

n=(1:80)';
x = high_bw_signal(n);

y1 = fdfir(x);
y2 = vfd(x,i0var+FD);

plot_signal_comparison(n,x,y1,y2,h,hvar,i0fixed,i0var,FD);

Figure contains 4 axes. Axes 1 with title designFracDelayFIR contains 4 objects of type line, scatter. These objects represent True Input Signal (delayed), Filtered signal. Axes 2 contains 2 objects of type line. These objects represent Signal, FIR Gain. Axes 3 with title dsp.VariableFractionalDelay contains 4 objects of type line, scatter. These objects represent True Input Signal (delayed), Filtered signal. Axes 4 contains 2 objects of type line. These objects represent Signal, FIR Gain.

Который должен использоваться: dsp.VariableFractionalDelay или designFracDelayFIR ?

Это решение в основном основано на требованиях фильтра и целевой платформе.

  • Для высокой пропускной способности и точного ответа групповой задержки, используйте designFracDelayFIR функция. Следует иметь в виду, что этот процесс проектирования более в вычислительном отношении интенсивен. Поэтому он лучше подходящий быть развернутым на оборудовании более высокого качества, особенно если настройка в реальном времени значения задержки желаема. Это также подходит для аппаратного развертывания более низкого уровня, если значение задержки фиксируется, и проект может быть сделан оффлайн.

  • Для изменяющихся во времени фильтров задержки, нацеленных на низкую эффективность вычислительный аппарат, используйте dsp.VariableFractionalDelay с 'Farrow' режим.

Для просмотра документации необходимо авторизоваться на сайте