LINUX.ORG.RU

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

 ,


0

1

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

Есть ряд значений - 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 значений частот?

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

>получить 125 значений частот?

А откуда появилось число 125? Это количество максимумов спектральных линий? Или количество участков, заданных на частотной оси? Условие задачи надо формулировать подробнее.

>толковый мануал

А. Оппенгейм, Р. Шафер Цифровая обработка сигналов (DJVU, 12,1 Мб).

quickquest ★★★★★
()
Ответ на: комментарий от quickquest

Спасибо за книгу, а то я видимо окончательно запутался. Вопрос снимается, буду курить маны.

Mr_Brightside
() автор топика
Вы не можете добавлять комментарии в эту тему. Тема перемещена в архив.