П.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

Литература

1. Гадзиковский В.И. Основы теории и проектирования цифровых фильтров: Учебное пособие для радиотехнических специальностей вузов. — М.: высшая школа, 1996. — 256 с.: ил.

2. Гадзиковский В.И.Методы цифрового моделирования радиотехнических систем: Учебное пособие. — Екатеринбург: УГТУ, 1995. — 182 с.: ил.

3. Изучение методов проектирования избирательных рекурсивных цифровых фильтров: Методические указания к лабораторной работе по курсу цифровой обработки сигналов / В.И. Гадзиковский. — Екатеринбург: Уральский филиал СибГУТИ, 2001. — 20 с.

1 Они отличаются от соответствующих функций проектирования аналоговых фильтров синтаксисом.

48