КИХ полуполосы фильтрует проект
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) разрабатывает lowpass фильтр полуполосы FIR порядка n с equiripple характеристикой. n должен быть ровным целым числом. fp определяет частоту ребра полосы пропускания, и должно удовлетворить 0 <fp <1/2, где 1/2 соответствует π/2 рад/выборка.
b = firhalfband(n,win) разрабатывает lowpass фильтр N-го порядка с помощью усеченного, метода ответа оконного импульса вместо equiripple метода. win является вектором длины n+1. Идеальный импульсный ответ является усеченным к длине, n + 1, и затем умноженный детально с окном задал in win.
b = firhalfband(n,dev,'dev') разрабатывает lowpass полуленточный фильтр N-го порядка с equiripple характеристикой. Входной параметр dev устанавливает значение для максимальной полосы пропускания и позволенной пульсации полосы задерживания.
b = firhalfband('minorder',fp,dev) разрабатывает lowpass фильтр минимального заказа, с ребром полосы пропускания fp. Пиковая пульсация ограничивается скалярным dev. Этот проект использует equiripple метод.
b = firhalfband('minorder',fp,dev,'kaiser') разрабатывает lowpass фильтр минимального заказа, с ребром полосы пропускания fp. Пиковая пульсация ограничивается скалярным dev. Этот проект использует метод окна Kaiser.
b = firhalfband(...,'high') возвращает highpass фильтр полуполосы FIR.
b = firhalfband(...,'minphase') разрабатывает КИХ-фильтр минимальной фазы, таким образом, что фильтр является спектральным фактором фильтра полуполосы (вспомните, что h = conv(b,fliplr(b)) является фильтром полуполосы). Это может быть полезно для разработки совершенной реконструкции, двухканальных КИХ-наборов фильтров. Опция minphase для firhalfband не доступна для проектов полуленточного фильтра оконных — b = firhalfband(n,win) и b = firhalfband('minorder',fp,dev,'kaiser').
В минимальных случаях фазы порядок фильтра должен быть нечетным.
halfbandЭтот пример разрабатывает полуленточный фильтр минимального заказа с заданной максимальной пульсацией.
b = firhalfband('minorder',.45,0.0001);
impz(b)
Вы видите, что импульсный ответ является нулем для каждой альтернативной выборки.
[1] Saramaki, T, “конечный импульсный проект фильтра ответа”, руководство для цифровой обработки сигналов. С.К. Митра и Дж.Ф. Wiley-межнаука редакторов кайзера, Нью-Йорк, 1993, глава 4.