Discussion :

[mathieu3392]

Bonjour à tous,
je vous écris car je dois dans mon application développer un Vumetre qui permet de visualiser une fréquence spécifique extraite d'un son en cours de lecture.

j'utilise DirectSound pour lire le son et je possède un buffer circulaire qui simule le comportement du buffer circulaire DirectSound me permettant d'obtenir chaque octet en cours de lecture.

j'ai donc à l'heure actuelle un vumetre qui permet de visualiser le signal global de mon son en cours de lecture. J'ai cherché du côté de Fourier avec la FFT rapide de Fourier mais je ne comprend pas comment l'appliquer à mon cas!

Pourriez-vous m'aider ? J'ai trouvé l'algorithme FFT mais je ne vois pas comment on l'utilise. Par où fournis-t-on les octets et par où fournis-t-on la fréquence à isoler désirée ? Puis-je directement utiliser les octets en entrée?

Merci de votre aide, je ne trouve pas grand chose en documentation. Merci

[pseudocode]

si ca peut t'aider:

http://www.captain.at/howto-fftw-spectrograph.php

[mathieu3392]

Merci c'est très gentil, mais ca ne m'aide pas trop trop;
J'ai trouvé cela : http://www.unilim.fr/pages_perso/jea...t.htm#section2

Mais c'est pareil, c'est l'utilisation que je ne sais pas faire.

[droggo]

Qia,

Pourriez-vous m'aider ? J'ai trouvé l'algorithme FFT mais je ne vois pas comment on l'utilise. Par où fournis-t-on les octets et par où fournis-t-on la fréquence à isoler désirée ? Puis-je directement utiliser les octets en entrée?

On ne "fournit" pas des octets à un algorithme, mais à son implémentation.

Pas besoin de chercher bien longtemps pour trouver une fft implémentée dans n'importe quel langage que tu désires, en général avec un petit mode d'emploi, même si parfois c'est très bref.

[mathieu3392]

oui bien sûr j'en trouve, mais en C avec plein d'autres fonctions comme l'équalisation. Il m'est donc difficile d'isoler cki correspond à l'équaliser et au vumetre. C'est là que je demande votre aide : c'est très théorique et je ne sais pas comment m'y prendre:

Fourier sert à décomposer une signal complexe en sinusoides d'après ce que j'ai compris : 1 sinusoide = 1 fréquence !(?) Ca voudrais donc dire que mon vumetre pour une certaine fréquence aura des évolutions de hauteur linéaires ! ou faudra-t-il faire la moyenne des quelques fréquences avant et aprés la fréquence voulue? Vous voyez ce que je veux dire?

[pseudocode]

Fourier sert à décomposer une signal complexe en sinusoides d'après ce que j'ai compris

En fait, ca fait un changement de repère: Volume(temps) -> Puissance(frequence).

Après la FFT, tu as un tableau S[f] qui te donne l' "importance" de chaque frequence "f".

1 sinusoide = 1 fréquence !(?)

Oui.

Ca voudrais donc dire que mon vumetre pour une certaine fréquence aura des évolutions de hauteur linéaires ! ou faudra-t-il faire la moyenne des quelques fréquences avant et aprés la fréquence voulue? Vous voyez ce que je veux dire?

oui, il faudra comptabiliser aussi les frequences voisines dans ton vuemetre.

[mathieu3392]

merci c'est déjà plus clair.
En gros après la "magique" fftw je récupère un tableau de fréquences à partir de mon octet fourni !

C'est assez simple comme cela, il faudra juste ensuite que je récupère les valeurs proches de ma valeur de fréquence désirée; que j'en fasse la moyenne, j'obtiendrai ainsi une puissance; laquelle varie pour un 8 bits de -127 à +128 ou de 0 à 255 en non signé.. et que je fasse ca à chaque octet de mon son en temps réel (ou groupe d'octets).

Ma dernière question (après je pense pouvoir me débrouiller) :quelle unité FFTW dois-je prendre pour Delphi 2006 ? J'en ai trouvé plusieurs versions sur le web.

Et enfin, pourquoi une notion de nombre complexes pour du son? Comment je fais pr prendre mon octet et aboutir à mon tableau de puissances en fonction de la fréquence ?

(je ne vous demande pas la solution précodée, juste des pistes; j'ai passé cette semaine sur cela, j'avoue avoir tout vu sans avoir trop avancé.. )

Merci à ceux qui m'ont déjà aidé et aux prochains !

[pseudocode]

merci c'est déjà plus clair.
En gros après la "magique" fftw je récupère un tableau de fréquences à partir de mon octet fourni !

Non, pas a partir d'un seul octet. Tu dois lui passer ton buffer de valeurs. Plus il y a d'octet, meilleur sera ton estimation des frequences.

Ma dernière question (après je pense pouvoir me débrouiller) :quelle unité FFTW dois-je prendre pour Delphi 2006 ? J'en ai trouvé plusieurs versions sur le web.

aucune idée.

Et enfin, pourquoi une notion de nombre complexes pour du son? Comment je fais pr prendre mon octet et aboutir à mon tableau de puissances en fonction de la fréquence ?

Ah... toujours la bonne vieille meme question. En fait, l'algo de la FFT est un calcul de DFT complexe.

Et c'est plus simple/rapide de faire une FFT (donc réel+complexe) puis virer la partie complexe que de calculer seulement une DFT réelle.