clear,clear all; warning off;
Fe=48000;
duration=120; % durée du sweep 
cd('C:\Documents and Settings\JMLC\Mes documents\Audio\cellules');
repertoire=uigetdir;
cd (repertoire);
[nomfichier,nomchemin]=uigetfile('*.wav','choix d''un fichier .wav');
nomchemin;nomfichier;
[signal,Fe,Nbits]=wavread([nomchemin nomfichier]);
tic;
Fe; Nbits;
M=length(signal);  %signal=zeros(1,M);
signal=signal';
maxsignal=max(abs(signal));
Pow2=nextpow2(M);N=2^Pow2;
Fft1=fft(signal,N);
%#################################################
% lecture du chirp initial
[balayage,Fs,Nbits]=wavread(['C:\Documents and Settings\JMLC\Bureau\balayage_.wav']);
balayage=balayage';
Fft2=fft(balayage,N);
i3=1+round(20* N/48000);
 for i = 1:i3       %                        boucle 3  d
    Fft1(i)=0.0;
    Fft1(N-i+1)=0.0;
end 
coef=(19900/24000);
imid=1+(N/2.0);
i3=1+round(coef* N/2.0);
i4=(2*imid)-i3;
 for i = i3:i4       %                        boucle 3  d
    Fft1(i)=0.0;
end 
Fft2=Fft1./Fft2;
Fft1=ifft(Fft2);
pulse = real(Fft1);
pulse = fftshift(pulse);
if abs(max(pulse))< abs(min(pulse));pulse=-pulse;end
pmax=max(pulse);
imax=find(pulse==pmax);
Lp=2^15;
pulc = pulse(imax-Lp : imax+Lp-1);
pmax=max(abs(pulc));
pulc=pulc/pmax;
figure(1);plot(pulc);
pmax=max(abs(pulse));
pulse=pulse/pmax;
wavwrite(pulse,48000,16,'C:\Documents and Settings\JMLC\Bureau\pulse.wav')
toc