m77so · May 6, 2023 15:51
diff --git a/a.py b/a.py
 import librosa
 import librosa.display
 import numpy as np
 import matplotlib.pyplot as plt
 import scipy

 # def librosa.lpc(y, order): return a

 def lpc_formant(a, fs):
  poles = np.roots(a)
  intns = np.abs(poles)
  bandw = np.log(intns) * fs / (2*np.pi) / (-2)
  ff = np.angle(poles) * fs / 2.0 / np.pi
  return ff[(ff>10) & (ff < fs / 2 - 10) & (bandw < 400)]

 def lpc_spectrogram(a, fs, worN= 513):
  freqs, h  =scipy.signal.freqz(1.0, a, worN=worN)
  return freqs * fs / 2 / np.pi, np.abs(h)
 
diff --git a/memo2.md b/memo2.md
diff --git a/voice-s.md b/voice-s.md
	import librosa
	import librosa.display
	import numpy as np
	import matplotlib.pyplot as plt
	import scipy

	# def librosa.lpc(y, order): return a

	def lpc_formant(a, fs):
	poles = np.roots(a)
	intns = np.abs(poles)
	bandw = np.log(intns) * fs / (2*np.pi) / (-2)
	ff = np.angle(poles) * fs / 2.0 / np.pi
	return ff[(ff>10) & (ff < fs / 2 - 10) & (bandw < 400)]

	def lpc_spectrogram(a, fs, worN= 513):
	freqs, h =scipy.signal.freqz(1.0, a, worN=worN)
	return freqs * fs / 2 / np.pi, np.abs(h)