[an error occurred while processing this directive] [an error occurred while processing this directive] This document last modified Monday, 17-Jun-2019 01:50:10 MSK [an error occurred while processing this directive]

УДК 577.3

Биофизика, том 46, вып. 2, стр. 319-329, 2001

Количественный анализ вариабельности электрокардиограмм, типичных для полиморфных аритмий

А.В. Москаленко*), Н.И. Кукушкин*), Ч.Ф. Стармер***), А.А. Деев,*) К.Н. Кукушкина**), А.Б. Медвинский*)

*)Институт теоретической и экспериментальной биофизики РАН,
Пущино, Россия

**)Институт биофизики клетки РАН,
Пущино, Россия

***)Medical University of South Carolina,
171 Ashley Ave., Charleston, SC 29425, USA

В работе развит новый подход к анализу вариабельности электрокардиограмм (ЭКГ), типичных для полиморфных аритмий. Для таких ЭКГ часто невозможно идентифицировать отдельные QRS-комплексы, что делает неприменимыми развитые ранее методы количественного анализа ЭКГ. Предлагаемый в данной работе подход базируется на численной оценке вариабельности соседних участков ЭКГ. При этом отпадает необходимость идентификации отдельных QRS-комплексов. На основе этого подхода развивается анализ нормированной изменчивости ЭКГ (АНИ-метод), в рамках которого каждому участку и/или всей ЭКГ в целом сопоставляются два индекса, характеризующие вариабельность осцилляций и её изменения во времени. Эти индекс позволяют как численно оценивать полиморфизм отдельной ЭКГ, так и сравнивать между собой различные электрокардиограммы. Предлагаемый подход к анализу вариабельности ЭКГ может оказаться полезным для исследований механизмов и диагностики полиморфных аритмий.

Ключевые слова: полиморфные аритмии; ЭКГ; количественный анализ.


Скачать:

Электрокардиография в настоящее время является наиболее широко используемым на практике методом дифференциальной диагностики аритмий сердца. Более того, многие из аритмий могут быть обнаружены лишь путем регистрации электрокардиограмм (ЭКГ). Есть все основания считать, что в дальнейшем роль электрокардиографии в дифференциальной диагностике нарушений сердечного ритма будет возрастать. Вместе с тем остается необходимость постоянного совершенствования методов дифференциальной диагностики сердечных аритмий с целью увеличения их чувствительности и специфичности, что требует развития новых алгоритмов анализа ЭКГ и применения более точных диагностических критериев.

В качестве параметров, позволяющих характеризовать ЭКГ с хорошо выраженными QRS-комплексами, кардиологи традиционно используют амплитуду и длительность зубцов QRS-комплекса, длительность внутрицикловых и межцикловых сегментов (например, RR-интервалов) и изменение их величины вдоль записи (отметим, что при этом важное значение придается форме как отдельных зубцов, так и QRS-комплекса в целом, а также - феномену отсутствия каких-либо зубцов, либо - появлению дополнительных). Все эти особенности ЭКГ отражают особенности проведения электрического возбуждения по миокарду.

В последние десятилетия при анализе ЭКГ все чаще применяются математические методы. В частности, используются классические методы анализа временных рядов, например, вычисление статистических параметров последовательности RR-интервалов или спектрального распределения энергии [1, 2]. Развитие теории хаоса ввело набор совершенно новых средств анализа ЭКГ: сечение Пуанкаре, реконструкция фазового портрета, вычисление показателя Ляпунова и энтропии Колмогорова и др. [3,4,5,6]. В качестве примера применения математических методов анализа ЭКГ можно указать работу [7], посвященную использованию так называемого DFA-метода для анализа последовательностей RR-интервалов при длительной (многочасовой) регистрации ЭКГ.

