%% =================================================================
%% Operador promedio simple
%% =================================================================

clear all; %  Limpiar todas las variables definidas

% Definimos el vector frecuencia angular digital w:
w = (0:pi/100:2*pi) ;

%___________________________________________________________________

% RESPUESTA EN FRECUENCIA

H = 1/3 + (2/3)*cos(w) ;

% Respuesta en frecuencia
figure(1); clf;
subplot(3,1,1)
plot(w,H,'b')
hold on
cero(length(w))=0; % graficar el eje horizontal para observar el cambio de signo
plot(w,cero,'k')
hold off
axis([0 2*pi -0.5 1])
title('Respuesta en frecuencia')
ylabel('H(\omega)');
xt      = [0; pi; 2*pi];
xtlabel = [ '0' ; '{/Symbol p}' ; '2{/Symbol p}' ];
set(gca,'xtick',xt,'xticklabel',xtlabel);

% Espectro de amplitud
subplot(3,1,2)
plot(w,abs(H),'r')
title('Espectro de amplitud')
ylabel('|H(\omega)|');
axis([0 2*pi 0 1])
x1 = acos(-0.5);
x2 = 2*pi - x1;
x1tick = num2str(x1,'%4.2f');
x2tick = num2str(x2,'%4.2f');
xt      = [ 0:pi:2*pi x1 x2];
xtlabel = ['0' ; '{/Symbol p}' ; '2{/Symbol p}' ; x1tick ; x2tick ];
set(gca,'xtick',xt,'xticklabel',xtlabel);

% Espectro de fase
subplot(3,1,3)
plot(w,angle(H),'g')
title('Espectro de fase')
xlabel('\omega');
ylabel('\Phi(\omega)');
axis([0 2*pi -pi*0.5 pi*1.5]);
xt      = [ 0:pi:2*pi x1 x2];
xtlabel = ['0' ; '{/Symbol p}' ; '2{/Symbol p}' ; x1tick ; x2tick ];
yt      = 0:pi:pi;
ytlabel = [ '0' ; '{/Symbol p}' ];
set(gca,'xtick',xt,'xticklabel',xtlabel,'ytick',yt,'yticklabel',ytlabel);

% Guardar grafico en un postscript
print('promedio_simple.eps','-depsc'); % -depsc: device eps color

% Convertimos a pdf llamando el sistema
system('epstopdf promedio_simple.eps');
