islinphase
Определите, имеет ли фильтр линейную фазу
Синтаксис
flag = islinphase(b,a)
flag = islinphase(sos)
flag = islinphase(d)
flag = islinphase(...,tol)
Описание
flag = islinphase(b,a)flag, равняйтесь true если коэффициенты фильтра в b и a задайте линейный фильтр фазы. flag равно false если фильтр не имеет линейной фазы.
flag = islinphase(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 = islinphase(d)true если цифровой фильтр, d, имеет линейную фазу. Использование designfilt сгенерировать d на основе технических требований частотной характеристики.
flag = islinphase(...,tol)tol, определить, когда два числа достаточно близки, чтобы быть рассмотренными равными. Если не заданный, tol, значения по умолчанию к eps^(2/3).
Примеры
Линейная и нелинейная фаза
Используйте метод окна, чтобы спроектировать десятый порядок КИХ-фильтр lowpass с нормированной частотой среза 0.55. Проверьте, что фильтр имеет линейную фазу.
d = designfilt('lowpassfir','DesignMethod','window', ... 'FilterOrder',10,'CutoffFrequency',0.55); flag = islinphase(d)
flag = logical
1
[phs,w] = phasez(d); plot(w/pi,phs) xlabel('Frequency \omega/\pi') ylabel('Phase')
БИХ-фильтры в целом не имеют линейной фазы. Проверьте оператор путем построения восьмого порядка Баттерворт, Чебышев, и эллиптические фильтры с подобными техническими требованиями.
ord = 8; Wcut = 0.35; atten = 20; rippl = 1; [zb,pb,kb] = butter(ord,Wcut); sosb = zp2sos(zb,pb,kb); [zc,pc,kc] = cheby1(ord,rippl,Wcut); sosc = zp2sos(zc,pc,kc); [zd,pd,kd] = cheby2(ord,atten,Wcut); sosd = zp2sos(zd,pd,kd); [ze,pe,ke] = ellip(ord,rippl,atten,Wcut); sose = zp2sos(ze,pe,ke);
Постройте фазовые отклики фильтров. Определите, есть ли у них линейная фаза.
fv = fvtool(sosb,sosc,sosd,sose,'Analysis','phase'); legend(fv,'Butterworth','Chebyshev I','Chebyshev II','Elliptic')
phs = [islinphase(sosb) islinphase(sosc) ...
islinphase(sosd) islinphase(sose)]phs = 1x4 logical array
0 0 0 0
Смотрите также
designfilt | digitalFilter | isallpass | ismaxphase | isminphase | isstable