Захват и удержание показателей вариабельности сердечного ритма метрономизированной частотой дыхания у здоровых добровольцев Текст научной статьи по специальности «Клиническая медицина»

5. Диагностическое и прогностическое значение маркера системного воспаления С-реактивного протеина у больных с острыми коронарными синдромами / А.Н. Пархоменко, Я.М. Лутай, Г.В. Пономарева [и др.] //Укр. Бан. Сеть. - 2001.

6. Malarstig A. Raised interleukin-10 is an indicator of poor outcome and enhanced systemic inflammation in patients with acute coronary syndrome / A. Malarstig, P. Eriksson, A. Hamsten [et al.] // Heart. - 2008.

- Vol.94. - P.724-729.

7. TH 1/TH2 imbalance, measured by circulating and intracytoplasmic inflammatory cytokines-immuno-logical alterations in acute coronary syndrome and stable coronary artery disease / P. Szodoray, O. Timar, K. Veres [et al.] // Scand J Immunol - 2006. - № 64 (3). - P. 336-344.

8. Predictive value of CRP and troponin T in patients with unstable angina: a comparative analysis. CAPTURE Investigators. / C. Heeschen, C.W. Hamm, J. Bruemmer [et al.] // JACC. - 2000. - № 35. -P. 1535-1542.

9. Динамика воспалительного процесса у больных с острым коронарным синдромом и больных со стабильной стенокардией. Сообщение 1. Биохимические и иммунологические аспекты / В.В. Кухарчук, К.А. Зыков, В.П. Масенко [и др.] // Кардиологический вестник. - 2007 - Том 02, - № 2.

© Копиця М.П., Литвін О.І., 2009

УДК: 612.213

ЗАХВАТ И УДЕРЖАНИЕ ПОКАЗАТЕЛЕЙ ВАРИАБЕЛЬНОСТИ СЕРДЕЧНОГО РИТМА МЕТРОНОМИЗИРОВАННОЙ ЧАСТОТОЙ ДЫХАНИЯ У ЗДОРОВЫХ ДОБРОВОЛЬЦЕВ

А.Л. Кулик1, А. К. Задерихин2, В.И. Шульгин2, А.В. Мартыненко1, Н.И. Яблучанский1 'Харьковский национальный университет имени В.Н. Каразина, Украина

Национальный Аэрокосмический Университет имени Н.Е. Жуковского «Харьковский авиационный институт», Украина

На 4 здоровых добровольцах в возрасте от 22 до 29 лет (2 мужчины и 2 женщины) с индексом массы тела от 21 до 23 кг/м2 изучены закономерности частотной зависимости спектральных показателей вариабельности сердечного ритма (ВСР) при различных паттернах дыхания. Показано, что при одно-направленых изменениях (рост, снижение) частоты метрономизированного дыхания происходит захват и удержание контура регуляции, реализующиеся через изменение спектральных показателей ВСР. Захвата и удержания ЧСС не происходит. Метрономизированное дыхание на частотах 6-8 дыханий в минуту приводит к усилению ЬБ спектра ВСР, а на частотах 10-12 дыханий в минуту - НЕ7 спектра ВСР; дыхание на частоте 9 дыханий в минуту приводит к активации обеих областей спектра ВСР. Для получения картины захвата и удержания контура регуляции достаточно использования длительности шага и емкости буфера в 1 минуту со скоростью изменения частоты дыхания в 1 дыхание в минуту. Применение паттернов дыхания с этими параметрами может быть удобным инструментом быстрого и эффективного подбора оптимальной частоты дыхания для «тонкой подстройки» работы регуляторных систем.

КЛЮЧЕВЫЕ СЛОВА: вариабельность сердечного ритма, биообратная связь, частота дыхания

ЗАХОПЛЕННЯ ТА УТРИМАННЯ ПОКАЗНИКІВ ВАРІАБЕЛЬНОСТІ СЕРЦЕВОГО РИТМУ МЕТРОНОМІЗІРОВАНОЮ ЧАСТОТОЮ ДИХАННЯ У ЗДОРОВИХ ДОБРОВОЛЬЦІВ

