code: def mfcc_filter_banks(
      sampling_rate,
      num_fft,
      lowfreq=133.33,
      linc=200 / 3,
      logsc=1.0711703,
      num_lin_filt=13,
      num_log_filt=27,
  ):
      """
      Computes the triangular filterbank for MFCC computation
      (used in the stFeatureExtraction function before the stMFCC function call)
      This function is taken from the scikits.talkbox library (MIT Licence):
      https://pypi.python.org/pypi/scikits.talkbox
      """
  
      if sampling_rate < 8000:
          nlogfil = 5
  
      # Total number of filters
      num_filt_total = num_lin_filt + num_log_filt
  
      # Compute frequency points of the triangle:
      frequencies = np.zeros(num_filt_total + 2)
      frequencies[:num_lin_filt] = (
          lowfreq + np.arange(num_lin_filt) * linc
      )
      frequencies[num_lin_filt:] = frequencies[
          num_lin_filt - 1
      ] * logsc ** np.arange(1, num_log_filt + 3)
      heights = 2.0 / (frequencies[2:] - frequencies[0:-2])
  
      # Compute filterbank coeff (in fft domain, in bins)
      fbank = np.zeros((num_filt_total, num_fft))
      nfreqs = np.arange(num_fft) / (1.0 * num_fft) * sampling_rate
  
      for i in range(num_filt_total):
          low_freqs = frequencies[i]
          cent_freqs = frequencies[i + 1]
          high_freqs = frequencies[i + 2]
  
          lid = np.arange(
              np.floor(low_freqs * num_fft / sampling_rate) + 1,
              np.floor(cent_freqs * num_fft / sampling_rate) + 1,
              dtype=int,
          )
          lslope = heights[i] / (cent_freqs - low_freqs)
          rid = np.arange(
              np.floor(cent_freqs * num_fft / sampling_rate) + 1,
              np.floor(high_freqs * num_fft / sampling_rate) + 1,
              dtype=int,
          )
          rslope = heights[i] / (high_freqs - cent_freqs)
          fbank[i][lid] = lslope * (nfreqs[lid] - low_freqs)
          fbank[i][rid] = rslope * (high_freqs - nfreqs[rid])
  
      return fbank, frequencies
  
imports: from __future__ import print_function
  import math
  import numpy
  import sys
  from scipy.fftpack import fft
  import matplotlib.pyplot
  from scipy.signal import lfilter
  from scipy.fftpack.realtransforms import dct
  from tqdm import tqdm
label_desc: Use image node_id 10 for calling a global function with example usage: mfcc_filter_banks(sampling_rate, num_fft, lowfreq, linc, logsc, num_lin_filt, num_log_filt) and returns: fbank,  frequencies