2.3 Энтропия цифровых рядов амплитудной и фазовой составляющих ЭКГ
Количество информации, недостающее для описания ярусной диаграммы в форме частного сообщения определяется известным из теории связи уравнением 1 [105]:
Рассмотрение точечной RR-интервалограммы и R-граммы в форме макросостояния «замороженного» системного ритма с набором микросостояний фазовых или ярусных значений RR-интервалов и R-амплитуд позволяет применить к ней известную трактовку энтропии для идеального газа и формулу Больцмана:
где Г - термодинамическая вероятность, определяемая суммарным количеством возможных равновероятных микропереходов.
В нашем случае величина Г определяет число возможных комбинаций значений RR-интервалов и R-амплитуд на ярусной диаграмме ритма, а р= 1/Г - вероятность одной комбинации. Тогда
и обе формулы становятся адекватными с точностью до постоянного множителя [94].
Для оценки информации, содержащейся в полном сообщении, обычно используется Шенноновский подход к расчёту информационной энтропии [105]:
где pi- вероятность исхода того или иного значения выборки, определяемая по /-тому столбцу гистограммы [105]. Используя количество ярусов k = ΔX∕Δx, получаем, что на один столбец гистограммы приходится один ярус диаграммы. Преобразовав(3) с использованием формулы Стирлинга [106], получаем, что исходная форма записи может быть представлена в виде [17]:
48
Как было показано в работе [107], при практической реализации расчёта информационной энтропии по подходам Шеннона (2) и Больцмана (3) для виртуальных рядов, сгенерированных по нормальному закону, разница получаемых оценок информационной энтропии стремится к нулю при объеме выборки n—→∞.
С ростом nколичество ярусов к растет в логарифмической зависимости. При этом величина Bрастет экспоненциально, как показано на рисунке 2.3 а, а величина абсолютной погрешности ΔI* = B/nмонотонно падает к нулевому значению при n → ∞ (рисунок 2.3 б).
Величина I* монотонно растет с ростом п, но при этом скорость роста также монотонно падает [94]. Результаты анализа коротких записей ЭКГ условноздоровых людей в состоянии покоя (σ = 40 - 70 мс) показывают, что расчетная относительная погрешность ΔI*∕I* составляет 3 - 5 % при п= 500 и не превышает 2 - 3 % для двадцатиминутной записи при п = 1300 - 2000. Полученные оценки погрешности включают в себя также отклонения функции экспериментального распределения от нормального закона. Для оценки информационной
энтропии ритма сердца в интервале Δσ величиной абсолютной погрешности ΔI*можно пренебречь при n > 2000 [61]. С учетом принятой точности
На рисунке 2.4 приведены данные точечной диаграммы зависимости выборочных расчетных значений Iςот объема nвыборок данных 366 регистраций ЭКГ больных и здоровых людей. Каждая точка на диаграмме адекватно соответствует одной из 366 ритмограмм [94]. При сопоставлении данных распределения I∑(π)с данными анамнеза обнаружено [100], что общее функциональное состояние организма человека адекватно величине тангенса угла наклона ΔI√Δπ.При этом значение 8,3 бит на единицу отсчета характерно для здоровых молодых людей и в рамках проведенных исследований является максимально возможным. Минимальное значение 1,2 бит на единицу отсчета отмечались в работах [94] для тяжелых больных в состоянии комы.
Рисунок 2.4 - Точечная диаграмма зависимости выборочных расчетных значений Iςот объема nвыборок данных по 366 регистрациям ЭКГ больных и здоровых людей [94]
Точечные данные распределения I∑(n) и по указанному признаку позволяют оценивать текущий уровень функционального состояния любого здорово
50 го или больного человека.
Распределение T∑(n) приобретает качество «диаграммы функциональных состояний», разделенной на верхнюю «запрещенную» и нижнюю «разрешенную» - функциональную область. Границей раздела является линейная функция, тангенс угла наклона которой равен ΔI∑∕Δnmax, определяющий «математическое ожидание» энтропии
системы X, состояния которой распределены по нормальному закону. Таким образом, для ритма сердца здоровых молодых людей при n → ∞ функции I*(σ) и H(X) становятся идентичными [108].
При формировании RR-интервалограммы и R-граммы количество информации набирается дискретно: от систолы к систоле, поэтому принципиально невозможно использовать понятие скорости набора информации, являющееся основным параметром для технических устройств связи. В этом случае масштабной единицей становится переменный интервал времени события (RR-интервал и амплитуда R). RR-интервалограммы и R- граммы представляют собой номерные ряды последовательностей n таких событий.
Величина I*, определенная с точностью до величины ABmax/n, определяет среднее количество информации, недостающее до полного описания одного отсчета [94,109].
Информационная энтропия I* обладает основными свойствами физической энтропии — при фиксированных внешних условиях растет с ростом n, принимая максимальное совместимое с внешними условиями значение. Также она определяется отношением количества информации i∑, недостающего для полного описания диаграммы к объему выборки п, и по смыслу является средним количеством информации, недостающим для описания одного микроперехода на диаграмме [98]. По сравнению с другими параметрами вариабельности сердечного ритма, информационная энтропия I* не теряет адекватного физического смысла для многомодального распределения
и для RR-интервалограммы имеет постоянную, четко выраженную «правую границу условного здоровья» [108]. Очевидно, что разные ФСО могут дать один и тот же результат по величине I* для RR-интервалограммы и R-граммы [110], что требует разработки системных методов, уточняющих оценку функционального состояния ССС по информационным характеристикам интервальных рядов.
Применение энтропии в анализе ритма сердца в форме RR- интервалограмм и R-грамм выгодно отличается от других параметров исследования тем, что существует граничное значение энтропии при нормальном и равномерном законе распределения при n→∞. Поэтому возникает возможность сравнивать не только энтропийные показатели с эталонным значением, но сами формы распределений [111-113].
2.4
Еще по теме 2.3 Энтропия цифровых рядов амплитудной и фазовой составляющих ЭКГ:
- 2.2 Особенности представления амплитудной и фазовой составляющих ЭКГ-сигнала
- Физическая модель амплитудно-фазового сопряжения ЭКГ- составляющих у условно-здоровых молодых людей
- Разработка метода оценки амплитудно-фазового сопряжения ЭКГ-составляющих на основе информационного подхода
- Разработка метода оценки амплитудно-фазового сопряжения ЭКГ-составляющих на основе статистического подхода
- Оценка корреляции параметров амплитудной и фазовой составляющих ЭКГ-сигнала
- Исследование функций распределения параметров амплитудной и фазовой составляющих ЭКГ-сигнала в ансамбле регистраций экспериментальной базы данных
- Исследование сопряжения амплитудной и фазовой составляющих ЭКГ-сигнала в форме циркуляционной кривой
- Экспериментальное программно-алгоритмическое обеспечение для анализа параметров амплитудно-фазового сопряжения ЭКГ- составляющих и управления диагностическим процессом
- 4.3. Исследование сопряжения амплитудной и фазовой составляющих ЭКГ-сигнала в форме двупараметрической гистограммы зависимости длительности RR-интервала от амплитуды зубца R
- Исследование сопряжения амплитудной и фазовой составляющих ЭКГ-сигнала в форме параметрической диаграммы сопряжения