Практическое введение в создание цифровых фильтров

В этом примере показано, как спроектировать КИХ, и БИХ фильтрует на основе технических требований частотной характеристики с помощью designfilt функция в продукте Signal Processing Toolbox®. Концентраты в качестве примера на фильтрах lowpass, но большинстве результатов применяются к другим типам ответа также.

Этот пример фокусируется на проекте цифровых фильтров, а не на их приложениях. Если вы хотите узнать больше о приложениях цифрового фильтра, смотрите Практическое Введение в Цифровую Фильтрацию.

КИХ-создание фильтра

Технические требования фильтра lowpass

Идеальный фильтр lowpass является тем, который оставляет все частотные составляющие без изменений сигнала ниже обозначенной частоты среза, ωc, и отклонения все компоненты выше ωc. Поскольку импульсная характеристика, требуемая реализовать идеальный фильтр lowpass, бесконечно длинна, невозможно спроектировать идеального КИХ фильтр lowpass. Конечные приближения длины к идеальной импульсной характеристике приводят к присутствию пульсаций в обоих полосу пропускания (ω<ωc) и полоса задерживания (ω>ωc) из фильтра, а также к ненулевой ширине перехода между полосой пропускания и полосой задерживания.

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

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

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

Треугольник используется, чтобы изучить степени свободы, доступные при выборе технических требований проекта. Поскольку сумма углов фиксируется, каждый может в большей части выбора значения двух из технических требований. Третья спецификация будет определена конкретным алгоритмом проекта. Кроме того, как с углами в треугольнике, если мы делаем одни из технических требований большими/меньшими, они повлияют один или оба из других технических требований.

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

КИХ-проекты минимального порядка

Проекты минимального порядка получены путем определения полосы пропускания и частот полосы задерживания, а также неравномерности в полосе пропускания и затухания в полосе задерживания. Алгоритм проекта затем выбирает минимальную длину фильтра, которая выполняет технические требования.

Спроектируйте минимальный порядок КИХ-фильтр lowpass с частотой полосы пропускания 0.37*pi рад/отсчет, частота полосы задерживания 0.43*pi рад/отсчет (следовательно, ширина перехода равняется 0.06*pi рад/отсчет), неравномерность в полосе пропускания 1 дБ и затухание в полосе задерживания 30 дБ.

Fpass = 0.37; 
Fstop = 0.43;
Ap = 1;
Ast = 30;

d = designfilt('lowpassfir','PassbandFrequency',Fpass,...
  'StopbandFrequency',Fstop,'PassbandRipple',Ap,'StopbandAttenuation',Ast);

hfvt = fvtool(d);

Figure Filter Visualization Tool - Magnitude Response (dB) contains an axes and other objects of type uitoolbar, uimenu. The axes with title Magnitude Response (dB) contains 2 objects of type line.

Получившийся порядок фильтра может быть запрошен с помощью filtord функция.

N = filtord(d)
N = 39

Можно использовать info функция, чтобы получить информацию о параметрах раньше проектировала фильтр

info(d)
ans = 17x44 char array
    'FIR Digital Filter (real)                   '
    '-------------------------                   '
    'Filter Length  : 40                         '
    'Stable         : Yes                        '
    'Linear Phase   : Yes (Type 2)               '
    '                                            '
    'Design Method Information                   '
    'Design Algorithm : Equiripple               '
    '                                            '
    'Design Specifications                       '
    'Sample Rate     : N/A (normalized frequency)'
    'Response        : Lowpass                   '
    'Specification   : Fp,Fst,Ap,Ast             '
    'Passband Ripple : 1 dB                      '
    'Stopband Atten. : 30 dB                     '
    'Passband Edge   : 0.37                      '
    'Stopband Edge   : 0.43                      '

  • designfilt функция выбирает алгоритм проекта equiripple по умолчанию. Линейная фаза equiripple фильтры желательна, потому что для данного распоряжения у них есть самое маленькое максимальное отклонение от идеального фильтра.

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

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

dk = designfilt('lowpassfir','PassbandFrequency',Fpass,...
  'StopbandFrequency',Fstop,'PassbandRipple',Ap,...
  'StopbandAttenuation',Ast, 'DesignMethod', 'kaiserwin');

addfilter(hfvt,dk);
legend(hfvt,'Equiripple design', 'Kaiser window design')

