Один раздел или расположенный каскадом фильтр allpass
dsp.AllpassFilter возразите фильтрует каждый канал входного использования реализации фильтра Allpass. Чтобы импортировать этот объект в Simulink®, используйте блок MATLAB® System.
Примечание
Поддержка массива ячеек AllpassCoefficients, WDFCoefficients, и LatticeCoefficients был удален. Используйте N-by-1 или N-by-2 числовой массив вместо этого. Для получения дополнительной информации см. Вопросы совместимости.
Отфильтровать каждый канал входа:
Создайте dsp.AllpassFilter объект и набор его свойства.
Вызовите объект с аргументами, как будто это была функция.
Чтобы узнать больше, как Системные объекты работают, смотрите то, Что Системные объекты?.
Allpass = dsp.AllpassFilterAllpass, это фильтрует каждый канал входного сигнала независимо с помощью фильтра allpass со структурой по умолчанию и коэффициентами.
Allpass = dsp.AllpassFilter(Name,Value)Allpass, с каждым набором свойств к заданному значению.
Чтобы использовать объектную функцию, задайте Системный объект как первый входной параметр. Например, чтобы выпустить системные ресурсы Системного объекта под названием obj, используйте этот синтаксис:
release(obj)
Примечание: Этот пример запускается только в R2016b или позже. Если вы используете более ранний релиз, заменяете каждый вызов функции с эквивалентным step синтаксис. Например, myObject (x) становится шагом (myObject, x).
Создайте фильтры Allpass
Fs = 48000; % in Hz FL = 1024; APF1 = dsp.AllpassFilter('AllpassCoefficients',... [-0.710525516540603 0.208818210000029]); APF2 = dsp.AllpassFilter('AllpassCoefficients',... [-0.940456403667957 0.6;... -0.324919696232907 0],... 'TrailingFirstOrderSection',true);
Создайте Средство оценки Передаточной функции, чтобы оценить, что передаточная функция между случайным входом и Allpass отфильтровала выход
TFE = dsp.TransferFunctionEstimator('FrequencyRange',... 'onesided','SpectralAverages',2);
Создайте ArrayPlot, чтобы построить ответ величины
AP = dsp.ArrayPlot('PlotType','Line','YLimits', [-80 5],... 'YLabel','Magnitude (dB)','SampleIncrement', Fs/FL,... 'XLabel','Frequency (Hz)','Title','Magnitude Response',... 'ShowLegend', true,'ChannelNames',{'Magnitude Response'});
Отфильтруйте Вход и покажите ответ величины предполагаемой передаточной функции между входом и отфильтрованным выходом
tic; while toc < 5 in = randn(FL,1); out = 0.5.*(APF1(in) + APF2(in)); A = TFE(in, out); AP(db(A)); end
Передаточной функцией фильтра allpass дают
.
c является allpass полиномиальным содействующим вектором. Порядок, n, передаточной функции являются длиной векторного c.
В минимальной форме множителя и цифровой форме волны, фильтр allpass реализуется как каскад или второго порядка (biquad) разделы или разделы первого порядка. Когда коэффициенты заданы как N-by-2 матрица, каждая строка матрицы задает коэффициенты фильтра второго порядка. Последний элемент последней строки может быть проигнорирован на основе запаздывающей установки первого порядка. Когда коэффициенты заданы как N-by-1 матрица, каждый элемент в матрице задает коэффициент фильтра первого порядка. Каскад всех разделов фильтра формирует фильтр allpass.
В форме решетки коэффициенты заданы как вектор.
Эти структуры в вычислительном отношении более экономичны и структурно более устойчивы по сравнению с типовыми БИХ-фильтрами, таковы как df1, df1t, df2, df2t. Для всех структур фильтр allpass может быть одно разделом или несколькими-разделами (расположенный каскадом) фильтр. Различные разделы могут иметь различные порядки, но они все реализованы согласно той же структуре.
Эта структура понимает фильтр allpass с минимальным количеством необходимых множителей, равных порядку n. Это также использует 2n задержите модули и 2n сумматоры. Множители используют заданные коэффициенты, которые равны полиномиальному векторному c в allpass передаточной функции. В этой секции второго порядка минимальной структуры множителя содействующий вектор, c, равен [0.1 -0.7].
Эта структура использует n множители, но только n задержите модули, за счет требования 3n сумматоры. Чтобы использовать эту структуру, задайте коэффициенты в форме цифрового фильтра волны (WDF). Получите эквивалент WDF обычного allpass содействующего использования allpass2wdf(allpass_coefficients). Чтобы преобразовать коэффициенты WDF в эквивалентную allpass полиномиальную форму, использовать wdf2allpass(WDF coefficients). В этой секции второго порядка структуры WDF содействующий вектор w равен allpass2wdf([0.1 -0.7]).
Эта структура решетки использует 2n множители, n задержите модули и 2n сумматоры. Чтобы использовать эту структуру, задайте коэффициенты как вектор.
Можно получить решетку, эквивалентную из обычного allpass содействующего использования transpose(tf2latc(1, [1 allpass_coefficients])). В следующей секции второго порядка структуры решетки содействующий вектор вычисляется с помощью transpose(tf2latc(1,
[1 0.1 -0.7])). Используйте эти коэффициенты для фильтра, который функционально эквивалентен минимальной структуре множителя с коэффициентами [0.1 - 0.7].
[1] Регалии, Филип А. и Митра Сэнджит К. и Вэйдьянэзэн, P. P. (1988) “Цифровой Фильтр все-Передачи: Базовый блок Обработки сигналов AVersatile”. Продолжения IEEE, Издания 76, № 1, 1988, стр 19–37
[2] М. Лутовэк, Д. Тошич, Б. Эванс, создание фильтра для обработки сигналов Используя MATLAB и Mathematica. Верхний Сэддл-Ривер, NJ: Prentice Hall, 2001.
