Конструкция полутонового фильтра FIR
b = firhalfband(n,fp)
b = firhalfband(n,win)
b = firhalfband(n,dev,'dev')
b = firhalfband('minorder',fp,dev)
b = firhalfband('minorder',fp,dev,'kaiser')
b = firhalfband(...,'high')
b = firhalfband(...,'minphase')
b = firhalfband(n,fp) конструирует низкочастотный полуполосный FIR-фильтр порядка n с чёткой характеристикой. n должно быть четным целым числом. fp определяет граничную частоту полосы пропускания, и она должна удовлетворять 0 < fp < 1/2, где 1/2 соответствует δ/2 рад/образец.
b = firhalfband(n,win) создает фильтр N-го порядка нижних частот с использованием метода усеченной оконной импульсной характеристики вместо метода equiripple. win является n+1 вектор длины. Идеальная импульсная характеристика усечена до длины n + 1, а затем умножены точка за точкой на окно, указанное в win.
b = firhalfband(n,dev,'dev') конструирует фильтр нижних частот N-го порядка с характеристикой equiripple. Входной аргумент dev устанавливает значение для максимальной разрешенной полосы пропускания и пульсации полосы останова.
b = firhalfband('minorder',fp,dev) конструирует низкочастотный фильтр минимального порядка с краем полосы пропускания fp. Пиковая пульсация ограничена скаляром dev. В этой конструкции используется метод equiripple.
b = firhalfband('minorder',fp,dev,'kaiser') конструирует низкочастотный фильтр минимального порядка с краем полосы пропускания fp. Пиковая пульсация ограничена скаляром dev. В этой конструкции используется метод окна Кайзера.
b = firhalfband(...,'high') возвращает фильтр FIR high pass half band.
b = firhalfband(...,'minphase') конструирует минимально-фазовый КИХ-фильтр таким образом, что фильтр является спектральным коэффициентом полуполосного фильтра (напомним, что h = conv(b,fliplr(b)) - полуполосный фильтр). Это может быть полезно для проектирования идеальной реконструкции, двухканальных банков фильтров FIR. minphase опция для firhalfband недоступен для конструкций фильтров полуполосы на основе окон - b = firhalfband(n,win) и b = firhalfband('minorder',fp,dev,'kaiser') .
В случаях минимальной фазы порядок фильтрации должен быть нечетным.
halfbandВ этом примере создается полуполосный фильтр минимального порядка с указанной максимальной пульсацией.
b = firhalfband('minorder',.45,0.0001);
impz(b)
Видно, что импульсная характеристика равна нулю для каждого альтернативного образца.
[1] Сарамаки, Т, «Проектирование фильтра с конечной импульсной характеристикой», Руководство по цифровой обработке сигналов. С.К. Митра и Джей Ф. Кайзер Эдс. Wiley-Interscience, N.Y., 1993, глава 4.