import time
import itertools
import math
from collections import OrderedDict
from synthplayer.sample import Sample
from synthplayer.synth import WaveSynth, key_freq, octave_notes, major_chord_keys
from synthplayer.synth import Sine, Triangle, Pulse, Square, SquareH, Sawtooth, SawtoothH, WhiteNoise, Linear, Harmonics
from synthplayer.synth import FastSine, FastPulse, FastSawtooth, FastSquare, FastTriangle
from synthplayer.synth import EchoFilter, EnvelopeFilter, AbsFilter, ClipFilter, DelayFilter
from synthplayer.playback import Output
from synthplayer import params
# some note frequencies for octaves 1 to 7
notes = [
None,
OrderedDict((note, key_freq(4+i)) for i, note in enumerate(octave_notes)),
OrderedDict((note, key_freq(16+i)) for i, note in enumerate(octave_notes)),
OrderedDict((note, key_freq(28+i)) for i, note in enumerate(octave_notes)),
OrderedDict((note, key_freq(40+i)) for i, note in enumerate(octave_notes)),
OrderedDict((note, key_freq(52+i)) for i, note in enumerate(octave_notes)),
OrderedDict((note, key_freq(64+i)) for i, note in enumerate(octave_notes)),
OrderedDict((note, key_freq(76+i)) for i, note in enumerate(octave_notes))
]
def demo_tones():
synth = WaveSynth()
with Output(nchannels=1, mixing="sequential", queue_size=2) as out:
for wave in [synth.square_h, synth.square, synth.sine, synth.triangle, synth.sawtooth, synth.sawtooth_h]:
print(wave.__name__)
for note, freq in list(notes[4].items())[6:]:
print(" {:f} hz".format(freq))
sample = wave(freq, duration=0.4).fadein(0.02).fadeout(0.1)
out.play_sample(sample)
print("pulse")
for note, freq in list(notes[4].items())[6:]:
print(" {:f} hz".format(freq))
sample = synth.pulse(freq, duration=0.4, pulsewidth=0.1).fadein(0.02).fadeout(0.1)
out.play_sample(sample)
print("harmonics (only even)")
for note, freq in list(notes[3].items())[6:]:
print(" {:f} hz".format(freq))
harmonics = [(n, 1/n) for n in range(1, 5*2, 2)]
sample = synth.harmonics(freq, 0.4, harmonics).fadein(0.02).fadeout(0.1)
out.play_sample(sample)
print("noise")
sample = synth.white_noise(frequency=440, duration=1.5).fadein(0.1).fadeout(0.1)
out.play_sample(sample)
out.wait_all_played()
def demo_song():
synth = WaveSynth()
notes = {note: key_freq(49+i) for i, note in enumerate(['A', 'A#', 'B', 'C', 'C#', 'D', 'D#', 'E', 'F', 'F#', 'G', 'G#'])}
tempo = 0.3
def synth_sample(freq, duration):
harmonics = [(1, 1), (2, 1/2), (4, 1/4), (6, 1/6)]
a = synth.harmonics(freq, duration, harmonics)
return a.envelope(0.05, 0.2, 0.8, 0.5)
silence = Sample.from_array([0]*int(synth.samplerate*tempo*2), synth.samplerate, numchannels=1)
song = "A A B. A D. C#.. ; A A B. A E. D.. ; A A A. F#.. D C#.. B ; G G F#.. D E D ; ; "\
"A A B. A D C#.. ; A A B. A E D. ; A A A. F#.. D C#.. B ; G G F#.. D E D ; ; "
with Output(synth.samplerate, synth.samplewidth, 1, mixing="sequential", queue_size=50) as out:
for note in song.split():
if note == ";":
print()
out.play_sample(silence)
continue
print(note, end=" ", flush=True)
if note.endswith(".."):
sample = synth_sample(notes[note[:-2]], tempo*4)
elif note.endswith("."):
sample = synth_sample(notes[note[:-1]], tempo*2)
else:
sample = synth_sample(notes[note], tempo)
out.play_sample(sample)
print()
out.wait_all_played()
def demo_plot():
from matplotlib import pyplot as plot
plot.title("Various waveforms")
synth = WaveSynth(samplerate=1000)
freq = 4
s = synth.sawtooth(freq, duration=1)
plot.plot(s.get_frame_array())
s = synth.sine(freq, duration=1)
plot.plot(s.get_frame_array())
s = synth.triangle(freq, duration=1)
plot.plot(s.get_frame_array())
s = synth.square(freq, duration=1)
plot.plot(s.get_frame_array())
s = synth.square_h(freq, duration=1)
plot.plot(s.get_frame_array())
s = synth.pulse(freq, duration=1, pulsewidth=0.2)
plot.plot(s.get_frame_array())
plot.show()
def modulate_amp():
from matplotlib import pyplot as plot
synth = WaveSynth()
freq = 220
s1 = synth.triangle(freq, duration=2)
m = synth.sine(2, duration=2, amplitude=0.4, bias=0.5)
s1.modulate_amp(m)
plot.title("Amplitude modulation by another waveform")
plot.plot(s1.get_frame_array())
plot.show()
with Output(nchannels=1) as out:
out.play_sample(s1)
out.wait_all_played()
s1 = synth.triangle(freq, duration=2)
m = Sine(3, amplitude=0.4, bias=0.5)
s1.modulate_amp(m)
plot.title("Amplitude modulation by an oscillator")
plot.plot(s1.get_frame_array())
plot.show()
with Output(nchannels=1) as out:
out.play_sample(s1)
out.wait_all_played()
def envelope():
from matplotlib import pyplot as plot
synth = WaveSynth()
freq = 440
s = synth.triangle(freq, duration=1)
s.envelope(0.05, 0.1, 0.6, 0.4)
plot.title("ADSR envelope")
plot.plot(s.get_frame_array())
plot.show()
with Output(nchannels=1) as out:
out.play_sample(s)
out.wait_all_played()
def fm():
synth = WaveSynth(samplerate=8000)
from matplotlib import pyplot as plot
freq = 2000
lfo1 = Sine(1, amplitude=0.4, samplerate=synth.samplerate)
s1 = synth.sine(freq, duration=3, fm_lfo=lfo1)
plot.title("Spectrogram")
plot.ylabel("Freq")
plot.xlabel("Time")
plot.specgram(s1.get_frame_array(), Fs=synth.samplerate, noverlap=90, cmap=plot.cm.gist_heat)
plot.show()
with Output(nchannels=1, mixing="sequential") as out:
synth = WaveSynth()
freq = 440
lfo1 = Linear(5, samplerate=synth.samplerate)
lfo1 = EnvelopeFilter(lfo1, 1, 0.5, 0.5, 0.5, 1)
s1 = synth.sine(freq, duration=3, fm_lfo=lfo1)
s_all = s1.copy()
out.play_sample(s1)
lfo1 = Sine(1, amplitude=0.2, samplerate=synth.samplerate)
s1 = synth.sine(freq, duration=2, fm_lfo=lfo1)
s_all.join(s1)
out.play_sample(s1)
lfo1 = Sine(freq/17, amplitude=0.5, samplerate=synth.samplerate)
s1 = synth.sine(freq, duration=2, fm_lfo=lfo1)
s_all.join(s1)
out.play_sample(s1)
lfo1 = Sine(freq/6, amplitude=0.5, samplerate=synth.samplerate)
s1 = synth.sine(freq, duration=2, fm_lfo=lfo1)
s_all.join(s1)
out.play_sample(s1)
lfo1 = Sine(1, amplitude=0.4, samplerate=synth.samplerate)
s1 = synth.triangle(freq, duration=2, fm_lfo=lfo1)
s_all.join(s1)
out.play_sample(s1)
freq = 440*2
lfo1 = Sine(freq/80, amplitude=0.4, samplerate=synth.samplerate)
s1 = synth.triangle(freq, duration=2, fm_lfo=lfo1)
s_all.join(s1)
out.play_sample(s1)
# s_all.write_wav("fmtestall.wav")
out.wait_all_played()
def pwm():
from matplotlib import pyplot as plot
synth = WaveSynth(samplerate=1000)
pwm_lfo = Sine(0.05, amplitude=0.49, bias=0.5, samplerate=synth.samplerate)
s1 = synth.pulse(4, amplitude=0.6, duration=20, pwm_lfo=pwm_lfo)
plot.figure(figsize=(16, 4))
plot.title("Pulse width modulation")
plot.ylim([-35000, 35000])
plot.plot(s1.get_frame_array())
plot.show()
with Output(nchannels=1) as out:
synth = WaveSynth()
lfo2 = Sine(0.2, amplitude=0.48, bias=0.5)
s1 = synth.pulse(440/6, amplitude=0.5, duration=6, fm_lfo=None, pwm_lfo=lfo2)
out.play_sample(s1)
# s1.write_wav("pwmtest.wav")
out.wait_all_played()
def oscillator():
from matplotlib import pyplot as plot
l2 = SquareH(4, samplerate=1000).blocks()
plot.subplot(2, 1, 1)
plot.title("Square from harmonics")
blocks = list(itertools.islice(l2, 2))
values = sum(blocks, [])
plot.plot(values)
harmonics = [(1, 1)]
harmonics.extend([(n, 1/n) for n in range(2, 8*2, 2)])
l3 = Harmonics(4, harmonics, samplerate=1000).blocks()
plot.subplot(2, 1, 2)
plot.title("Even harmonics")
blocks = list(itertools.islice(l3, 2))
values = sum(blocks, [])
plot.plot(values)
plot.show()
def bias():
from matplotlib import pyplot as plot
synth = WaveSynth(samplerate=1000)
waves = [synth.sine(2, 4, 0.02, bias=0.1),
synth.triangle(2, 4, 0.02, bias=0.2),
synth.pulse(2, 4, 0.02, bias=0.3, pulsewidth=0.45),
synth.harmonics(2, 4, [(n, 1 / n) for n in range(1, 8)], 0.02, bias=0.4),
synth.sawtooth(2, 4, 0.02, bias=0.5),
synth.sawtooth_h(2, 4, 7, 0.02, bias=0.6),
synth.square(2, 4, 0.02, bias=0.7),
synth.square_h(2, 4, 7, 0.02, bias=0.8),
synth.white_noise(frequency=100, duration=4, amplitude=0.02, bias=0.9)]
for wave in waves:
plot.plot(wave.get_frame_array())
plot.title("All waveforms biased to levels above zero")
plot.show()
def lfo_envelope():
synth = WaveSynth(samplerate=100)
lfo = Linear(1000, samplerate=synth.samplerate)
lfo = EnvelopeFilter(lfo, 2, 1, 4, 0.3, 2, stop_at_end=True).blocks()
samples = sum(list(lfo), [])
from matplotlib import pyplot as plot
plot.title("LFO Envelope")
plot.plot(samples)
plot.show()
def a440():
synth = WaveSynth(samplerate=44100, samplewidth=4)
a440 = synth.sine(440, duration=2)
with Output.for_sample(a440) as out:
out.play_sample(a440)
out.wait_all_played()
def echo_sample():
synth = WaveSynth(samplerate=44100)
lfo = Linear(1, -0.0001, min_value=-99999)
s = synth.pulse(220, .5, fm_lfo=lfo).fadeout(.2)
with Output(s.samplerate, s.samplewidth, s.nchannels) as out:
e = s.copy().echo(1, 4, 0.5, 0.4) # echo
out.play_sample(e)
e = s.copy().echo(1, 30, 0.15, 0.5) # simple "reverberation" (simulated using fast echos)
out.play_sample(e)
out.wait_all_played()
def echo_lfo():
synth = WaveSynth()
s = Sine(440, amplitude=25000, samplerate=synth.samplerate)
s = EnvelopeFilter(s, .2, .2, 0, 0, 1.5, stop_at_end=True)
s = EchoFilter(s, .15, 5, 0.3, 0.6)
s = ClipFilter(s, -32000, 32000)
blocks = [[int(v) for v in block] for block in s.blocks()]
frames = sum(blocks, [])
import matplotlib.pyplot as plot
plot.plot(frames)
plot.show()
samp = Sample.from_array(frames, synth.samplerate, 1)
with Output.for_sample(samp) as out:
out.play_sample(samp)
out.wait_all_played()
def lfo_func():
rate = 1000
s = Sine(1, amplitude=100, bias=40, samplerate=rate)
s = AbsFilter(s)
s = ClipFilter(s, minimum=20, maximum=80)
s = DelayFilter(s, 0.5)
s = list(itertools.islice(s.blocks(), math.ceil(rate*2/params.norm_osc_blocksize)))
import matplotlib.pyplot as plot
a = plot.subplot(111)
a.set_ylim([-50, 100])
a.plot(sum(s, []))
plot.show()
def bells():
def makebell(freq):
synth = WaveSynth()
duration = 2
divider = 2.2823535
fm = Triangle(freq/divider, amplitude=0.5)
s = synth.sine(freq, duration, fm_lfo=fm)
# apply ADSR envelope that resembles bell amp curve, see http://www.hibberts.co.uk/make.htm
s.envelope(0, duration*0.25, .5, duration*0.75)
s.echo(2, 5, 0.06, 0.6)
return s.make_32bit(False)
b_l1 = makebell(key_freq(56))
b_l2 = makebell(key_freq(60))
b_h1 = makebell(key_freq(78)).amplify(0.7)
b_h2 = makebell(key_freq(82)).amplify(0.7)
b_h3 = makebell(key_freq(84)).amplify(0.7)
bells = b_l1.mix_at(1.0, b_h1)
bells.mix_at(1.5, b_h2)
bells.mix_at(2, b_h3)
bells.mix_at(3, b_l2)
bells.mix_at(4, b_h2)
bells.mix_at(4.5, b_h3)
bells.mix_at(5, b_h1)
bells.make_16bit()
with Output.for_sample(bells) as out:
out.play_sample(bells)
out.wait_all_played()
def stereo_pan():
synth = WaveSynth()
# panning a stereo source:
wave = Sample("samples/SOS 020.wav").clip(6, 12).normalize().fadein(0.5).fadeout(0.5).lock()
osc = Sine(0.4)
panning = wave.copy().pan(lfo=osc).fadeout(0.2)
with Output.for_sample(panning) as out:
out.play_sample(panning)
out.wait_all_played()
# panning a generated mono source:
fm = Sine(0.5, 0.1999, bias=0.2)
wave = synth.triangle(220, 5, fm_lfo=fm).lock()
osc = Sine(0.4)
panning = wave.copy().pan(lfo=osc).fadeout(0.2)
with Output.for_sample(panning) as out:
out.play_sample(panning)
out.wait_all_played()
def osc_bench():
rate = 44100
duration = 2.0
def get_values(osc):
values = list(itertools.islice(osc.blocks(), int(rate*duration/params.norm_osc_blocksize)))
fm = FastSine(220)
print("GENERATING {:g} SECONDS SAMPLE DATA {:d} HZ USING LFO.".format(duration, rate))
print(" WAVEFORM: with-FM / no-FM / optimized")
# sine
print(" Sine: ", end="")
start = time.time()
get_values(Sine(440, fm_lfo=fm))
duration1 = time.time()-start
start = time.time()
get_values(Sine(440))
duration2 = time.time()-start
start = time.time()
get_values(FastSine(440))
duration3 = time.time()-start
print("{:.3f} / {:.3f} / {:.3f}".format(duration1, duration2, duration3))
# triangle
print(" Triangle: ", end="")
start = time.time()
get_values(Triangle(440, fm_lfo=fm))
duration1 = time.time()-start
start = time.time()
get_values(Triangle(440))
duration2 = time.time()-start
start = time.time()
get_values(FastTriangle(440))
duration3 = time.time()-start
print("{:.3f} / {:.3f} / {:.3f}".format(duration1, duration2, duration3))
# square
print(" Square: ", end="")
start = time.time()
get_values(Square(440, fm_lfo=fm))
duration1 = time.time()-start
start = time.time()
get_values(Square(440))
duration2 = time.time()-start
start = time.time()
get_values(FastSquare(440))
duration3 = time.time()-start
print("{:.3f} / {:.3f} / {:.3f}".format(duration1, duration2, duration3))
# sawtooth
print(" Sawtooth: ", end="")
start = time.time()
get_values(Sawtooth(440, fm_lfo=fm))
duration1 = time.time()-start
start = time.time()
get_values(Sawtooth(440))
duration2 = time.time()-start
start = time.time()
get_values(FastSawtooth(440))
duration3 = time.time()-start
print("{:.3f} / {:.3f} / {:.3f}".format(duration1, duration2, duration3))
# pulse
print(" Pulse: ", end="")
start = time.time()
get_values(Pulse(440, fm_lfo=fm))
duration1 = time.time()-start
start = time.time()
get_values(Pulse(440))
duration2 = time.time()-start
start = time.time()
get_values(FastPulse(440))
duration3 = time.time()-start
print("{:.3f} / {:.3f} / {:.3f}".format(duration1, duration2, duration3))
# square_h
print(" Square_H: ", end="")
start = time.time()
get_values(SquareH(440, fm_lfo=fm))
duration1 = time.time()-start
start = time.time()
get_values(SquareH(440))
duration2 = time.time()-start
print("{:.3f} / {:.3f}".format(duration1, duration2))
print("Sawtooth_H: ", end="")
start = time.time()
get_values(SawtoothH(440, fm_lfo=fm))
duration1 = time.time()-start
start = time.time()
get_values(SawtoothH(440))
duration2 = time.time()-start
print("{:.3f} / {:.3f}".format(duration1, duration2))
print(" Noise: ", end="")
start = time.time()
get_values(WhiteNoise(5000))
duration1 = time.time()-start
print(" {:.3f}".format(duration1))
print(" Linear: ", end="")
start = time.time()
get_values(Linear(0, 0.0001))
duration1 = time.time()-start
print(" {:.3f}".format(duration1))
def vibrato():
synth = WaveSynth()
duration = 3
def make_sample(freq):
fmfm = Linear(0, 0.002, max_value=99999)
fm = Sine(0.05, amplitude=0.5, fm_lfo=fmfm)
s1 = synth.sawtooth(freq, duration, amplitude=0.6, fm_lfo=fm)
s1.envelope(0.01, 0.1, 0.6, 2)
return s1
with Output(synth.samplerate, nchannels=1, mixing="sequential") as out:
for f in [220, 330, 440]:
sample = make_sample(f)
out.play_sample(sample)
out.wait_all_played()
def harmonics():
synth = WaveSynth()
freq = 1500
num_harmonics = 6
h_all = synth.harmonics(freq, 1, [(n, 1/n) for n in range(1, num_harmonics+1)])
even_harmonics = [(1, 1)] # always include fundamental tone harmonic
even_harmonics.extend([(n, 1/n) for n in range(2, num_harmonics*2, 2)])
h_even = synth.harmonics(freq, 1, even_harmonics)
h_odd = synth.harmonics(freq, 1, [(n, 1/n) for n in range(1, num_harmonics*2, 2)])
h_all.join(h_even).join(h_odd)
import matplotlib.pyplot as plot
plot.title("Spectrogram")
plot.ylabel("Freq")
plot.xlabel("Time")
plot.specgram(h_all.get_frame_array(), Fs=synth.samplerate, noverlap=90, cmap=plot.cm.gist_heat)
plot.show()
def chords():
synth = WaveSynth()
with Output(nchannels=1, mixing="sequential", queue_size=1) as out:
for rootnote in octave_notes:
chord_keys = major_chord_keys(rootnote, 4)
print("chord", rootnote, ["{0} {1}".format(note, octave) for note, octave in chord_keys])
freqs = [notes[octave][key] for key, octave in chord_keys]
for i in range(1, len(freqs)):
assert freqs[i] > freqs[i-1]
samples = [synth.sine(freq, 1.5, amplitude=0.333) for freq in freqs]
s = samples[0].mix(samples[1]).mix(samples[2]).fadein(0.1).fadeout(0.1)
out.play_sample(s)
out.wait_all_played()
if __name__ == "__main__":
harmonics()
osc_bench()
lfo_func()
bells()
echo_sample()
echo_lfo()
demo_plot()
a440()
demo_tones()
demo_song()
modulate_amp()
envelope()
lfo_envelope()
pwm()
fm()
oscillator()
bias()
stereo_pan()
vibrato()
chords()
