holger@0: function numBinsPerSemitone = getNumBinsPerSemitoneFromFreqVec(freqsInHz)
holger@0: % numBinsPerSemitone = getNumBinsPerSemitoneFromFreqVec(freqsInHz)
holger@0: %
holger@0: % infers the number of bins per semitone from a given constant-Q frequency
holger@0: % vector in Hz
holger@0: 
holger@0: ratios = freqsInHz(2:end) ./ freqsInHz(1:end-1);
holger@0: avgRatio =  mean(ratios);
holger@0: numBinsPerSemitone = round( 1/(12*log2(avgRatio)) );