%load aa_delay;
%load period_delay;
%aa=aa_delay;
%period=period_delay;
y=y1;
conv_fft_abc;
y1=y;
convolution1=wave;
y=y2;
conv_fft_abc;
y2=y;
convolution2=wave;
integration_essai;
%phase=unwrap(phase);
delay=(t.*phase/(2*pi));
for i=1:n1;
   xx(i)=i+delay(i);
%xx(i)=i-delay(i);
end;
%xxx=[1:n1];
%yyy=interp1(xx,y1,xxx);
%yyy=yyy';
%iii=round(delay(1)+2);
%rrr=corrcoef(yyy(iii:n1),y2(iii:n1));
%rrr=rrr(1,2)




  



