MORE PYTHON

Changes

2 changed files (+62/-11)

• example.wav
• music.py
  • @@ -21,13 +21,57 @@ (1, 3, 6),
    (2, 3, 10), (1 / 2, 3, 9, 1 / 2), (1 / 2, 3, 8, 1 / 2), - (6, 3, 7, 1/2), + (6, 3, 7, 1 / 2), + ] + + melody = [ + (3, 3, 7), + (3, 3, 3), + (3, 3, 8), + (2, 3, 7), + (2, 3, 5), + (3, 3, 7), + (2, 3, 3), + (2, 2, 10), + (2, 3, 0), + (6, 3, 3), + ] + + bass = [ + (1, 0, 7), + (1, 2, 7), + (1, 2, 7), + (1, 0, 7), + (1, 2, 3), + (1, 2, 3), + (1, 0, 7), + (1, 2, 8), + (1, 2, 8), + (1, 0, 7), + (1, 2, 7), + (1, 0, 7), + (1, 2, 5), + (1, 0, 7), + (1, 2, 7), + (1, 2, 7), + (1, 0, 7), + (1, 2, 3), + (1, 0, 7), + (1, 1, 10), + (1, 0, 7), + (1, 2, 0), + (1, 0, 7), + (1, 2, 3), + (1, 2, 3), + (1, 0, 7), + (1, 2, 3), + (1, 2, 3), ] music = [] - def process(notes, start, mult): + def process(notes, start, speed=1, gain=1): """ Adds a list of notes to the music list """
  • @@ -36,12 +80,14 @@ for note in notes:
    vol = 1 if len(note) == 4: vol = note[3] - music.append((t, t + note[0] * mult, note[1], note[2], vol)) - t += note[0] * mult + music.append((t, t + note[0] / speed, note[1], note[2], vol * gain)) + t += note[0] / speed # Process all lists of notes - process(intro, 0, 1 / 3) + process(intro, 0, 3) + process(melody, 8, 3) + process(bass, 8, 3, 1 / 2) # And sort them music.sort()
  • @@ -57,6 +103,8 @@ def tone(f, t):
    """ Returns the intensity of a tone of frequency f sampled at time t """ + # https://dsp.stackexchange.com/questions/46598/mathematical-equation-for-the-sound-wave-that-a-piano-makes + # https://youtu.be/ogFAHvYatWs?t=254 # return int(2**13*(1+square(t, 440*2**(math.floor(5*t)/12)))) # Y = sum([math.sin(2 * i * math.pi * t * f) * math.exp(-0.0004 * 2 * math.pi * t * f) / 2**i for i in range(1, 4)]) # Y += Y * Y * Y
  • @@ -73,14 +121,17 @@ def at(t):
    """ Returns the total intensity of music sampled at time t """ - i = bisect.bisect(music, (t, 10**9)) - 1 + i = bisect.bisect(music, (t, 10**9)) + # This is actually pretty efficient ngl + # Because people usually don't have that many overlapping notes ret = 0 - while i >= 0 and music[i][0] + music[i][1] >= t: - ret += music[i][4] * tone(freq(music[i][2], music[i][3]), t - music[i][0]) - i -= 1 - return int(2**13 * ret) + for j in range(max(i - 20, 0), i): + m = music[j] + if m[0] + m[1] >= t: + ret += m[4] * tone(freq(m[2], m[3]), t - m[0]) + return int(2**12 * ret) # Print out music encoded in s16 to standard output - for i in range(0, 10 * bitrate): + for i in range(0 * bitrate, 16 * bitrate): sys.stdout.buffer.write(struct.pack("h", at(i / bitrate)))