%skripta koja usporeduje originalni zig i umetnuti zig
%zig    = originalni zig
%zigR   = izvuceni zig

%alpha  = 0.1;
%ucitavaje zigova......


maxcorrelation = -10000.0;

%2 razine razlaganja -> 6 detalja
sig_n   = size(zig); %broj elemenata originalnog ziga
in_n    = sig_n;

td = 1; %brojac
for i=1:6
    n=zigR(td,1);
    s1 = 0.0;
    s2 = 0.0;
    s3 = 0.0;
    for j=1:n
        in_x=0; sig_x=0;
        sig_x=zig(mod(j-1, sig_n(1))+1 , 1);
        
        in_x=zigR(j+1,1);
        
        s1 = s1 + sig_x * in_x;
        s2 = s2 + in_x * in_x;
        s3 = s3 + sig_x * sig_x;
        
        td = td + 1;
        
    end
    
    td = td + 1;
    
    correlation = s1 / sqrt (s2*s3);
    
    if (correlation > maxcorrelation)
        maxcorrelation = correlation;
    end
    
    rezultat_korelacija(i) = correlation;
   
   
end
    rezultat_korelacija
    maxcorrelation