Особую проблему представляет анализ полиморфных желудочковых тахикардий (ПЖТ) – высокочастотных аритмий, характеризующихся высокой вариабельностью ЭКГ. Трудность решения этой проблемы связана с тем, что на ЭКГ ПЖТ практически невозможно различить отдельные QRS-комплексы, что делает непригодными традиционные способы анализа ЭКГ. ПЖТ часто отождествляют с желудочковой фибрилляцией, предшественницей внезапной сердечной смерти. Классическое определение фибрилляции состоит в следующем: “Фибрилляция желудочков сердца - патологическое (хаотическое, разрозненное и разновременное) сокращение отдельных мышечных волокон желудочков сердца, неспособное поддерживать их эффективную работу. (…) ЭКГ в этих случаях характеризуется стадийностью, отражающей степень гипоксии и метаболических изменений в миокарде. В характеристике каждой стадии имеет значение вид и соотношение фибриллярных осцилляций - ритмичных изоморфных и аритмичных полиморфных.” [10]. Отсюда ясно, что проблема оценки вариабельности ПЖТ представляет существенный интерес и существует объективная потребность в развитии новых методов количественного анализа ЭКГ. Данная работа посвящена поиску путей решения этой проблемы.

Ранее в работе [11], мы, используя рекомендации, предложенные в [12], предприняли попытку ввести количественные критерии различий полиморфных, мономорфных и квазимономорфных ЭКГ в зависимости от изменения амплитуды и частоты регистрируемых сигналов. При этом вариабельность полиморфных аритмий численно не оценивалась.

В данной работе основное внимание уделяется поиску параметров вариабельности ЭКГ, способных дать количественную оценку полиморфности ЭКГ при ПЖТ. Для решения этой задачи нами развит новый подход к анализу аритмий, базирующийся на численной оценке вариабельности соседних участков ЭКГ (при этом отпадает необходимость идентификации отдельных QRS-комплексов). Каждому участку и/или всей ЭКГ мы ставим в соответствие два индекса, которые характеризуют вариабельность осцилляций и её изменения во времени. Изменения эти индексов позволяют численно как оценивать полиморфизм отдельной ЭКГ, так и сравнивать между собой различные электрокардиограммы. Предлагаемый подход к анализу вариабельности ЭКГ может оказаться полезным для исследований механизмов, а также для диагностики полиморфных аритмий.

АНАЛИЗ НОРМИРОВАННОЙ ИЗМЕНЧИВОСТИ ЭКГ ПРИ ПЖТ

1. Критерии, которым должны удовлетворять методы анализа нормированной изменчивости ЭКГ при ПЖТ. Мы полагаем, что характеристика вариабельности (изменчивости) ЭКГ при ПЖТ должна:
1) представлять собой число или набор чисел, которые дают возможность сравнивать два участка ЭКГ или две ЭКГ произвольной длительности;
2) быть нормированной таким образом, чтобы не зависеть от конкретной формы и амплитуды QRS-комплексов (на тех участках, где эти комплексы различимы);
3) позволять отслеживать динамику полиморфизма в ходе аритмии даже в тех случаях, когда не удаётся различать отдельные QRS-комплексы;
4) давать возможность получения новой информации о характере ЭКГ при желудочковых тахиаритмиях, не противоречащую результатам, получаемым с применением традиционных качественных методов анализа полиморфных ЭКГ.

Семейство методов, удовлетворяющее указанным выше критериям, будем в дальнейшем называть анализом нормированной изменчивости (АНИ) ЭКГ при ПЖТ. В данной работе представлена одна из возможных реализаций такого анализа, или АНИ-метод.

2. Экспертная оценка вариабельности простых ЭКГ. В качестве шага, предшествующего разработке АНИ-метода, и с целью тестирования этого метода мы формулируем правила, позволяющие оценивать вариабельность ЭКГ в тех простых случаях, когда удаётся различать отдельные QRS-комплексы. Этот шаг необходим для проверки работы развиваемого нами подхода для тех случаев, когда, характеризуя полиморфизм ЭКГ, удаётся сравнить амплитуду и ширину соседних QRS-комплексов и получить количественную оценку. Для таких простых случаев легко затем сопоставить полученные подобным образом оценки вариабельности ЭКГ с теми оценками, которые получаются в результате применения АНИ-метода. Качественное совпадение оценок будет свидетельствовать о выполнении критерия 4 (см выше). Далее мы предполагаем, что метод, адекватно (т.е. в соответствии со сформулированными правилами оценки вариабельности простых случаев ПЖТ) классифицирующий простые формы полиморфной ЭКГ, даёт приемлемые результаты и в тех случаях, когда сравнение невозможно, а именно, когда отдельные QRS-комплексы практически неразличимы.

