Оркестр, где каждый инструмент вступает по своему расписанию ✎

Оригинал: "Adaptive Bayesian Power Spectrum Analysis of Multivariate Nonstationary Time Series"
· Zeda Li, Robert T. Krafty
arXiv:1706.05661 · 2017-06-18 · arXiv non-exclusive · собственный разбор · 2 мин
Новый метод делит сложный сигнал на куски, позволяя разным частям меняться независимо, и находит моменты перемен с высокой точностью.

Оркестр играет сложную пьесу, и каждая группа инструментов вступает по своему расписанию. Скрипки могут замедлиться, пока духовые держат прежний темп, а ударные вдруг взорвутся дробью. Дирижёр слышит всё сразу, но определить, в какой именно момент изменился ритм у флейт, а не у виолончелей, почти невозможно. Именно так выглядят сигналы, которые учёные снимают с мозга во время сна или с поверхности океана в годы Эль-Ниньо: много голосов, и каждый живёт своей жизнью.

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

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

На записи мозговых волн во время сна метод нашёл переход от медленной фазы к быстрой с вероятностью 94 процента, и граница легла в четверти секунды от той, что отметил врач вручную.

Для данных Эль-Ниньо метод предположил, что характер колебаний изменился в 1980-х годах, и дал этому событию вероятность 65 процентов. Это не окончательный приговор, а скорее подсказка климатологам, где стоит искать причины перемен. В тестах на искусственных сигналах с известными точками разрыва метод ошибался реже, чем его конкуренты, и работал быстрее одного из самых точных прежних алгоритмов.

Такой инструмент пригодится везде, где важна динамика связей между несколькими сигналами: в медицине сна для диагностики нарушений, в нейронауке для анализа активности разных зон мозга, в климатологии для мониторинга индексов океана и атмосферы. Оркестр сложен, но теперь у нас есть способ услышать, когда именно каждая группа инструментов меняет свой темп.

В анализе ЭЭГ метод автоматически определил переход от медленного сна к быстрому с точностью до четверти секунды, что сравнимо с работой опытного врача.

Идея адаптивного разбиения сигнала на сегменты с разной динамикой напоминает анализ мозговых волн в фильме «Начало», где герои отслеживают изменения состояний сознания по паттернам активности.

f^{-1}(u,\omega) = \Theta(u,\omega) \Psi(u,\omega)^{-1} \Theta(u,\omega)^*
Обратная спектральная матрица выражается через нижнюю треугольную матрицу Θ и положительную диагональную матрицу Ψ, что гарантирует положительную определённость.
Понятия
нелинейный спектральный анализ Лапласиана вейвлет-анализ вариационный метод байесовская иерархическая модель вероятностная модель вейвлет-преобразование неинвазивный мониторинг мозга гауссовский процесс алгоритм Метрополиса-Хастингса параметр апостериорное распределение спектр мощности обращение времени перестановочный тест вложенная выборка уравнение Больцмана байесовский вывод стохастические процессы спектральный анализ с помощью ИИ эмпирическая ортогональная функция