script pour une représentation fréquentielle en utilisant le matlab par FFT

[IOMP]

salut tout le monde
j'ai un problème au niveau du spectre du signal .
la méthode que j'ai utiliser c'est celle la, étant x le signal
X=fft(x)
pxx=X.*conj(X)/N ;%pour éliminer la partie imaginaire
f=fe*((0:N)/N)
plot((f,pxx*((0:N)/N)

merci pour votre aide
IOMP
-----

[IOMP]

resalut
est ce que le scrip est juste ou non ?
merci

[invitea29b3af3]

salut

Les parenthèses ne sont pas équilibrées dans ton code. Et N c'est la longueur du signal ?
En fait je sais pas exactement ce que t'essaies de faire, mais si c'est la FFT de x que tu veux, avec des Hertz sur l'axe fréquentiel, tu peux utiliser ça:

% "s" est ton signal, Fe sa fréquence d'échantillonnage
N=length(s);
f=-Fe/2:Fe/N:Fe/2-Fe/N; % tu divises l'espace fréquenciel [-Fe/2, Fe/2] en N parties
S=fft(s)/N;
figure,plot(f,fftshift(abs(S))) % version avec ampli linéaire
figure,plot(f,fftshift(20*log10(abs(S)))) % version avec magnitude