Для экспертной оценки вариабельности таких ЭКГ, в которых удаётся идентифицировать отдельные QRS-комплексы, мы полагаем, что вариабельность, характеризующая переход от i-го комплекса к последующему, (i+1)-ому, равна

(1)

где Ai, Ai+1 - размах осцилляций в двух соседних комплексах, а Ti и Ti+1 - ширина этих комплексов; wA = 1 при Ai = Ai+1 и wT = 1 при Ti = Ti+1. Таким образом, вариабельность комплексов, размах осцилляций и ширина которых одинаковы, равна нулю, т.е. при Ai = Ai+1 и Ti = Ti+1 wi = 0.

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

(2)

где N - число переходов.

С целью облегчения процедуры проверки работы АНИ-метода для случаев простых ПЖТ нами использовались так называемые “искусственные ЭКГ”, т.е. кривые с заранее заданными значениями V, которые были сконструированы из элементов псевдо-ЭКГ, полученных нами на основе обработки результатов экспериментов по многоэлектродному картированию распространения возбуждения миокарда во время искусственно вызванных желудочковых аритмий (метод подробно описан в [11]).

На рис. 1-I показаны два типа элементов, наиболее характерных для псевдо-ЭКГ при ПЖТ, которые использовались в качестве базовых для построения искусственных ЭКГ: двухкомпонентный (а) и трёхкомпонентный (b). Ширина двухкомпонентного элемента Ta = 71 мс , а ширина трёхкомпонентного элемента Tb = 182 мс, так что wT = 2.56 и wA = 1,56.

На рис. 1-II показан пример искусственной ЭКГ, состоящей из двух участков, внутри каждого из которых переход от одного комплекса к другому характеризуется величиной
wi = 0, в то время как для перехода от одного участка к другому, т.е. от элемента b к элементу a, w5 = 3. Следовательно, в соответствии с формулой (2), вариабельность ЭКГ в целом, V, в данном случае равна 0,16.

На рис. 1-III представлен пример более сложной искусственной ЭКГ. Она получена из элементов a и b при вариации амплитуды одного из этих базовых комплексов. В данном случае, в соответствии с (2), V = 2,8.

На рис. 2-I и 3-I представлены искусственные ЭКГ, отобранные для демонстрации АНИ-метода и последующего сравнения результатов, полученных с его помощью, с соответствующими значениями V.

3. Определение мгновенной вариабельности ЭКГ с использованием АНИ-метода. С целью определения локальных характеристик вариабельности ЭКГ мы предлагаем сравнивать некоторый оцифрованный участка ЭКГ с другим оцифрованным участком, который рассматривается в качестве образца для сравнения. Его длительность (окно образца, ОО) мы выбрали равной минимально возможной длительности QRS-комплекса, характерной для данного биологического вида в отсутствие аритмии. Процедура сравнения осуществляется для каждого момента времени. Получаемая в результате сравнения локальная характеристика вариабельности ЭКГ (индекс мгновенной вариабельности, Ii) позволяет численно оценивать различие между i-ым участком ЭКГ, соответствующим текущему положению ОО, и последующим ближайшим во времени j-ым участком ЭКГ, аналогичным образцу.

Формально процедура сравнения реализована следующим образом. Каждому из участков ЭКГ соответствует вектор

Fk = (fk, fk+h, ... ,fk+(p-1)h),

где fn - величина сигнала в момент n, n изменяется в интервале от k до k+(p-1)h; p - размерность так называемого пространства вложения, а шаг h называется лагом вложения (см., например, [4,5]). Расстояние между двумя векторами определяется нормой

rij = |Fi - Fj|. (3)

