Среднее значение во временной области БПФ и среднее значение для частотного бина

12

У меня есть несколько испытаний физиологических данных. Я делаю частотный анализ для анализа мощности (амплитуды) на определенных частотах, представляющих интерес. Является ли усреднение нескольких испытаний одинаковой длины с последующим взятием одного БПФ усредненного сигнала по сравнению с вычислением БПФ для каждого испытания, а затем усреднение частотных бинов одинаковым? На практике я считаю, что это не так.

В частности, сигнал, естественно, имеет сильный 1 / f компонент, и это подчеркивается, если я вычисляю БПФ каждого отдельного испытания, а затем усредняю ​​амплитуды (действительную часть) каждого частотного бина. Эти два эквивалента? Есть ли правильный способ сделать вещи? или при каких принципиальных условиях следует выбирать между усреднением по временной области и усреднением по частотным бинам?


источник

Ответы:

8

Позвольте мне уточнить.

  • Преобразование Фурье не представляет гистограмму сигнала. Преобразование Фурье - это линейное преобразование, которое переносит сигнал из временной области (комплексная функция) в частотную область (еще одна комплексная функция). Требуется сложная функция для другой сложной функции.
  • Преобразование Фурье является линейным, как указывалось выше.
  • Этап в ваших образцах имеет значение, как указано выше. Если данные по пробам изменяются по фазе, то вы не хотите усреднять перед выполнением преобразования Фурье, но вы также не хотите усреднять после преобразования Фурье. Вы хотите усреднить после преобразования Фурье и норма. Ниже я подробно остановлюсь на том, что именно нужно сделать.

Основная проблема заключается в том, что вопрос поставлен неправильно. Это не «я должен взять преобразование Фурье перед усреднением или после усреднения». Потому что это не имеет значения из-за линейности преобразования Фурье.

Правильный вопрос, который нужно задать: «Должен ли я взять амплитуду преобразования Фурье до усреднения или после усреднения». На этот вопрос ответ раньше.

Вот подробности.

Предположим, что ваши выборочные данные представлены последовательностями:

d1=d1[n1],d1[n2],...d1[nN]

d2=d2[n1],d2[n2],...d2[nN]

d3=d3[n1],d3[n2],...d3[nN]

...

dM=dM[n1],dM[n2],...dM[nN]

где - данные из M испытаний, а - выборочные моменты времени, тогда:d1,...dMn1,...nN

F1=j=1M|F{dj}||F{j=1Mdj}|=F2

Таким образом, хотя преобразование является линейным,не является.F|F|

Кроме того, хотя является реальным для всех , нет, ноявляется.dj[ni]i,jF{dj}|F{dj}|

Что касается того, что вы должны сделать, вы должны взять преобразование Фурье отдельных испытаний (через БПФ), получить амплитуду отдельных испытаний и усреднить их вместе.

Наконец, что такое . - это краткий термин для частотного спектра «естественных» сигналов (обычно люди думают об изображениях).1/f1/f

Когда люди говорят, что есть большой компонент , это означает, что амплитуда как функция частоты выглядит как . Это полностью волнистые ... вероятно, от биолога: p1/f1/f

Обратное преобразование Фурье является некоторой функцией знака, но это бесполезно. Это мнимая функция знака! Вещественные функции генерируют симметричное преобразование Фурье.1/f

Фактически, говоря, что спектр , говорит вам кое-что о сигнале, но он не позволяет вам восстановить сигнал. Все, что вы знаете, это то, что, Это не позволяет однозначно определить потому что вся информация о фазе исчезла, и мы знаем, что структура сигнала сильно зависит от его фазы .| F { x ( t ) } | = | 1 / ф | х ( т )1/f|F{x(t)}|=|1/f|x(t)

Что говорит ? Просто он содержит много низких частот и немного высоких частот.1/f

Не менее важный вопрос: что дает вам усреднение? и что более важно, как интерпретировать результат? Настройтесь на завтра для более углубленного обсуждения: p

Тханг
источник
1
+1 Спасибо за разъяснения. Я чувствую, что я неправильно понял основную проблему, стоящую за этим вопросом, и я думаю, что это становится ближе к сути.
Jstarr
2

Во-первых, БПФ - это алгоритм. Преобразование называется преобразованием Фурье! Он представляет гистограмму сигналов. В дискретном случае высокое показание в частотных областях означает много энергии на этой частоте.

Не следует усреднять данные до БПФ, поскольку информация о фазе приведет к значительным изменениям в данных.

Представьте себе 2 образца, каждый из которых состоит из чистого косинуса. В реальном мире вы никогда не захватите этот косинус в одной и той же начальной точке. Один косинус будет смещен относительно другого (или оба имеют разные сдвиги относительно начала. Математически это означает, что y1 = cos (wt-A) y2 = cos (wt-B), где A & B - сдвиги. В вашей модели эти два лучше показать как одно и то же. С небольшой математикой я могу выбрать эти значения так, чтобы y2-y1 = 0. Среднее значение нуля равно нулю и совсем не то, что вы хотите. Это фазовая проблема.

Если ваша цель - найти средний спектр, который вы должны усреднить по спектрам, не усредняйте сигналы!

Михаил
источник
Спасибо за это. Конечно, было беспокойство, что, если есть шум с правильной разностью фаз от испытания к испытанию, это отменит все или сигналы на интересующих частотах. Я до сих пор не понимаю, почему 1 / f усиливается усреднением по частотным бинам, а не усреднением по временной области, за которым следует БПФ.
@ user1487551 Что означает сильный 1 / f компонент? Обратное преобразование Фурье, равное 1 / f, является функцией знака и может указывать на то, что ваши данные содержат большой кусок, когда система стабилизируется. Вы, вероятно, должны показать сюжет или некоторые данные.
Михаил
1
Обратите внимание, что вы можете выполнять усреднение по времени, если у вас есть фазовое задание, с которым вы можете синхронизировать захват данных, и это может быть очень эффективным для снижения минимального уровня шума, но в противном случае вам действительно необходимо выполнить усреднение по ансамблю в частотной области.
Пол Р
2

Если я не полностью от основания или неправильно понять ваш вопрос, ответ да : По линейности ДПФА, усреднение сигналов во время , а затем принимать ДПФ среднего эквивалентен усреднением ДПФА сигналов.

Чтобы показать это, давайте определим некоторые переменные:

  • xn[] : образец пробной временной области в момент времениthn
  • Xk[] : пробная выборка в частотной области на частотеthk

«Средний» сигнал во временной области задается как . Взяв его DFT, мы имеем1L=0Lxn[]

n=0N11LLxn[]ei2πkn/N .

Переключая порядок суммирования, мы можем написать

1L=0Ln=0N1xn[]ei2πkn/N,

но это так же, как

1L=0LXk[l]

что аналогично усреднению ДПФ каждого тривала. Это то, что мы хотели показать.

jstarr
источник
Для усреднения во временной области вам нужна какая-то фазовая ссылка, с которой вы можете синхронизироваться. Когда это возможно, хотя иногда вы можете получить лучшие результаты, чем при усреднении по ансамблю в частотной области, главным образом из-за уменьшения эффектов квантования и т. Д.
Пол Р.