isstable

Определите, является ли фильтр стабильным

Синтаксис

flag = isstable(b,a) flag = isstable(sos) flag = isstable(d)

Описание

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

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

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

Примеры

свернуть все

Стабильность фильтра

Открыть Live Script

Разработайте фильтр Butterworth БИХ шестого порядка с использованием разделов второго порядка. Задайте нормированную частоту 3-dB 0,7. Определите, является ли фильтр стабильным.

[z,p,k] = butter(6,0.7,'high');
SOS = zp2sos(z,p,k);    
flag = isstable(SOS)

flag = logical
   1

zplane(z,p)

Figure contains an axes. The axes contains 4 objects of type line, text.

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

d = designfilt('highpassiir','DesignMethod','butter','FilterOrder',6, ...
               'HalfPowerFrequency',0.7);
dflg = isstable(d)

dflg = 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 = [1 -0.5];
a = [1 -0.999999999];
act_flag1 = isstable(b,a)

act_flag1 = logical
   1

act_flag2 = isstable(single(b),single(a))

act_flag2 = logical
   0

См. также

Введенный в R2013a

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