CC BY
36
эффективного медико-технологического процесса. Эти проекции могут быть найдены с использованием современной тепловизионной техники.
Во-вторых, такие системы должны быть снабжены блоками для объективной оценки эффекта аналгезии при выбранном режиме стимуляции.
В-третьих, подсистема стимуляции должна быть самонастраивающейся по сигналам от блока объективного контроля уровня боли.
И, наконец, в-четвертых, система коррекции боли должна включать в себя подсистему биологической обратной связи, которая будет модулировать режимы стимуляции проекций разных отделов головного мозга, связанных с болью. Это необходимо для индивидуализации режимов стимуляции в связи с тем, что восприятие боли может зависеть от психогенных, этнических и культуральных особенностей людей.
УДК 612.821
Е.Н. Винарская, Р.А. Кууз, Г.И. Фирсов
ПРОБЛЕМЫ ОЦЕНКИ ДИСКРИМИНАНТНО-ПРОГНОСТИЧЕСКОЙ СПОСОБНОСТИ КОЛИЧЕСТВЕННЫХ ХАРАКТЕРИСТИК НЕЛИНЕЙНОЙ ДИНАМИКИ ПРОЦЕССОВ РЕГУЛЯЦИИ ПОЗЫ В ЗАДАЧАХ КЛИНИЧЕСКОЙ НЕВРОЛОГИИ И СПОРТА
Основные пути совершенствования метода компьютерной стабилографии, направленные на непосредственное его применение в физиологических, психологических и медицинских исследованиях, зависят, в частности, от разработки способов оценки стабило-грамм с позиции теории сигналов и систем. Такой подход требует привлечения достаточно хорошо формализованных способов описания целостной структуры стабилограммы. В настоящее время известно два основных подхода к анализу стабилограмм - либо отдельно в сагиттальном и фронтальном направлениях, либо совместно. Первый подход достаточно традиционен, его широко используют при оценке характера стабилографической кривой и ее изменениях при различных патологиях или функциональных нагрузках [1,2]. Второй подход связан как с рассмотрением совместно стабилограмм в сагиттальном и фронтальном направлениях, так и анализом траектории проекции общего центра давления на плоскость опоры. В последние годы появились работы, в которых к анализу стабилографиче-ской информации применяются методы нелинейной динамики, в частности, расчеты корреляционной размерности, ляпуновских показателей, энтропии Шеннона-Колмогорова-Климонтовича [3-8].
В литературе, посвященной анализу процессов регуляции вертикальной позы, неоднократно высказывалось мнение об отсутствии корреляционных связей между сагиттальной и фронтальной составляющими стабилограммы. Обосновывалось это обычно малыми значениями коэффициента корреляции r. Наши эксперименты показали, что подобное заключение не всегда справедливо, связь между ортогональными составляющими стаби-лограммы может быть нелинейной, что не может быть установлено с помощью коэффициента корреляции. Кроме того, корреляционная связь между фронтальной и сагиттальной стабилограммами зависит от функционального состояния человека.
Поэтому при обработке стабилограмм помимо вычисления значений коэффициента корреляции должен выполняться расчет взаимнокорреляционной функции, дисперсионного отношения и взаимной дисперсионной функции [9].
Дисперсионное (корреляционное) отношение случайной величины Y относительно X определяется как
2 БЫ (Х\У)
пу\* - б
где
БЫ(7|7) -Ы[Ы(7|Х)-Ы(7)]2 - | | уф(у\х)йу-т
д>(х)ёх - диспер-
сия условного математического ожидания, характеризующая ту часть флуктуаций переменной У, которая вызвана влиянием переменной X. Дисперсионное отношение пх(У) определяет, в какой мере соблюдается функциональная зависимость между переменными х и у. При пх(у) « 1 зависимость практически функциональная; чем ближе г|х(у) к нулю, тем она больше нарушается. В случае независимости г|х(у) = 0. В общем случае 0 < пх(у) < 1. Дисперсионное отношение пх(у) можно интерпретировать как количественную характеристику меры определенности случайной величины у по значениям случайной величины х. При этом, дисперсионное отношение всегда больше или равно коэффициенту корреляции и не является симметричным, т.е. пх(у) ^ пу(х).
Использование дисперсионного отношения позволяет выявить наличие функциональной связи между колебаниями тела в сагиттальном и фронтальном направлениях, особенно у больных с органическими поражениями ЦНС. Так, для больного А.Б. (диагноз - рассеянный склероз), коэффициент корреляции составил 0,03, а дисперсионное отношение было больше чем на порядок (п = 0,28). Вместе с тем, при некоторых заболеваниях отличие величины дисперсионного отношения от коэффициента корреляции было сравнительно небольшим, так у больной С.И. (паркинсонизм) при значении коэффициента корреляции 0,72 величина дисперсионного отношения (как во фронтальном направлении по отношению к сагиттальному, так и наоборот, составил примерно 0,73). Такие большие значения г и п свидетельствуют о возможном наличии единого мощного источника колебаний по обеим плоскостям.
Дисперсионное отношение п и коэффициент корреляции г достаточно полно описывают общую форму статистической безынерционной взаимосвязи двух случайных процессов, подразумевающей отсутствие сдвига во времени (по фазе) между значениями двух связанных процессов. Рассмотрение же динамической системы регуляции позы делает целесообразным обращение к взаимному дисперсионному анализу [9]. Взаимная дисперсионная функция для каждой пары значений її, ^ равна дисперсии условного математического ожидания сечения одной функции х относительно сечения другой функции у сдвинутой на интервал ^, - Ъ
вух (/,, /2) - ы [ы (7 їх,) - М7ї]2 -
- | | у^Су ;«, |х2 ;/і)ауч - {у,уу(у«, 'АЩ
(х
где ^(У1',^ \х,г^2) - условная плотность вероятности У(^) относительно Х(г2);
фу ( у,;А) и (х2;^) --одномерные плотности вероятности случайных функций У(^) и Х(12).
Как и взаимная корреляционная функция, взаимная дисперсионная функция позволяет оценить величину инерционной статистической связи процессов во времени. Количественной мерой величины этой связи служат коэффициенты максимальной корреляции (и, соответственно, максимального дисперсионного отношения), равные максимальным значениям соответственно взаимной корреляционной или дисперсионной функции. Указанные коэффициенты позволяют оценить степень связанности двух процессов, даже если между ними имеются фазовые сдвиги.
Так, для случая функционального левостороннего гемипареза у больного В.В. при значении коэффициента корреляции при стоянии с закрытыми глазами 0,0677 максимальный коэффициент корреляции составил 0,4448, что говорит о наличии определенной линейной инерционной статистической связи колебаний во фронтальной и сагиттальной плоскостях. При этом у здоровых обследуемых различие в величинах дисперсионного отношения максимального дисперсионного отношения очень мало, не более 10-15%.
При некоторых видах неврологических патологий наблюдается возбуждение взаимосвязанных колебаний центра тяжести, в том числе при неврозах и истерии [10]. Об этом может свидетельствовать замеченная нами высокая степень когерентности колебаний на стабилограммах в саггитальной и фронтальной плоскостях в достаточно широких частотных диапазонах (от 4 до 8 Гц).
Большое значение функции когерентности свидетельствует о наличии сильной линейной инерционной статистической связи колебаний в двух плоскостях и может быть объяснено либо существованием единого мощного источника колебаний в ЦНС, либо синхронизацией колебаний в саггитальной и фронтальной плоскостях. Выявленный феномен может быть проявлением нарушения программирования регуляции равновесия вертикальной позы и заинтересованности систем, отвечающих за это программирование, в частности, лобных отделов. В частности, возможным механизмом, обусловливающим появление площадок функции когерентности, выступает электрическая деятельность мозга, проявляющаяся в тета-ритме ЭЭГ с частотой 4-7 Гц и амплитудой в десятки мкВ, при этом выраженность тета-ритма зависит от степени эмоционального, умственного напряжения, фона основной активности и возраста. Существующая взаимосвязь между активностью тета-ритма и умственным напряжением, эффективностью деятельности объясняется тем, что тета-ритм отражает активность срединно-стволовых образований головного мозга и является электрофизиологическим коррелятором механизма, квантующим поток извлекаемых из памяти энграмм.
Для оценки связанности колебаний в сагиттальной и фронтальной плоскостях весьма перспективно использование взаимной обобщенной информации, подробно рассмотренное в работе [5].
Для количественной оценки степени различия стабилограмм, в которых проявляется системная организация целенаправленной деятельности человека-оператора, нами использовались некоторые показатели степени хаотичности процессов перемещения центра тяжести человека, зафиксированные в стабилограммах. В частности, рассматривалось применение физической энтропии Больцмана-Шеннона [6], которая может быть вычислена по одной экспериментальной реализации, для выявления изменений в состоянии человека. Энтропия является относительной характеристикой, поэтому она не может быть использована при диагностике состояния различных людей. Однако ее удобно использовать для слежения за изменением состояния человека при предъявлении ему различных тестовых задач, т.е. для решения задачи мониторинга.
Кроме шенноновской энтропии эффективным представляется использование энтропийного коэффициента плотности распределения вероятности KH = 0,5ехр[/ш (п)]ст где 1ш (n) - энтропия плотности распределения вероятностей (информация по Шеннону), а ст - стандартное отклонение. Зная гистограмму плотности распределения вероятности мгновенных значений стабилограммы, энтропийный коэффициент может быть определен bN _!/N d
по формуле кн = “—10 2 ni lg nt , где b - ширина столбца гистограммы, N - объем
2ст 1=1
выборки, ст - СКО, d - число столбцов гистограммы, n - число наблюдений в i-м столбце. При этом для любых законов распределения величина KH лежит в пределах 0 -г- 2,066, причем максимальное значение KH - 2,066 имеет гауссово распределение, а, например,
арксинусоидальному распределению соответствует значение KH = 1,11.
Другим подходом к оценке гистограммных характеристик стабилограмм является использование как коэффициента эксцесса кэ, изменяющегося от 1 до œ, так и контрэксцесса k = k31, значение которого может меняться в пределах от 0 до 1. Рассеяние оценки контрэксцесса к приближенно аппроксимируется формулой S(k) = с{к)/k « 4\s2 -1)3 /V29Ñ
где е - эксцесс распределения, а рассеяние оценок энтропийного коэффициента KH может быть найдено из соотношений ст( Kh ) = 0,9/(kKH( KhN )05) и S(Kh) = v(Kh)KH = 0,9/(kK_2 (KhN)0’5).
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. ГурфинкельВ.С., Коц Я.М., ШикМ.Л. Регуляция позы человека. - М.: Наука, 1965. - 256 с.
2. Скворцов Д.В. Клинический анализ движений. Стабилометрия. - М.: АОЗТ "Антидор", 2000. - 192 с.
3. Розенблюм М.Г., Фирсов Г.И. Стохастические автоколебания в системе регуляции вертикальной позы тела человека. I. Стратегия управления позой и динамическая модель // Биомеханика (София). Т. 24. - 1992. - С.34-41. II. Модельный и натурный эксперимент // Биомеханика (София). Т. 25. - 1992. - С.37-43.
4. Firsov G.I., Rosenblum M.G., Landa P.S. Deterministic 1/f fluctuations in biomechanical system // Noise in physical systems and 1/f fluctuations. - New York: AIP Press, 1993. - Р. 716719.
5. Rosenblum M.G., Firsov G.I., Kuuz R.A., Pompe B. Human Postural Control - Force Plate Experiments and Modelling// Nonlinear Analysis of Physiological Data. / Ed. by H.Kantz, J.Kurths and G.Mayer-Kress. - Berlin: Springer, 1998. - P. 283-306.
6. Винарская Е.Н., Кууз Р.А., Розенблюм М.Г., Фирсов Г.И. Методы нелинейной динамики в задачах оценки функционального состояния человека // Известия ТРТУ. Медицинские информационные системы. - Таганрог: Изд-во ТРТУ, 2000, № 4 (18). - С. 7 - 10.
7. Кононов А.Ф., Слива С.С. Применение методов нелинейной динамики для анализа системы поддержания вертикальной позы человека // Российский журнал биомеханики.
2000. - № 2. - С. 70.
8. Колесников А.А., Слива С.С., Кононов А.Ф. Аналитическое конструирование агрегиро-ваннных регуляторов: синергетическое управление биомеханическими системами // Управление и информационные технологии. Том. 1. - СПб.: СПбГЭТУ “ЛЭТИ”, 2003. -С. 42-46.
9. Райбман Н.С., Капитоненко В.В., Овсепян Ф.А., Варлаки П.М. Дисперсионная идентификация / Под ред. Н.С. Райбмана. - М.: Наука, 1981. - 336 с.
10. Кууз Р.А., Фирсов Г.И. Применение методов компьютерной стабилометрии для решения задач функциональной диагностики в неврологии // Биомедицинская радиоэлектроника.
2001. - № 5-6. - С. 2433.
УДК 612.821
И.Н. Статников, Г.И. Фирсов
ДИСПЕРСИОННЫЙ АНАЛИЗ РЕЗУЛЬТАТОВ КОМПЛЕКСНОГО ЭКСПЕРИМЕНТАЛЬНОГО ИССЛЕДОВАНИЯ ПРОЦЕССОВ ВОСПРИЯТИЯ И ПЕРЕРАБОТКИ ВИБРОТАКТИЛЬНОЙ ИНФОРМАЦИИ ЧЕЛОВЕКОМ-ОПЕРАТОРОМ В ЭКСТРЕМАЛЬНЫХ УСЛОВИЯХ
Исследования сложных биомеханических систем класса человек-оператор, находя-
