2009-12-17 3 views
1

Я выполнил fft (быстрое преобразование Фурье) в форме временного ряда в Matlab, но у меня, кажется, есть странная волна на самом деле в графике преобразования Фурье, хотя есть шипы, эта волна выглядит как то, что я ожидаю увидеть только во временной области. Есть ли какие-либо соображения программирования, почему это может произойти?Странный артефакт в моем преобразовании Фурье

+1

вы могли бы быть более конкретным? Картинки? – monksy

+0

Как мне добавить картинку? – Jon

+1

Загрузите его на сайт для размещения изображений и ссылку на него. –

ответ

7

Преобразование Фурье очень похоже на инверсное преобразование Фурье. Шип в одном - это волна в другом. Следовательно, если в вашей серии имеется один пакет данных outlier, у вас будет волновой компонент в частотной области.

Возможная проблема, связанная с программированием, может быть неинициализированной точкой данных, например. обеспечивая 1023 точек данных для БПФ с разрешением 1024 точки.

+1

, но Matlab должен заполнить fft нулями, если это неправильная длина? – Jon

+0

Если данные представляют собой, скажем, косинусоидальную волну, которая равна 1 с обоих концов, то дополнительный ноль в конце будет выбросом. –

+0

Как я могу удалить выбросы, тогда они, кажется, разлагают его. – Jon

7

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

2

В настоящее время Стив проводит прекрасную серию преобразований Фурье on his blog. Он специально говорит о двумерных преобразованиях, но вы можете найти его обсуждение окна полезным.

6

Обратите внимание, что преобразование Фурье непрерывного времени сигнала конечной длины может иметь такие вещи, которые выглядят как «всплески» в частотной области. См. Графики в this post непрерывного преобразования Фурье за ​​один период косинусного сигнала и десяти периодов косинусного сигнала.

Например, сигнал косинуса с бесконечной степенью сжатия имеет простое преобразование Фурье, которое представляет собой пару импульсов при +/- частоте косинуса. Но если у вас есть только десять периодов сигнала косинус преобразование Фурье выглядит следующим образом:

alt text

+0

Точно. Звон. –