ESEMPI DI CALCOLO DELLA DFT Uso della DFT per ampliare la finestra di osservazione in frequenza e nel tempo 2)
Consideriamo lo stesso segnale precedente, calcoliamo la DFT,
inseriamo dei campioni nulli nel tempo ricalcoliamo la DFT: quello che
dobbiamo ottenere è un ampliamento della finestra di osservazione nel dominio
delle frequenze.Viceversa se gli zeri sono inseriti nello spettro e calcoliamo
la IDFT. k=[-201:1:201]; f0=1000; fc=5000; x=k/fc; y=sin(pi*f0*x)./(pi*f0*x); y(202)=1; plot(x*1000,y) xlabel('ms'); ff=fft(y); figure freq=[0:5000/403:5000*(1-1/403)]; plot(freq,abs(ff)) xlabel('Hz'); zoom j=1; for
i=1:403
yi(j)=y(i);
yi(j+1)=0;
j=j+2; end ffi=fft(yi); m=length(yi); freq=[0:10000/m:10000*(1-1/m)]; figure plot(freq,abs(ffi)) xlabel('Hz'); j=1; for
i=1:403
ffii(j)=ff(i);
ffii(j+1)=0;
j=j+2; end yii=ifft(ffii); n=length(yii); xi=[2*x(1):(x(403)-x(1))/402:2*x(403)+(x(403)-x(1))/402]; figure plot(xi*1000,real(yii))
|