from scipy.io import wavfile import numpy as np import matplotlib.pyplot as plt srate, data = wavfile.read("nasty.wav") data = data.flatten() for i in range(1, len(data)): data[i] = data[i-1] - data[i] plt.hist(data, range=(-2**23, 2**23), bins=128) plt.show()