isminphase

Определите, является ли фильтр минимальной фазой

Синтаксис

flag = isminphase(b,a)
flag = isminphase(sos)
flag = isminphase(d)
flag = isminphase(...,tol)
flag = isminphase(h)

Описание

Фильтр является минимальной фазой, когда все нули ее передаточной функции включены или в модульном кругу, или числитель является скаляром. Эквивалентным определением для минимального фильтра фазы является причинная и устойчивая система с причинной и устойчивой инверсией.

flag = isminphase(b,a) возвращает логический выходной параметр, flag, равняйтесь true если фильтр, заданный коэффициентами числителя, b, и коэффициенты знаменателя, a, минимальный фильтр фазы.

flag = isminphase(sos) возвращает true если фильтр, заданный вторым порядком, разделяет матрицу, sos, минимальная фаза. sos K-by-6 матрица, где количество разделов, K, должно быть больше или быть равно 2. Каждая строка sos соответствует коэффициентам второго порядка (biquad) фильтр. i th строка sos матрица соответствует [bi(1) bi(2) bi(3) ai(1) ai(2) ai(3)].

flag = isminphase(d) возвращает true если цифровой фильтр, d, имеет минимальную фазу. Используйте designfilt сгенерировать d на основе технических требований частотной характеристики.

flag = isminphase(...,tol) использует допуск, tol, определить, когда два числа достаточно близки, чтобы быть рассмотренными равными. Если не заданный, tol, значения по умолчанию к eps^(2/3).

 Детали для вычислений с фиксированной точкой

flag = isminphase(h) определяет если dfilt объект фильтра h минимальная фаза.

Примеры

свернуть все

Спроектируйте шестой порядок lowpass БИХ-фильтр Баттерворта с помощью разделов второго порядка. Задайте нормированную частоту на 3 дБ 0,15. Проверяйте, имеет ли фильтр минимальную фазу.

[z,p,k] = butter(6,0.15);
SOS = zp2sos(z,p,k);            
min_flag = isminphase(SOS)
min_flag = logical
   1

Перепроектируйте фильтр с помощью designfilt. Проверяйте, что нули и полюса передаточной функции включены или в модульном кругу.

d = designfilt('lowpassiir','DesignMethod','butter','FilterOrder',6, ...
               'HalfPowerFrequency',0.25);
d_flag = isminphase(d)
d_flag = logical
   1

zplane(d)

Учитывая фильтр, заданный с набором числителя с одинарной точностью и коэффициентов знаменателя, проверяйте, имеет ли он минимальную фазу для различных значений допуска.

b = single([1 1.00001]);  
a = single([1 0.45]);               
min_flag1 = isminphase(b,a)        
min_flag1 = logical
   0

min_flag2 = isminphase(b,a,1e-3)
min_flag2 = logical
   1

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

| | | | |

Введенный в R2013a

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