t = -80 nH = 30 minf0 = 400 maxf0 = 600 f0et = 5 maxhd = 0.2 stocf = 0.5 y, yh, ys = hpsModel(x, fs, w, N, t, nH, minf0, maxf0, f0et, maxhd, stocf) print "time taken for computation " + str(time.time()-str_time) if __name__ == '__main__': (fs, x) = wp.wavread(os.path.join(os.path.dirname(os.path.realpath(__file__)), '../../sounds/oboe-A4.wav')) w = np.blackman(801) N = 1024 t = -80 nH = 30 minf0 = 400 maxf0 = 600 f0et = 5 maxhd = 0.2 stocf = 0.5 y, yh, ys = hpsModel(x[0:44100], fs, w, N, t, nH, minf0, maxf0, f0et, maxhd, stocf) wp.play(y, fs) wp.play(yh, fs) wp.play(ys, fs)
return y def defaultTest(): str_time = time.time() (fs, x) = wp.wavread(os.path.join(os.path.dirname(os.path.realpath(__file__)), '../../sounds/sax-phrase-short.wav')) w = np.blackman(701) N = 1024 t = -80 nH = 30 minf0 = 400 maxf0 = 700 f0et = 5 maxhd = 0.2 y = harmonicModel(x, fs, w, N, t, nH, minf0, maxf0, f0et, maxhd) print "time taken for computation " + str(time.time()-str_time) if __name__ == '__main__': (fs, x) = wp.wavread('../../sounds/sax-phrase-short.wav') w = np.blackman(901) N = 1024 t = -90 nH = 40 minf0 = 350 maxf0 = 700 f0et = 10 maxhd = 0.2 maxnpeaksTwm = 5 y = harmonicModel(x, fs, w, N, t, nH, minf0, maxf0, f0et, maxhd, maxnpeaksTwm) wp.play(y, fs)
ys[ri:ri+Ns] += sw*ysw # overlap-add for sines yst[ri:ri+Ns] += sw*ystw # overlap-add for residual pin += H # advance sound pointer y = ys+yst # sum of sinusoidal and residual components return y, ys, yst def defaultTest(): str_time = time.time() (fs, x) = wp.wavread(os.path.join(os.path.dirname(os.path.realpath(__file__)), '../../sounds/mridangam.wav')) w = np.blackman(601) N = 2048 t = -70 stocf = 0.2 y, ys, yst = spsModel(x, fs, w, N, t, stocf) print "time taken for computation " + str(time.time()-str_time) if __name__ == '__main__': (fs, x) = wp.wavread(os.path.join(os.path.dirname(os.path.realpath(__file__)), '../../sounds/mridangam.wav')) w = np.blackman(601) N = 2048 t = -70 stocf = 0.2 y, ys, yst = spsModel(x, fs, w, N, t, stocf) wp.play(y, fs) wp.play(ys, fs) wp.play(yst, fs)
def playYsClicked(self): wp.play(self.ys, self.fs)
def defaultTest(): str_time = time.time() (fs, x) = wp.wavread(os.path.join(os.path.dirname(os.path.realpath(__file__)), '../../sounds/sax-phrase-short.wav')) w = np.blackman(801) N = 1024 t = -90 nH = 50 minf0 = 350 maxf0 = 700 f0et = 10 maxhd = 0.2 y, yh, yr = hprModel(x, fs, w, N, t, nH, minf0, maxf0, f0et, maxhd) print "time taken for computation " + str(time.time()-str_time) if __name__ == '__main__': (fs, x) = wp.wavread(os.path.join(os.path.dirname(os.path.realpath(__file__)), '../../sounds/sax-phrase-short.wav')) w = np.blackman(801) N = 1024 t = -90 nH = 50 minf0 = 350 maxf0 = 700 f0et = 10 maxhd = 0.2 maxnpeaksTwm = 5 y, yh, yr = hprModel(x, fs, w, N, t, nH, minf0, maxf0, f0et, maxhd, maxnpeaksTwm) wp.play(y, fs) wp.play(yh, fs) wp.play(yr, fs)
def playYClicked(self): wp.play(self.y, self.fs)
def playYhClicked(self): wp.play(self.yh, self.fs)
w = np.blackman(801) N = 1024 t = -90 nH = 50 minf0 = 350 maxf0 = 700 f0et = 10 maxhd = 0.2 stocf = 0.2 y, yh, yst = hpsModel(x, fs, w, N, t, nH, minf0, maxf0, f0et, maxhd, stocf) print "time taken for computation " + str(time.time()-str_time) if __name__ == '__main__': (fs, x) = wp.wavread(os.path.join(os.path.dirname(os.path.realpath(__file__)), '../../sounds/sax-phrase-short.wav')) w = np.blackman(801) N = 1024 t = -90 nH = 50 minf0 = 350 maxf0 = 700 f0et = 10 maxhd = 0.2 stocf = 0.5 maxnpeaksTwm = 5 y, yh, yst = hpsModel(x, fs, w, N, t, nH, minf0, maxf0, f0et, maxhd, stocf, maxnpeaksTwm) wp.play(y, fs) wp.play(yh, fs) wp.play(yst, fs)
def playXClicked(self): wp.play(self.x, self.fs)
def playYsClicked(self, ys, fs): wp.play(ys, fs)
def playYhClicked(self, yh, fs): wp.play(yh, fs)
def playYClicked(self, y, fs): wp.play(y, fs)
def playYClicked(self, x, fs): wp.play(x, fs)
yrw[:hNs-1] = fftbuffer[hNs+1:] # undo zero-phase window yrw[hNs-1:] = fftbuffer[:hNs+1] ys[ri:ri+Ns] += sw*ysw # overlap-add for sines yr[ri:ri+Ns] += sw*yrw # overlap-add for residual pin += H # advance sound pointer y = ys+yr # sum of sinusoidal and residual components return y, ys, yr def defaultTest(): str_time = time.time() (fs, x) = wp.wavread(os.path.join(os.path.dirname(os.path.realpath(__file__)), '../../sounds/mridangam.wav')) w = np.blackman(901) N = 1024 t = -70 y, ys, yr = sprModel(x, fs, w, N, t) print "time taken for computation " + str(time.time()-str_time) if __name__ == '__main__': (fs, x) = wp.wavread(os.path.join(os.path.dirname(os.path.realpath(__file__)), '../../sounds/mridangam.wav')) w = np.blackman(601) N = 2048 t = -70 y, ys, yr = sprModel(x, fs, w, N, t) wp.play(y, fs) wp.play(ys, fs) wp.play(yr, fs)