Changeset b507607 for src

Timestamp:

Sep 17, 2009, 4:02:04 AM (15 years ago)

Author:

Paul Brossier <piem@piem.org>

Branches:

feature/autosink, feature/cnn, feature/cnn_org, feature/constantq, feature/crepe, feature/crepe_org, feature/pitchshift, feature/pydocstrings, feature/timestretch, fix/ffmpeg5, master, pitchshift, sampler, timestretch, yinfft+

Children:

4e0fbe6

Parents:

addc9ec

Message:

src/spectral/filterbank_mel.c: correct computation of coefficients, use float samplerate, use POW, make sure we don't try to set unexisting filters

Location:

src/spectral

Files:

• filterbank_mel.c (modified) (6 diffs)
• filterbank_mel.h (modified) (2 diffs)

src/spectral/filterbank_mel.c

@@ -27 +27 @@
 
 void
-aubio_filterbank_set_mel_coeffs (aubio_filterbank_t * fb, uint_t samplerate,
+aubio_filterbank_set_mel_coeffs (aubio_filterbank_t * fb, smpl_t samplerate,
     smpl_t freq_min, smpl_t freq_max)
 {
@@ -42 +42 @@
   uint_t logFilters = 27;
   uint_t allFilters = linearFilters + logFilters;
+
+  if (allFilters > n_filters) {
+    AUBIO_WRN("not enough Mel filters, got %d but %d needed\n",
+        n_filters, allFilters);
+  }
 
   //buffers for computing filter frequencies
@@ -65 +70 @@
   for (filter_cnt = 0; filter_cnt < logFilters + 2; filter_cnt++) {
     freqs->data[0][filter_cnt + linearFilters] =
-        lastlinearCF * (pow (logSpacing, filter_cnt + 1));
+        lastlinearCF * (POW(logSpacing, filter_cnt + 1));
   }
 
@@ -84 +89 @@
   }
 
-  //AUBIO_DBG("filter tables frequencies\n");
-  //for(filter_cnt=0; filter_cnt<allFilters; filter_cnt++)
-  //  AUBIO_DBG("filter n. %d %f %f %f %f\n",
-  //    filter_cnt, lower_freqs->data[0][filter_cnt], 
-  //    center_freqs->data[0][filter_cnt], upper_freqs->data[0][filter_cnt], 
-  //    triangle_heights->data[0][filter_cnt]);
-
   //filling the fft_freqs lookup table, which assigns the frequency in hz to each bin
   for (bin_cnt = 0; bin_cnt < win_s; bin_cnt++) {
@@ -96 +94 @@
   }
 
-  //building each filter table
-  for (filter_cnt = 0; filter_cnt < allFilters; filter_cnt++) {
 
-    //TODO:check special case : lower freq =0
-    //calculating rise increment in mag/Hz
-    smpl_t riseInc =
-        triangle_heights->data[0][filter_cnt] /
-        (center_freqs->data[0][filter_cnt] - lower_freqs->data[0][filter_cnt]);
+  /* zeroing begining of filter */
+  fvec_zeros(filters);
 
-    //zeroing begining of filter
+  /* building each filter table */
+  for (filter_cnt = 0; filter_cnt < n_filters; filter_cnt++) {
+
+    /* skip first elements */
     for (bin_cnt = 0; bin_cnt < win_s - 1; bin_cnt++) {
-      filters->data[filter_cnt][bin_cnt] = 0.0;
       if (fft_freqs->data[0][bin_cnt] <= lower_freqs->data[0][filter_cnt] &&
           fft_freqs->data[0][bin_cnt + 1] > lower_freqs->data[0][filter_cnt]) {
@@ -115 +110 @@
     bin_cnt++;
 
-    //positive slope
+    /* compute positive slope step size */
+    smpl_t riseInc =
+        triangle_heights->data[0][filter_cnt] /
+        (center_freqs->data[0][filter_cnt] - lower_freqs->data[0][filter_cnt]);
+
+    /* compute coefficients in positive slope */
     for (; bin_cnt < win_s - 1; bin_cnt++) {
       filters->data[filter_cnt][bin_cnt] =
           (fft_freqs->data[0][bin_cnt] -
           lower_freqs->data[0][filter_cnt]) * riseInc;
-      //if(fft_freqs->data[0][bin_cnt]<= center_freqs->data[0][filter_cnt] && fft_freqs->data[0][bin_cnt+1]> center_freqs->data[0][filter_cnt])
+      
       if (fft_freqs->data[0][bin_cnt + 1] > center_freqs->data[0][filter_cnt])
         break;
     }
-    //bin_cnt++;
+    bin_cnt++;
 
-    //negative slope
+    /* compute negative slope step size */
+    smpl_t downInc =
+        triangle_heights->data[0][filter_cnt] /
+        (upper_freqs->data[0][filter_cnt] - center_freqs->data[0][filter_cnt]);
+
+    /* compute coefficents in negative slope */
     for (; bin_cnt < win_s - 1; bin_cnt++) {
+      filters->data[filter_cnt][bin_cnt] += 
+          (upper_freqs->data[0][filter_cnt] -
+          fft_freqs->data[0][bin_cnt]) * downInc;
 
-      //checking whether last value is less than 0...
-      smpl_t val =
-          triangle_heights->data[0][filter_cnt] - (fft_freqs->data[0][bin_cnt] -
-          center_freqs->data[0][filter_cnt]) * riseInc;
-      if (val >= 0)
-        filters->data[filter_cnt][bin_cnt] = val;
-      else
-        filters->data[filter_cnt][bin_cnt] = 0.0;
-
-      //if(fft_freqs->data[0][bin_cnt]<= upper_freqs->data[0][bin_cnt] && fft_freqs->data[0][bin_cnt+1]> upper_freqs->data[0][filter_cnt])
-      //TODO: CHECK whether bugfix correct
       if (fft_freqs->data[0][bin_cnt + 1] > upper_freqs->data[0][filter_cnt])
         break;
     }
-    //bin_cnt++;
-
-    //zeroing tail
-    for (; bin_cnt < win_s; bin_cnt++)
-      filters->data[filter_cnt][bin_cnt] = 0.f;
+    /* nothing else to do */
 
   }

src/spectral/filterbank_mel.h

@@ -46 +46 @@
   \param n_filters number of filters
   \param win_s window size
-  \param samplerate
+  \param samplerate audio sampling rate
   \param freq_min lowest filter frequency
   \param freq_max highest filter frequency
@@ -52 +52 @@
 */
 void aubio_filterbank_set_mel_coeffs (aubio_filterbank_t * fb,
-    uint_t samplerate, smpl_t freq_min, smpl_t freq_max);
+    smpl_t samplerate, smpl_t freq_min, smpl_t freq_max);
 
 #ifdef __cplusplus