Отметим, что для периодического сигнала rij = 0, когда |i-j| = mT, m = 0,1,2,3,... Для апериодического сигнала rij тем меньше, чем меньше отличаются друг от друга i-й и j-й участки кривой.

Ближайший во времени аналогичный образцу j-ый участок ЭКГ (т.е. такой участок, который менее других отличается от образца) находится на некотором временном расстоянии t i = j-i, j>i от соответствующего ему i-го участка-образца. Процедура поиска j-го участка и определения t i предполагает, что j однозначно зависит от i.

Если величина ti каким-то образом определена, то для оценки подобия двух участков можно воспользоваться нормой (3) в p-мерном пространстве вложения. Таким образом, мы постулируем:

(4)

где момент i определяет локализацию участка-образца, а Si - размах амплитуды ЭКГ на этом участке.

Для автоматизированного поиска участка ЭКГ, аналогичного образцу и ближайшего к нему во времени, мы последовательно с шагом h сканируем ЭКГ в некотором окне шириной L (окне сканирования, ОС), сопоставляя образец с каждым последующим участком внутри ОС. Сопоставление может проводиться с помощью различных методов, например, с применением автокорреляционной функции или нормы в пространстве вложения. В простых случаях, когда QRS-комплексы различимы, значение t i, независимо от применяемого метода, равно межкомплексному интервалу. Для сложных случаев, связанных с анализом реальных высокочастотных полиморфных аритмий, мы формально приравниваем t i положению минимума нормы (3); иными словами, мы принимаем, что минимальное значение rij в ОС определяет локализацию j-го участка, который менее других отличается от образца (аналогичен ему). При этом формула (4) приобретает вид:

(5)

Такой подход позволяет значительно сократить затраты машинного времени на вычисление Ii, поскольку объединяются процедуры поиска участка-аналога и сравнения исследуемого участка ЭКГ с участком-образцом. Ширина ОС L выбрана равной максимально возможной (для данного биологического вида) длительности QRS-комплекса.

Необходимо отметить то обстоятельство, что поиск участков-аналогов является самостоятельной процедурой, в общем случае совершенно независимой от процедуры сравнения исследуемого участка с участком-образцом: поиск может осуществляться на основе анализа одного набора признаков, необязательно связанных с полиморфными свойствами ЭКГ, в то время как определение вариабельности основывается на наборе признаков, непосредственно связанным с полиморфизмом.

4. Исследование мгновенной вариабельности искусственных ЭКГ. Результаты вычисления функции Ii с применением АНИ-метода к искусственным ЭКГ представлены на рис. 2 (для простых искусственных ЭКГ с единственным переходом, для которого значение wi отлично от нуля) и на рис. 3 ( для более сложных полиморфных искусственных ЭКГ). Во всех представленных на рис. 2 и 3 случаях ширина базовых комплексов, использованных для построения ЭКГ, та же, что на рис. 1-I; меняется только амплитуда комплекса b.

На рис. 3-II можно видеть тенденцию роста средних по времени значений Ii с возрастанием значений V (табл. 1), характеризующих вариабельность соответствующих кривых (рис. 3-I).

Отметим, что на участках, соответствующих переходам, которые характеризуются одинаковыми значениями вариабельности wi обнаруживаются некоторые различия вида функции Ii, хотя можно было бы ожидать, что функция Ii. будет иметь одинаковый вид для одинаковых wi. Так на рис. 2-II представлены графики функции Ii для соответствующих простых искусственных ЭКГ (рис. 2-I); заметим, что переход в паре a и b характеризуется изменением периода колебаний, в то время как в случаях c и d период не изменяется. Видно, что для каждой из симметричных пар ЭКГ (a-b и c-d соответственно) вид функции Ii не сохраняется неизменным, несмотря на то, что соответствующие значения вариабельности перехода wi (1) для этих ЭКГ одинаковы. На наш взгляд, такое различие обусловлено, с одной стороны, не скомпенсированной ошибкой распознавания положения участков-аналогов.