ОЛ. Кулик1, О.К Задеріхін2, ВЛ. Шульгін2, О.В. Мартиненко1, МЛ. Яблучанський1 'Харківський національний університет імені В.Н. Каразіна, Україна

Національний Аерокосмічний Університет імені М.Є. Жуковського «Харківський авіаційний інститут», Україна

На 4 здорових добровольцях у віці від 22 до 29 років (2 чоловіки і 2 жінки) з індексом маси тіла від 21 до 23 кг/м2 вивчено закономірності частотної залежності спектральних показників варіабельності серцевого ритму (ВСР) при різних патернах дихання. Показано, що при однонаправлених змінах (зростання, зниження) частоти метрономізірованого дихання відбувається захоплення і утримання контуру регуляції, що реалізуються через зміну спектральних показників ВСР. Захоплення і утримання ЧСС не відбувається. Метрономізіроване дихання на частотах 6-8 дихань в хвилину приводить до посилення и спектру ВСР, а на частотах 10-12 дихань в хвилину - НЕ7 спектру ВСР; дихання на частоті 9 дихань в хвилину приводить до активації обох областей спектру ВСР. Для отримання картини захоплення і утримання контура регуляції достатньо використання тривалості кроку і місткості буфера в 1 хвилину із швидкістю зміни частоти дихання в 1 дихання в хвилину. Застосування патернів дихання з цими параметрами може бути зручним інструментом швидкого і ефективного підбору оптимальної частоти дихання для «тонкого підстроювання» роботи регуляторних систем.

КЛЮЧОВІ СЛОВА: варіабельність серцевого ритму, біофідбек, частота дихання

CAPTURE AND WITHHOLDING OF HEART RATE VARIABILITY INDEXES WITH PACED BREATHING IN HEALTHY VOLUNTEERS

A.L. Kulik1, O.K. Zaderykhin2, VI. Shulgin 2, A. V. Martynenko1, M.I. Iabluchanskyi1

'V.N. Karazin Kharkov National University, Ukraine

2National Aerospace University, named by N. E. Zhukovskii «Kharkov Aviation Institute», Ukraine

At 4 healthy volunteers in age from 22 to 29 years (2 men and 2 women) with the body mass index from 21 to 23 kg/m2 frequency dependence of spectral indexes of heart rate variability (HRV) on at different paced breathing patterns were studied. It is shown that at the unidirectional breath rate changes (growth, decline) a capture and withholding of regulation contours occurs, realized through the change of spectral indexes of HRV. A capture and withholding of heart rate, however, did not occur. Paced breathing with 6-8 breathings per minute results in strengthening of LF domain of HRV, and 10-12 breathings per minute - HF domain of HRV; paced breathing with 9 breathings per minute result in activation of both HRV domains. To recieve adequate pattern of capture and withholding of regulation contour usage of duration of breathing step and capacity of HRV buffer - 1 minute and a speed of change of breathing frequency - 1 breathing per minute is enough. Application of breathing patterns with these parameters can be the handy instrument of rapid and effective selection of optimum frequency of breathing for «fine tuning» of regulatory systems activity.

KEY WORDS: HRV biofeedback, breathing frequency

Состояние здоровья человека во многом зависит от качества и режимов функционирования регуляторных систем [1, 2]. Одним из распространенных неинвазивных методов их оценки является вариабельность сердечного ритма (ВСР) [3, 4].

Среди различных методов вмешательства в состояние регуляторных систем одним из эффективных является управление дыханием [5, 6].

С помощью метрономизированного дыхания на соответствующих частотах можно избирательно вызывать повышение отдельных звеньев спектра ВСР, чем существенным образом влиять на регуляторные процессы [1, 3, 7].

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

Притом, что зависимость ВСР от частоты дыхания изучалась [3-7], в выполненных работах внимание концентрируется на зависимости параметров ВСР от частоты дыхания на изолированных значениях. Возможности же изменений параметров ВСР под влиянием однонаправленных изменений дыхания не изучались.

Исследование выполнено в рамках НИР ХНУ имени В.Н. Каразина «Разработка и исследование системы автоматического управления вариабельностью сердечного ритма», № госрегистрации 010911000622.

