% =================================================================
% DIPOLO ELEMENTAL CON UN UNICO CERO
% =================================================================

clear;
clf;

rho = 2; % 0.5         % defino el módulo del cero
Wo  = -pi;             % defino la fase del cero 
Zo  = rho*exp(-1i*Wo); % defino el cero

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




for i = 1:N
  H(i) = Z(i)-Zo; % defino la respuesta de frecuencia
end

AMP =  abs(H);   % espectro de amplitud
FAS = -angle(H); % espectro de fase [-pi,+pi]

%%FAS = unwrap(FAS); %  espectro de fase  en [0,2pi]



% Preparo los tics para los gráficos
xt      = -pi:pi/4:pi;
xtlabel = ['{/Symbol -p}';'{/Symbol -3p/4}';'{/Symbol -p/2}';'{/Symbol -p/4}';...
           '0';'{/Symbol p/4}';'{/Symbol p/2}';'{/Symbol 3p/4}';'{/Symbol p}'];
yt      = [-2*pi -pi -1 0 1 pi 2*pi];
ytlabel = ['-{/Symbol 2p}';'-{/Symbol p}';'-1';'0';'1';'{/Symbol p}';'{/Symbol 2p}'];

plot(W,AMP,'r-;Amplitud;',"linewidth",2); hold on
plot(W,FAS,'b:;Fase;',"linewidth",2); hold off
title('Filtro con un solo cero',"fontsize", 14)
xlabel('Frecuencia {/Symbol w}',"fontsize",14)
axis([-pi pi -2*pi-0.5 2*pi+0.5])
set(gca,"xtick",xt,"xticklabel",xtlabel,"ytick",yt,"yticklabel",ytlabel);

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