//Autor: Pedro Ledoux

//Gabarito:
//Filtro passa altas Chebyshev
//Inicio da banda de transição: 180kHz
//Inicio da banda de passagem: 450kHz
//Máxima atenuação na banda de passagem: -0.5dB
//Minima atenuação na banda de rejeição: -60dB

s=poly(0,'s'); //Declaração da incognita s da qual podemos escrever polinômios
z=poly(0,'z'); //Declaração da incognita z.

a=11.177965;//epsulon multiĺicado por 2 elevado a n-1=5, termo o qual aparece na função de Chebyshev no denominador.

s=2*%pi*450000/z;  //Declarando s em função de z de forma que a transformação de frequência é realizada

//Polinômio em s da função de transferência do passa baixas normatizado. No abbiente scilab esse polinômio será mostrado de forma que s será substituido por z.
p1=1/( a*(s^2+s*0.1553+1.0230228)*(s^2+s*0.4242879+0.5900)*(s^2+s*0.5795881+0.1569974))

//O Scilab possui a função de plotagem bode() que nos permite plotar o grafico do polinômio em diagrama de Bode disponibilizando simultaneamente o de modulo e o de fase.
p=syslin('c',p1) ;
bode(p,1000,100000000);
pp1=1/(s^2+s*0.1553+1.02302)
pp2=1/(s^2+s*0.4242879+0.5900)
pp3=1/(a*(s^2+s*0.5795881+0.15699))

pp=pp1*pp2*pp3