Лабы / L 4 / L4_1 / blase

N=6; Vs=0.5;

[b1,a1]= butter(N,Vs);

[b2,a2]=cheby1(N,0.5,Vs);

[b3,a3]=ellip(N,0.5,30,Vs);

h1=abs(freqz(b1,a1,256));

h2=abs(freqz(b2,a2,256));

h3=abs(freqz(b3,a3,256));

subplot(311),plot(h1);

title('butter() фильтр Баттерворда','FontName','courier');

subplot(312),plot(h2);

title('cheby1() филтр Чебышева','FontName','courier');

subplot(313),plot(h3);

title('ellip() эллиптический филтр','FontName','courier');

pause; clg;

an=angle(freqz(b1,a1,256));

subplot(311),plot(angle(freqz(b1,a1,256)));

title('Фазовый сдвиг передаточной характеристики фильтра

Баттерворда','FontName','courier');

subplot(312),plot(angle(freqz(b2,a2,256)));

title('Фазовый сдвиг передаточной характеристики филтра

Чебышева','FontName','courier');

subplot(313),plot(angle(freqz(b3,a3,256)));

title('Фазовый сдвиг передаточной характеристики эллиптического

филтра','FontName','courier');

N=50; Vs=0.5;

b1=FIR1(N,Vs,boxcar(N+1));

b2=FIR1(N,Vs,bartlett(N+1));

b3=FIR1(N,Vs,hamming(N+1));

a=[]; a=[1 a];

h1=abs(freqz(b1,a,256));

h2=abs(freqz(b2,a,256));

h3=abs(freqz(b3,a,256));

subplot(311),plot(h1);

title('Прямоугольное окно (Boxcar)','FontName','courier');

subplot(312),plot(h2);

title('Треугольное окно (Bartlett)','FontName','courier');

subplot(313),plot(h3);

title('Косинусное окно (Hamming)','FontName','courier');

pause; clg;

plot(h1,'r');

hold;

plot(h2,'g');

plot(h3,'b');