polystab

Стабилизируйте полином

Синтаксис

b = polystab(a)

Описание

polystab стабилизирует полином относительно модульного круга; это отражает корни с величинами, больше, чем 1 внутренняя часть модульный круг.

b = polystab(a) возвращает вектор-строку b содержа стабилизированный полином. a вектор с полиномиальными коэффициентами, обычно в z-области:

A(z)=a(1)+a(2)z1++a(m+1)zm.

Примеры

свернуть все

Используйте метод окна, чтобы спроектировать КИХ-фильтр 25-го Одера с нормированной частотой среза 0.4π рад/отсчет. Проверьте, что это имеет линейную фазу, но не минимальную фазу.

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')

Figure Filter Visualization Tool - Phase Response contains an axes and other objects of type uitoolbar, uimenu. The axes with title Phase Response contains 2 objects of type line. These objects represent h, hmin.

Проверьте, что два фильтра имеют идентичные ответы величины.

hfvt = fvtool(h,1,hmin,1);
legend(hfvt,'h','hmin')

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. These objects represent 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);

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

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