3)  Campioniamo un segnale qualsiasi considerando un numero abbastanza elevato di campioni e calcoliamo la DFT.   Dopo di che cambiamo alternativamente il segno del vettore dei campioni nel dominio del tempo ed effettuiamo la DFT nuovamente : lo spettro periodico ora dovrà essere traslato in frequenza di una quantità pari alla frequenza di Nyquist.

Consideriamo come segnale un rettangolo periodicizzato:

clear all

tempo=42e-3;

T=1/50000;

campioni=round((tempo/T));

k=(1:1:campioni);

xf=(k-1)*T;

yf=[ones(1,campioni/70) zeros(1,(34/35)*campioni) ones(1,campioni/70)];

figure

plot(xf,yf);

xlabel('sec');

title('grafico del segnale nel tempo');

zoom

ff=fft(yf);

figure

freq=[0:1/(T*campioni):(1/T)*(1-1/campioni)];

plot(freq,abs(ff));

xlabel('Hz');

title('spettro del segnale');

zoom

for i=1:length(ff)

   yi(i)=yf(i)*cos(pi*(i-1));

end

ffi=fft(yi);

figure

plot(freq,abs(ffi));

xlabel('Hz');

title('spettro del segnale elaborato');

zoom