def test_square(): signal = SquareSignal(200) duration = signal.period * 3 segment = signal.make_wave(duration, framerate=10000) segment.plot() decorate(xlabel='Time (s)') wave = signal.make_wave(duration=0.5, framerate=10000) wave.apodize() spectrum = wave.make_spectrum() spectrum.plot() decorate(xlabel='Frequency (Hz)')
import sys sys.path.insert(1, 'dsp-modulo') from thinkdsp import SinSignal from thinkdsp import SquareSignal from thinkdsp import decorate import thinkplot import numpy senalRectangular = SquareSignal(freq=2, amp=1.0, offset=numpy.pi) waveRectangular = senalRectangular.make_wave(duration=1, start=0, framerate=44100) senalSenoidal = SinSignal(freq=2, amp=1, offset=0) harmonico1 = SinSignal(freq=senalSenoidal.freq * 3, amp=0.3, offset=0) harmonico2 = SinSignal(freq=senalSenoidal.freq * 5, amp=0.25, offset=0) harmonico3 = SinSignal(freq=senalSenoidal.freq * 7, amp=0.20, offset=0) harmonico4 = SinSignal(freq=senalSenoidal.freq * 9, amp=0.15, offset=0) harmonico5 = SinSignal(freq=senalSenoidal.freq * 11, amp=0.1, offset=0) harmonico6 = SinSignal(freq=senalSenoidal.freq * 13, amp=0.05, offset=0) waveSenoidal = senalSenoidal.make_wave(duration=1, start=0, framerate=44100) mezcla = senalSenoidal + harmonico1 + harmonico2 + harmonico3 + harmonico4 + harmonico5 + harmonico6 waveMezcla = mezcla.make_wave(duration=1, start=0, framerate=44100) waveRectangular.plot() waveMezcla.plot() thinkplot.show() espectroRectangular = waveRectangular.make_spectrum() espectroMezcla = waveMezcla.make_spectrum()
from thinkdsp import SawtoothSignal import matplotlib.pyplot as plt from thinkdsp import decorate from thinkdsp import SquareSignal signal = SquareSignal(200) duration = signal.period * 3 segment = signal.make_wave(duration, framerate=10000) segment.plot() decorate(xlabel='Time (s)') wave = signal.make_wave(duration=0.5, framerate=10000) wave.apodize() wave.make_audio() wave.write(filename='6.9.2.wav') spectrum = wave.make_spectrum() spectrum.plot() decorate(xlabel='Frequency (Hz)') plt.show()
from thinkdsp import TriangleSignal, SquareSignal from thinkdsp import decorate import matplotlib.pyplot as plt import numpy as np square = SquareSignal(1100) cysquare = square.make_wave(duration=2, framerate=10000) cysquare.play('temp11.wav') square = square.make_wave(duration=0.5, framerate=10000) square.make_spectrum().plot() decorate(xlabel='Frequency (Hz)') from thinkdsp import SinSignal sin = SinSignal(500).make_wave(duration=2, framerate=10000) sin.play('temp22.wav') plt.show()