Teoria dei segnali: DFT di un coseno finestrato

[BrianMay85]

Ciao a tutti.
Ho un problema con la trasformata discreta di un coseno moltiplicato per un rect in matlab.
So che dovrebbero venire 2 sinc centrati in -f0 e f0 con f0=1/(periodo cos).
A me vengono 2 metà sinc quasi alla fine dell'asse delle frequenze, è normale che siano a metà?!?!
Ho usato il comando fft e fftshift per "adattare" il segnale.
Grazie!
Ciao

God save the QUEEN

[goblyn]

T=0.7; % periodo coseno
P=1; % potenza sinusoide
Np=10; % numero di periodi del coseno da finestrare
N=100; % numero di campioni del segnale
Ttot=T*Np; % tempo totale di finestratura
Tc=Ttot/N; % tempo di campionamento
t=[0:Tc:(N-1)*Tc]; % asse dei tempi
fc=1/Tc; % frequenza di campionamento
fn=fc/2; % frequenza di Nyquist
f=[-fn:fc/N:fn-fc/N]; % asse delle frequenze

x=sqrt(2*P)*cos(2*pi/T*t);
X=abs(fft(x));

plot(f/fc,fftshift(X/N))
grid