% Mostrar la realción entre frecuencias de los dominios analógico
% y discreto en la transformación bilineal

clear;close;
N=500;
w=(-N/2:N/2);
color = ['r' 'b' 'm' 'g'];
i=1;
for fm=[1 10 50 250]
  W=2*atan(w/(2*fm));
  plot(w,W,color(i));hold on;
  gtext({['fm='num2str(fm) 'Hz']},'FontName','Lucida','Fontsize', ...
	10);
  i=i+1;
end
grid;
print -depsc2 warping;