Лабораторная работа: Расчет цифровых фильтров с бесконечными импульсными характеристиками
[b,a]=zp2tf(z,p,k);
figure;
zplane(z,p);
figure;
freqs(z,p);
%%%%%%%%%%%%Фильтр Чебышева 2 рода%%%%%%%%%%
[n3, Wn3] = cheb2ord(Wp, Ws, Rp, Rs);
[z,p,k] = cheby2(n3,Rs,Wn3,'High');
[b,a]=zp2tf(z,p,k);
figure;
zplane(z,p);
figure;
freqs(z,p);
%%%%%%%%%%%Фильтр Эллиптический%%%%%%%%%%%%%
[n4, Wn4] = ellipord(Wp, Ws, Rp, Rs);
[z,p,k] = ellip(n4,Rp,Rs,Wn4,'high');
[b,a]=zp2tf(z,p,k);
figure;
zplane(z,p);
figure;
freqs(z,p);
sprintf( 'chebysheva1 - %d',n2)
sprintf( 'chebysheva2 - %d',n3)
sprintf( 'elipticheskiy - %d',n4)
Рисунок 7– Диаграмма нулей и полюсов фильтра Баттерворта
.
Рисунок 8– АЧХ , ФЧХ фильтра Баттерворта