Приложение 2

П.2.1 Листинг программы синтеза цифрового ФНЧ Баттерворта

в системе MATLAB

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

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

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

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

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

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

>> Fd=5000;

>> fp=1000;

>> fs=1500;

>> Rp=3;

>> Rs=40;

>> [n,Wn]=buttord(2*fp/Fd,2*fs/Fd,Rp,Rs)

>> [a,b]=butter(n,Wn);

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

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

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

>> gridon

П.2.2 Листинг программы синтеза цифрового ФВЧ Чебышева типа 1

в системе MATLAB

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

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

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

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

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

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

>> Fd=5000;

>> fp=1000;

>> fs=667;

>> Rp=1.5;

>> Rs=40;

>> [n,Wn]=cheb1ord(2*fp/Fd,2*fs/Fd,Rp,Rs)

>> [a,b]=cheby1(n,Rp,Wn,’high’);

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

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

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

>> gridon

П.2.3 Листинг программы синтеза цифрового полосового фильтра

Чебышева типа 2 в системе MATLAB

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

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

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

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

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

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

>> Fd=70000;

>> fp1=20000;

>> fp2=22000;

>> fs1=19300;

>> fs2=22700;

>> Rp=1.5;

>> Rs=40;

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

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

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

>> [a,b]=cheby2(n,Rs,Wn);

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

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

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

>> gridon