Проверьте, что Системный объект фильтра в дискретном времени является стабильным
Спроектируйте фильтр БИХ типа I Чебышева и определите, имеет ли фильтр минимальную фазу и является ли он стабильным.
Использование fdesign.lowpass и design функции, проектирование БИХ-фильтра типа I Чебышева с неравномерностью в полосе пропускания 0,5 дБ и частотой отсечения 3 дБ при 9600 Гц.
Fs = 48000; % Sampling frequency of input signal d = fdesign.lowpass('N,F3dB,Ap', 10, 9600, .5, Fs); filt = design(d,'cheby1','Systemobject',true)
filt =
dsp.BiquadFilter with properties:
Structure: 'Direct form II'
SOSMatrixSource: 'Property'
SOSMatrix: [5x6 double]
ScaleValues: [6x1 double]
InitialConditions: 0
OptimizeUnityScaleValues: true
Show all properties
Использование isminphase function, определить, имеет ли фильтр минимальную фазу.
isminphase(filt)
ans = logical
1
Проверьте местоположение полюсов и нули передаточной функции фильтра на плоскости z. По определению, полюса и нули минимальной фазы фильтра должны быть на или внутри модуля круга.
zplane(filt)
Все минимальные фильтры фазы стабильны. Чтобы проверить, является ли разработанный фильтр стабильным, используйте isstable функция.
isstable(filt)
ans = logical
1