Short Time Fourier Transform как сделать в Mathematica ?

[qweras]

Subj.

[Vince Diesel]

Чего бы не по определению? Из того, что написано на вики, следует, что это преобразование Фурье от некоторой функции.

[qweras]

Vince Diesel

Смысл STFT в анализе изменений спектра со временем. Интересует сама функция которая выделит частоты, которые наиболее часто встречаются в анализируемом потоке данных. Никто не подскажет как сделать STFT в Matematica ?

[arseniiv]

qweras в сообщении #922786 писал(а):

Интересует сама функция которая выделит частоты, которые наиболее часто встречаются в анализируемом потоке данных.

Так это уже не STFT, а функция от него. Можно банально брать взвешенное среднее в каждый момент времени. Можно делать что-то другое, если вам нужно несколько частот.

Вам STFT нужно, обработчик или всё и сразу? Если STFT в списке, то по совету Vince Diesel присмотритесь к формуле преобразования. Функция умножается на сдвинутую оконную, а потом берётся преобразование Фурье. От выбора оконной функции зависит пригодность того, что получится, вашим целям.

Ещё есть такая вещь как мгновенная частота. Возможно, вам пригодится и она.

[qweras]

{x} = {data[[;; , 8]]}
{y} = Fourier[{x}, FourierParameters -> {0, 1}]
Mg := Abs[y]
ListLinePlot[Mg]

Почему-то график получается симметричный, а как сделать что-бы график был убывающией и справа после минимума опять не поднимался вверх ? И почему-то, на симметричном графике в начале и в конце максимумы похоже не помещаются на график - как исправить ?

Как на этом изображении:
http://habrastorage.org/storage3/39d/a6e/71b/39da6e71bcb260da8146236b02d5a384.png

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

[arseniiv]

qweras в сообщении #929637 писал(а):

Почему-то график получается симметричный

Вы про частоту Найквиста знаете?

А вкратце, возьмите первую половину массива.

[qweras]

arseniiv в сообщении #929640 писал(а):

А вкратце, возьмите первую половину массива.

Если "не вкратце" - как правильно сделать, что-бы убывало до нуля или спектр уходил ниже ноля ? (Как на втором изображении из прошлого поста).

[arseniiv]

Ах это… Тогда вопрос будет про децибелы. Логарифмируйте! (Только правильно логарифмируйте.)