П.2.4 Листинг программы синтеза эллиптического цифрового режекторного фильтра в системе matlab

Исходные данные:

• частота дискретизации кГц;

• граничные частоты полосы пропускания кГц,кГц;

• максимальное затухание в полосе пропускания дБ;

• граничные частоты полос задерживания Гц,кГц;

• минимальное затухание в полосах задерживания дБ.

>> Fd=70000;

>> fp1=20000;

>> fp2=22000;

>> fs1=20250;

>> fs2=21750;

>> Rp=1.5;

>> Rs=40;

>> Wp=[2*fp1/Fd 2*fp2/Fd];

>> Ws=[2*fs1/Fd 2*fs2/Fd];

>> [n,Wn]=ellipord(Wp,Ws,Rp,Rs)

>> [a,b]=ellip(n,Rp,Rs,Wn,’stop’);

>> [sos,g]=tf2sos(a,b,’up’,’none’)

>> [h,w]=freqz(a,b);

>> plot(w/(2*pi),20*log10(abs(h)))

>> gridon