%% =================================================================
%% RAÍCES DE ONDÍCULAS EQUIVALENTES A PARTIR DE LA AUTOCORRELACIÓN
%% =================================================================

clear all; close;

% Autocorrelación y raíces
x = [ 0.125, -0.84375, 1.578125, -0.84375, 0.125 ];
r = roots(x)

% Fase mínima
a  = conv([1,-1/r(1)],[1,-1/r(2)]);
ka = sqrt(x(1)/xcorr(a)(1)); % factor de normalización
a  = a*ka                    % normalizo

% Fase mixta 
b  = conv([1,-1/r(1)],[1,-1/r(3)]);
kb = sqrt(x(1)/xcorr(b)(1)); % factor de normalización
b  = b*kb                    % normalizo

% Fase mixta
c  = conv([1,-1/r(4)],[1,-1/r(2)]);
kc = sqrt(x(1)/xcorr(c)(1)); % factor de normalización
c  = c*kc                    % normalizo

% Fase máxima 
d  = conv([1,-1/r(4)],[1,-1/r(3)]);
kd = sqrt(x(1)/xcorr(d)(1)); % factor de normalización
d  = d*kd                    % normalizo

% Autocorrelaciones
aa = xcorr(a)
bb = xcorr(b)
cc = xcorr(c)
dd = xcorr(d)
