code: def phormants(x, sampling_rate):
      N = len(x)
      w = np.hamming(N)
  
      # Apply window and high pass filter.
      x1 = x * w
      x1 = lfilter([1], [1.0, 0.63], x1)
  
      # Get LPC.
      ncoeff = 2 + sampling_rate / 1000
      A, e, k = lpc(x1, ncoeff)
      # A, e, k = lpc(x1, 8)
  
      # Get roots.
      rts = np.roots(A)
      rts = [r for r in rts if np.imag(r) >= 0]
  
      # Get angles.
      angz = np.arctan2(np.imag(rts), np.real(rts))
  
      # Get frequencies.
      frqs = sorted(angz * (sampling_rate / (2 * math.pi)))
  
      return frqs
  
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 17 for calling a global function with example usage: phormants(x, sampling_rate) and returns: frqs