Цель исследования. Установить индивидуальные особенности и общие закономерности частотной зависимости ВСР у здоровых добровольцев при различных программах (паттернах) дыхания для создания эффективных техник биообратной связи в за-

дачах повышения качества регуляторных систем человека.

МАТЕРИАЛЫ И МЕТОДЫ

Обследовано 4 здоровых добровольца в возрасте от 22 до 29 лет (2 мужчины и 2 женщины), с индексом массы тела от 21 до 23 кг/м2. Частота сердечных сокращений варьировала от 81 до 85 в минуту, артериальное давление находилось в диапазоне -100/60 и 120/80 мм рт.ст.

У всех испытуемых с помощью компьютерного диагностического комплекса «Саг-сПоЬаЬ 2009» («ХАИ-Медика») с частотой дискретизации сигнала 1000 Гц проводилось по 8 мониторных записей ЭКГ на семиминутных интервалах. Записи производились в положении сидя при двух вариантах емкости буфера - 256 интервалов и 1 минута, и при различных параметрах изменения частоты дыхания (ЧД) - возрастании с 6 до 12 дыханий в минуту и убывании с 12 до 6 в минуту. При этом при составлении программы изменения ЧД использавалось два варианта шага изменения ЧД - 0,5 и 1 дыхание в минуту с длительностью каждого шага - 30 секунд и 1 минута, соответственно. Ритм дыхания задавался метрономом, встроенным в программу «СагсПоЬаЬ 2009».

Полученные ритмограммы при помощи быстрого преобразования Фурье раскладывали на три типа волн: медленные (УЬБ) -частотой от 0,0033 до 0,05 Гц - преимущественно связаны с терморегуляцией, гуморальной (калликреинкинировая, ренинангио-тензивая, гормональные, иные) регуляцией и симпатическим звеном вегетативной нервной системы; средние (ЬГ) - от 0,05 до 0,15 Гц - с симпатическим звеном вегетативной нервной системы и быстрые (НГ) - от 0,15 Гц до 0,40 Гц - с парасимпатическим звеном вегетативной нервной системы.

На каждом интервале записей изучались: общяя мощность (ТР, мс2) спектра ВСР, мощности спектров доменов низких (УЬР, мс2), средних (ЬБ, мс2) и высоких (НБ, мс2) частот, соотношение ЬР/НР как мера симпа-товагального балланса и соотношение УЬР/(ЬР+НР) как мера гуморальновегетативного баланса.

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

а) ТР ВСР, буфер 256 RR, шаг 1 дыхание в минуту

с) ТР ВСР, буфер 256 шаг 0,5 дыханий в минуту

Рис. 1. Реакция ТР

На рис. 2 представлена реакция УЬБ ВСР на изменение ЧД. На графике 2А отражено изменение УЬР ВСР при увеличении и уменьшении ЧД с шагом 1 дыхание в минуту при емкости буфера 256 ЯЯ-интср-валов и на графике 2В - при емкости буфера 1 минута. Графики 2С и 2Б, соответственно, отображают изменение УЬР при увеличении и уменьшении ЧД с шагом 0,5 дыхания в

го дыхания.

В программе Microsoft Excel 2003 по каждому добровольцу проводилось построение графиков показателей ВСР в зависимости от частоты дыхания. Графики анализировались с учетом индивидуальных и общих закономерностей в изменениях исследуемых показателей.

РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ

На рис. 1 представлена реакция ТР ВСР на изменение ЧД. На графике 1А отражено изменение ТР при увеличении и уменьшении ЧД с шагом 1 дыхание в минуту при емкости буфера 256 RR-интервалов и на графике 1В -при емкости буфера 1 минута. Соответственно, графики 1C и 1D отображают изменение ТР при увеличении и уменьшении ЧД с шагом 0,5 дыхания в минуту с длительностью шагов в 30 секунд при тех же емкостях буфера. Из графиков видно, что при обоих шагах ЧД с ее увеличением наблюдается снижение ТР ВСР, а при снижении - увеличение. Степень увеличения, однако, не достигает уровня значений, зафиксированных до начала увеличения ЧД. Также следует отметить, что при буфере емкостью 1 минута регистрируются более низкие значения ТР ВСР (~ в 2 раза), чем при буфере емкостью 256 RR-интервалов.

