[N,W]=cheb1ord(FC/10000,FR/10000,R,A);
[B,A]=cheby1(N,R,W);
[H,F]=freqz(B,A,512,20000);
modulodB=20*log10(abs(H));
fase=unwrap(angle(H));
Fc=W*10000;