Такие ошибки распознавания можно обнаружить по виду функции t i. Это иллюстрируется рис. 2-III. Так, например, для случая d межкомплексный интервал не меняется со временем, что должно было бы обусловливать неизменность величины t i. Видно, что вблизи перехода от меньших амплитуд к большим в данном случае наблюдается резкое изменение t i, что свидетельствует об ошибке в определении локализации участка-аналога. Такая ошибка возникает при очень резком различии характеристик двух соседних комплексов (ср. со случаем a, где различие амплитуд существенно меньше). Столь резкие изменения характеристик при переходе от одного комплекса к другому сравнительно редко встречаются в экспериментальных и клинических ЭКГ.

Предложенный нами алгоритм поиска участка, аналогичного образцу, может приводить к ошибкам двух видов: 1) в случае, когда участок-аналог оказывается вне ОС, для вычисления Ii берется значение нормы на одном из концов ОС, и величина Ii оказывается завышенной; 2) в случае, когда в ОС попадают несколько аналогичных образцу участков, причем ближайший из них отличается от образца более, чем какой-либо из последующих, Ii оказывается заниженным. Как показывает наша практика, эти ошибки частично взаимно компенсируются.

Поскольку полиморфизм ЭКГ принято оценивать по достаточно продолжительному участку ЭКГ, то АНИ-метод оценки вариабельности ЭКГ не должен существенно зависеть от локальных характеристик сигнала и ошибок, связанных с их определением. Поэтому, чтобы снизить влияние указанных недостатков примененного алгоритма, мы вводим усреднённые во времени характеристики вариабельности ЭКГ.

5. Получение усреднённых характеристик вариабельности ЭКГ с использованием АНИ-метода. Для дальнейших исследований участков ЭКГ произвольной длины и с целью получения интегральных характеристик вариабельности всего исследуемого участка ЭКГ нами взяты за основу: (a) среднее значение функции Ii (в дальнейшем - индекс вариабельности ЭКГ, V1) и (b) соответствующее ему значение коэффициента изменчивости [13] функции Ii (в дальнейшем - индекс изменчивости вариабельности ЭКГ, V2). Из табл. 1 видна тенденция роста V1 с переходом к более вариабельным кривым, т.е. к таким искусственным ЭКГ, для которых характерен рост величины V. Это свидетельствует о чувствительности предложенного метода к вариабельности сигнала. Напомним, что вычисление V1, в отличие от V, не привязано к наличию выраженных QRS-комплексов. Возможно, это обстоятельство обусловливает некоторое рассогласование в направлении изменения этих величин при переходе от кривой f к кривой i (табл. 1).

6. Использование индексов вариабельности для слежения за динамикой ЭКГ. С целью анализа изменений ЭКГ во времени вычисляются индексы вариабельности (V1i и V2i) в некотором временном окне фиксированной ширины (окне усреднения - ОУ). В результате, при сдвиге ОУ вдоль оси времени соответствующая траектория в пространстве индексов (V1i, V2i) отображает динамику изменений ЭКГ. При этом следует учитывать то обстоятельство, что свойства конца записи ЭКГ длиной в это окно оказываются неопределенными, в связи с чем нет смысла делать ОУ слишком широким. Ширину ОУ в соответствии с принятым определением желудочковой тахикардии [9] мы выбрали равной примерно шестикратной длительности QRS-комплекса, что для исследуемых нами образцов ЭКГ составляет 400 мс. В дальнейшем будет предполагаться, что, в отличие от V1i и V2i, индексы V1 и V2 (пункт 5) характеризуют участок ЭКГ, длина которого значительно превышает ОУ.

На рис. 4 показана сложная искусственная ЭКГ и временные зависимости величин Ii, V1i и V2i,, характеризующие эту ЭКГ. Искусственная ЭКГ (рис. 4-I) получена путём последовательной сшивки четырёх фрагментов i, h, g и f (в порядке убывания интегральной экспертной характеристики вариабельности V; см. табл. 1), показанных на рис. 3-I. Можно видеть, что как функция Ii (рис.4-II), так и функции V1i и V2i, (рис.4-IV и рис.4-V, соответственно) удовлетворительно отражают изменения на ЭКГ (рис.4-I).