Figure Filter Visualization Tool - Magnitude Response (dB) contains an axes and other objects of type uitoolbar, uimenu. The axes with title Magnitude Response (dB) contains 3 objects of type line. These objects represent Equiripple design, Kaiser window design.

N = filtord(dk)
N = 52

Определение параметров частоты в герц

Если вы знаете частоту дискретизации, на уровне которой будет действовать фильтр, можно задать частоту дискретизации и частоты в герц. Перепроектируйте минимальный порядок equiripple фильтр для частоты дискретизации 2 кГц.

Fpass = 370;
Fstop = 430;
Ap = 1;
Ast = 30;
Fs = 2000; 

d = designfilt('lowpassfir','PassbandFrequency',Fpass,...
  'StopbandFrequency',Fstop,'PassbandRipple',Ap,...
  'StopbandAttenuation',Ast,'SampleRate',Fs);
  
hfvt = fvtool(d);

Figure Filter Visualization Tool - Magnitude Response (dB) contains an axes and other objects of type uitoolbar, uimenu. The axes with title Magnitude Response (dB) contains 2 objects of type line.

Фиксированный порядок, фиксированная ширина перехода

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

Считайте 30-й порядок КИХ-фильтром lowpass с частотой полосы пропускания 370 Гц, частотой полосы задерживания 430 Гц и частотой дискретизации 2 кГц. Существует два метода разработки, доступные для этого конкретного набора технических требований: equiripple и наименьшие квадраты. Давайте спроектируем один фильтр для каждого метода и сравним результаты.

N = 30;
Fpass = 370;
Fstop = 430;
Fs = 2000; 

% Design method defaults to 'equiripple' when omitted
deq = designfilt('lowpassfir','FilterOrder',N,'PassbandFrequency',Fpass,...
  'StopbandFrequency',Fstop,'SampleRate',Fs);

dls = designfilt('lowpassfir','FilterOrder',N,'PassbandFrequency',Fpass,...
  'StopbandFrequency',Fstop,'SampleRate',Fs,'DesignMethod','ls');

hfvt = fvtool(deq,dls);
legend(hfvt,'Equiripple design', 'Least-squares design')

Figure Filter Visualization Tool - Magnitude Response (dB) contains an axes and other objects of type uitoolbar, uimenu. The axes with title Magnitude Response (dB) contains 3 objects of type line. These objects represent Equiripple design, Least-squares design.

Фильтры Equiripple идеально подходят для приложений, в которых определенному допуску нужно соответствовать, такие как разработка фильтра с данным минимальным затуханием в полосе задерживания или данной максимальной неравномерностью в полосе пропускания. С другой стороны, эти проекты не могут быть желательными, если мы хотим минимизировать энергию ошибки (между идеальным и фактическим фильтром) в полосе пропускания/полосе задерживания.

  • Если вы хотите уменьшать энергию сигнала как можно больше в определенном диапазоне частот, используйте проект наименьших квадратов.

В примерах выше, спроектированные фильтры имели ту же пульсацию в полосе пропускания и в полосе задерживания. Мы можем использовать веса, чтобы уменьшать пульсацию в одной из полос при хранении порядка фильтра зафиксированным. Например, если вы хотите, чтобы пульсация полосы задерживания была одной десятой из этого в полосе пропускания, необходимо дать полосе задерживания десять раз вес полосы пропускания. Перепроектируйте фильтр equiripple с помощью того факта.

deqw = designfilt('lowpassfir','FilterOrder',N,'PassbandFrequency',Fpass,...
  'StopbandFrequency',Fstop,'SampleRate',Fs,...
  'PassbandWeight',1,'StopbandWeight',10);

hfvt = fvtool(deq,deqw);
legend(hfvt,'Equiripple design', 'Equiripple design with weighted stopband')

Figure Filter Visualization Tool - Magnitude Response (dB) contains an axes and other objects of type uitoolbar, uimenu. The axes with title Magnitude Response (dB) contains 3 objects of type line. These objects represent Equiripple design, Equiripple design with weighted stopband.

Фиксированный порядок, фиксированная частота среза

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

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

Например, считайте 100-й порядок КИХ-фильтром lowpass с частотой среза 60 Гц и частоты дискретизации 1 кГц. Сравните проекты, которые следуют из использования Окна Хэмминга и Окна Чебышева с 90 дБ затухания бокового лепестка.

dhamming = designfilt('lowpassfir','FilterOrder',100,'CutoffFrequency',60,...
  'SampleRate',1000,'Window','hamming');

