clear;
% Especificaciones
wp = 2*pi*[1500 3500];
ws = 2*pi*[500 5000];
Ap = 1;As = 60;

% Pasar las especificaciones a un filtro pasobajo
% Tablas 1 y 2,
Bw = wp(2)-wp(1);
wx2 = wp(1)*wp(2);
if ws(1)*ws(2)<wx2,
	ws(1) = wx2/ws(2);
else
	ws(2) = wx2/ws(1);
end
wx = sqrt(wx2);

% Las especificaciones del prototipo de filtro pasobajo son
wplpp = Bw;
wslpp = ws(2)-ws(1);

% Determinación del orden del filtro, de e2 y de v3
e2 = 10^(Ap*0.1)-1;
n = ceil((log10(sqrt((10^(0.1*As)-1)/e2)))/log10(wslpp/wplpp));
v3 =  1/(e2^(1/(2*n)));

% Función de transferencia del prototipo de filtro pasobajo
[Z,P,K] = buttap(n);
[Blpp,Alpp] = zp2tf(Z,P,K);

% Pasar a filtro pasabanda con los parámetros wx2 y Bw
% El valor de Bw hay que multiplicarlo por v3 para desnormalizar
[B,A] = lp2bp(Blpp,Alpp,wx,Bw*v3);

[H,W] = freqs(B,A,500);
Hdb=20*log10(abs(H));
plot(W/(2*pi),Hdb);hold;
plot([wp ws]/(2*pi),-[Ap Ap As As],'*r');zoom;grid;hold off;