Utilisation de la FFT sous R

Fermé
Yann Béli - 20 juin 2012 à 14:25
Bonjour,

Actuellement en stage dans une banque, et débutant avec Fourier, je cherche à vérifier quelques transformées de base avec le logiciel R. À commencer par la transformée de la loi normale, qui donne normalement une loi normale... J'ai donc écrit le code suivant dans l'éditeur :


a <- 0
for (i in 1:2048) # On choisit 2048 = 2^11
{
a[i] <- dnorm((i-1024)/256)
}
b <- fft(a)
windows()
plot(a)
windows()
plot(1/Mod(b))

Ce code me donne donc deux graphes : le premier, une loi normale, celle que l'on veut transformer ; le deuxième, 1/Mod(b) , une fonction ressemblant à une loi normale mais qui n'en est pas une (on s'en convainc en écrivant plot(b), qui retourne une figure complexe étrange).
Encore une fois, ceci a pour but de me familiariser avec l'utilisation de la FFT sous R. Que puis-je faire pour corriger mon code ?


A voir également: