Discussion :

[Jeanmi7]

Bonjour,

J'ai pu voir sur votre forum des sujets en lien avec ma question mais je préfère ouvrir un nouveau sujet pour plus de clarté. Après pas mal de recherches sur le web je n'ai pas trouvé de solutions à mon problème.

Je voudrais réaliser la FFT d'un fichier wav, de façon à avoir son spectre de fréquence, cependant je ne maîtrise pas assez le concept mathématiques et Python pour le faire moi même. Ainsi pour débuter j'ai trouvé sur internet ce code https://sites.google.com/site/haskel...on-numpy-scipy, cependant quand je le fais tourner sur mes fichiers wav j'ai des erreurs de ce style :

Traceback (most recent call last):
File "<tmp 2>", line 24, in <module>
y=data[:,1]
IndexError: too many indices

Je pense que cela peut être du au fait que mon fichier ne soit pas au bon format. J'ai donc essayé de modifié la fréquence d’échantillonnage sur Audacity mais sans succès. Pour information le fichier a les caractéristiques suivantes: Mono, 44100Hz, 32 bits flottant.

Merci de votre aide.

[tyrtamos]

Bonjour,

Mes souvenir de fft sont trop loin pour que je puisse beaucoup t'aider. Cependant, par curiosité, j'avais fait il y a quelques temps un petit code .
Voilà mon petit code de test. Il charge un fichier wav, et affiche son spectre dans une fenêtre :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
import matplotlib.pyplot as plt
from scipy.io import wavfile
from scipy.fftpack import rfft
 
fs, data = wavfile.read('ringin.wav')  # chargement du fichier wav.
# fs est la fréquence d'échantillonnage / s => 44100 points/seconde
 
print(fs, len(data))  # 44100 49728
 
a = data.T[0]  # 1er canal. type de data et a: class 'numpy.ndarray'
 
print(len(a))
 
b = [(elem / (2 ** 8)) * 2 - 1 for elem in a]  # changement d'échelle des valeurs
 
c = rfft(b)  # liste de complexes. type de c: class 'numpy.complex128'
 
d = len(c) // 2  # on n'a besoin que de la moitié de la liste fft
 
plt.plot(abs(c[0:(d - 1)]), 'r')  # affichage
plt.show()

[Jeanmi7]

Merci pour ton aide mais j'ai encore une erreur :

Traceback (most recent call last):
File "<tmp 3>", line 12, in <module>
print(len(a))
TypeError: object of type 'numpy.int16' has no len()


Je ne vois pas trop comment la résoudre.

EDIT: c'est bon j'ai trouvé un autre fichier son qui fonctionne. Cependant il n'y a pas d'échelle en abscisse, l'unité utilisée est-elle arbitraire ?