зультатам появится возможность проанализировать перспективы дальнейшего развития кластерного анализа. В работе осуществлены все этапы обработки больших данных текстовой информации, начиная с этапа предварительной обработки и до момента графического отображения, что также ново.
Практическая значимость результатов работы актуальна из-за стремительного развития электронного обучения во всем мире. Определение наиболее подходящих методов анализа, фреймворков, стеммера и других инструментов позволит обработать эти данные быстрее и уменьшить число ошибок.
Применение технологии больших данных даст возможность педагогам и сотрудникам образовательных учреждений получать актуальную информацию об обучающихся и реагировать на изменения процесса обучения, а в результате - прогнозировать популярность определенного вида деятельности, вносить изменения в учебный контент.
УДК 519.62, 57.034
В.В. Сельвинский, С.А. Журавлёв
МОДЕЛИРОВАНИЕ РАБОТЫ СЕРДЦА И НЕРВНЫХ ИМПУЛЬСОВ
В данной статье рассматриваются известные модели Зимана биений сердца и распространения нервных импульсов. Обе модели имеют как общие свойства, так и существенные различия. В статье исследуется влияние параметров указанных биологических процессов, таких как чувствительность, величина электрохимического импульса, на форму и частоту сердцебиений, амплитуду удлинения мышечного волокна. Исследование моделей проводится в пакете MathCad.
SIMULATION OF HEART AND NERVOUS PULSES
This article discusses the known models of Zeeman heart beat and propagation of nerve impulses. Both models have both common properties and significant differences. The paper explores the effect of parameters of said biological processes, such as sensitivity> magnitude of electrochemical pulse> on heart rate> amplitude of muscle fiber elongation and electrochemical potential. Model research is carried out in a MathCad package.
DOI: 10/22250/jasu.8
Сердце - важнейшая составная часть любого теплокровного организма, оно обеспечивает его жизнедеятельность достаточно продолжительное время - от рождения и до умирания. Вопросам функционирования сердечно-сосудистой системы посвящено огромное количество работ отечественных и зарубежных авторов, специалистов в области медицины, химии, биологии.
В данной работе рассматриваются чисто механические аспекты сердечных сокращений, а также формирования нервных импульсов, хотя мотивы для конструирования соответствующих математических моделей, безусловно, опираются на происходящие при этом химические и биологические процессы.
Известно, что рабочий цикл сердечной мышцы состоит из быстрого сокращения (систола), инерционной паузы на переключение, затем следует быстрое расслабление (диастола), опять инерционная пауза на переключение, далее опять следует быстрое сокращение, и циклы повторяются. Данный процесс описывается двумерной моделью Зимана, представленной в работе [1]:
£х = х-х3-Ь, Ь = х-х(), (1)
где х - удлинение мышечного волокна; Ь - характеристика электрохимического стимулирования сердечных сокращений (потенциал); s «1 - малый параметр реакции мышечного волокна на электрохимическое воздействие (всюду по ходу вычислений полагаем s - 0.001).
Фазовый портрет системы (1) содержит одну особую точку Е с координатами
X = Xq , Ъ — 6Q = Xq — XQ .
Наличие малого параметра s в первом уравнении системы (1) придает переменной х способность быстро изменяться по сравнению с переменной Ъ. Это свойство проявляется тем сильнее, чем дальше находится
изображающая точка (х,Ь) от кривой Ъ = х-х3 на фазовой
плоскости (рис. 1). Поэтому в удалении от кривой Ь = х-х3 фазовые траектории почти параллельны оси х; в окрестности этой кривой они изгибаются в сторону возрастания Ъ при х > х0, в сторону убывания Ъ - при х < х0. Также в правой
верхней части фазового портрета - там, где Ь> х-х3, изображающая точка (х,Ь) будет перемещаться преимущественно вниз вдоль фазовых траекторий, в левой верхней части, где
b < х-х - вверх.
Сделанные замечания допускают возможность существования цикла типа EABCDE при условии, что изображающая точка (х,Ь), попадая в некоторую 5-окрестность особой точки Е(х0,Ь0), далее получает импульс АЬ > 0, достаточный для ее выхода в окрестность точки А, а далее вдоль фазовой траектории через положения, близкие к В, С, D и опять Е. В данном случае размер 8 > 0 можно трактовать как чувствительность организма на состояние замирания мышечного волокна, в том смысле, что чем меньше 5, тем выше чувствительность. Таким образом, основными параметрами исследуемого процесса являются чувствительность 8 и электрохимический импульс
J 1 и \\ 1 и 1 1 1 )f 1 А V 1 X Е |\\ 1 ^ А i i
> 1 i с ь1 V \\ А V 1 A i f ^ В \\ 1 t [ 1 b t \\
Рис. 7. Фазовый портрет динамической системы (1).
АЪ, а основными характеристиками процесса могут быть частота к и амплитуда Ах мышечных сокращений, амплитуда изменения электрохимического потенциала Аъ, а также форма этих колебаний.
Главный интерес представляет адекватность данной модели особенностям реального процесса. С этой точки зрения основное внимание следует уделять относительным значениям величин, характеризующих исследуемый процесс; абсолютным значениям можно не придавать конкретного смысла.
На рис. 2-4 представлены результаты численных исследований указанных зависимостей в пакете МаЙ1сас1. В расчетах использовались встроенная подпрограмма интегрирования дифференциальных уравнений гкйхес1, а также графические возможности пакета МаШсас!.
На рис. 2 представлены характерные формы сердечного цикла в зависимости от положения особой точки Е(х0,Ь0) при условии, что импульс АЬ достаточно велик - настолько, что обеспечивает вылет траектории цикла в начале (пунктирная линия) и далее за гребень волны кубической параболы
Ъ=х-х .
На рис. 2а, 26 особая точка Е является устойчивой (это нетрудно проверить после линеаризации системы (1) при х > 1/д/з ), изображающая точка, попадая в ¿-окрестность точки Е, получает импульс АЬ > 0 и попадает в положение А, далее по циклу ВСБЕ. Циклы 2а, 26 отличаются участками БЕ, где происходит медленное сокращение 2а или расслабление 26 мышечного волокна, а также сравнительно медленный рост 2а или падение 2б электрохимического потенциала.
На рис. 2в особая точка Е является неустойчивой (|х| < 1/д/з ). В этом случае изображающая
точка, покидая окрестность точки Е, далее перемещается по циклу АВСВ, без дополнительных импульсов, за счет размеренного периодического изменения электрохимического потенциала. В случае
окрестность точки Е, за счет импульса АЬ > 0 выбрасывается в положение А, далее образуя неполный цикл АВЕ. Такой цикл естественно ассоциировать с понятием «конвульсия».
■0.6 -0.3 0 0.3 0.6 -г 0.8 -0.4 0 0.4 01 я) х=1; ¿=0.1; ¿Ъ=0.435 6) х =1.25; ¿=0.1; = 1.188
■0.6 -0.3 0 0.3 0.6 -0.6 -0.3 0 0.3 0.6 в) х =0; £=0.1; ¿£=0.485 г) х=0.62; ¿=0.1; ^?=0.867 Рис. 2. Влияние положения особой точки на форму сердечного цикла.
Явление «конвульсии» можно наблюдать также в окрестности устойчивой особой точки Е (рис.За), когда порог чувствительности достаточно высок (при малых 8 > 0), а электрохимический импульс ЛЬ мал. Если порог чувствительности достаточно низок (при больших 8 > 0), а электрохимический импульс АЬ также мал, то в окрестности устойчивой особой точки Е могут возникнуть дополнительные систолы типа А\\В\\Е\\ (рис.36) при полном сердечном цикле ЕАВСБ.
-0.6 -0.3 0 0.3 0.6 0.6 -0.3 0 0.3 0.6
а) х. =1.1; [5= 0.1; Л=0.51б б) X, =1; ¿5=0.1; ^=0.365
Рис. 3. Явление «конвульсии» и дополнительной систолы в окрестности
устойчивой особой точки.
Важными характеристиками процесса сердцебиения являются также частота (пульс), амплитуда колебаний мышечного волокна. В результате обработки семейства решений системы (1) на рис. 4а-б представлены зависимости частоты к и амплитуды Ах от порога чувствительности а на рис. 4в - амплитудо-частотная характеристика сердцебиения. На основании анализа данных можно сделать вывод, что с понижением порога чувствительности (8 возрастает) возрастает частота и уменьшается амплитуда колебаний мышечного волокна. Это согласуется с тем, что в основной массе у пожилых людей с возрастом повышается в среднем давление и слабеет сердцебиение.
х =1; ДЬ= 0.485
О 35 34 33 32 31
Рис. 4. Графики частоты и амплитуды колебаний мышечного волокна.
При передаче нервного импульса основную роль играет аксон - часть нервной клетки, химические процессы в которой вызывают формирование электрохимического импульса. Параметрами процесса являются показатели натриевой (-х) и калиевой (а) проницаемости поверхности аксона, а также разность потенциалов Ъ на внешней и внутренней сторонах клеточной мембраны. При увеличении разности потенциалов на АЬ показатель натриевой проницаемости (-лг) резко возрастает (х убывает), формируется импульс, затем изменение калиевой проницаемости а сравнительно медленно возвращает аксон в нейтральное состояние.
Такой процесс описывается трехмерной моделью Зимана:
£х = -(х3 + ах + Ь), а- -2х - 2а, Ь - -а -1, (2)
где в > 0 - малая постоянная.
Система (2) имеет единственную устойчивую особую точку (х = 1, а--\\ ¿ = 0), в окрестности которой аксон приобретает спокойное состояние. Передача нервного импульса к сердцу также носит
циклический характер и может быть изображена фазовой траекторией в трехмерном фазовом пространстве (х,а,Ь) (рис. 5). Фазовая траектория изображающей точки состоит: из участка ЕА, на котором происходит резкое изменение разности потенциалов ЛЬ; из участка АВ, на котором происходит резкое изменение показателя натриевой проницаемости (-х); из участка ВСЕ, на котором происходит медленное изменение показателя калиевой проницаемости а, а также восстановление всех характеристик аксона в спокойном состоянии.
Как и для модели сердцебиения (1), здесь можно представить зависимости основных характеристик (частота к, амплитуды изменений Ах, Аа, Аь, формы) от чувствительности 5, амплитудо-частотных характеристик, однако разумное толкование полученных результатов должно опираться на глубокое понимание электрохимических процессов.
На рис. 6 демонстрируются некоторые из указанных зависимостей.
Рис. 6. Зависимости основных характеристик передачи нервного импульса.
Резкие скачки на графиках пока еще не получили своего объяснения, - возможно, здесь влияет специфика вычислительного процесса.
Рис. 5. Фазовая траектория цикла нервного импульса.