DSP/TP1/DiogoEliseuHugo_TP1_13.m

%% Exercício 13
function [Vmax, Lmax] = DiogoEliseuHugo_TP1_13(X, Y, GRAF)
    n1 = 1:size(X, 2);
    n2 = 1:size(Y, 2);
    n3 = 1:(size(X, 2) * 2) - 1;
    n4 = 1:(size(Y, 2) * 2) - 1;
    n5 = 1:(max(size(X, 2), size(Y, 2)) * 2) - 1;
    
    rxx = xcorr(X);
    ryy = xcorr(Y);
    X0 = max(rxx);
    Y0 = max(ryy);
    R0xx = rxx / X0;
    R0yy = ryy / Y0;
    R0xy = xcorr(X, Y, size(n5, 1)) / sqrt(X0 * Y0);
    
    n5 = size(R0xy, 1);
    
    if (GRAF == 1)
        figure
        subplot(5,1,1),plot(n1, X)
        xlabel('n'),title('$x[n]$', 'Interpreter', 'latex', 'FontSize', 14) 
        subplot(5,1,2),plot(n2, Y)
        xlabel('n'),title('$y[n]$', 'Interpreter', 'latex', 'FontSize', 14) 
        subplot(5,1,3),plot(n3, R0xx)
        xlabel('n'),title('$\rho_{xx} [n]$', 'Interpreter', 'latex', 'FontSize', 14) 
        subplot(5,1,4),plot(n4, R0yy)
        xlabel('n'),title('$\rho_{yy} [n]$', 'Interpreter', 'latex', 'FontSize', 14) 
        subplot(5,1,5),plot(n5, R0xy)
        xlabel('n'),title('$\rho_{xy} [n]$', 'Interpreter', 'latex', 'FontSize', 14) 
    end
    
    Vmax = max(R0xy);
    Lmax = find(R0xy == Vmax, 1);
end