Eu tenho um arquivo ASCII contendo pares de valores 2048 x e y. Eu só quero saber como traçar o fft de y no MATLAB. Estou escrevendo seguindo o código do MATLAB, mas não consegui encontrar o resultado apropriado.

Como posso fazer isso? Isso é o que eu tentei:

I = load('data1.asc');

for i = 1:2048
    y = I(:,2);
end

plot(x)

Fs = 40000;                    
T = 1/Fs;                   
L = 2000;     
NFFT = 2^nextpow2(L);
Y = abs(fft(y,NFFT))/L;
f = Fs/2*linspace(0,1,NFFT/2+1);

figure, plot(f,2*abs(Y(1:NFFT/2+1))) 
axis([0 40000 0 40])
xlabel('Frequency (Hz)')
ylabel('|Y(f)|')