Определите, стабилен ли фильтр
flag = isstable(b,a)
flag = isstable(sos)
flag = isstable(d)
flag =
isstable(hs)
flag =
isstable(hs,'Arithmetic',arithtype)
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 th строка матрицы sos соответствует [bi(1) bi(2) bi(3) ai(1) ai(2) ai(3)].
flag = isstable(d)true, если цифровой фильтр, d, стабилен. Используйте designfilt, чтобы сгенерировать d на основе спецификаций частотной характеристики.
flag =
isstable(hs)true, если Система фильтра object™ hs стабильна. У вас должно быть программное обеспечение DSP System Toolbox™, чтобы использовать этот синтаксис.
flag =
isstable(hs,'Arithmetic',arithtype)hs на основе заданного arithtype. arithtype может быть 'double', 'single' или 'fixed'. Когда вы задаете 'double' или 'single', функция выполняет дважды - или анализ с одинарной точностью. Когда вы задаете 'fixed', арифметические изменения в зависимости от установки свойства CoefficientDataType и заблокирован ли Системный объект или разблокирован. У вас должно быть программное обеспечение DSP System Toolbox, чтобы использовать этот синтаксис.
Детали для вычислений с фиксированной точкой
Когда вы не задаете арифметику для non-CIC структур, функция использует арифметику с двойной точностью, если Системный объект фильтра находится в разблокированном состоянии. Если Системный объект заблокирован, функция выполняет анализ на основе заблокированного типа входных данных. Структуры CIC только поддерживают вычисления с фиксированной точкой.
Разработайте шестой порядок Баттерворт highpass БИХ-фильтр с помощью разделов второго порядка. Задайте нормированную частоту на 3 дБ 0,7. Определите, стабилен ли фильтр.
[z,p,k] = butter(6,0.7,'high');
SOS = zp2sos(z,p,k);
flag = isstable(SOS) flag = logical
1
zplane(z,p)
Перепроектируйте фильтр с помощью designfilt и проверяйте его на устойчивость.
d = designfilt('highpassiir','DesignMethod','butter','FilterOrder',6, ... 'HalfPowerFrequency',0.7); dflg = isstable(d)
dflg = logical
1
zplane(d)
Создайте фильтр и определите его устойчивость в двойной и одинарной точности.
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
designfilt | digitalFilter | isallpass | islinphase | ismaxphase | isminphase | zplane