Стабилизировать полином
b = polystab(a)
polystab стабилизирует полином относительно единичной окружности; она отражает корни с величинами более 1 внутри единичной окружности.
b = polystab(a) возвращает вектор строки b содержащий стабилизированный полином. a - вектор полиномиальных коэффициентов, обычно в z-области:
.. + a (m + 1) z − m.
Используйте оконный метод для проектирования КИХ-фильтра 25-го порядка с нормированной частотой отсечения рад/выборка. Убедитесь, что она имеет линейную фазу, но не минимальную фазу.
h = fir1(25,0.4); h_linphase = islinphase(h)
h_linphase = logical
1
h_minphase = isminphase(h)
h_minphase = logical
0
Использовать polystab преобразование фильтра линейной фазы в фильтр минимальной фазы. Постройте график фазовых откликов фильтров.
hmin = polystab(h)*norm(h)/norm(polystab(h)); hmin_linphase = islinphase(hmin)
hmin_linphase = logical
0
hmin_minphase = isminphase(hmin)
hmin_minphase = logical
1
hfvt = fvtool(h,1,hmin,1,'Analysis','phase'); legend(hfvt,'h','hmin')
Убедитесь, что два фильтра имеют одинаковую амплитудную характеристику.
hfvt = fvtool(h,1,hmin,1); legend(hfvt,'h','hmin')
polystab находит корни многочлена и сопоставляет те корни, которые находятся вне единичной окружности, с внутренней частью единичной окружности:
v = roots(a); vs = 0.5*(sign(abs(v)-1)+1); v = (1-vs).*v + vs./conj(v); b = a(1)*poly(v);