[source,fs] = wavread(uigetfile("*.wav","","Select Source"));
target = wavread(uigetfile("*.wav","","Select Target"));
nfft = 2^nextpow2(fs);
target_spectrum=pspect(nfft/2,nfft,"kr",target,8);
source_spectrum=pspect(nfft/2,nfft,"kr",source,8);
matching_spectrum = sqrt(target_spectrum./source_spectrum);
cepst=ifft(log(abs(matching_spectrum)));
cepst(2:nfft/2)=2*cepst(2:nfft/2);
cepst(nfft/2+2:nfft)=0;
mpimp=real(ifft(exp(fft(cepst))));
mpimp=mpimp/max(abs(mpimp));
wavwrite(mpimp,fs,32,uiputfile("*.wav","","Select Output"));