Ь) ТР ВСР, буфер 1 мин, шаг 1 дыхание в минуту

«I) ТР ВСР, буфер 1 мин, шаг 0,5 дыханий в минуту ВСР на изменение ЧД

минуту с длительностью шагов в 30 секунд при тех же емкостях буфера. Из графиков видно, что, несмотря на колебательный характер, изменение УЬР ВСР следует таким же тенденциям, что и изменение ТР - повышение с ростом ЧД и снижение при уменьшении ЧД. При обоих шагах изменения ЧД с ее увеличением наблюдается снижение УЬР ВСР, а со снижением - увеличение. При

этом степень увеличения также не достигает уровня значений, зафиксированных до начала увеличения ЧД, а значения УЬБ, регист-

a) VLF ВСР, буфер 256 RR, шаг 1 дыхание в минуту

с) УЬБ ВСР, буфер 256 RR, шаг 0,5 дыханий в минуту

Рис. 2. Реакция УЬБ

На рис. 3 представлена реакция ЬБ и НБ ВСР на изменение ЧД. На графике ЗА отражено изменение и и НБ ВСР при увеличении и уменьшении ЧД с шагом 1 дыхание в минуту при емкости буфера 256 ЯЯ-интср-валов и на графике ЗВ - при емкости буфера 1 минута. Соответственно, графики ЗС и ЗБ, отображают изменение ЬБ и НБ при увеличении и уменьшении ЧД с шагом 0,5 дыхания в минуту с длительностью шагов в 30 секунд при тех же емкостях буфера. На графиках можно проследить ряд закономерностей - при увеличении ЧД вначале наблюдаются наибольшие значения ЬБ ВСР и наи-

рируемые при при буфере емкостью 1 минута почти в два раза ниже, чем при использовании буфера емкостью 256 RR-интервалов.

b) VLF ВСР, буфер 1 мин, шаг 1 дыхание в минуту

d) VLF ВСР, буфер 1 мин, шаг 0,5 дыхание в минуту ВСР на изменение ЧД

меньшие HF ВСР. По мере повышения ЧД значения LF ВСР снижаются, пока не достигают минимального уровня, на котором сохраняются. Значения HF ВСР на некотором промежутке изменения ЧД остаются низкими, возрастая на ЧД, соответствующей достижению минимума значениями LF ВСР, после чего плавно снижаются. При уменьшении ЧД наблюдается обратная картина. Как и в отношении предыдущих показателей ВСР, в буфере емкостью 1 минута имели место низкие значения ТР ВСР, чем в буфере емкостью 256 RR-интервалов.

с) ЬБ и ОТ ВСР, буфер 256 RR, шаг 0,5 дыханий в минуту сі) ЬБ и ОТ ВСР, буфер 1 мин, шаг 0,5 дыхание в минуту

Рис. 3. Реакция ЬБ и ОТ ВСР на изменение ЧД

На рис. 4 представлена реакция соотношения ЬР/НР на изменение ЧД. На графике 4А отражено изменение соотношения ЬР/НР при увеличении и уменьшении ЧД с шагом 1 дыхание в минуту при емкости буфера 256 ЯЯ-интсрвалов. а на графике 4В - при емкости буфера 1 минута. Графики 4С и 4Б, соответственно, отображают изменение соот-

ношения ЬР/НР при увеличении и уменьшении ЧД с шагом 0,5 дыхания в минуту с длительностью шагов в 30 секунд при тех же емкостях буфера. Полученные при различных параметрах графики соотношения ЬР/ НР по своей форме повторяют поведение и ВСР в соответствующих условиях.

а) Соотношение ЬР/НР, буфер 256 ЯЯ, шаг 1 дыхание

с) Соотношение ЬР/ОТ буфер 256 шаг 0,5 дыханий в минуту

Ь) Соотношение ЬР/ОТ, буфер 1 мин, шаг 1 дыхание в минуту

с1) Соотношение ЬР/ОТ, буфер 1 мин, шаг 0,5 дыханий в минуту

Рис. 4. Реакция соотношения ЬР/НР ВСР на изменение ЧД

