Вариабельность сердечного ритма
Вариабельность сердечного ритма (ВСР), также Вариабельность ритма сердца (англ. Heart rate variability, HRV), — физиологическое явление, проявляющееся в изменении интервала между началами двух соседних сердечных циклов. Оценивается по изменениям интервала времени между соседними сердечными сокращениями (сердцебиениями).
Реже используются следующие синонимичные термины: «вариабельность длины цикла», «вариабельность RR» (где R — точка, соответствующая пику комплекса QRS волны электрокардиограммы (ЭКГ) и RR — интервал между последовательными R) и «вариабельность периода сердца». Wiki
Методы математического анализа
Статистические методы временного домена
- SDNN - стандартное отклонение интервалов NN. Часто рассчитывается в течение 24 часов.
\[\overline{X} = \frac{1}{n} \sum_{i=1}^n x_{i}\]
\[SDNN = \sqrt{ \frac{1}{n} \big(x_i - \overline{X} \big)^2 }\] - RMSSD - «среднеквадратичное значение последовательных различий», квадратный корень среднего значения квадратов конечных разностей соседних NN.
\[diff_{RR} = diff(RR)\]
\[RMSSD = \sqrt{ \frac{1}{n} diff^2_{RR_i} }\] - SDSD - «стандартное отклонение последовательных различий», стандартное отклонение последовательных различий между соседними NN.
\[S(X) = \sqrt{ \frac{1}{n} \big(x_i - \overline{X}\big)^2 }\]
\[SDSD = S(diff(RR))\] - NN50 - количество пар последовательных NN, которые отличаются более чем на 50 мс.
\[NN50 = \sum_{i=1}^n diff_{RR_i} > 50\] - pNN50 - доля NN50, деленная на общее количество NN.
\[pNN50 = 100 \frac{NN50}{length(RR)}\] - NN20 - количество пар последовательных NN, которые отличаются более чем на 20 мс.
\[NN20 = \sum_{i=1}^n diff_{RR_i} > 20\] - pNN20 - доля NN20, деленная на общее количество NN.
\[pNN20 = 100 \frac{NN20}{length(RR)}\]
Методы частотного домена
При анализе ВСР под спектральной плотностью мощности (СПМ) временной последовательности понимают СПМ стационарного (в широком смысле) случайного процесса, реализацией которого является эта последовательность. Следует иметь в виду, что любой из применяемых спектральных методов — это метод оценки СПМ, а не точного её построения. Если целью исследования является выявление межгрупповых различий в спектральных характеристиках ВСР обследуемых из разных групп, то оценки СПМ ВСР всех больных всех групп надо проводить при помощи одного и того же метода. Стандартная длина участка хронокардиограммы для спектрального анализа составляет 256 кардиоциклов, что соответствует промежутку времени 3,5—5 минут; жёсткие требования стационарности делают невыполнимым спектральный анализ 24-часовой хронокардиограммы.
Методы в частотной области используются для подсчитывания количество NN-интервалов, которые соответствуют каждой полосе частот. Стандартами рекомендовано различать следующие полосы (компоненты) частот:
- высокой частоты (ВЧ, HF) от 0,15 до 0,4 Гц,
- низкой частоты (НЧ, LF) от 0,04 до 0,15 Гц и
- очень низкой частоты (ОНЧ, VLF) от 0,003 до 0,04 Гц и
- ультранизкочастотный (УНЧ, ULF) менее 0,003 Гц.
Измерение мощности ОНЧ, НЧ, ВЧ обычно осуществляется в абсолютных единицах мощности (мс2), но НЧ и ВЧ могут быть дополнительно выражены в нормализованных единицах, которые показывают относительный вклад каждого из компонентов в пропорции к общей мощности за вычетом ОНЧ-компонента. УНЧ будет получен лишь при использовании спектрального анализа для всего 24-х часовой периода наблюдения.
Расчет данных показателей производится с помощью либо чистого БПФ, с последующим получением спектральной плотности мощности (СПМ) руками, либо с применением уже готовых методов расчета СПМ. Их много, но на мой вкус стоит использовать 2 из них, т.к. они уже реализованы в пакете NumPy: метод Уэлча и Lomb method (спектральный анализ методом наименьших квадратов).
Использование указанных методов объясняется тем, что спектр реальных измеряемых сигналов (в нашем случае RR-интервалов) имеет стохастический характер и в большей степени состоит просто из шума. Чтобы избавится от шума, необходимо просто как-то усреднить сигнал. Это и делают указанные методы. Например, метод Уэлча разбивает оригинальный сигнал на несколько частей (таким образом, несколько увеличивает размер вашей выборки сигналов) и усредняет их спектры. Сокращая сигналы, к которым вы применяете БПФ, и усредняя, спектры передаются по низким частотам, и поэтому пики не так уж и узки: это цена, которую нужно заплатить за избавление от шума, т.е. мы теряем часть полезной информации о сигнале, но сжираем огромное количество шума, содержащегося в этом сигнале.
Далее, полученный спектр банально интегрируется в нужных диапазонах частот (см. выше) с помощью любого численного метода. Например, метод трапеций, который так же реализован уже в NumPy.