dchebwin = designfilt('lowpassfir','FilterOrder',100,'CutoffFrequency',60,...
  'SampleRate',1000,'Window',{'chebwin',90});

hfvt = fvtool(dhamming,dchebwin);
legend(hfvt,'Hamming window', 'Chebyshev window')

Figure Filter Visualization Tool - Magnitude Response (dB) contains an axes and other objects of type uitoolbar, uimenu. The axes with title Magnitude Response (dB) contains 3 objects of type line. These objects represent Hamming window, Chebyshev window.

Существуют другие пути, которыми можно задать фильтр с фиксированным порядком: фиксированная частота среза, неравномерность в полосе пропускания и затухание в полосе задерживания; фиксированная ширина перехода; и зафиксированная частота (на 3 дБ) на уровне половинной мощности.

БИХ-создание фильтра

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

  • БИХ-фильтры полезны, когда вычислительные ресурсы в большом почете. Однако устойчивые, причинные БИХ-фильтры не могут иметь совершенно линейной фазы. Избегайте БИХ-проектов в случаях, где линейность фазы является требованием.

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

Фильтры Баттерворта

Фильтры Баттерворта являются максимально плоскими БИХ-фильтрами. Плоскость в полосе пропускания и полосе задерживания заставляет полосу перехода быть очень широкой. Крупные заказы требуются, чтобы получать фильтры с узкими ширинами перехода.

Спроектируйте Фильтр Баттерворта минимального порядка с частотой полосы пропускания 100 Гц, частота полосы задерживания 300 Гц, максимальная неравномерность в полосе пропускания 1 дБ и затухание в полосе задерживания на 60 дБ. Частота дискретизации составляет 2 кГц.

Fp = 100;
Fst = 300;
Ap = 1;
Ast = 60;
Fs = 2e3;

dbutter = designfilt('lowpassiir','PassbandFrequency',Fp,...
  'StopbandFrequency',Fst,'PassbandRipple',Ap,...
  'StopbandAttenuation',Ast,'SampleRate',Fs,'DesignMethod','butter');

Фильтры Чебышевский Тип 1

Фильтры Чебышевский Тип 1 достигают меньших ширин перехода, чем Фильтры Баттерворта того же порядка путем обеспечения неравномерности в полосе пропускания.

  • У Баттерворта и фильтров Чебышевский Тип 1 оба есть максимально плоские полосы задерживания. Для данного порядка фильтра компромисс между шириной перехода и неравномерностью в полосе пропускания.

Спроектируйте фильтр Чебышевский Тип 1 с теми же техническими требованиями как Фильтр Баттерворта выше.

dcheby1 = designfilt('lowpassiir','PassbandFrequency',Fp,...
  'StopbandFrequency',Fst,'PassbandRipple',Ap,...
  'StopbandAttenuation',Ast,'SampleRate',Fs,'DesignMethod','cheby1');

Чебышевские фильтры типа II

  • Чебышевские фильтры Типа II имеют максимально плоские полосы пропускания и equiripple полосы задерживания.

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

Спроектируйте минимальный порядок фильтр Типа II Чебышева с теми же техническими требованиями как в предыдущих примерах.

dcheby2 = designfilt('lowpassiir','PassbandFrequency',Fp,...
  'StopbandFrequency',Fst,'PassbandRipple',Ap,...
  'StopbandAttenuation',Ast,'SampleRate',Fs,'DesignMethod','cheby2');

Эллиптические фильтры

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

  • Эллиптические фильтры достигают данной ширины перехода с наименьшим порядком.

dellip = designfilt('lowpassiir','PassbandFrequency',Fp,...
  'StopbandFrequency',Fst,'PassbandRipple',Ap,...
  'StopbandAttenuation',Ast,'SampleRate',Fs,'DesignMethod','ellip');

Сравните ответ и порядок четырех БИХ-фильтров.

  • Для тех же ограничений спецификации метод Баттерворта дает к самому высокому порядку, и эллиптический метод дает к самому маленькому.

FilterOrders = [filtord(dbutter) filtord(dcheby1) filtord(dcheby2) filtord(dellip)]
FilterOrders = 1×4

     7     5     5     4

hfvt = fvtool(dbutter,dcheby1,dcheby2,dellip);
axis([0 1e3 -80 2]);
legend(hfvt,'Butterworth', 'Chebyshev Type I',...
  'Chebyshev Type II','Elliptic')

