%%%%%%%%%%%%%%%%%ALIGN SIGNALS%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%zamik v soundi med ušesoma: 0.00058309037s

[s1,Fs] = audioread('Sinus.wav');

[s2,Fs2] = audioread('Sinus2.wav');

s1(numel(s2)) = 0; %dodamo ničle da je enako vzorcev
len1=size(s1);
len2=size(s2);

F1=fft(s1);  %hitra fourierjeva transformacija
F2=fft(s2);


%P=atan2(imag(F1),real(F1)); %faze za pregled
%plot(P);

%P2=atan2(imag(F2),real(F2));
%plot(P2);


amplituda=abs(F1); %amplituda
amplituda2=abs(F2);

[M,I]=max(amplituda); %index vzorca z najvišjo amplitudo



%temp=(Fs/len1(1))*I;

faza=atan2(imag(F1(I)),real(F1(I))); %faza na vzorcu z najvišjo amplitudo
faza2=atan2(imag(F2(I)),real(F2(I)));%faza na drugem vzorcu na enakem položaju
n=faza-faza2; % razlika med fazama na vzorcu

FsPerioda=Fs/200;

temp=n*FsPerioda; %križni račun : zamik v radianih -> zamik v vzorcih

zamik=temp/(2*pi);