ameliorer resolution fft

Bonjour à tous
j'essais de faire la transformée de fourier d'un signal sinusoïdal avec la fonction fft de matlab .
voici mon code
n etant la longueur du signal que je veux
fs : la fréquence échantillonnage
nb:le nombre de poit que je veux représenter de ma fft

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13

function [ftp]=tft(signal,n,fs,nb) figure(1); signal=signal(1:n); subplot(2,1,1) plot(signal); %longueur signal fn=hanning(n); s=signal'.*fn;% le fenetrage w = 0:fs/(n):(fs/2); ft=fftshift(abs(fft(s,n))); ftp=ft(fix(length(ft)/2)+1:(length(ft))); subplot(2,1,2); plot(w(1:nb),ftp(1:nb));

---

et voici ce que j'obtient quand j'applique cette fonction sur un signal sinusoïdal de fréquence 1 hz avec une fréquence échantillonnage de 1khz
http://www.heberger-image.fr/data/images/41069_fft.png
mon problème c'est que je veux obtenir un dirac normalement et non pas cette sorte de cone en la frequence fondamental de mon signal.
qu'est ce que je dois modifier dans mon code pour obtenir celà
Merci

Bonjour.

C'est normal. Vous avez 'apodisé' (fenêtre de Hanning), le signal que vous analysez n'est plus une sinusoïde mais une sinusoïde modulée par la fonction de Hanning. Son spectre n'est pas un Dirac. Vous aurez un Dirac si et seulement si 1/ vous ne fenêtrez pas; 2/ la période est un sous-multiple de la fenêtre (ex: 16 points pour une FFT de 256).

Ce que vous donne la FFT est le spectre de la fenêtre répétée à l'infini. Si vous ne remplissez pas la condition 2, la sinusoïde répétée est "cassée" (la phase n'est pas continue) et cette discontinuté a une signature spectrale. On utilise justement les fenêtres pour minimiser cette pollution, mais elles élargissent les pics spectraux.

Par ailleurs, utiliser la fonction stem à la place de plot. Elle ne relie pas les points.