LINUX.ORG.RU
ФорумTalks

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

 


0

2

Доброй ночи. Кому не сложно, сделайте преобразование Фурье вот этому файлу: http://ifolder.ru/31375930

Буду безгранично благодарен, и напишу его ник на самом высоком дереве на Ямале, фотку выложу сюда.

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

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

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

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

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

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

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

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

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

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

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

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

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 туго, подсказать мало чем могу

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

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

vk.com

вступите

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

dikiy ★★☆☆☆ ()

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

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

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

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

Eddy_Em ☆☆☆☆☆ ()
Ответ на: комментарий от dikiy

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

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

Eddy_Em ☆☆☆☆☆ ()
Ответ на: комментарий от AITap

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

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