% =================================================================
% DIPOLO ELEMENTAL CON UN UNICO POLO
% =================================================================

clear;
clf;

rho = 1.1;             % defino el módulo del polo
Wp  = pi/4;            % defino la fase del polo 

Zp  = rho*exp(-1i*Wp); % defino el polo

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)   = 1./(Z(i)-Zp);
  POT(i) = abs(H(i))**2;
end

% 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}'];
max     = (1/(rho-1))**2;
yt      = 0:max/10:max;

plot(W,POT,'b-')
title('Filtro con un solo polo',"fontsize", 14)
xlabel('Frecuencia {/Symbol w}')
ylabel('Espectro de potencia')
axis([-pi pi 0 max])
set(gca,"xtick",xt,"xticklabel",xtlabel,"ytick",yt);

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