Figure Filter Visualization Tool - Magnitude Response (dB) contains an axes and other objects of type uitoolbar, uimenu. The axes with title Magnitude Response (dB) contains 5 objects of type line. These objects represent Butterworth, Chebyshev Type I, Chebyshev Type II, Elliptic.

Масштабируйте в полосу пропускания, чтобы видеть различия в пульсации.

axis([0 150 -3 2]);

Figure Filter Visualization Tool - Magnitude Response (dB) contains an axes and other objects of type uitoolbar, uimenu. The axes with title Magnitude Response (dB) contains 5 objects of type line. These objects represent Butterworth, Chebyshev Type I, Chebyshev Type II, Elliptic.

Соответствие точно с техническими требованиями полосы пропускания или полосы задерживания

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

  • Используйте 'MatchExactly' параметр, чтобы ограничить алгоритм проекта совпадать с одной полосой точно. Другая полоса превышает свою спецификацию.

По умолчанию Чебышевский Тип, который я проектирую, совпадает с полосой пропускания, Баттерворт и Чебышевский Тип II совпадают с полосой задерживания, и эллиптические проекты совпадают и с полосой пропускания и с полосой задерживания (в то время как частота ребра полосы задерживания превышена):

dellip1 = designfilt('lowpassiir','PassbandFrequency',Fp,...
  'StopbandFrequency',Fst,'PassbandRipple',Ap,...
  'StopbandAttenuation',Ast,'SampleRate',Fs,'DesignMethod','ellip',...
  'MatchExactly','passband');

dellip2 = designfilt('lowpassiir','PassbandFrequency',Fp,...
  'StopbandFrequency',Fst,'PassbandRipple',Ap,...
  'StopbandAttenuation',Ast,'SampleRate',Fs,'DesignMethod','ellip',...
  'MatchExactly','stopband');

hfvt = fvtool(dellip, dellip1, dellip2);
legend(hfvt,'Matched passband and stopband','Matched passband',...
  'Matched stopband');
axis([0 1e3 -80 2]);

Figure Filter Visualization Tool - Magnitude Response (dB) contains an axes and other objects of type uitoolbar, uimenu. The axes with title Magnitude Response (dB) contains 4 objects of type line. These objects represent Matched passband and stopband, Matched passband, Matched stopband.

Совпадающая полоса пропускания и соответствующий - оба проекта имеет пульсацию точно 1 дБ в значении частоты полосы пропускания 100 Гц.

Сравнение групповой задержки

С БИХ-фильтрами мы должны рассмотреть не только компромисс ширины пульсации/перехода, но также и степень искажения фазы. Мы знаем, что невозможно иметь линейную фазу в целом интервале Найквиста. Таким образом мы можем хотеть, чтобы видеть, как совсем не линейный фазовый отклик. Хороший способ сделать это должно выглядеть на (идеально постоянным) групповая задержка и видеть, насколько плоский это.

Сравните групповую задержку четырех БИХ-фильтров, спроектированных выше.

  • Если фаза является проблемой, имейте в виду, что Баттерворт и Чебышевские проекты Типа II имеют самую плоскую групповую задержку и таким образом вводят наименьшее количество искажения.

hfvt = fvtool(dbutter,dcheby1,dcheby2,dellip,'Analysis','grpdelay');
legend(hfvt,'Butterworth', 'Chebyshev Type I',...
  'Chebyshev Type II','Elliptic')

Figure Filter Visualization Tool - Group delay contains an axes and other objects of type uitoolbar, uimenu. The axes with title Group delay contains 4 objects of type line. These objects represent Butterworth, Chebyshev Type I, Chebyshev Type II, Elliptic.

Заключения

В этом примере вы изучили, как использовать designfilt получить разнообразие КИХ lowpass и БИХ фильтрует с различными ограничениями и методами разработки. designfilt может также использоваться, чтобы получить highpass, полосу пропускания, bandstop, произвольную величину, дифференциатор и Гильбертовы проекты. Узнать больше обо всех доступных параметрах, галерее Filter Design.

Дополнительные материалы для чтения

Для получения дополнительной информации о создании фильтра и анализе, см. документацию программного обеспечения Signal Processing Toolbox®. Для получения дополнительной информации о фильтре приложения видят Практическое Введение в Цифровую Фильтрацию.

Смотрите также

| |

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