На рис. 5 в пространстве индексов (V1i, V2i) представлены траектории (рис.4-IV и рис.4-V), отражающие динамику вариабельности ЭКГ (изображённой на рис. 4-I) в пределах каждого из участков i, h, g и f. Отметим, что каждый из участков компактно проецируется в существенно различные области пространства индексов. Важно указать, что эти траектории существенно различаются даже в тех случаях, когда соответствующие значения индексов V1 и V2 изменяются не очень сильно (ср., например, участки i и g на рис. 5) и по этому признаку различие вариабельности двух участков не может быть установлено. Другими словами, построение траекторий в пространстве индексов дает дополнительную информацию относительно вариабельности ЭКГ.

Подчеркнём, что использование АНИ-метода для слежения за динамикой искусственных ЭКГ (рис. 5) не связано с распознаванием отдельных пиков; это даёт основание предполагать, что анализ изменений каждого из индексов: V1i и V2i, - может оказаться полезным и для реальных полиморфных ЭКГ, даже тех, где QRS-комплексы не поддаются идентификации. В данной работе в качестве первого шага к такому анализу предлагаемый нами метод используется для исследования вариабельности псевдо-ЭКГ, полученных нами экспериментально.

7. Исследования вариабельности экспериментальных псевдо-ЭКГ. Исследуемые в данной работе псевдо-ЭКГ представляют собой временные ряды, полученные путём взвешенного суммирования отдельных электрограмм, зарегистрированных в наших экспериментах в ходе двустороннего многоэлектродного картирования процессов распространения возбуждения во время искусственно вызванных аритмий в изолированном препарате - стенке правого желудочка сердца суслика. Используемая экспериментальная модель и метод получения псевдо-ЭКГ подробно изложены в работе [11]. В данном разделе представлены результаты применения АНИ-метода для анализа вариабельности псевдо-ЭКГ, характерных для полиморфных аритмий.

На рис. 6 показаны три примера подобных псевдо-ЭКГ и соответствующие им траектории в пространстве (V1i, V2i). Заметно, что при видимом увеличении вариабельности псевдо-ЭКГ (от рис. 6а и к рис. 6b) траектория попадает в область больших значений индексов V1i и V2i. Индексы V1 и V2, характеризующие вариабельность ЭКГ в целом (на рисунке не показаны), при этом также изменяются, причём аналогичным образом. Из рис. 6с видно, что переход от полиморфного участка ЭКГ (сегмент 1) к сравнительно более регулярным осцилляциям (сегмент 2) сопровождается переходом в пространстве индексов (V1i, V2i) из одной области этого пространства в другую, так что полиморфному участку соответствуют в среднем большие значения каждого из этих индексов.

На рис. 7 представлены два примера псевдо-ЭКГ, которые зрительно практически неотличимы друг от друга по вариабельности сигнала во времени. Однако, изменения в характере и локализации соответствующих траекторий (а значит, и в величине индексов V1 и V2) в пространстве (V1i, V2i) позволяют сделать вывод о заметном различии в характере этих аритмий: кривая a более вариабельна, чем кривая b.

ОБСУЖДЕНИЕ

Нами предложен метод, который позволяет проводить численный анализ вариабельности ЭКГ.

Другим существенным результатом нашего исследования является, что, кроме вариабельности V1i, ЭКГ характеризуется индексом V2i, описывающим изменчивость самой вариабельности ЭКГ. Это позволяет различать аритмии с одинаковой вариабельностью, но с различным характером ее изменчивости. Таким образом появляется возможность различать “более устойчивый” и “менее устойчивый” полиморфизм. Отметим вместе с тем, что изменения вариабельности ЭКГ, оцениваемые по изменению индекса V1, обычно неплохо коррелируют с изменениями V2. Исследование природы этого факта может открыть новый взгляд на механизм полиморфных аритмий. Важным представляется также исследование корреляции изменений характера картин возбуждения миокарда и изменений в характере траекторий в пространстве (V1i, V2i).

