FFT scilab

[invite6f166b31]

Bonjour à Tous! alors voilà mon petit problème, j'aimerais trouvé une explication. Je fais un progamme en scilab, dans lequel je dois analyser un signal temporel. Je veux en faire la FFT. je change mon échelle de temps en échelle de fréquence, et je fais fft(signal). Seulemnt je ne comprends pas "l'unité" du résultat. J'ai bien un seul pic qui sort à la fréquence de mon signal, mais faut il que je fasse une normalisation sur mon siganl avant de faire sa FFT? quelle est l'unité du résultat, des mW, des dB????
Si quelqu'un peut m'expliquer un peu comment cela fonctionne je lui en serait reconnaissant.
Merci
-----

[lou_ibmix_xi]

En général c'est plutôt la densité spectrale de puissance que l'on regarde, soit "fft(x) .* conj(fft(x))", auquel cas l'unité du résultat est l'unité de départ au carré. Par contre, suivant l'implémentation de la fft, il faut effectivement normaliser le résultat par le nb de point utilisé par la fft

[invite6f166b31]

Merci beaucoup. Donc si je comprends bien, lorsque j'ai mon signal temporelle, dont je veux faire l'étude en fréquentielle, je lui applique la fft, puis je multiplie par son conjugué, je divise le tout par le nombre d'élément de mon signal de départ, et j'obtiens donc mon signal en fréquentielle. Mon signal temporel est en volt, j'obtiens donc des volts au carré. Je fais 20log10 du tout et mon signal est en dB. Maintenant mon autre problème concerne le passage du temps en fréquence....j'ai défini mon échelle de temps en seconde comme ça : t=TeSec*(0:nbPoint-1);
je ne peux pas faire 1/t pour obtenir mon échelle en fréquence, cela me semble étrange.... non?

[lou_ibmix_xi]

je lui applique la fft, puis je multiplie par son conjugué

Si c'est la densité spectrale de puissance qui t'intéresse: la répartition de la puissance en f° de la fréquence.

je divise le tout par le nombre d'élément

Pas tout à fait, tu divises chaque résultat de FFT par le nb de points utilisés pour calculer la FFT, donc pour la PSD c'est le nb de points au carré. C'est comme ça pour MATLAB, je ne sais pas pour SCILAB.

je ne peux pas faire 1/t pour obtenir mon échelle en fréquence, cela me semble étrange.... non?

Je ne comprends pas trop ce que tu veux dire, mais l'axe des fréquences est données par

F = (0:nbPoint-1) / (TeSec * nbPoint) - 1/(2*TeSec)

Mais là encore attention, suivant l'implémentation de la FFT, et c'est le cas sous MATLAB, les fréquences négatives peuvent se retrouver "après" les fréquences positives.

[A voir en vidéo sur Futura]

[invite6f166b31]

Excusez moi je m'exprime mal, mais vous avez tout à fait compris mon problème pour l'échelle de fréquence. Mais lorsque je trace ma FFT en fonction de l'échelel de fréquence que vous m'avez indiqué j'obtiens un dirac à une fréquence négative, alors que si j'implémente mon échelle de fréquence comme ceci

F = (0:nbPoint-1) / (TeSec * nbPoint)

mon dirac sort à 2MHz qui est la fréquence que j'attends justement. ma méthode est elle quand même correct? J'ai beacoup de mal avec la notion de passage du domaine temporelle à fréquentielle... et avec les fréquence négative