На рис. 5 представлена реакция соотношения УЬР/(ЬР+НР) на изменение ЧД. На графике 5А отражено изменение соотношения УЬР/(ЬР+НР) при увеличении и уменьшении ЧД с шагом 1 дыхание в минуту при емкости буфера 256 ЯЯ-интсрвалов и на графике 5В - при емкости буфера 1 минута. Соответственно, графики 5С и 5Б, отображают изменение соотношения УЬР/(ЬР+НР)

при увеличении и уменьшении ЧД с шагом

0,5 дыхания в минуту с длительностью шагов в 30 секунд при тех же емкостях буфера. Представленные изменения соотношения УЬР/(ЬР+НР) подобны изменениям УЬР ВСР и носят колебательный характер с тенденцией роста с увеличением ЧД и снижением соотношения при уменьшении ЧД.

А) Соотношение УЬР/(ЬР+НР), буфер 256 Ш1, шаг 1 дыхание

С) Соотношение УЫ7(ЬР+НР), буфер 256 шаг 0,5 дыхание в минуту

В) Соотношение УЬР/(ЬР+НР), буфер 1 мин, шаг 1 дыхание в минуту

О) Соотношение УЬР/(ЬР+НР), буфер 1 мин, шаг 0,5 дыханий в минуту

Рис. 5. Реакция соотношения УЬР/(ЬР+НР) ВСР на изменение ЧД

Испытуемый 1

Испытуемая 2

Испытуемый 3

Испытуемый 1

Испытуемый d

испытуемая 4

В целом, как при шаге 1 дыхание в минуту, так и шаге 0,5 дыхания в минуту, независимо от направления изменения ЧД и емкости буфера происходит захват и удержание контура регуляции метрономизированным дыханием.

Использование более «плавного» (при шаге приращения метрономизированного дыхания на 0,1 дыхание в минуту при длительности каждого шага приращения - 6 се-

кунд) изменения ЧД недостаточно для захвата и удержания контура регуляции (рис. 6), поэтому такие параметры программы изменения дыхания были отброшены.

Что касается ЧСС, ее захвата и удержания не происходит ни при каких емкостях буфера и параметрах изменения ЧД. Для примера на рис. 7 показана реакция ЧСС на увеличение ЧД при буфере 1 мин и шаге 1 дыхание в минуту.

Частота дыхания, дых/мин

-Испытуемый 1

- Испытуемая 2 и Испытуемый 3

- Испытуемая 4

Рис. 6. Реакция ТР ВСР на увеличение ЧД (буфер 1 мин, шаг 0,1 дыхание в минуту, длительность шага - 6 секунд)

160

140

120

тмтиштгагажпіптіїшгашітмітіггаиііішпшїїпштшішштгаштиїїїїтштїїштїїтшпіпшшптіпіішіпііміпіїшіш

1 31 61 91 121 151 181 211 241 271 301 331 361 391 421 451 481 511

— Испытуемый 1 —'— Испытуемая 2 ------Испытуемый 3 -----Испытуемая 4

0

Рис. 7. Реакция ЧСС на увеличение ЧД (буфер 1 мин, шаг 1 дыхание в минуту)

Несмотря на то, что при емкости буфера 256 ЯЯ-интсрвалов регистрируются большие значения показателей ВСР, чем при емкости в 1 минуту, следует отметить, что в таком случае данные при этом получаются с задержкой в 2-2,5 минуты, независимо от направления изменения ЧД. Это наиболее заметно при изучении графиков, отражающих

реакцию показателей, претерпевающих резкое изменение на изучаемом диапазоне ЧД (Ы7, НБ и ЬР/НБ). Результаты находят объяснение в том, что буфер емкостью 256 ЯЯ-интервалов соответствует примерно 3 минутам, охватывая не только текущую ЧД, но и несколько предыдущих (2 - при шаге 1 дыхание в минуту и 5 - при шаге 0,5 дыхания в

минуту с длительностью шага в 30 секунд), несколько скрадывая изменения картины ВСР. В идеале, перекрест значений ЬБ и НБ должен происходить на частоте 9 дыханий в минуту. Данная ЧД соответствует частоте сердечных волн 0,15 Гц (граница между Ы7 и НБ доменами спектра ВСР).

