n=[0:1:9999]';
u=cos(2*pi*0.01*n)+cos(2*pi*0.02*n)+...
    cos(2*pi*0.03*n); % entrada
%=======================================
% Respuesta en frecuencia
% H1(w)= 2*exp(-j*pi/4) en banda de paso
% Fase constante en banda de paso
%=======================================
y1=2*cos(2*pi*0.01*n-pi/4)+2*cos(2*pi*0.02*n-pi/4)+...
    2*cos(2*pi*0.03*n-pi/4); % salida
figure(1)
subplot(211)
plot(n,u)
axis([4000 4500 -4 4])
title('Entrada')
subplot(212)
plot(n,y1/2)
axis([4000 4500 -4 4])
title('Salida')
%=======================================
% Respuesta en frecuencia del filtro
% H2(w)= 2*exp(-j*w*5) en banda de paso
% Fase proporcional a w, en banda de paso
%=======================================
y2=2*cos(2*pi*0.01*(n-5))+2*cos(2*pi*0.02*(n-5))+...
    2*cos(2*pi*0.03*(n-5)); % salida
figure(2)
subplot(211)
plot(n,u)
axis([4000 4500 -4 4])
title('Entrada')
subplot(212)
plot(n,y2/2)
axis([4000 4500 -4 4])
title('Salida')