In [1]:
%load_ext autoreload
%autoreload 2
%matplotlib inline
import matplotlib.pyplot as plt
from IPython.display import Audio, display
import numpy as np
# from pysas import World, waveread # https://github.com/shunsukeaihara/pysas
In [2]:
import soundfile
def pad2(i):
return '{0:02d}'.format(i)
# read phoneme WAVs from file
datas = []
path = "/Users/pi/Dev/pythonStuff/VowelGen/4. Analyze my own vowels/vowels/me/"
for i in range(12):
file = path + pad2(i) + '.aiff'
print(file)
data, samplerate = soundfile.read( file )
data = data[:,0]
datas.append( data )
In [3]:
from pysas import World
world = World(samplerate)
signal = datas[0]
f0, spec_mat, aperiod_mat = world.analyze(signal)
plt.plot(f0)
Out[3]:
In [4]:
spec = np.sum( spec_mat[40:80], axis=0 )
plt.plot(np.log(spec))
Out[4]:
