import numpy as np import scipy.io.wavfile as wf def cosinewave(f): return lambda t: np.iinfo(np.int16).max*np.sin(2*np.pi*f*t) fs = 44100 tn = np.linspace(0, 0.5, int(fs*0.5)) def getTone(n): f = 440*(2**((n-49)/12)) return cosinewave(f)(tn) notes = [52, 52, 59, 59, 61, 61, 59, 59, 57, 57, 56, 56, 54, 54, 56, 52, 59, 59, 57, 57, 56, 56, 54, 54] output = [] for n in notes: output.append(getTone(n)) output = np.array(output) wf.write("song.wav", fs, output.astype(np.int16))