У меня есть файл ascii, содержащий 2048 пар значений x и y. Я просто хочу знать, как построить fft of y в MATLAB. Я пишу следующий код MATLAB, но не смог найти подходящий результат.

Как я могу это сделать? Вот что я попробовал:

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)|')