В работах [1, 3, 7] проводилось изучение спектральных показателей ВСР на отдельных изолированных ЧД, и результаты, полученные нами при росте ЧД, в целом совпадают с их данными. Что касается изменений спектральных показателей ВСР при снижении ЧД, таких данных в литературе нет.

Несмотря на имеющиеся работы о связи спектральных показателей ВСР с ЧД при метрономизированном дыхании [1-8], они не трактовались как захват контура регуляции.

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

ВЫВОДЫ

1. При использовании адекватных программ (паттернов) изменения частоты метрономизированного дыхания можно добиться захвата и удержания контура регуляции, что подтверждается соответствующими изменениями спектральных

ЛИТЕРАТУРА

показателей ВСР.

2. Захват и удержание контура регуляции можно получить при увеличении и уменьшении частоты метрономизированного дыхания.

3. При изменениях частоты метрономизированного дыхания происходит захват и удержание спектральных показателей ВСР, но не ЧСС.

4. Метрономизированное дыхание у здоровых добровольцев на частотах 6-8 дыханий в минуту приводит к усилению ЬБ спектра ВСР, на частотах 10-12 дыханий в минуту - НБ спектра ВСР. Метрономизированное дыхание на частоте 9 дыханий в минуту приводит к активации обеих областей спектра ВСР.

5. Для получения захвата и удержания контура регуляции требуется использование достаточной длительности каждого шага, величины приращения частоты метрономизированного дыхания и емкости буфера, например, длительности шага и емкость буфера в 1 минуту с, приращением частоты дыхания в 1 дыхание в минуту.

6. При использовании «плавного» изменения частоты дыхания с шагом 0,1 дыхание в минуту и длительностью каждого шага - 6 секунд захвата и удержания контура регуляции не происходит.

7. Задачи биообратной связи, основанные на технологии метрономизированного дыхания и ВСР, должны решаться при использовании достаточных параметров изменения частоты дыхания.

1. Kobayashi H. Does paced breathing improve the reproducibility of heart rate variability measurements? /

H. Kobayashi // J. Physiol Anthropol. - 2009. - № 28(5). - P. 225-230. '

2. Ng J. Autonomic effects on the spectral analysis of heart rate variability after exercise. / J.Ng, S. Sundaram, A.H. Kadish [et al.] // Am J. Physiol Heart Circ Physiol. - 2009. - № 297(4). - P. 1421-1428.

3. Fang Y. Effect of different breathing patterns on nonlinearity of heart rate variability. / Y. Fang, J.T. Sun, C. Li [et al.] //Conf Proc IEEE Eng Med Biol Soc. - 2008. - P. 3220-3223.

4. Shields R.W. Jr. Heart rate variability with deep breathing as a clinical test of cardiovagal function. / R.W. Jr. Shields // Cleve Clin J Med. - 2009,- № 76. - Suppl 2. - P. 37-40.

5. Perakakis P. Breathing frequency bias in fractal analysis of heart rate variability. / P. Perakakis, M.

Taylor, E. Martinez-Nieto [et al.] '// Biol Psychol. - 2009. - № 82(1). - P. 82-88.

6. Van de Louw A. Breathing cardiovascular variability and baroreflex in mechanically ventilated patients. /

A. Van de Louw, C. Medigue, Y. Papelier [et al.] //Am J Physiol Regul Integr Comp Physiol. - 2008. -№295(6).-P. 1934-1940. ' '

7. Lehrer P. M. Resonant frequency biofeedback training to increase cardiac variability: rationale and

manual for training. / P. M. Lehrer, E. Vaschillo, B. Vaschillo // Appl Psychophysiol Biofeedback. -2000.-№25.-P. 177-191. ' '

8. Buchner T. On the nature of heart rate variability in a breathing normal subject: a stochastic process

analysis. / T. Buchner, M. Petelczyc, J.J. Zebrowski [et al.] // Chaos. - 2009. - № 19(2). - P. 028504-1028504-5. '

© Кулик О.Л., Задеріхін O.K., Шульгін В.1., Мартыненко О.В., Яблучансъкий М.1., 2010