Вместе с тем следует отметить, что предложенная в данной работе реализация (5) более универсального алгоритма (4) анализа вариабельности полиморфных аритмий не свободна от недостатков. Так, например, выбранный в ходе применения нашего метода участок-аналог может иногда не соответствовать таковому, определяемому экспертом визуально. В таких случаях (промахах распознавания) получаемые значения Ii оказываются или заниженными, или завышенными (см. пункт 4). Отметим ещё два типичных артефакта. Во-первых, это - “феномен начала”, т.е. резкое падение значения Ii почти до нуля на отрезке, предшествующем первым осцилляциям. Во-вторых, это - “феномен окончания”, т.е. резкое увеличение Ii и V1i при сравнительно небольшом колебании V2i, которые в сочетании прорисовывают в пространстве (V1i, V2i) характерную “улыбку”, не отражающую действительных изменений вариабельности ЭКГ (рис. 8). Оба артефакта возникают как результат выбранной нами нормировки на размах колебаний (5). Разработка более совершенных модификаций предложенного нами подхода к анализу полиморфизма может явиться предметом будущих исследований.

АНИ-метод и его модификации могут явиться полезным инструментом для классификации аритмий, поскольку вместо расплывчатого, качественного понятия “полиморфизм” появляется количественная оценка изменчивости ЭКГ как для каждого момента времени, так и для участков ЭКГ произвольной длительности.

Работа поддержана грантом Медицинского Университета Южной Каролины RB0-676.

СПИСОК ЛИТЕРАТУРЫ

1.The European Society of Cardiology and the North American Society of Pacing and Electrophysiology. Heart Rate Variability // Circulation. 1996. V. 93. P. 1043-1065.

2.Кнеппо П. // Теоретические основы электрокардиологии. Москва: Медицина. 1979. С. 410-432.

3.Babloyantz A., Destexhe A. // Biol. Cybern. 1988. V. 58. P.203-211.

4.Kaplan D., Glass L. // Understanding Nonlinear Dynamics. New-York: Springer, 1995.

5.Kantz H., Schreiber T. // Nonlinear tine series analysis. Cambridge: Cambridge University, 1997.

6.Skinner J.E., Carpeggiani C., Landisman C.E., Fulton K.W.// Circ. Res. 1991. V. 68. P. 966-976.

7.Absil P.-A., Sepulchre R., Bilge A., Gerard P. // Physica A. 1999. V. 272. P.235-244.

8.Robles de Medina E.O., Fisch Ch., Bernard R., Coumel Ph., Damato A.N., Krikler D., Mazur N.A., Meijler F.L., Mogensen L., Moret P., Pisa Z., Rosenbaum M.B., Wellens H.J. // Кардиология. 1981. т. 22. стр. 86-89.

9.Waldo A.L., Akhtar M., Brugada P., Henthorn R.W., Scheinman M.M., Ward D.E., Wellens H.J.J. // JACC. 1985. V. 6. P. 1174-1177.

10.Семёнов В.Н. // Большая Медицинская Энциклопедия. Т. 26. Москва: Советская энциклопедия. 1985. С. 275.

11.Кукушкин Н.И., Сидоров В.Ю., Медвинский А.Б., Ромашко Д.Н., Бурашников А.Ю., Стармер Ч.Ф., Саранча Д.Ю., Баум О.В. // Биофизика. 1998. Т. 43. С. 1043-1059.

12.Coumel Ph. // Antiarrhythmic Drugs. Berlin, Heidelberg: Springer. 1989. P. 86.

13.Войтинский Е.Я., Лившиц М.Е., Ромм Б.И., Рыжиков В.С. // Анализ биопотенциалов на цифровой адаптивной системе. Ленинград: Наука. 1972.

РИСУНКИ


