comparison Example VamPy plugins/PyMFCC.py @ 68:44d56a3d16b7

Various fixes to the example plugins

 		maxHz: higher bound of warping (default = Nyquist frequency)
 		'''
 		self.sampleRate = sampleRate
 		self.NqHz = sampleRate / 2.0
 		self.minHz = minHz
-		if maxHz is None : maxHz = self.NqHz
+		if maxHz is None : maxHz = 11025
 		self.maxHz = maxHz
 		self.inputSize = inputSize
 		self.numBands = numBands
 		self.valid = False
 		self.updated = False
+		
+	def reset(self):
+		# reset any initial conditions
+                self.updated = False
+		return None
 		
 	def update(self): 
 		# make sure this will run only once 
 		# if called from a vamp process
 		if self.updated: return self.valid
 		self.updated = True
 		self.valid = False
-#		print 'Updating parameters and recalculating filters: '
-#		print 'Nyquist: ',self.NqHz
-		
-		if self.maxHz > self.NqHz : 
-			raise Exception('Maximum frequency must be smaller than the Nyquist frequency')
-		
-		self.maxMel = 1000*log(1+self.maxHz/700.0)/log(1+1000.0/700.0)
-		self.minMel = 1000*log(1+self.minHz/700.0)/log(1+1000.0/700.0)
-#		print 'minHz:%s\nmaxHz:%s\nminMel:%s\nmaxMel:%s\n' \
-#		%(self.minHz,self.maxHz,self.minMel,self.maxMel)
+                # print 'Updating parameters and recalculating filters: '
+                # print 'Nyquist: ',self.NqHz
+                maxHz = self.maxHz
+		if maxHz > self.NqHz : maxHz = self.NqHz
+		minHz = self.minHz
+                if minHz > self.NqHz : minHz = self.NqHz
+		self.maxMel = 1000*log(1+maxHz/700.0)/log(1+1000.0/700.0)
+		self.minMel = 1000*log(1+minHz/700.0)/log(1+1000.0/700.0)
+                # print 'minHz:%s\nmaxHz:%s\nminMel:%s\nmaxMel:%s\n' \
+                # %(self.minHz,self.maxHz,self.minMel,self.maxMel)
 		self.filterMatrix = self.getFilterMatrix(self.inputSize,self.numBands)
 		self.DCTMatrix = self.getDCTMatrix(self.numBands)
 		self.filterIter = self.filterMatrix.__iter__()
 		self.valid = True
 		return self.valid

 		'''Compute DCT of input matrix.'''
 		return numpy.dot(self.DCTMatrix,data_matrix)
 		
 	def getMFCCs(self,warpedSpectrum,cn=True):
 		'''Compute MFCC coefficients from Mel warped magnitude spectrum.'''
-		mfccs=self.dct(numpy.log(warpedSpectrum))
-		if cn is False : mfccs[0] = 0.0
-		return mfccs
+                eps = 1e-8
+                mfccs=self.dct(numpy.log(warpedSpectrum + eps))
+                if cn is False : mfccs[0] = 0.0
+                return mfccs
 	
 
 class PyMFCC(melScaling): 
 	
 	def __init__(self,inputSampleRate):

 		
 
 	def setParameter(self,paramid,newval):
 		self.valid = False
 		if paramid == 'minHz' :
-			if newval < self.maxHz and newval < self.NqHz :
-				self.minHz = float(newval)
+                        self.minHz = float(newval)
 		if paramid == 'maxHz' :
-			if newval < self.NqHz and newval > self.minHz+1000 :
-				self.maxHz = float(newval)
-			else :
-				self.maxHz = self.NqHz
+                        self.maxHz = float(newval)
 		if paramid == 'cnull' :
 			self.cnull = int(not int(newval))
 		if paramid == 'melbands' :
 			self.numBands = int(newval)
 		if paramid == 'two_ch' :