You could try to plot a spectrogram of both signals like so:
spectrogram(x, hann(nfft), nfft/2, nfft, fsample)
When you calculate a PSD with pwelch
, you only get the average spectrum and you lose all temporal information: it will not show you if a signal is louder at the beginning or at the end. Calculating a spectrogram is also pretty similar to what your ears perceive, they measure a sort of spectrum as a function of time. So if you can hear the difference, it should show up as a difference in the spectrogram with enough SNR.
When you make a spectrogram, you should try to play with the frequency and temporal resolution, since they follow some sort of Heisenberg law. Making the FFT window short (by chosing a small nfft
) will give you good temporal resolution, but worse frequency resolution, and the other way around.