close all;
clear;

%1) phi(t) je ortogonalne na svoje posuny

max=10;
k=-(floor(max));
for j=1:(2*max)+1
    m=round(k*pocet/dlzka);
    for i=1:pocet
        
        if (i-m)>0 & (i-m)<=pocet
            sucin(i)=phi(i)*phi(i-m);
        else
            sucin(i)=0;
        end  
    end
    autokorelacia(j)=integral(sucin,t);
    K(j)=k;
    if(k==0)
%        plot(t-k,phi,'red');
    else
 %       plot(t-k,phi,'blue');
    end
    k=k+1;
    clear sucin;

end

figure('name','Ortogonalita phi na svoje celociselne posuny (autokorelacia)');
stem(K,autokorelacia,'-ob');

