Для извлечения фазовой характеристики фильтра доступны функции MATLAB ®. Учитывая частотный отклик, функцияabs возвращает величину и angle возвращает фазовый угол в радианах. Просмотр величины и фазы фильтра Баттерворта с помощью fvtool:
d = designfilt('lowpassiir','FilterOrder',9, ... 'HalfPowerFrequency',400,'SampleRate',2000); fvtool(d,'Analysis','freq')
Также можно нажать кнопку «Величина и фазовый отклик» на панели инструментов или выбрать «Анализ» > «Величина и фазовый отклик» для отображения графика.
unwrap функция также полезна в частотном анализе. unwrap разворачивает фазу, чтобы сделать ее непрерывной при разрывах фаз на 360 ° путем добавления кратных ± 360 °, при необходимости. Посмотреть, какunwrap полезно разработать фильтр FIR нижних частот 25-го порядка:
h = fir1(25,0.4);
Получение частотной характеристики с помощью freqz и постройте график фазы в градусах:
[H,f] = freqz(h,1,512,2); plot(f,angle(H)*180/pi) grid
Трудно различить скачки 360 ° (артефакт арктангенс-функции внутриangle) от скачков на 180 °, которые означают нули в частотной характеристике.
unwrap устраняет скачки на 360 °:
plot(f,unwrap(angle(H))*180/pi)
Кроме того, можно использовать phasez для просмотра развернутой фазы:
phasez(h,1)
