Введение
Анализ связей проводится очень давно во многих областях науки. Это позволяет глубже взглянуть на структуру и механизмы взаимодействия, возникающие между отдельными элементами, включёнными в сложный комплекс, функционирующий как единое целое, например, сердечно-сосудистую систему, являющуюся объектом интереса со стороны многих современных исследователей [1–9]. Анализ направленных связей — полезный и чувствительный метод, который может дать ценную информацию, которая позволяет улучшить понимание устройства систем регуляции сердечно-сосудистой системы, что важно для решения таких прикладных задач, как диагностика патологий и реабилитация после перенесения заболевания.
Но его адаптация, применение и интерпретация результатов при анализе таких сложных систем, как системы биологической природы, по их зашумленным сигналам требует проведения специализированных исследований.
Поскольку детальная информация о процессах регуляции сердечно-сосудистой системы изучена недостаточно, то тестирование метода анализа её сигналов на её сигналах не представляется целесообразной тактикой. В связи с этим, ранее была предложена математическая модель сердечно-сосудистой системы, которая показала свои преимущества по сравнению с другими известными моделями [3]. На временных рядах данной математической модели были изучены особенности применения методов диагностики направленных связей, поскольку внутреннее устройство исследуемого объекта известно. После чего будет возможно провести сопоставление полученных на модели результатов с результатами, полученными при анализе экспериментальных данных.
Материал и методы
Методы
Для моделирования фазовой динамики из исходных сигналов требуется выделить их фазы [10–12], для чего в работе было использовано преобразование Гильберта. Сравнительно простой, не вносящий дополнительных искажений метод, позволяющий точно выделить фазу, особенно хорошо работающий для сигналов, фильтрованных в относительно узкой полосе частот.
Для оценки направленных связей в работе использовался метод, основанный на моделировании фаз временных рядов взаимодействующих систем, подробно рассмотренный в работах [13, 14]. Суть метода состоит в оценке меры воздействия по улучшению прогноза, производимого использованием совместной модели фазовой динамики исследуемых систем по сравнению с прогнозом, полученным с использованием индивидуальной модели. Превышение оценки меры воздействия (в дальнейшем обозначаемые G) в направлении воздействия со стороны одной из систем на другую может быть интерпретировано как направление преобладающего влияния в этом случае, положение максимума на зависимости оценки меры воздействия может быть использовано для оценки задержки в связи, соответствующей этому направлению воздействия.
Используемые для тестирования метода временные ряды были получены при помощи математической модели сердечно-сосудистой системы человека, воспроизводящей основной сердечный ритм, регуляторные воздействия, артериальное давление, а также учитывающей влияние дыхания на процессы регуляции, предложенная в работах [15, 16].
Данные
Экспериментальные данные математической модели, включающей четыре дифференциальных уравнения первого порядка, представляли собой модельные временные ряды электрокардиограмм (ЭКГ) и систолического артериального давления (АД) длительностью 60 минут при частоте дискретизации 250 Гц, перевыбранные в дальнейшем до 10 Гц, интегрированные методом Эйлера с шагом 0,004.
В качестве данных для сравнения были использованы 12 записей ЭКГ и фотопетизмограмм (ФПГ) длительностью 60 минут, полученные со здоровых субъектов, находившихся в состоянии лёжа. Регистрация сигналов осуществлялась с помощью серийного энцефалографа-анализатора ЭЭГА-21/26 «Энцефалан-131-03» (модель 10) с частотой дискретизации 250 Гц, перевыбранные в дальнейшем до 10 Гц. Полоса пропускания аналогового сигнала составляла 0,05–100 Гц.
Экспериментальные и модельные сигналы подвергались предварительной обработке, включавшей выделение эквидистантной кардиоинтервалограмм (КИГ) из ЭКГ и полосовую фильтрацию в диапазоне 0,05–0,15 Гц, проводимую в соответствии с методическими указаниями [17, 18].
Результаты
Были рассчитаны оценки меры воздействия одного сигнала на другой, полученные при различных значениях аддитивных равномерно распределённых фазовых шумов. На Рисунке 1. приведены оценки меры воздействия при отсутствии шума. 5% уровень значимости, оцененный с помощью суррогатных данных, полученных из исходных сигналов путем случайного задания фаз, отложен штриховой линией.
На графике мы можем наблюдать, что детектирована значимая двунаправленная связь, причём для направления воздействия КИГ→АД уровень воздействия выше, чем в обратном направлении.
На Рисунке 2 приведены зависимости максимальных оценок мер воздействия, полученных с помощью метода моделирования фазовой динамики от уровня фазового шума. Видно, что абсолютные значения оценок уменьшаются с ростом уровня фазового шума. Усреднённые по ансамблю значения оценок уровней воздействия, полученные по экспериментальным записям здоровых субъектов принимают значения GКИГ→АД=43,3±1,8×10-5 и GАД→КИГ=42,5±5,6×10-5 (цифры приведены в формате «среднее ± стандартное отклонение»). Диапазон, соответствующий минимальному и максимальному с учётом стандартного отклонения уровню оценок мер воздействия, приведён на Рисунке 2 красной и синей линиями. Уровень фазового шума шума, для которых оценки мер воздействия по модельным данным принимает значения из диапазона, составляет 10–13%. Зависимости самих оценок воздействия, полученных для экспериментальных данных и модельных данных с 10% фазовым шумом для сравнения приведены на Рисунках 3 и 4.
Оценки экспериментальные и с шумом 10% приведены на Рисунке 3.
На Рисунке 5. приведены зависимости восстановленных оценок задержек в воздействии GКИГ→АД и GАД→КИГ от уровня фазового шума. Из графика видно, что данные оценки приобретают случайный характер для уровней фазовых шумов свыше 17%. Кроме того, максимальные значения оценок мер воздействия (Рисунок 2) для шумов этой интенсивности принимают отрицательные значения, что не позволяет детектировать связь.
Рисунок 1. Зависимости оценок мер воздействия GКИГ→АД и GАД→КИГ от пробного времени задержки d для модельных данных при отсутствии фазового шума. Уровень значимости отложен штриховой линией.
Рисунок 2. Зависимости максимальных значений оценок мер воздействия GКИГ→АД и GАД→КИГ от уровня фазового шума.
Рисунок 3. Типичные зависимости оценок мер воздействия GКИГ→АД и GАД→КИГ от пробного времени задержки d для экспериментальных данных одного из испытуемых.
Рисунок 4. Зависимости оценок мер воздействия GКИГ→АД и GАД→КИГ от пробного времени задержки d для модельных данных в присутствии 10% фазового шума.
Рисунок 5. Зависимости восстановленных оценок задержек в воздействии GКИГ→АД и GАД→КИГ от уровня фазового шума.
Обсуждение
Определение направления связи является сложной, но решаемой задачей. Основной проблемой при анализе нестационарных и зашумленных сигналов сердечно-сосудистой системы, глубинные свойства которых до конца не изучены, является интерпретация результатов использования сложных методов. Результаты исследований простых эталонных систем напрямую сопоставлять с результатами анализа экспериментальных данных невозможно из-за сложной природы и богатой динамики, присущей объектам биологического происхождения. Именно использование развитых и как можно более близких к исследуемым объектам математических моделей призвано помочь с решением этой проблемы. Конечно, любая математическая модель является упрощением по-сравнению с реальной сложной системой биологической природы, но если модель охватывает важные качественные особенности, то она вполне пригодна для тестирования методов. Используемая модель включает в себя много существенных особенностей, по сравнению с другими известными (особенно, в плане соответствия спектральных и фазовых характеристик). Поэтому мы считаем, что она позволяет использование для тестирования метода, основанного на моделировании фаз временных рядов взаимодействующих систем, предназначенного для анализа связей между взаимодействующими автоколебательными системами, для которых наблюдались явления фазового и частотного захвата [19].
Уровень шума модельных данных, при котором результаты похожи на экспериментальные составляет примерно 10–12%, возможно это оценка типичного измерительного шума в экспериментах, обусловленного несовершенством измерительной аппаратуры либо недостаточно изученными процессами биологического происхождения [20].
Заключение
Метод, основанный на моделировании фаз временных рядов взаимодействующих систем, предназначенный для анализа сигналов, был протестирован при анализе реализаций модели, воспроизводящей основной сердечный ритм, регуляторные воздействия, артериальное давление, а также учитывающий влияние дыхания на процессы регуляции при разных уровнях фазовых шумов. Выяснилось, что при уровне шума порядка 10% значения оценок мер воздействия снижаются до уровня, наблюдаемого в экспериментальных данных. При уровне шума до 15% результаты позволяют оценить направленные связи. При значениях шума свыше 17% определение задержки и даже наличия воздействия по модельным данным не представляется возможным.
Конфликт интересов
Работа выполнена при поддержке гранта РФФИ № 17-02-00307.
- Stefanovska A., Luchinsky D.G., McClintock P.V.E. Modelling couplings among the oscillators of the cardiovascular system. Phys Meas 2001; 22(3): 551–564. http://dx.doi.org/10.1088/0967-3334/22/3/312.
- Боровкова Е.И., Караваев А.С., Пономаренко В.И. и др. Сопоставление методов диагностики фазовой синхронизованности по тестовым данным, моделирующим нестационарные сигналы биологической природы. Известия Саратовского Университета. Новая серия. Серия Физика 2015; 15(3): 36–42. http://dx.doi.org/10.18500/1817-3020-2015-15-336-42.
- Ишбулатов Ю.М., Караваев А.С., Пономаренко В.И. и др. Модель системы автономной регуляции сердечно-сосудистой системы с контуром барорефлекторного контроля среднего артериального давления в виде автогенератора с запаздыванием. Известия Саратовского Университета. Новая серия. Серия Физика 2015; 15(2): 32–38. http://dx.doi.org/10.18500/1817-3020-2015-15-2-32-38.
- Хорев В.С., Киселев А.Р., Шварц В.А. и др. Исследование запаздывания в связи между контурами регуляции сердечно-сосудистой системы у здорового человека методом моделирования фазовой динамики. Известия Саратовского Университета. Новая серия. Серия Физика 2016; 16(4): 227–237. http://dx.doi.org/10.18500/1817-3020-2016-16-4-227-237.
- Kralemann B., Frühwirth M., Pikovsky A., et al. In vivo cardiac phase response curve elucidates human respiratory heart rate variability. Nat Commun 2013; 4: 2418. http://dx.doi.org/10.1038/ncomms3418.
- Burgess D.E., Hundley J.C., Brown D.R., et al. First-order differential-delay equation for the baroreflex predicts the 0.4-Hz blood pressure rhythm in rats. Am J Physiol 1997; 273: 1878–1884. http://dx.doi.org/10.1152/ajpregu.1997.273.6.R1878.
- Красников Г.В., Танканаг А.В., Коняева Т.Н. и др. Оценка изменений в системах регуляции кровотока в коже человека при локальном нагреве. Российский физиологический журнал им. И.М. Сеченова 2007; 93(4): 394–401. http://www.ncbi.nlm.nih.gov/pubmed/17654864.
- Тихонова И.В., Танканаг А.В., Чемерис Н.К. Возрастные особенности динамики амплитуд колебаний кровотока кожи в процессе постокклюзионной реактивной гиперемии. Физиология человека 2010; 36(2): 114–120. http://www.ncbi.nlm.nih.gov/pubmed/20432699.
- Allen J., Di Maria C., Mizeva I., et al. Finger microvascular responses to deep inspiratory gasp assessed and quantified using wavelet analysis. Phys Meas 2013; 34(7): 769. http://dx.doi.org/10.1088/0967-3334/34/7/769.
- Pikovsky A., Rosenblum M., Kurths J. Phase synchronization in regular and chaotic systems. Int J Bif Chaos 2000; 10(10): 2291–2305. http://dx.doi.org/10.1142/S0218127400001481.
- Кузнецов А.П., Станкевич Н.В. Синхронизация генераторов квазипериодических колебаний. Нелинейная динамика 2013; 9(3): 409–419. http://dx.doi.org/10.20537/nd1303002.
- Kiselev A.R., Karavaev A.S., Gridnev V.I., et al. Method of estimation of synchronization strength between low-frequency oscillations in heart rate variability and photoplethysmographic waveform variability. Russ Open Med J 2016; 5(1): e0101. http://dx.doi.org/10.15275/rusomj.2016.0103.
- Rosenblum M., Pikovsky A. Detecting direction of coupling in interacting oscillators. Phys Rev E 2001; 64: 045202(R). http://dx.doi.org/10.1103/PhysRevE.64.045202.
- Smirnov D.A., Bezruchko B.P. Revealing mutual influence of oscillatory systems from the observation data. Radiophysics and Quantum Electronics 2013; 55: 662–675. http://dx.doi.org/10.1007/s11141-013-9404-6.
- Karavaev A.S., Ishbulatov J.M., Ponomarenko V.I., et al. Model of human cardiovascular system with a loop of autonomic regulation of the mean arterial pressure. J Am Soc Hypertens 2016; 10(3): 235–243. http://dx.doi.org/10.1016/j.jash.2015.12.014.
- Караваев А.С., Ишбулатов Ю.М., Киселев А.Р. и др. Модель сердечно-сосудистой системы человека с автономным контуром регуляции среднего артериального давления. Физиология человека 2017; 43(1): 70–80. http://dx.doi.org/10.7868/S0131164616060096.
- Heart rate variability: Standards of measurement, physiological interpretation and clinical use. Task Force of the European Society of Cardiology and the North American Society of Pacing Electrophysiology. Circulation 1996; 93(5): 1043–1065. https://www.ncbi.nlm.nih.gov/pubmed/8598068.
- Баевский Р.М., Иванов Г.Г., Чирейкин Л.В. и др. Анализ вариабельности сердечного ритма при использовании различных электрокардиографических систем. Вестник аритмологии 2002; (24): 65–86. https://elibrary.ru/item.asp?id=9166688.
- Караваев А.С., Киселев А.Р., Гриднев В.И. и др. Фазовый и частотный захват 0,1 Гц-колебаний в ритме сердца и барорефлекторной регуляции артериального давления дыханием с линейно меняющейся частотой у здоровых лиц. Физиология человека 2013; 39(4): 93–104. http://www.ncbi.nlm.nih.gov/pubmed/25486835.
- Klimovitch G.V. Near-carrier oscillator spectrum due to flicker and white noise. Proceedings of ISCAS 2000; 1: 703–706. http://dx.doi.org/10.1109/ISCAS.2000.857192.
Поступила в редакцию 2 марта 2018. Принята в печать 30 марта 2018.
© 2018, Хорев В.С., Ишбулатов Ю.М., Караваев А.С., Попова Ю.В., Киселев А.Р., Безручко Б.П.
© 2018, Кардио-ИТ
Ответственный автор: Киселев Антон Робертович. Адрес для переписки: НИИ кардиологии, Саратовский ГМУ им. В.И. Разумовского, 112, ул. Большая Казачья, г. Саратов, 410012, Россия. E-mail: kiselev@cardio-it.ru.