figure(6)
hold on;

%WsNew=373.92;
%WpNew=3450.15;
Wp1 = [511 3233];
Ws1 = [740 1731];


[N2,Wn2]=buttord(Wp1,Ws1,Rp,Rs,'s')
[Ns,Ds]=butter(N2,Wn2,'s')

printsys(Ns,Ds)%принтира Ns и Ds

[Ts,ws]=freqs(Ns,Ds);%изчислява аналоговата предавателна функция
ma1=-20*log10(abs(Ts));%затихването в децибели

plot(ws,ma1,'g--'); %крива на затихване на аналогов ЛФ
plot(w*Fs/(2*pi),ma,'r-.');%крива на затихване на цифров ЛФ

ylabel('a,dB');
xlabel('f,Hz');
title('Криви на затихването на аналогов ЛФ и цифров ЛФ');
legend('Крива на затихване на аналогов ЛФ','Крива на затихване на цифров ЛФ');

plot([0 500],[25 25]);
plot([500 500],[0 25]);
plot([707 1414],[3 3]);
plot([2000 2000],[0 25]);
plot([2000 Fn],[25 25]);
axis([0 Fn 0 45]);
xlabel('f,Hz');
ylabel('a,dB');
title('Габарит на ЛФ');