CH02_SEC04_2_SpectrogramBeethoven.ipynb

Kernel: Python 3 (ipykernel)

import scipy
import scipy.io
import IPython.display
import numpy as np
import librosa

import matplotlib.pyplot as plt

#Example to load mat file (first 40 seconds)
#file = scipy.io.loadmat('../DATA/beethoven_40sec.mat')
#sound = file["y"][0] # data
#T = int(file["T"][0][0]) # length in seconds
#FS = int(file["FS"][0][0]) #sample frequency
#samples = T*FS
#t = np.arange(samples)

#original mp3 file
FS = 24000
y, sr = librosa.load('../DATA/beethoven.mp3', sr=FS)

plt.figure(figsize=(20,5))
plt.plot(y)

#for "mat" file use command
#plt.plot(t, sound)

IPython.display.Audio(y, rate=FS)

#For mat file use command 
#IPython.display.Audio(sound, rate=FS)

plt.rcParams["figure.figsize"] = (20,5)
plt.specgram(y, NFFT=8192, Fs=FS, noverlap=2048)
plt.colorbar()