isminphase

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

Синтаксис

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

Описание

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

flag = isminphase(sos) прибыль true если фильтр задан матрицей секций второго порядка, sos, является минимальной фазой. sos является K-by-6 матрицей, где число сечений, K, должно быть больше или равно 2. Каждая строка sos соответствует коэффициентам фильтра второго порядка (биквад). I-й ряд 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).

Примеры

свернуть все

Фильтры минимальной фазы

Открыть сценарий в реальном времени

Проектирование низкочастотного БИХ-фильтра Butterworth шестого порядка с использованием секций второго порядка. Укажите нормированную частоту 3-dB, равную 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)

Figure Filter Visualization Tool - Pole-Zero Plot contains an axes and other objects of type uitoolbar, uimenu. The axes with title Pole-Zero Plot contains 4 objects of type line, text.

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

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

Документация