График нулей и полюсов для систем в дискретном времени
zplane( строит графики нулей, заданных в вектор-столбец z,p)z и полюса, указанные в вектор-столбец p в текущую фигуру окне. Символ 'o' представляет нуль и символ 'x' представляет собой полюс. График включает модуль круг для ссылки.
Если z и p являются матрицами, тогда zplane строит графики полюсов и нулей в столбцах z и p в разных цветах.
zplane( находит нули и полюса передаточной функции, представленной цифровым фильтром, d)d. Использовать designfilt чтобы сгенерировать d на основе спецификаций частотной характеристики. Эта диаграмма нулей и полюсов отображается в FVTool.
Для данных, дискретизированных с частотой 1000 Гц, постройте графики полюсов и нулей эллиптического lowpass цифрового фильтра 4-го порядка с частотой отсечения 200 Гц, 3 дБ пульсации в полосе пропускания и 30 дБ ослабления в полосе пропускания.
[z,p,k] = ellip(4,3,30,200/500);
zplane(z,p)
grid
title('4th-Order Elliptic Lowpass Digital Filter')
Создайте тот же фильтр с помощью designfilt. Использование zplane для построения графика полюсов и нулей. Обратите внимание, что этот синтаксис zplane вызывает fvtool.
d = designfilt('lowpassiir','FilterOrder',4,'PassbandFrequency',200, ... 'PassbandRipple',3,'StopbandAttenuation',30, ... 'DesignMethod','ellip','SampleRate',1000); zplane(d)
Проектируйте полосно-пропускной фильтр Чебышевского типа II 8-го порядка с затуханием в полосе задерживания 20 дБ. Задайте частоты ребра полосы упора как рад/образец и рад/образец.
[b,a] = cheby2(8/2,20,[1 5]/8);
Использование zplane для построения графика полюсов и нулей передаточной функции.
zplane(b,a)
Визуализируйте нулевую фазовую характеристику фильтра. Наложите модуль круг, полюс и нулевые положения.
[hw,fw] = zerophase(b,a,1024,"whole"); z = roots(b); p = roots(a); plot3(cos(fw),sin(fw),hw) hold on plot3(cos(fw),sin(fw),zeros(size(fw)),'--') plot3(real(z),imag(z),zeros(size(z)),'o') plot3(real(p),imag(p),zeros(size(p)),'x') hold off xlabel("Real") ylabel("Imaginary") view(35,40) grid
Вы можете переопределить автоматическое масштабирование zplane использование
axis([xmin xmax ymin ymax])
после вызова zplane. Это масштабирование полезно, когда один или несколько нулей или полюсов имеют настолько большую величину, что другие тесно сгруппированы вокруг источника и их трудно различить.