CC BY
55
УДК 57.087.1
Российский
Журнал
Биомеханики
www.biomech.ru
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ВЕСТИБУЛЯРНОГО НИСТАГМА. ЧАСТЬ II. МЕХАНИКО-ИНФОРМАЦИОННАЯ МОДЕЛЬ
ОАО «ВымпелКом», 127083, Россия, Москва, ул. Восьмого марта, 10, стр. 14, e-mail: [email protected]
Аннотация. Предложена механико-информационная модель вращательного горизонтального вестибулярного нистагма. Эта модель в зависимости от вестибулярной стимуляции горизонтальных полукружных каналов и отолитовых рецепторов способна описать ответное нистагменное движение глаз. Параметры модели определены на основе литературных данных, а также по результатам лабораторных исследований. Модель носит стохастический характер и включает в себя моделирование медленной фазы нистагма, логику переключения между медленной и быстрой фазами, моделирование работы саккадического механизма как марковского процесса и построение управлений для глазодвигательных мышц. Для полученных в результате моделирования траекторий были вычислены значения таких параметров нистагма, как длительности и амплитуды быстрых и медленных фаз нистагма. Построены гистограммы для всех четырех параметров нистагма. Для удобства сравнения результатов моделирования с результатами эксперимента гистограммы для всех параметров были объединены в одну. Проведено сравнение распределений параметров нистагма, полученных для механико-информационной модели и статистической модели, с распределениями параметров нистагма, записанного в ходе исследования здорового испытуемого в лабораторных условиях. Кроме того, проведено сравнение распределений параметров нистагма, полученных для механико-информационной модели, с распределениями параметров нистагма, полученных для статистической модели. Для сравнения распределений параметров использовался индекс Чекановского-Съеренсена. Результаты всех трех сравнений позволяют утверждать, что распределения параметров модельных нистагмов, полученных с помощью механико-информационной и статистической моделей, достаточно хорошо соответствуют как друг другу, так и распределению параметров нистагма, записанного в лабораторных условиях.
Ключевые слова: вестибулярный нистагм, цепь Маркова, модель Голдберга-Фернандеса, тест Лиллиефорса, индекс Чекановского-Съеренсена.
Вестибулярный нистагм является одним из древнейших эволюционно закрепившихся рефлексов, играющих важную роль в жизни человека.
Моделирование этой вестибуло-глазодвигательной реакции, как и других, является насущной задачей. Можно указать несколько прикладных областей, в которых имеется актуальная потребность в таких моделях. Прежде всего следует упомянуть медицинские приложения, например [2, 4, 18], в которых использование математических моделей позволяет оценивать тяжесть состояния больного и следить за динамикой заболевания. Показано также, что мониторинг вестибуло-окулярных реакций может быть полезен для оценки функционального состояния человека-оператора, подверженного
© Сучалкина А.Ф., 2014
Сучалкина Анна Федоровна, старший эксперт ОАО «ВымпелКом», Москва
А.Ф. Сучалкина
Введение
высокому уровню напряженности труда и профессионального стресса [3]. Моделирование вестибуло-окулярных реакций, несомненно, необходимо при решении весьма актуальной и интересной задачи создания вестибулярного протеза [8]. Наконец, модели движения глаз человека могут быть полезны при создании или проектировании биомехатронных систем, опирающихся на зрение человека [11], или при разработке антропоморфных роботов, обладающих системой технического зрения [14].
Обзор литературных источников, посвященных вопросам моделирования вестибулярного нистагма, приведен в работе [1]. Там же приводится описание использованного экспериментального оборудования и сценариев проведения лабораторных исследований.
При математическом описании модели будут использованы системы координат и обозначения промежутков времени и переменных, введенные в работе [1].
Структура модели
В работе рассматривается моделирование движения глаз испытуемого, сидящего прямо на вестибулярном кресле в условиях темноты, для случая пассивного вращения вокруг вертикальной оси. На рис. 1 приведена общая схема механико-информационной модели вестибулярного нистагма. Модель получила такое название, потому что в ней воспроизводятся как механические, так и информационные процессы, протекающие при реализации нистагма.
В структуре модели можно условно выделить четыре основных блока:
1. Оценка движения головы и формирование медленной фазы нистагма. Входным звеном вестибуло-окулярных реакций является механическое движение купул полукружных каналов и мембран отолитов, вызываемое воздействием сил инерции, возникающих при ускоренном движении головы. Связанный с этим движением изгиб волосковых пучков волосковых клеток изменяет ионные токи, протекающие в вестибулярном рецепторе. В модели в качестве носителя полезного сигнала рассматривается изменение частоты афферентной импульсации в первичном вестибулярном нейроне. После этого моделируются оценки углового ускорения головы, которые строятся вестибулярными ядрами по информации, поступающей от полукружных каналов и отолитов. На основе этих оценок, взятых с соответствующими весовыми коэффициентами, вычисляется относительное угловое движение глаз, которое после наложения на него быстрых фаз будет составлять медленную фазу нистагма.
2. Логика переключения фаз. Переключение от медленной фазы к быстрой происходит в те моменты времени, когда некоторая бинарная случайная величина, подчиняющаяся марковскому закону, принимает значение, равное единице, а также одновременно выполнены следующие четыре условия:
1) абсолютная величина относительного угла поворота глаза превысила заданное пороговое значение;
2) амплитуда угла поворота глаза за текущую медленную фазу превысила пороговое значение;
3) величина угловой скорости глаза превысила пороговое значение;
4) от момента окончания предыдущей быстрой фазы прошло достаточное время, которое исключает наложение быстрых фаз.
3. Моделирование быстрой фазы. Для задания быстрой фазы нистагма определяются ее амплитуда и длительность как случайные величины с известным законом распределения, зависящим от текущего значения абсолютной угловой скорости головы. Траектория движения глаза во время быстрой фазы представляет собой косинусоиду, заданную на половине периода.
Рис. 1. Общая схема модели нистагма
4. Задание движения глаз. При движении глаза в горизонтальной плоскости поворот осуществляется за счет согласованного сокращения и расслабления внутренней и наружной прямых мышц глаза. В данном блоке формируются команды для этих мышц.
Формирование медленной фазы нистагма
Для моделирования движения купулы полукружного канала в условиях ускоренного углового движения головы используется первая часть модели Голдберга-Фернандеса [13], записанная в виде передаточной функции:
Х1Х2
*5,/С0 = (- l№0>w(» = (- 1)г + + ^ sa(s), i =1; 2,
где xsj описывает угловое смещение купулы; sa - угловое ускорение головы; xj и т2 ~ большая и малая постоянные времени полукружного канала; i = \;2 обозначают левый и правый полукружные каналы соответственно.
На выходе получается величина vsi, которая пропорциональна изменению афферентной импульсации, возникающей в вестибулярных механорецепторах:
TaS
vs,iО) = Я20)х(» = kQ—^--(1 + Tis)xsi(s), i= \;2,
(l+T^)
, спайк/с
где Kq - константа, имеющая размерность Град .
Поскольку вестибулярные рецепторы являются эволюционно одними из самых древних, предполагается, что в центральной нервной системе сформировалась внутренняя модель вестибулярного аппарата [17], которая позволяет получить внутреннюю оценку углового ускорения головы. Таким образом, величина
sïï^s) = (- \ун10)-я20)-^0) =
Kq
= f_ , V (1 +Т!5)(1 +Т2^) 1 (1 + ( л . = 1 9
моделирует внутренние оценки абсолютного углового ускорения головы, построенные в центральной нервной системе по информации, поступающей от полукружных каналов.
Линейное ускорение, возникающее во время вращения вокруг вертикальной оси, можно разложить на две составляющие: тангенциальную и центростремительную. Соответственно, на мембраны отолитов действуют силы инерции. Сила, направленная тангенциально, воздействует на утрикулярную и саккулярную мембраны, а направленная центробежно - только на утрикулярную. Вследствие воздействия сил инерции происходит смещение мембран отолитового аппарата.
Согласно модели, предложенной в работе [17], для случая вращения вокруг вертикальной оси можно написать следующие уравнения, позволяющие для утрикулюсов получить величины vui j(s) и vu2 ¡(s) и для саккулюсов - величины vsac i(.s), которые пропорциональны изменению частоты афферентной импульсации (спайк/с), порождаемой механорецепторами отолитовых органов:
(1+X,s)
VMl,/0) = (- 1)'#о(»юС0р = (- 1УК^—----Т5ю(5)р,
(1 +X0iJj(l +1o2S)
VU2,i0) = (- l^OVCOp = (- 1 УК , (1 ^--ю20)p,
(1 +Xols)(l +Xo2s)
(1+X,5)
(1 +X0iJj(l +1o2S)
где К, x0i, x02, ts - известные константы; i = 1; 2 - для левого и правого утрикулюсов и саккулюсов; р - расстояние от центра масс утрикулюса и саккулюса до оси вращения.
В центральной нервной системе образуются внутренние оценки абсолютного углового ускорения головы, построенные по тангенциальным составляющим ускорения:
1 — 1 • 1 С1 + хо1^)С1 + хо2 5)
•*юИ1,/С0 = -(-1 )гя 0КМ0) = -(-1 У--——--
' р 0 ' р К (1 + т^) '
_ 1 — 1 . 1 (1 +хо15)(1 +То2 я)
= "НЖОК^О) = -(- 1 У--—---
р Р К (1 + ъя) '
Рассмотрим утрикуляриый сигнал о центростремительном ускорении, для которого можно записать приближенное равенство:
=/(ю2) -/(®о) ~ 2Ю0/(Ю^)АЮ.
Таким образом, по информации, поступающей от утрикулюса, можно построить оценку абсолютной угловой скорости головы
~ ( \ ("О' тт ( \ ( \ ("О' (1 + Т01^)(1 + Х02^) ( .
®м2,г) =- •/ 7\ Н0{$)\и2Дд) =-—1—рГ -/1 ч- ^Д*) -
2ш0/(со§)р 2^соо/(соо)р +
Для получения осредненной оценки абсолютного углового ускорения головы выписывается следующее соотношение:
srö(s) = |
^uu(s) + SO)uU(s) S(äu2J(■S') + ) SGW,! (,?) + S(0sac2 (s)
2 2
+
2
Далее в модель включается уравнение, связывающее оценку абсолютного углового ускорения головы и командное относительное угловое ускорение глаза. Благодаря вестибуло-окулярному рефлексу происходит компенсаторный поворот глазных яблок с ускорением, противоположным по направлению абсолютному ускорению головы. Как известно [6], отношение среднего относительного углового ускорения глаза на протяжении медленной фазы к среднему абсолютному угловому ускорению головы на том же промежутке называют коэффициентом усиления к5. Коэффициент усиления не является постоянным и изменяется от цикла к циклу, поэтому при моделировании медленной фазы он считается случайной величиной с известным законом распределения. Таким образом, в модели определяется
52фДя) = -кц (и)^ ю(я),
где - командное относительное угловое ускорение глаза; п - номер текущей медленной фазы. Поскольку в начальный момент времени требуется выполнение условий 5фД0) = ю(0) = 0 и фДО) = 0(0) = 0, получается
Ф.О) (л)0О),
где ф^ - командный угол поворота глаза на протяжении медленной фазы; 0 - оценка абсолютного углового положения головы.
Переключение между медленной и быстрой фазами
Саккадический механизм, использующий обратную связь от мышц глаз [16], обеспечивает переключение фаз нистагма. Для моделирования работы саккадического механизма рассматривается бинарная случайная величина г, принимающая значения 0
или 1 по заданному вероятностному закону (рис. 2). Матрица вероятностей переходов от одного значения к другому имеет вид
где р и д - вероятности того, что переменная остается в состоянии 1 и О соответственно. В дискретные моменты времени с шагом А? = 0,001 с разыгрываются значения переменной г.
Считается, что пока величина г имеет значение 0, переключение от медленной фазы к быстрой невозможно, а как только случайная величина принимает значение 1, проверяются условия возможности переключения к быстрой фазе.
Текущая медленная фаза нистагма с номером п в момент времени прерывается быстрой фазой, только если одновременно выполняются следующие пять условий:
1) случайная величина г, заданная марковским законом с приведенной выше матрицей вероятностей переходов, в момент времени 4 приняла значение, равное 1;
2) командный угол поворота глаза превысил пороговое значение фт:
3) амплитуда поворота глаза за текущую медленную фазу превысила пороговое значение фа:
где t¿ - некоторая заданная константа.
Неравенства 2), 3) и 4) задают минимальные угол поворота глаза, амплитуду поворота глаза и скорость поворота глаза, при которых медленная фаза нистагма прерывается быстрой фазой. Условие 4) не позволяет прерывать быстрыми фазами медленные фазы непосредственно после того, как произошла смена знака угловой скорости головы.
P = [/ W
|фХ4 )| > Ф™;
|ф*(4) -Ф(^)| > Фа;
4) относительная угловая скорость глаза превысила пороговое значение ют:
5) чтобы не произошло наложения быстрых фаз, проверяется условие
и,
q
Р
Рис. 2. Состояния и вероятности переходов
Таким образом, когда угол поворота глаза превышает значение фш, амплитуда медленной фазы нистагма превышает значение фа, скорость поворота глаза превышает значение юш и не происходит наложения быстрых фаз, как только цепь Маркова генерирует значение г=\, медленная фаза нистагма прерывается быстрой фазой.
Моделирование быстрой фазы
Моделирование быстрого саккадического движения для каждой конкретной саккады осуществляется определением значений двух параметров: амплитуды А и длительности Т саккады, как это было предложено в работе [10].
Для задания быстрой фазы нистагма используется часть косинусоиды, причем при определении направления саккадического скачка учитывается, является на момент переключения функция фДО, описывающая медленную фазу нистагма, возрастающей или убывающей:
' V
Здесь Ау и Ту - амплитуда и длительность текущей быстрой фазы, случайные величины, формирующиеся следующим образом. Длительность Ту генерируется как случайная величина, имеющая усеченное нормальное распределение на отрезке \Т\; Т2] (табл. 1). После этого генерируется случайная величина Ау, имеющая усеченное нормальное распределение на отрезке \А\,Апричем дополнительно проверяется
А/
условие, что отношение — попадает в интервал [vi; у2] (см. табл. 1).
Т/
Таким образом, полностью сформирован закон движения глаза на протяжении текущего цикла нистагма
Фс= 1 С(0 - + Фс(//), / 6 ^ Ту], ^ + Т"
Здесь фс - команда, подающаяся на глазодвигательные мышцы.
Управление мышцами глаз
На основании результатов исследований можно сделать вывод, что у здорового человека значения параметров нистагма для правого и левого глаз практически совпадают, поэтому далее строится модель движения одного глаза. Для примера рассматривается правый глаз.
Задаются управления для внутренней и наружной прямых мышц глаза, при этом используется композиционный способ построения управления [12]:
V/ (0 = (0 +
где Vj - управлениеу-й мышцы, имеющее следующие составляющие: - тоническую, Vsj - статическую, - кинематическую, Ущ - динамическую; управления имеют размерность силы, у = 1; 2.
Составляющие задаются следующим образом. Тонические компоненты постоянны и представляют собой долю ц максимальной силы ушах, развиваемой мышцей:
Vn (О = VT2{t) = pVmax.
Статические компоненты определяются упругими свойствами мышц: vsl{t)=kRqc(0, vS2(t)=-kRqß),
где к - коэффициент упругости мышцы; Я - радиус глаза. Кинематические компоненты характеризуют вязкие свойства мышц:
где g - коэффициент вязкого трения. Наконец, динамические компоненты управления определяются инерционными характеристиками глаза:
где М- масса глаза.
Обозначив силы, развиваемые мышцами глаза, через ^ и запишем уравнения
где - текущая длина мышцы; /0 - длина мышцы в ненапряженном состоянии, поэтому 1\ (?) =/0-^фс(0 и /2(0 =/о +^фс(0. Величина силы, развиваемой мышцей, считается положительной, если она является тянущей. Завершает модель уравнение движения глаза
Результаты моделирования. Сравнение с результатами исследований
В табл. 1 приведены значения констант модели и указаны законы распределения параметров, являющихся случайными величинами.
Среди перечисленных выше параметров обсудим дополнительно коэффициент усиления ks, который не обсуждался в работе [1]. Его значения были вычислены как отношения амплитуды поворота глаза к амплитуде поворота головы для всех медленных фаз нистагма. Распределение значений ks показано на гистограмме (рис. 3). С помощью теста Лиллиефорса в приложении MATLAB проверено, что значения ks этой записи имеют усеченное нормальное распределение [15]. Уровень значимости Р выбран равным 0,05. Для выбранного массива числовых значений коэффициента усиления ks в приложении MATLAB были вычислены его математическое ожидание и среднеквадратическое отклонение, они внесены в табл. 1.
На основании выписанных выше уравнений и данных табл. 1 в среде Simulink приложения MATLAB написана программа, реализующая механико-информационную модель нистагма. На вход программы подается числовой массив, моделирующий экспериментальный стимул - угловое положение головы 0(i). Движение головы после начального поворота представляет собой синусоидальное вращение с периодом 8 с и амплитудой 60°. Детали эксперимента описаны в работе [1]. В результате выполнения программы получается модельная траектория, отвечающая вестибулярному нистагму, который возникает у здорового человека в ответ на входной стимул (рис. 4).
v*i(0 =gR^cit), vK2{t)=-gR<vc{t),
Fj it) = k(lj it) - /о) + gij (t) + V ,- (0, j = 1; 2,
Таблица 1
Параметры модели
Параметр Значение Размерность Источник
5,7 с Из литературы
0,003 с Из литературы
ко 29,24 спайк/(град-с) Подбирается по результатам испытаний для обеспечения качественного совпадения
80 с Из литературы
Т/ 0,049 с Из литературы
К 33,3 - Из литературы
^ 10 с Из литературы
5 с Из литературы
0,016 с Из литературы
Р 0 - Не изменяется
Ч 0,98 - Подбирается
аг- ц, = 0,08, а = 0,11 град Подбирается
Ь 7 - Подбирается
Т1 ц=37,с = 40 мс Подбирается
Тх 35 мс Подбирается
т2 160 мс Подбирается
ах 1,54 град Подбирается
а2 33 град Подбирается
VI 0,03 град/с Подбирается
0,1015 град/с Подбирается
120 мс Подбирается
Фа 1 град Из собственных исследований
Ф 1,1 град Из собственных исследований
1 град/с Из собственных исследований
К ц,= 0,48, а = 0,2 - Из собственных исследований
V- 0,05 - Из литературы
"^тах 0,0245 Н Из литературы
К 0,0856 Н/м Из литературы
К 0,0127 м Из литературы
с 0,0342 Н-с/м Из литературы
¿о 0,03 м Из литературы
м 0,0107 кг Из литературы
Рис. 3. Гистограмма коэффициента усиления
Рис. 4. Фрагмент модельной траектории
Было проведено сравнение распределений параметров нистагма, полученных в результате моделирования, с распределением параметров нистагма, записанного при исследовании одного из испытуемых, а также с распределением параметров траектории, полученной с помощью статистической модели [1]. Сравнивались распределения четырех параметров: длительностей Ту и Тх и амплитуд Ау иАх быстрых и медленных фаз нистагма. Для сравнения были построены гистограммы этих параметров (рис. 5). Величина экспериментальной выборки (количество циклов нистагма) п\ = 53, величины выборок, полученных для механико-информационной и статистической моделей, одинаковы и равны я2 = ^з = 63. По формуле к — 5^(п), где п = шт(п1, п2, п3), вычислено количество интервалов для гистограмм к = 9 [9].
Для попарного сравнения распределений использовался индекс /чс Чекановского-Съеренсена [7]. Поскольку для сравнения выборок параметров играет роль только количество частот, попавших в каждый из интервалов, гистограммы для всех четырех параметров были объединены в одну гистограмму. Таким образом, объединенная гистограмма имеет 36 столбцов. По формуле 2к-2 число степеней свободы получается равным 70 [5]. Статистический уровень значимости Р выбран равным 0,05. В табл. 2 приведены значения индекса /чс, а также величина t, которая показывает статистическую значимость отличия дополнения (1 - /чс) от нуля, и критическое значение, взятое из таблицы Стьюдента.
Рис. 5. Гистограммы параметров нистагма (эксперимент (слева), механико-информационная (в середине) и статистическая (справа) модели): а - длительности быстрых фаз; б - длительности медленных фаз; в - амплитуды быстрых фаз; г -
амплитуды медленных фаз
Таблица 2
Сравнение распределений параметров экспериментального и модельных нистагмов с помощью индекса Чекановского-Съеренсена
Сравнение ^чс t ¿Cm
Эксперимент - механико-информационная модель 0,825 1,483 1,994
Эксперимент - статистическая модель 0,884 0,956
Механико-информационная модель - статистическая модель 0,833 1,475
В результате всех трех сравнений оказалось, что параметр t имеет значения, существенно меньшие критического значения . Этот факт позволяет сделать вывод о том, что в этих трех случаях нулевая гипотеза о принадлежности выборок одной генеральной совокупности статистически подтверждается, т.е. распределения параметров модельных нистагмов, полученных с помощью механпко--информационной и статистической моделей, достаточно хорошо соответствуют как друг другу, так и распределению параметров нистагма, записанного в лабораторных условиях.
Заключение
На базе вестибулометрического комплекса проведено изучение вестибулярного нистагма, возникающего при вращении с ускорением. На основе обработки экспериментальных записей нистагменных циклов создана механико-информационная модель вестибулярного нистагма. Для этой модели подобран такой набор параметров, что модельная траектория глаз имеет те же характеристики, что и запись нистагма здорового испытуемого. Распределения параметров нистагма, записанного в эксперименте, сравнены с помощью индекса Чекановского-Съеренсена с аналогичными распределениями параметров модельных траекторий, полученных с помощью механико-информационной и статистической моделей. При выбранном уровне значимости можно говорить об отсутствии достоверного различия между выборками. Таким образом, можно заключить, что была успешно построена механико-информационная модель вестибулярного нистагма и подобраны значения ее параметров для одного из испытуемых. Предложенная модель является достаточно компактной, при этом в ней нашли отражения основные механические и информационные процессы, протекающие во время реализации вестибулярного нистагма.
список литературы
1. Боков Т.Ю., Сучалкина А.Ф., Якушева Е.В., Якушев А.Г. Математическое моделирование вестибулярного нистагма. Часть I. Статистическая модель II Российский журнал биомеханики. -2014. - Т. 18, № 1. - С. 48-64.
2. Доценко В.И., Егорова Е.А., Каспранская Г.Р., Муратова Е.А., Сучалкина А.Ф., Якушев А.Г. Экспериментальное изучение и математическое моделирование вестибулярного нистагма у здоровых людей и больных ДЦП II Фундаментальная и прикладная математика. - 2005. - Т. 11, № 8. -С. 195-204.
3. Доценко В.И., Штефанова О.Ю., Якушев А.Г., Боков Т.Ю., Якушев A.A. Новые алгоритмы оценки морфофункциональных нарушений работы головного мозга и глубины экспериментальных патологических состояний организма с использованием показателей глазодвижений II Фундаментальные, клинические и гигиенические основы и аппаратно-методическое обеспечение системы медико-психологической реабилитации пациентов, подверженных высокому уровню напряженности труда и профессионального стресса: сб. тр. всерос. науч. шк. для молодежи. -Ростов н/Д: Изд-во Юж. федер. ун-та, 2011. - С. 48-59.
4. Доценко В.И., Штефанова О.Ю., Якушев А.Г., Каспранская Г.Р., Боков Т.Ю. Отображение степени поражения головного мозга и динамики восстановительных процессов у неврологических больных с использованием новых показателей глазодвижений II Вестник восстановительной медицины. -2011, №5,-С. 6-10.
5. Орлов А.И. О проверке однородности двух независимых выборок II Заводская лаборатория. - 2003. -Т. 69, № 1. - С. 55-60.
6. Орлов И.В. Вестибулярная функция. - СПб.: Наука, 1998. - 248 с.
7. Песенко Ю.А. Принципы и методы количественного анализа в фаунистических исследованиях. - М.: Наука, 1982,- 287 с.
8. Садовничий В.А., Александров В.В., Александрова Т.Е., Вега Р., Сидоренко Г.Ю., Сото Е., Шуленина Н.Э. Динамическая имитация стабилизации и потери вертикальной позы и тестирование прототипов вестибулярного протеза II Современные проблемы математики и механики. Т. 1. Прикладные исследования. - 2009. - Вып. 1. - С. 154-164.
9. Шторм Р. Теория вероятностей. Математическая статистика. Статистический контроль качества. -М.: Мир, 1970. - 368 с.
10. Якушев А.Г. Композиционный способ программирования саккадического скачка глаза II Вестник Моск. ун-та. Сер. 1. Математика. Механика. - 2011. - № 1. - С. 70-73.
11. Якушев А.Г. О возможности визуального (с помощью взора) управления II Вестник Моск. ун-та. Сер. 1. Математика. Механика. - 2012. - № 4. - C. 32-36.
12. Якушев А.Г., Напалков Д.А., Ратманова П.О., Кручинина А.П., Штефанова О.Ю., Гинзбург Е.А. Композиционный способ определения управления глазодвигательными мышцами при саккаде II Росссийский журнал биомеханики. - 2011. - Т. 15, № 1. - С. 99-109.
13. Fernandez C., Goldberg J. Physiology of peripheral neurons innervating semicircular canals of the squirrel monkey. II. Response to sinusoidal stimulation and dynamics of peripheral vestibular system // J. Neurophysiol. - 1971. - Vol. 34. - P. 661-675.
14. Laschi C. Vision and eye movements in humans and robots [Электронный ресурс]. - URL: www.robot.uji.es/documents/iurs06/documents/Laschi_i.pdf (дата обращения: 20.02.2014).
15. Lilliefors H. On the Kolmogorov-Smirnov test for normality with mean and variance unknown // J. Amer. Stat. Assoc. - 1967. - Vol. 62, № 318. - P. 399-402.
16. Schmid R.M. System analysis of the vestibule-ocular system // Fifth symposium on the role of the vestibular organs in space exploration / Naval aerospace medical institute, Naval aerospace medical centre. 19-21 August, 1970, Pensacola, Florida. - P. 237-249.
17. Telban R.J., Cardullo F.M., Guo L. Investigation of mathematical models of otolith organs for human centered motion cueing algorithms [Электронный ресурс]. - URL: http://arc.aiaa.org/doi/pdf/ 10.2514/6.2000-4291 (дата обращения: 25.02.2014).
18. Zupan L.H., Merfeld D.M., Darlot C. Using sensory weighting to model the influence of canal, otolith and visual cues on spatial orientation and eye movements // Biol. Cybern. - 2002. - Vol. 86. - P. 209-230.
MATHEMATICAL MODELLING OF VESTIBULAR NYSTAGMUS. PART II. MECHANICO-INFORMATIONAL MODEL
A.F. Suchalkina (Moscow, Russia)
In the paper, a mechanicoinformational model of vestibular nystagmus is proposed. Model parameters are determined on the basis of published data as well as on the results of laboratory investigations. The model is of stochastic nature and involves the following blocks: simulation of nystagmus slow phase, logic of switch between the slow and fast phases, simulation of saccadic mechanism functioning as a Markov process and arrangement of signals to control the extraocular muscles. For the trajectories by the model for each nystagmic cycle, a set of four parameters are calculated. This set four parameters: durations and amplitudes of fast and slow phases of nystagmus. Histograms are drawn for all four parameters of nystagmus. To compare the simulation results with the experimental results, all four histograms are merged into one. Distributions of nystagmus parameters obtained from the mechanicoinformational and statistical models are compared with distribution of nystagmus parameters of a healthy individual registered in a laboratory environment. Besides that, the distributions of nystagmus parameters obtained from the mechanicoinformational model are compared with distribution of nystagmus parameters obtained from the statistical model. Czekanowski-Sorensen index is used for comparison of the recordings. The results of all three comparisons allow us to conclude that the mechanicoinformational model of vestibular nystagmus is successfully constructed and its parameter values are chosen for one of the subjects. The proposed model is fairly compact, while it reflects the basic mechanical and information processes that occur during the realization of vestibular nystagmus.
Key words: vestibular nystagmus, Markov process, Goldberg-Fernandez model, Lilliefors test, Czekanowski-Sorensen index.
Получено 20 марта 2014
