Traitement du signal (FFT Transformée de Fourier)

[invite1894b156]

bonjour à tous,

je suis en train d'étudier un signal très bruité provenant d'un accéléromètre (fréquence 100Hz) et j'effectue un lissage à l'aide d'un passage dans le domaine fréquentiel grâce à une transformée de Fourier.
J'effectue une coupure des données (mise à 0) pour toutes les fréquences supérieures à 1Hz afin de pouvoir observer le comportement des variations d'accélération se produisant sur une longue période.

j'obtiens des résultats concluants et un lissage plutôt propre mais j'observe aussi l'apparition d'un offset sur toutes les valeurs du résultat de la transformée inverse.



légende :
- en rouge : le résultat de ma FFT inverse avec un filtrage à 1 Hz
- en noir : le résultat attendu

Je voulais donc savoir si quelqu'un connaissait l'origine de ce décalage ?

ps: je travaille sous javascript et utilise une bibliothèque pour effectuer mes transformées et transformées inverses.

merci à tous.
-----

[invite15928b85]

Bonsoir.

Intéressez-vous à la composante continue ( raie du spectre à fréquence nulle ).
N'y aurait-il pas un "démoyennage" caché dans votre calcul de FFT, ou quelque chose comme cela ?
Faites une simulation spécifique avec un signal constant par exemple.

[invite1894b156]

merci pour votre réponse.
j'avais déjà regardé la composante continue (f=0Hz) sans que cela ne change grand chose au résultat.

je vais étudier plus en profondeur la bibliothèque que j'utilise pour voir si le problème ne vient effectivement pas de la.

Sinon je vais tester, comme vous me le conseillez, un signal continu en entrée de ma FFT et observer le résultat.

je me lance de ce pas

merci encore

[invite1894b156]

re-bonjour,

j'ai testé mon programme avec un signal constant et je ne retrouve pas d'offset. La valeur en sortie est bien celle mise en entrée.

En supprimant la composante continue du signal (0Hz) On obtient un résultat intéressant mais qui ne convient toujours pas (en fait on inverse la valeur de l'offset)


le problème reste donc entier :/ il ne ne me reste plus qu'a me plonger dans les lignes de codes de la librairie de ma FFT ^^'.

merci encore

[A voir en vidéo sur Futura]

[invite9cecbb6f]

Avant de vous pencher dans le code de vôtre librairie vous pouvez essayer avec matlab/octave/numpy... pour voir si le problème vient de votre librairie

[invite1894b156]

alors j'ai effectué des recherches et il s'avère qu'au final il suffisait de diviser par 2 les données dont la fréquence était 0Hz :


au final je retrouve bien mon résultat

merci