from audiodisplay import Audio

%pylab inline

max_time = 3
f = 440.0
rate = 44100.0
times = np.arange(max_time*rate)/rate
signal = 2.0**16*np.sin(2*pi*f*times)/2.0
signal = signal.astype(np.int16)

plot(signal[0:500])

Audio(data=signal, rate=rate)

filename = 'data/Bach Cello Suite #3.wav'

Audio(filename=filename)

from scipy.io import wavfile

rate, bach = wavfile.read(filename)

fig, (ax1, ax2) = plt.subplots(1,2, figsize=(12,4))
ax1.plot(bach[:,0]); ax1.set_title('Raw audio signal');
ax1.locator_params(nbins=6)
ax2.specgram(bach[:,0]); ax2.set_title('Spectrogram');
ax2.locator_params(nbins=6)