%% =================================================================
%% DIPOLOS CON UN ÚNICO CERO DE FASE MÍNIMA Y MÁXIMA
%% =================================================================

clear;
clf;

W = -pi:pi/100:pi; % vector frec. ang. dig.
Z = exp(-1i*W);    % vector Z para evaluar H(Z)

% Z0 afuera del círculo unitario |Z|=1
Z0   = -2;
H1   =  Z - Z0;
AMP1 =  abs(H1);
FAS1 =  angle(H1); % exp(-iw) con w>0 en sentido horario

% Z0 adentro del círculo unitario |Z|=1
Z0   = -0.5;
H2   =  Z - Z0;
AMP2 =  abs(H2);
FAS2 =  angle(H2);

% Preparo los tics para los gráficos
xt         = -pi:pi/2:pi;
xtlabel    = ['{/Symbol -p}';'{/Symbol -p/2}';'0';'{/Symbol p/2}';'{/Symbol p}'];
ytAMP      = 0:4;
ytFAS      = -pi:pi:pi;
ytFASlabel = ['{/Symbol -p}';'0';'{/Symbol p}'];

subplot(2,2,1)
plot(W,AMP1,'b-')
title('Z_0 = -2')
ylabel('Espectro de amplitud')
axis([-pi pi 0 4])
set(gca,"xtick",xt,"xticklabel",xtlabel,"ytick",ytAMP);

subplot(2,2,2)
plot(W,AMP2,'r-')
title('Z_0 = -0.5')
ylabel('Espectro de amplitud')
axis([-pi pi 0 4])
set(gca,"xtick",xt,"xticklabel",xtlabel,"ytick",ytAMP);

subplot(2,2,3)
plot(W,FAS1,'b-')
xlabel('Frecuencia')
ylabel('Espectro de fase')
axis([-pi pi -pi pi])
set(gca,"xtick",xt,"xticklabel",xtlabel,"ytick",ytFAS,"yticklabel",ytFASlabel);

subplot(2,2,4)
plot(W,FAS2,'r-')
xlabel('Frecuencia')
ylabel('Espectro de Fase')
set(gca,"xtick",xt,"xticklabel",xtlabel,"ytick",ytFAS,"yticklabel",ytFASlabel);
axis([-pi pi -pi pi])

%print('dipolos.eps','-depsc'); % guardo el gráfico
