Python: Analyse de la Fréquence des Fichiers Audio

Je suis de générer des fichiers sonores qui jouent des tons de différentes fréquences, avec un certain nombre d'harmoniques.
En fin de compte, ces sons ne pourront être lus sur un appareil avec un petit haut-parleur.

J'ai la courbe de réponse en fréquence de l'enceinte et souhaitez effectuer les actions suivantes en Python:

  1. Tracer le spectre de fréquence d'un fichier audio. J'ai besoin d'un prendre la FFT du fichier et de l'intrigue avec gnuplot
  2. Appliquer une non-linéaire de la fonction de transfert basée sur la courbe de réponse en fréquence dans la feuille de données.
  3. Tracer le résultat une fois que la fonction est appliquée.

Personne ne sait :

  • Ce que le façon la plus simple à faire ce serait?
  • ou d'un Application (GNU/Linux) qui pourrait faire cela pour moi?
SciPy l'introduction comprend une FFT et de tracer le résultat: scipy.org/Getting_Started
Votre logo a l'air délicieux.

OriginalL'auteur waffleman | 2010-11-30

  1. 9

    Je sais que vous n'avez pas mentionné Pylab/Matplotlib, mais il fonctionne. Voici un exemple (en supposant l'unique signal du canal):

    x, fs, nbits = audiolab.wavread('schubert.wav')
audiolab.play(x, fs)
N = 4*fs    # four seconds of audio
X = scipy.fft(x[:N])
Xdb = 20*scipy.log10(scipy.absolute(X))
f = scipy.linspace(0, fs, N, endpoint=False)
pylab.plot(f, Xdb)
pylab.xlim(0, 5000)   # view up to 5 kHz

Y = X*H
y = scipy.real(scipy.ifft(Y))
    Qu'est-ce que H? Vouliez-vous dire N?
    Non, je voulais dire H. Désolé, j'ai juste sauté un tas d'étapes parce que je suis paresseux. Si h est votre réponse d'impulsion, puis H = scipy.fft(h). Dans le point 2 de votre question, vous avez une fonction de transfert en fréquence de la réponse directement, de sorte que serait H.

    OriginalL'auteur Steve Tjoa

  2. 3

    vous pouvez utiliser numpy et matPlotLib. Quelque chose comme le code ci-dessous:

    spectrum = numpy.fft.fft(signal)
frequencies = numpy.fft.fftfreq(len(spectrum))
pylab.plot(frequencies,spectrum)
pylab.show()

    Qui affiche un graphique du spectre fft.

    OriginalL'auteur Nemeth

  3. 0

    scipy a une FFT et les crochets bien dans gnuplot. Vous devriez être en mesure d'utiliser le signal module pour faire le calcul.

    OriginalL'auteur nmichaels