CC BY
36
УДК 519.876.5
DOI: 10.17586/0021-3454-2016-59-11-959-963
ПРОГРАММНЫЙ КОМПЛЕКС ИДЕНТИФИКАЦИИ ХАОТИЧЕСКИХ ПАРАМЕТРОВ МОЗГОВОГО КРОВООБРАЩЕНИЯ
А. А. Мусаев1'2, А. И. Загайнов3
1 Санкт-Петербургский государственный технологический институт (Технический университет),
190013, Санкт-Петербург, Россия E-mail: amusaev@technolog.edu.ru 2Санкт-Петербургский институт информатики и автоматизации РАН, 199178, Санкт-Петербург, Россия 3 Военно-космическая академия им. А. Ф. Можайского, 197198, Санкт-Петербург, Россия
Рассматриваются возможности разработанной автоматизированной системы, предназначенной для анализа параметров мозгового кровообращения. Изложены методы численной реализации фрактальной размерности и намечены пути ее вычисления в режиме реального времени. Показан интерфейс реализованного в среде Qt Creator программного обеспечения. Приведен пример построения фрактальных размерностей мозгового кровообращения пациента, рассчитанный в реальном времени с помощью доплерографа Multi-DopX.
Ключевые слова: восстановленный аттрактор, фазовое пространство, корреляционная размерность, линейная скорость кровотока, системное артериальное давление
Введение. В настоящее время развитие методов анализа биомедицинских сигналов имеет явную тенденцию, направленную на исследование возможностей нелинейных подходов к изучению их закономерностей. Особенно отчетливо эта тенденция наблюдается при построении нелинейных базисов для разложения исходного сигнала, учете нелинейных особенностей при разработке методов нелинейного регрессионного анализа и RS-статистики, а также при исследовании сигнала в конечномерных фазовых пространствах. Все чаще встречаются прикладные задачи, связанные с использованием таких категорий, как фрактал, хаос и, особенно, детерминированный хаос. Интерес к анализу фракталов основан на постулате о существовании конечномерной системы нелинейных уравнений для исходного отображения, что расширяет границы конструктивных исследований, в частности, для прогнозирования состояния нелинейных систем.
Методология фрактального анализа. Построение фрактала из исходного одномерного конечного сигнала связано с восстановлением его аттрактора. Данное положение вытекает из теории динамических систем и теории детерминированного хаоса. Этот факт отчетливо прослеживается в известной теореме Такенса, в которой предложен способ построения восстановленного аттрактора, принадлежащего гладкому многообразию. Положив в качестве координат вектора состояния ряд наблюдений, смещенный относительно себя на некоторое постоянное значение, получим
x (i ) = (a (i), a (i + x),..., a (i + i(n -1))) = (xb *2,..., xn ),
где a(i) — исходный временной ряд, т — временная задержка, вектор xi = (xb x2,..., xn) — координаты одной точки на восстановленном аттракторе, n — размерность пространства вложения, удовлетворяющая условиям теоремы Такенса: n > 2 [d^ ] +1, де dA — размерность восстановленного аттрактора.
Свойства построенного таким образом аттрактора метрически (и вероятностно) эквиваленты исходному аттрактору динамической системы. Однако для реальных временных рядов
эта теорема неприменима в силу их конечности. К сожалению, к настоящему моменту для конечных временных рядов подобные теоретические результаты не получены. Существуют лишь оценки длины временного ряда, необходимые для характеристики степени подобия аттрактора с конечным числом значений.
Другой параметр, требующий уточнения при численной реализации, — это параметр т. Самый приемлемый способ, рекомендуемый в литературе [1, 2], заключается в нахождении
1 N—т
первого нуля автокорреляционной функции: B(т) =-V (ak — a)(ak+т — a), где
N — т k=i
ak = a(kAt) — исходный временной ряд.
Применение этого метода связано с гипотезой некоррелированности координат точек аттрактора в силу ортогональности векторов базиса. Однако выбор параметра задержки таким способом не всегда является оптимальным. Альтернатива ему — построение функции средней взаимной информации по методу, предложенному в работах [2—4]. Интервал
[min ak,max ak ] делится на L равных частей. Как правило [2], значение L выбирается по фор-
k k
муле Старка: L = [log2 N ] +1. Событие „a(t) принадлежит i-му интервалу" обозначается как Ль
а событие „a(t+T) принадлежит j-му интервалу" — как Bj; P — вероятность соответствующего события. Функция средней взаимной информации определяется как
" P(ABj) " _ P( Л) P( Bj)
и в качестве оптимального параметра задержки выбирается первый локальный минимум построенной функции:
1 (т) = —VSPЩ )log2
i=1 j =1
тор1 = ш1п{1'(т,-) = 0, I'(тГ ) < 0, I'(тг+ ) > 0}.
т,
Функция средней взаимной информации является более точной мерой независимости точек аттрактора [2]; для некоторых тестовых данных (аттрактор системы уравнений Лоренца) получение значения оптимальной задержки этим способом более предпочтительно.
Фрактальные размерности и их численная реализация. Наиболее известными характеристиками аттрактора динамической системы являются вероятностные (фрактальные) размерности. Под вероятностью здесь понимается вероятность нахождения точки в определенной области самого аттрактора в фазовом пространстве. Вариант фрактальной размерности — размерность Реньи [5]:
( М ( 8) Л
1п I РЦ
Ва = —Нш-^-1, (1)
4 1 Г ц 8—>0 1п(1 / 8)
где М(е) — минимальное количество разбиений на кубические фрагменты со стороной 8, образующие область, полностью покрывающую аттрактор; р, — вероятность пересечения ,-м кубическим фрагментом фазовой траектории динамической системы; ц=1-р,.
Частными случаями размерности Реньи являются размерность Колмогорова, информационная размерность и корреляционная размерность, получаемые при ц, равном 0, 1, 2 соответственно. Для пространственно однородных аттракторов все эти размерности одинаковы. В общем случае, исходя из определения, размерность Реньи является монотонно убывающей функцией ц: ц < Ц2 ^ ^ Д2.
Следовательно, для рассматриваемых аттракторов параметр ц должен быть положительным. В настоящее время этот параметр принят равным двум (корреляционная размерность является оценкой информационной размерности, для ее вычисления разработан уни-
версальный алгоритм, из которого автоматически следует оценка соответствующей корреляционной или аппроксимированной энтропии). Согласно выражению (1) корреляционная размерность определяется как
(М (8) Л (М(8) Л
ln
1
Е Рг ln Е p
Dc = — lim—-J = lim—^-J. (2)
1 - 2 e^o ln(1 / s) s^o ln s
Выражение (2) удобно представить в асимптотической форме: limln(C(r)), где
r^0 ln r
1 m f1, a>0,
C (r)= lim—- £0(r — p(x-, Xj )) — корреляционный интеграл; 0(a)=« ' ' — функция
m^® m , j=1 [0,a<0
Хевисайда; p( xt, Xj ) — функция расстояния в и-мерном пространстве.
Вследствие фрактальности исследуемого объекта можно предположить, что C (r) ~ г°с, откуда следует ln(C(r)) ~ Dc ln r, и корреляционную размерность можно оценить по наклону прямой, полученной при логарифмировании корреляционного интеграла.
Программный комплекс фрактального анализа. Вычисление корреляционной размерности в режиме реального времени осложняется, прежде всего, ограниченными возможностями автоматизированного вычисления необходимых параметров (размерности пространства вложения, параметра задержки, скейлингового диапазона и пр.). Поэтому в предлагаемом комплексе фрактального анализа предусмотрена возможность изменения и автоматизированного вычисления этих параметров с помощью оригинальных алгоритмов, использующих функции вероятности нахождения точек на аттракторе. Скользящее окно (т.е. текущий рассматриваемый массив) перемещается по временному ряду (сигналу), тем самым изменяя вычисляемую характеристику во времени. При этом изменяется аттрактор системы, его размерность и параметры. Здесь необходимо подчеркнуть, что указанные параметры могут изменяться уже при незначительном перемещении скользящего окна, что при большом количестве значений временного ряда представляет определенную трудность, так как на каждом шаге необходимо выполнять пересчет всех точек аттрактора.
Программный комплекс был реализован в среде программирования Qt Creator. Основные возможности комплекса состоят в следующем:
— представлении исходных временных рядов;
— автоматическом изменении аттрактора сигнала в режиме реального времени и его перестроении при изменении параметра задержки, а также в визуализации соответствующего двухмерного аттрактора системы;
— расчете параметров (левого и правого) скейлингового диапазона;
— вычислении корреляционной размерности на скользящем окне;
— аппроксимации динамики изменения корреляционной размерности;
— вычислении и визуализации специальных параметров, служащих, в частности для характеристики деятельности системы мозгового кровотока (спектральная плотность в диапазоне M- и B-волн, фазовый сдвиг, кросс-амплитуда и др.).
В комплексе предусмотрены возможности изменения длины скользящего окна и его перемещения, изменения скорости воспроизведения наблюдаемого процесса (при считывании из файла), принудительного изменения указанных выше параметров аттрактора.
На рисунке приведен пример мультиграфической визуализации полученных результатов. Доплерограф Multi-DopX позволяет построить временные ряды линейной скорости кровотока и системного артериального давления. С помощью программного комплекса на одном графике могут быть построены сами временные ряды, сдвиг фаз между М-волнами, амплитуда
Р,2
V г =1 J
B-волн, корреляционные размерности сигналов. Для количественного анализа сигнала могут быть установлены пять возможных шкал — скорость V, давление P, фазовый сдвиг Дф, амплитуда А и корреляционная размерность Пс. Так, для исследуемого сигнала линейная скорость кровотока слева изменялась от 1,38 до 1,64, справа — от 1,38 до 1,68; как видно из построенных графиков (CD2L и CD2R) эти показатели подобны. При этом фазовые сдвиги M-волн между линейными скоростями и давлением (PSR и PSL) изменялись (в соответствующем масштабе шкал) от -0,8 до 0 и от -0,95 до -0,65 рад, а амплитуды в диапазоне B-волн (ABL и ABR) изменялись от 520 до 1320 мс и от 390 до 830 мс. Как следует из полученных результатов, построенная аппроксимация корреляционной размерности точнее отображает внутренний характер структуры обрабатываемого временного ряда, нежели рассмотренные спектральные показатели.
£>„ А Аф Р V
1,8 1,6 1,4 1,2 1
0,8 0,6 0,4 0,2 0
c
1400 0 , 4 180 180
1200 0, 2 0 160 1 40 160 1 40
1000 800 -0,2 -0,4 -0,6 120 100 120 100
600 -0,8 80 80
400 -1 60 60
200 -1 , 2 -1 , 4 40 20 40 20
0 -1 , 6 0
0
5000 10000 15000 20000 25000 30000 35000 t, c
Заключение. Предложен программный комплекс вычисления фрактальных размерностей линейной скорости кровотока и системного артериального давления системы мозгового кровотока в режиме реального времени. При программной реализации комплекса необходимые характеристики (размерность пространства вложения, параметр задержки, скейлинговый диапазон и пр.) могут быть изменены уже в процессе вычислений. На примере средней скорости мозгового кровотока и системного артериального давления пациента с помощью доп-лерографа МиШ-ОорХ получены временные ряды указанных показателей, вычислены их фрактальные размерности и спектральные компоненты, выявлены их диапазоны и особенности изменения.
список литературы
1. Махортых С. А., Сычев В. В. Алгоритмы вычисления характеристик стохастических сигналов и их применение к анализу электрофизиологических данных // Сб. тез.: Математическая и вычислительная биология, 4-я Пущинская конф. молодых ученых. 1999.
2. Головко В. А. Нейросетевые методы обработки хаотических процессов // Науч. сессия МИФИ-2005: VII Всерос. науч.-техн. конф. „Нейроинформатика-2005". М.: МИФИ, 2005. С. 43—91.
3. Меклер А. А. Применение аппарата нелинейного анализа динамических систем для обработки сигналов ЭЭГ // Актуальные проблемы современной математики: Ученые записки. 2004. Т. 13(2). С. 112—140.
4. Янсон Н. Б., Анищенко В. С. Моделирование динамических систем по экспериментальным данным // Изв. вузов. Прикладная нелинейная динамика. 1995. Т. 3, № 3. С. 112—121.
5. Перерва Л. М., Юдин В. В. Фрактальное моделирование: Учеб. пособие / Под общ. ред. В. Н. Гряника. Владивосток: Изд-во ВГУЭС, 2007.
Сведения об авторах
Александр Азерович Мусаев — д-р техн. наук, профессор; СПбГТИ(ТУ), кафедра системного анализа и
информационных технологий, зав. кафедрой; СПИИРАН, лаборатория информационных технологий в системном анализе и моделировании; E-mail: amusaev@technolog.edu.ru Артем Игоревич Загайнов — ВКА им. А. Ф. Можайского, кафедра математического и программно-
го обеспечения; ст. преподаватель; E-mail: zagainov239@gmail.com
Рекомендована СПИИРАН Поступила в редакцию
01.06.16 г.
Ссылка для цитирования: Мусаев А. А., Загайнов А. И. Программный комплекс идентификации хаотических параметров мозгового кровообращения // Изв. вузов. Приборостроение. 2016. Т. 59, № 11. С. 959—963.
SOFTWARE COMPLEX FOR IDENTIFICATION OF CHAOTIC PARAMETERS OF CEREBRAL CIRCULATION
A. A. Musaev 12, A. I. Zagaynov 3
1 St. Petersburg State Technological Institute (Technical University), 190013, St. Petersburg, Russia
2St. Petersburg Institute for Informatics and Automation of the Russian Academy of Sciences,
199178, St. Petersburg, Russia E-mail: amusaev@technolog.edu.ru
3 A. F. Mozhaysky Military Space Academy, 197198, St. Petersburg, Russia
Features of developed automated system for identification of chaotic parameters of cerebral circulation are considered. Methods for numerical implementation of the fractal dimension and implementation of the calculations in real time are described. The software interface developed in Qt Creator environment is demonstrated. An example of construction of fractal dimensions of cerebral circulation of a patient calculated in real time using Dopplerograph Multi-DopX is presented.
Keywords: reconstructed attractor, phase space, correlation dimension, linear velocity of blood flow, systemic arterial pressure
Data on authors
Alexander A. Musaev — Dr. Sci., Professor; St. Petersburg State Technological Institute (Technical University), Department of System Analysis and Information Technologies; Head of the Department; SPIIRAS, Laboratory of Information Technologies in System Analysis and Modeling; E-mail: amusaev@technolog.edu.ru Artem I. Zagaynov — A. F. Mozhaysky Military Space Academy, Department of Mathematics
and Software; Senior Lecturer; E-mail: zagainov239@gmail.com
For citation: Musaev A. A., Zagaynov A. I. Software complex for identification of chaotic parameters of cerebral circulation // Izv. vuzov. Priborostroenie. 2016. Vol. 59, N 11. P. 959—963 (in Russian).
DOI: 10.17586/0021-3454-2016-59-11-959-963
