[Scilab][R][БПФ] простой вопрос

Возник один глупый вопрос

Есть ряд значений - 8192 отчсчетов. Частотота дискретизации f=250Hz. Провожу БПФ в Scilab (или R - не важно), код элементарный:

    [Y]=read(some_file,8192,1);
    Fs=250;
    T=1/Fs;
    L=size(Y,1);
    t=(0:L-1)*T;
    NFFT=2^nextpow2(L);
    y=fft(Y(:,i));
    f=Fs/2*linspace(0,1,NFFT/2+1);
    plot(f,(abs(y(1:NFFT/2+1)))/L)
    title('Спектральная мощность сигнала','fontsize',4)
    xlabel('Частота, Гц','fontsize',2)
    ylabel('Мощность спектра, отн. ед.','fontsize',2)
    xgrid

Получаю чудесный график нормированного частотного спектра сигнала.

Вопрос: как проще всего из 4096 значений спектра (которые (abs(y(1:NFFT/2)))/L ) получить 125 значений частот?

Может кто сталкивался и посоветует толковый мануал (рус/англ), буду весьма признателен.