В этом примере показано, как спроектировать классические БИХ-фильтры. Начальное особое внимание находится на ситуации, для которой критический расчетный параметр является частотой среза, в которой степень фильтра затухает к половине (-3 дБ) номинальное значение полосы пропускания.
Пример иллюстрирует, как легкий это должно заменить проект Баттерворта или на Чебышева или на эллиптический фильтр того же порядка и получить более крутой спад за счет некоторой пульсации в полосе пропускания и/или полосе задерживания фильтра. После этого проекты минимального порядка исследуются.
Давайте спроектируем 8-й фильтр порядка с нормированной частотой среза 0.4pi. Во-первых, мы проектируем Фильтр Баттерворта, который является максимально плоским (никакая пульсация в полосе пропускания или в полосе задерживания):
N = 8; F3dB = .4; d = fdesign.lowpass('N,F3dB',N,F3dB); Hbutter = design(d,'butter','SystemObject',true)
Hbutter = dsp.BiquadFilter with properties: Structure: 'Direct form II' SOSMatrixSource: 'Property' SOSMatrix: [4x6 double] ScaleValues: [5x1 double] InitialConditions: 0 OptimizeUnityScaleValues: true Show all properties
Чебышевский Тип, который я проектирую, допускает управление пульсаций в полосе пропускания. В полосе задерживания нет все еще никаких пульсаций. Большие пульсации включают более крутой спад. Здесь, мы задаем пульсации от пика к пику 0.5 дБ:
Ap = .5; setspecs(d,'N,F3dB,Ap',N,F3dB,Ap); Hcheby1 = design(d,'cheby1','SystemObject',true); hfvt = fvtool(Hbutter,Hcheby1,'Color','white'); axis([0 .44 -5 .1]) legend(hfvt,'Butterworth','Chebyshev Type I');
Чебышевский проект Типа II допускает управление затухания в полосе задерживания. В полосе пропускания нет никаких пульсаций. Меньшее затухание в полосе задерживания включает более крутой спад. В этом примере мы задаем затухание в полосе задерживания 80 дБ:
Ast = 80; setspecs(d,'N,F3dB,Ast',N,F3dB,Ast); Hcheby2 = design(d,'cheby2','SystemObject',true); hfvt = fvtool(Hbutter,Hcheby2,'Color','white'); axis([0 1 -90 2]) legend(hfvt,'Butterworth','Chebyshev Type II');
Наконец, эллиптический фильтр может обеспечить более крутой спад по сравнению с предыдущими проектами путем разрешения пульсаций и в полосе задерживания и в полосе пропускания. Чтобы проиллюстрировать, что, мы снова используем ту же полосу пропускания и характеристику полосы задерживания как выше:
setspecs(d,'N,F3dB,Ap,Ast',N,F3dB,Ap,Ast); Hellip = design(d,'ellip','SystemObject',true); hfvt = fvtool(Hbutter,Hcheby1,Hcheby2,Hellip,'Color','white'); axis([0 1 -90 2]) legend(hfvt, ... 'Butterworth','Chebyshev Type I','Chebyshev Type II','Elliptic');
Путем увеличивания масштаб полоса пропускания мы проверяем, что все фильтры имеют ту же-3dB точку частоты и что только у Баттерворта и Чебышевских проектов Типа II есть совершенно плоская полоса пропускания:
axis([0 .44 -5 .1])
Если фаза является проблемой, полезно заметить, что Баттерворт и Чебышевские проекты Типа II - также те вводящие меньшее искажение (их групповая задержка является более плоской):
hfvt.Analysis = 'grpdelay';
В случаях, где частота среза на 3 дБ не представляет главного интереса, но вместо этого и полоса пропускания и полоса задерживания полностью заданы в терминах частот и суммы терпимых пульсаций, мы можем использовать метод проектирования минимального порядка:
Fp = .1; Fst = .3; Ap = 1; Ast = 60; setspecs(d,'Fp,Fst,Ap,Ast',Fp,Fst,Ap,Ast); Hbutter = design(d,'butter','SystemObject',true); Hcheby1 = design(d,'cheby1','SystemObject',true); Hcheby2 = design(d,'cheby2','SystemObject',true); Hellip = design(d,'ellip','SystemObject',true); hfvt = fvtool(Hbutter,Hcheby1,Hcheby2,Hellip, 'DesignMask', 'on',... 'Color','white'); axis([0 1 -70 2]) legend(hfvt, ... 'Butterworth','Chebyshev Type I','Chebyshev Type II','Elliptic');
7-й фильтр порядка необходим, чтобы выполнить спецификации с проектом Баттерворта, тогда как 5-й порядок достаточен с любой Чебышевские методы. Порядок фильтра может даже уменьшаться до 4 с эллиптическим проектом:
order(Hbutter)
ans = 7
order(Hcheby1)
ans = 5
order(Hcheby2)
ans = 5
order(Hellip)
ans = 4
С проектами минимального порядка идеальный порядок должен быть округлен до следующего целого числа. Этот дополнительный дробный порядок позволяет алгоритму на самом деле превышать технические требования. Мы можем использовать MatchExactly
отметьте к ограничению алгоритм проекта, чтобы совпадать точно с одной полосой. Другая полоса превысит свои технические требования. По умолчанию Чебышевский Тип, который я проектирую, совпадает с полосой пропускания, Баттерворт и Чебышевский Тип II совпадают с полосой задерживания, и затухания обеих полос являются соответствующими эллиптическим проектом (в то время как частота ребра полосы задерживания превышена):
Hellipmin1 = design(d, 'ellip', 'MatchExactly', 'passband',... 'SystemObject',true); Hellipmin2 = design(d, 'ellip', 'MatchExactly', 'stopband',... 'SystemObject',true); hfvt = fvtool(Hellip, Hellipmin1, Hellipmin2, 'DesignMask', 'on',... 'Color','white'); axis([0 1 -80 2]); legend(hfvt, 'Matched passband and stopband', ... 'Matched passband', 'Matched stopband', ... 'Location', 'Northeast')
Увеличение полосы пропускания, чтобы сравнить ребра полосы пропускания. Совпадающая полоса пропускания и соответствующий оба проекта имеет затухание точно 1 дБ в Fpass =.1:
axis([0 .11 -1.1 0.1]); legend(hfvt, 'Location', 'Southwest')
Мы проверяем, что получившийся порядок фильтров не изменился:
order(Hellip)
ans = 4
order(Hellipmin1)
ans = 4
order(Hellipmin2)
ans = 4
Результаты, представленные выше, могут быть расширены к highpass, полосе пропускания и bandstop типам ответа. Например, вот полосовые фильтры минимального порядка:
d = fdesign.bandpass('Fst1,Fp1,Fp2,Fst2,Ast1,Ap,Ast2', ... .35,.45,.55,.65,60,1,60);
Hbutter = design(d,'butter','SystemObject',true); Hcheby1 = design(d,'cheby1','SystemObject',true); Hcheby2 = design(d,'cheby2','SystemObject',true); Hellip = design(d,'ellip','SystemObject',true); hfvt = fvtool(Hbutter,Hcheby1,Hcheby2,Hellip, 'DesignMask', 'on',... 'Color','white'); axis([0 1 -70 2]) legend(hfvt, ... 'Butterworth','Chebyshev Type I','Chebyshev Type II','Elliptic',... 'Location', 'Northwest')