Спросить
Войти
Категория: Математика

МОДЕЛИРОВАНИЕ РАБОТЫ СЕРДЦА И НЕРВНЫХ ИМПУЛЬСОВ

Автор: Сельвинский В.В.

зультатам появится возможность проанализировать перспективы дальнейшего развития кластерного анализа. В работе осуществлены все этапы обработки больших данных текстовой информации, начиная с этапа предварительной обработки и до момента графического отображения, что также ново.

Практическая значимость результатов работы актуальна из-за стремительного развития электронного обучения во всем мире. Определение наиболее подходящих методов анализа, фреймворков, стеммера и других инструментов позволит обработать эти данные быстрее и уменьшить число ошибок.

Применение технологии больших данных даст возможность педагогам и сотрудникам образовательных учреждений получать актуальную информацию об обучающихся и реагировать на изменения процесса обучения, а в результате - прогнозировать популярность определенного вида деятельности, вносить изменения в учебный контент.

1. Благирев, A. Big data простым языком. - М.: ACT, 2019.
2. Изофатов, К.А. Кластерный и интеллектуальный анализ текстовой информации. Основные понятия и проблемы. - http://econf.rae.ru/pdf/201 l/03/120.pdf
3. Корнев, М.С. История понятия «Большие данные»(Е^ data): словари, научная и деловая периодика. -https://cyberleninka.ru/article/n/istoriya-ponyatiya-bolshie-damye-big-data-slovari-nauchnayaviewer
4. Мамедова, Г.А, Зайналова, JT.A, Меликова, Р.Т. Технологии больших данных в электронном образовании. -https://cyberleninka.rU/article/n/tehnologii-bolshih-damyh-v-elektronnom-obraz
5. Обзор алгоритмов кластеризации данных, -https://habr.com/ru/post/101338/
6. Тюрин, А.Г, Зуев, И.О. Кластерный анализ, методы и алгоритмы кластеризации. - https://rtj.mirea.ru/ upload/medialibrary/fba/09-tyurin.pdf

УДК 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/д/з ). В этом случае изображающая

точка, покидая окрестность точки Е, далее перемещается по циклу АВСВ, без дополнительных импульсов, за счет размеренного периодического изменения электрохимического потенциала. В случае

2г особая точка Е является устойчивой изображающая точка, периодически попадая в

окрестность точки Е, за счет импульса АЬ > 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) при полном сердечном цикле ЕАВСБ.

1.51"

-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

0
75 150 225 300
0 0.02 0.04 0.06 0.08 0.1 0 0.02 0.04 0.06 0.08 0.1

Рис. 4. Графики частоты и амплитуды колебаний мышечного волокна.

При передаче нервного импульса основную роль играет аксон - часть нервной клетки, химические процессы в которой вызывают формирование электрохимического импульса. Параметрами процесса являются показатели натриевой (-х) и калиевой (а) проницаемости поверхности аксона, а также разность потенциалов Ъ на внешней и внутренней сторонах клеточной мембраны. При увеличении разности потенциалов на АЬ показатель натриевой проницаемости (-лг) резко возрастает (х убывает), формируется импульс, затем изменение калиевой проницаемости а сравнительно медленно возвращает аксон в нейтральное состояние.

Такой процесс описывается трехмерной моделью Зимана:

£х = -(х3 + ах + Ь), а- -2х - 2а, Ь - -а -1, (2)

где в > 0 - малая постоянная.

Система (2) имеет единственную устойчивую особую точку (х = 1, а--\\ ¿ = 0), в окрестности которой аксон приобретает спокойное состояние. Передача нервного импульса к сердцу также носит

циклический характер и может быть изображена фазовой траекторией в трехмерном фазовом пространстве (х,а,Ь) (рис. 5). Фазовая траектория изображающей точки состоит: из участка ЕА, на котором происходит резкое изменение разности потенциалов ЛЬ; из участка АВ, на котором происходит резкое изменение показателя натриевой проницаемости (-х); из участка ВСЕ, на котором происходит медленное изменение показателя калиевой проницаемости а, а также восстановление всех характеристик аксона в спокойном состоянии.

Как и для модели сердцебиения (1), здесь можно представить зависимости основных характеристик (частота к, амплитуды изменений Ах, Аа, Аь, формы) от чувствительности 5, амплитудо-частотных характеристик, однако разумное толкование полученных результатов должно опираться на глубокое понимание электрохимических процессов.

На рис. 6 демонстрируются некоторые из указанных зависимостей.

501-— 2 461-— 2 4б1-0 0.02 0.04 0.06 0.08 0.1 & 0 0.02 0.04 0.06 0.08 0.1 & 50 100 150 200

Рис. 6. Зависимости основных характеристик передачи нервного импульса.

Резкие скачки на графиках пока еще не получили своего объяснения, - возможно, здесь влияет специфика вычислительного процесса.

1. Эрроусмит, Д., Плейс, К. Обыкновенные дифференциальные уравнения. Качественная теория с приложениями / пер. с англ. - М.: Мир, 1986. - 243 с.

Рис. 5. Фазовая траектория цикла нервного импульса.

БИЕНИЕ СЕРДЦА НЕРВНЫЙ ИМПУЛЬС МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ mathcad heart beat nerve impulse mathematical model differential equations
Другие работы в данной теме:
Контакты
Обратная связь
support@uchimsya.com
Учимся
Общая информация
Разделы
Тесты