Кто ни-будь, сделайте быстрое преобразование Фурье массиву в файле

#!/usr/bin/octave -f
x=dlmread(«graf.txt»);
x=x(:,2);
y=fft(x);
dlmwrite(«fft.txt»,y);

http://rghost.ru/38929312

Извиняюсь, если что-то спутал.

http://ompldr.org/vZWpuaA/t.png

оно?

ты вообще сколько частот ожидаешь? На спектре только две отчетливо видно. Но возможно, что некоторые потерялись в начале. Надо для этого езе оконную функцию наложить.

Ой, а разве это нормально что БПФ в минусы уходит? А Вы не могли бы сделать тоже, но без мнимой части?

Скорее всего нет, ожидал что то вроде: http://fotoifolder.ru/view_foto/18tdcyu_o42b/

Ой, а разве это нормально что БПФ в минусы уходит? А Вы не могли бы сделать тоже, но без мнимой части?

возьми просто abs(y).

Скорее всего нет, ожидал что то вроде: http://fotoifolder.ru/view_foto/18tdcyu_o42b/

вторую ж половину надо отрезать всегда :)

ну и рисовать на логарифмической шкале, тогда получишь, что и у меня вышло.

А можно поподробнее, если Вас не затруднит. Как читать логарифмическую шкалу? Я просто впервые сталкиваюсь с этим преобразованием, и честно говоря слабовато представляю себе что мы можем понять из этого графика...

А можно поподробнее, если Вас не затруднит. Как читать логарифмическую шкалу? Я просто впервые сталкиваюсь с этим преобразованием, и честно говоря слабовато представляю себе что мы можем понять из этого графика...

че там считать? На сетке же все пронумеровано.

и честно говоря слабовато представляю себе что мы можем понять из этого графика...

ну мы можем понять, что есть два частотных пика. А на каких частотах конкретно - зависит от длительности единиц отсчета.

данные тут http://ompldr.org/vZWpucA/t.txt

http://ru.wikipedia.org/wiki/Логарифмическая_шкала

#!/usr/bin/octave -f
# Octave 3.6.2, Fri Jun 29 00:53:25 2012 MSK <aitap@Tarkus>
x=dlmread('graf.txt');
x=x(:,2);
y=fft(x);
semilogy(abs(y));
print -dpng fft.png

http://rghost.ru/38929705

Так получается похожая на ту картинку, которую Вы ожидали.

// у меня с матчастью по FFT туго, подсказать мало чем могу

Огромное спасибо за помощь.

В этой группе: http://vk.com/sso.slaviane Через пару недель будут небольшие сюрпризы для Вас. Если интересно вступите :)

Огромное спасибо за помощь.

В этой группе: http://vk.com/sso.slaviane Через пару недель будут небольшие сюрпризы для Вас. Если интересно вступите :)

http://vk.com/sso.slaviane Через пару недель будут небольшие сюрпризы для Вас. Если интересно вступите :)

vk.com

вступите

Да упаси Патрег o_O

Поставь qtiplot или scidavis.

А нафига дублировать? Правая половинка после 600-с-чем-то не нужна.

Ну, а судя по графику, сигнала там практически нет - один шум.

Во-во, у меня так же получилось.

Явный побочный максимум где-то около двадцатки. Но дохловатый какой-то.

ТС, почитай мое пособие.

Хоть оно и для матлаба, в октаве отлично все работает. Там есть и про фурье-обработку данных.

У меня не открывается, просит залогиниться. Что там?

Пока ничего, скоро будет небольшая благодарность, не материальная конечно, так творческая благодарность.

Благодарю, вот это находка! Очень классно!