Рис.1. “Искусственные” ЭКГ и оценка их вариабельности. I - два типа комплексов, наиболее характерных для ПЖТ: двухкомпонентный (а) и трёхкомпонентный (b). Эти комплексы использовались в качестве элементов для построения искусственных ЭКГ. II - пример одной из простейших искусственно сконструированных ЭКГ. Кривая состоит из двух однородных (составленных из идентичных комплексов) участков. Место состыковки этих участков характеризуется не нулевой величиной вариабельности перехода. III - пример более сложной искусственной ЭКГ, полученной из элементов a и b путём вариации амплитуды одного из этих базовых комплексов.

Вернуться к тексту

Fig.1
Скачать для принтера(33Kb)

Рис.2. Приложение АНИ-метода к простым искусственным ЭКГ . I - четыре примера (a, b, c и d) простых искусственных ЭКГ (с одним и тем же масштабом по вертикальной оси). Образцы составлены таким образом, чтобы ЭКГ в каждой из пар a-b, и c-d были симметричны относительно места состыковки однородных участков. II - функции Ii для соответствующих простых искусственных ЭКГ. Масштаб по вертикальной оси для всех функций - один и тот же. III - функции t i для тех же ЭКГ. Временной масштаб, одинаковый для всех кривых, указан в верхней части рисунка.

Вернуться к тексту

Fig.2
Скачать для принтера(21Kb)

Рис.3. Приложение АНИ-метода к сложным искусственным ЭКГ . Столбец I - образцы искусственных ЭКГ (с одним и тем же масштабом по вертикальной оси). Столбец II - соответствующие им графики функции Ii. Масштаб по вертикальной оси для всех функций - один и тот же. Временной масштаб, одинаковый для всех кривых, указан в верхней части рисунка.

Вернуться к тексту

Fig.3
Скачать для принтера(32Kb)

Рис.4. Приложение АНИ-метода к искусственной ЭКГ, составленной из кривых i, h, g и f (рис. 3) в порядке убывания интегральной экспертной характеристики вариабельности V. I - искусственная ЭКГ; II - график функции Ii; III - функция t i; IV - функция V1i; V - функция V2i. Буквами i, h, g и fотмечены соответствующие участки кривых, каждый длительностью 2000 мс.

Вернуться к тексту

Fig.4
Скачать для принтера(51Kb)

Рис.5. Траектории i, h, g и f в пространстве индексов V1i, V2i, отражающие динамику соответствующих участков искусственной ЭКГ, показанной на рис. 4. Пунктиром отмечены величины индексов V1 и V2, они же приведены над каждой из траекторий.

Вернуться к тексту

Fig.5
Скачать для принтера(18Kb)

Рис.6. Приложение АНИ-метода к экспериментальным псевдо-ЭКГ. Столбец I - псевдо-ЭКГ, столбец II - соответствующие этим псевдо-ЭКГ траектории в пространстве индексов V1i и V2i.

Вернуться к тексту

Fig.6
Скачать для принтера(141Kb)

Рис.7. Приложение АНИ-метода к экспериментальным псевдо-ЭКГ, визуально сходным по вариабельности. Столбец I - псевдо-ЭКГ, столбец II - соответствующие им траектории в пространстве индексов V1i и V2i.

Вернуться к тексту

Fig.7
Скачать для принтера(38Kb)

Рис.8. “Феномен окончания” как один из типичных артефактов АНИ-метода. Слева - конечная стадия аритмии, справа - соответствующая ей траектория в пространстве индексов V1i и V2i.

Вернуться к тексту

Fig.8
Скачать для принтера(15Kb)

 


Таблицa 1

Численные характеристики вариабельности искусственных ЭКГ (рис. 2-I и рис. 3-I)

Вернуться к тексту

код кривой экспертная оценка вариабельности
V
индекс
вариабельности
V1´103
индекс изменчивости вариабельности
V2
a 0,2 0,4 0,27
b 0,2 1,1 0,36
c 0,3 1,2 0,35
d 0,3 0,5 0,32
e 0,6 3,0 0,73
f 1,0 3,7 0,91
g 1,8 6,1 1,52
h 2,0 7,5 2,43
i 2,8 6,2 1,78

 

[an error occurred while processing this directive]