CC BY
22
МЕТОДЫ СПЕКТРАЛЬНОГО АНАЛИЗА В ОЦЕНКЕ РЕЗУЛЬТАТОВ
МОНИТОРИНГА РИТМА СЕРДЦА ПРИ ЭМОЦИОНАЛЬНОЙ И
УМСТВЕННОЙ НАГРУЗКЕ
Осипенко Александр Викторович
соискатель, НИУ «МЭИ», РФ, г. Москва Е-mail: avo 77@yandex. ru
METHODS OF SPECTRAL ANALYSIS IN RECOGNITION OF THE HRV-MONITORING RESULTS DURING EMOTIONAL AND MENTAL STRESS
Alexander Osipenko
applicant, National Research University "MPEI", Russia, Moscow АННОТАЦИЯ
В статье анализируется применение методов спектрального анализа для обработки данных мониторинга ритма сердца, и предлагается методика предварительной обработки измерений для исключения физиологического шума, трендов и погрешностей сигнала в процессе длительной регистрации ВСР при эмоциональной и умственной нагрузке.
ABSTRACT
The article analyzes the application of methods of spectral analysis for data of HRV-monitoring, and proposes a methodology for pre-processing of the measurements to exclude physiological noise, trends and errors of the signal in the longtime-registration of HRV during emotional and mental stress
Ключевые слова: ВСР; методы и алгоритмы спектрального анализа; преобразования Фурье; корреляционная функция; спектральная плотность.
Keywords: HRV; methods and algorithms of spectral analysis; Fourier transform; correlation function; spectral density.
В последние годы показатели вариабельности сердечного ритма (ВСР) находят все более широкое применение, как в клинической практике, так и в оценке функциональных состояний в физиологии труда, в том числе и для оценки эмоциональной и умственной нагрузки. При этом для получения сопоставимых результатов необходимо обеспечить условие стационарности, что
^ created by free version of
S DociFreezer
практически невозможно при длительной регистрации особенно в условиях трудовой деятельности. Это обуславливает необходимость применения математических методов обработки данных, для исключения физиологического шума и погрешностей сигнала, исследованию этих методов и разработке методики оценки ВСР при длительной регистрации посвящена данная работа.
Как известно ВСР наблюдается даже в состоянии покоя, а непостоянство интервала между кардиоциклами находится в пределах некой средней величины, являющейся оптимальной для рассматриваемого функционального состояния организма. Причем при любом изменении (положения тела, позы и пр.) частота сердечных сокращений начинает подстраиваться под новый функциональный уровень, образуя тренды, сохраняющиеся и некоторое время после окончания двигательной активности, аналогичное влияние оказывает и суточная ритмика, выражающаяся проявлением периодической составляющей физиологических функций.
Вместе с тем, стандартами [2] предусмотрены:
• 5 минутные регистрации RR интервалов, для которых нет сложностей обеспечения покоя и стационарности,
• и 24 часовые, где включаются циркадные колебания, и позные тренды, поскольку обеспечить стационарность при длительной регистрации сложно.
Исключение названных факторов обычно обеспечивается применением спектрального анализа RR-интервалов и прежде всего преобразованиями Фурье, которые программно предусмотрены, как в обычных холтеровских мониторах, так и в специализированных стресс системах. Но различия в программной реализации данного метода в приборах разных производителей, могут вносить существенные расхождения в получаемые результаты.
В большинстве систем мониторинга ритма сердца результаты измерения RR интервалов, оцениваются как временные ряды, и к ним для исключения периодических колебаний, применяются методы спектрального анализа (преобразования Фурье). Между тем RR-интервал — это интервал времени от зубца R до зубца R или расстояние между ними на ЭКГ (по оси х). Ввиду этого,
created by free version of
DociFreezer
ряд интервалов RR на интервалограмме не вполне может рассматриваться как временной ряд, так как аргументом является не время, а порядковые номера кардиоцикла, длительность которых непостоянна, что приводит к неравномерности шага временного ряда. Поэтому предлагаются, геометрические методы или интерполяция, и хоть это не совсем корректно, но методы анализа временных рядов применяются к ЯЯ-интервалам достаточно широко.
Как и большинство других видов анализа, анализ временных рядов предполагает, что данные содержат систематическую составляющую и случайный шум, который затрудняет обнаружение регулярных компонент. Большинство методов исследования временных рядов включает различные способы фильтрации шума, позволяющие увидеть регулярную составляющую более отчетливо [1].
Разрабатываемая методика предполагает следующие способы фильтрации шума: 1) исключение погрешностей и ошибок измерений (в случае воздействия помех или пропуска сигналов из-за ненадежной фиксации датчика) и 2) исключение случайного шума и трендов, образуемых циркадной ритмикой, сменой позы и другой физической активностью.
Для исключения ошибок измерений, получаемых в процессе мониторинга ритма сердца, предусмотрена предварительная оценка результатов, и их корректировка. Из рядов данных исключаются измерения, находящиеся вне заданного интервала, границы которого назначаются как предельно возможные в условиях исследования значения измеряемых параметров. Также предварительная обработка предусматривает применение динамического спектра с возможность установки времени начала и конца стационарных участков и исключения из анализа части данных начального и конечного этапов измерений, и нестационарных участков, где имелась физическая двигательная активность, с возможностью замены некоторых измерений, средним значением ряда.
^ сгеа!ес1 Ьу ^ее уетоп
д РооРгеегег
В целях исключения периодическом составляющей циркаднои ритмики, и случайного шума и трендов, образуемых сменой позы, не выделенных предварительной оценкой измерений:
• обеспечиваем получение временного ряда с равномерным шагом применением преобразования Фурье не к интервалограмме, а к функции п^) (изменяя данные ряда с RR-интервалов на соответствующие значения пульса);
• применяем методы вычисления корреляционных функций и спектральных плотностей. Данную процедуру подробно рассмотрим ниже.
Главной задачей вычисления корреляционных функций и спектральных плотностей является исключение тренда и сглаживание высокочастотных компонентов. Для этого экспериментально полученная последовательность RR интервалов — хкоторую можно представить в виде суммы медленно меняющейся нестационарной составляющей хск и случайной составляющей 8где [ = 1, N ^ — число RR интервалов последовательности) подвергается спектральному анализу. Нестационарная составляющая практически тождественна рассмотренным выше медленным изменениям частоты пульса. Для получения информативного спектра случайной составляющей £ можно воспользоваться вычислениями нормированной корреляционной функцию R(К) относительно скользящего среднего, которая тождественна хск
В данном случае, скользящая средняя несет в себе ряд спектральных
составляющих, исключая их из общего спектра х Эти искажения можно легко
учесть, поскольку скользящее среднее последовательных отсчетов значений к
некоторой гармоники с периодом Т (выраженным числом RR интервалов, и
1
весовыми коэффициентами — (где т — интервал усреднения)) дает в хск
гармонику с тем же периодом Т, но с амплитудой в l раз меньше амплитуды первоначального ряда.
. пМ 1 sln-^г-
1=--5Т
т
^ сгеа!ес1 Ьу ^ее уетоп
д Роа^геегег
Спектр остатка е(1) при этом полностью сохраняет все гармоники с периодом Т, а амплитуда гармоник с т < Т < 1,7 т уменьшается не более, чем в два раза. Это обуславливает выбор интервала усреднения по длительности наибольшего периода Тмакс в процессе х (1), который интересен для исследования, из условия:
Т 1
. А макс ±
т >
1,7 1,7 Д
мин
Рмин — минимальная частота в спектре х (1). Скользящее среднее вычисляем по формуле:
*ск(1) = ^^-рХ^+Я
где Р = —_—, а Р+1< £ < N — Р,т.к.хск(1) не может вычисляться для
первых и последних членов ряда х (1).
Определяем нормированную корреляционную функция Я(К):
п_р_к
Н(К) = а2(Ы — т —К) X {[Х1 — — + к)]}
Ь=р
где а2 = — хск(0]2 — дисперсия процесса х(1) относительно
скользящего среднего, а К - порядковый номер значения Я(К). Число N (общее число значений нормированной корреляционной функции) обычно выбирают не менее (№2р)/ъ при ъ = 2+5, но не более 80—120, т. к. дальнейшее увеличение числа N мало влияет на точность измерений спектра. Вычисляем спектр 8(и) в дискретные моменты:
ык
1 г V"1 2пп
Б(п)=-[1 + 2 ) И (к) соБ——к]
п ^к к=1
^ сгеа!ес1 Ьу ^ее уетоп of
д Роа^геегег
я
где п — дискретная частота 0,1,2., связанная частотой Рп= — и круговой
ык
частотой ып =
2 пп
Ык
При данном способе выделения тренда, полученные спектры характеризуются большим число пиков, и появлением, так называемых боковых пиков, обусловленных незатухающим характером коррелограммы. Ввиду чего целесообразно сглаживание таких спектров, способом скользящей средней (или медианное сглаживание) с интервалом усреднения Ш8 = 3^5 или Р8 = 1^2.
В результате соответствующих преобразований выражение для вычисления спектра имеет вид:
ад^ 5(П + ;) для Р5 + 1 <п<М5 — РБ
Полученный спектр достаточно полно отражает частоты, для которых в коррелограмме укладываются не меньше двух периодов, т. е. частоты с периодом
"Т ^ и ^2
Тмакс< — ИЛИ ^мин > —.. 2
Требования к интервалу сглаживания и длительности коррелограммы находим из соотношений:
1 2 ™ ^ ——; ик =
мин
2! 2!
N > —--+2 V ~ —--+ т
1 мин 1 мин
Поскольку 8(п) является периодической функцией с периодом 2п. Т.к. cos[(2п — а)у] = cos а], спектр симметричен относительно 1 = п. Поэтому его достаточно найти только для таких п, при которых аргумент косинуса в формуле меняется от 0 до п, т. е. 0< п < Ык/2. Число N8=^/2 определяет наибольшее число точек, для которых необходимо вычислять спектр. Для оценки реальной интенсивности каждой гармоники полученного дискретного спектра достаточно ординаты умножить на а2, а в аргументе дискретную частоту п заменить на
где х — среднее арифметическое ряда х(1 /). Тогда аргумент выражен в Гц,
N кх
2
а ордината в а2.
Для наглядного представления веса каждой гармоники целесообразно вычислить нормированный спектр:
, 100$(п)
В этом случае каждая ордината спектра характеризует долю общей интенсивности спектра (в %), приходящуюся на данную гармонику.
Таким образом, предложенные подходы и методы спектрального анализа позволяют корректно выделять периодические циркадные компоненты и тренды, обусловленные изменениями положения тела, позы, двигательной активности. Что обеспечивается получением временного ряда с равномерным шагом, или получением размерности анализируемых периодов в герцах и в процентах спектральной плотности.
Применение методики на практике позволяет более точно учитывать динамику изменения отдельных составляющих спектра ЯЯ-интервалов, и соответственно учитывать их влияние при оценке функциональных состояний, в том числе и индивидуальных реакций на эмоциональную и умственную нагрузку.
Список литературы:
1. Андерсон Т. Статистический анализ временных рядов М. Мир 1976, — 756 с.
2. Task force report of European Society of Cardiology and the North American Society of Pacing and Electrophysiology. Heart rate variability: standards of measurement, physiological interpretation and clinical use. Circulation 1996: — р. 1043—1065.
^ created by free version of
S DociFreezer
