Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Методические указания к лабораторным работам 2,3,4,5 «Математическое моделирование биологических систем» (90

..pdf
Скачиваний:
3
Добавлен:
15.11.2022
Размер:
732.01 Кб
Скачать

Лабораторная работа №5 « Моделирование аритмии сердца». Цель работы: закрепление знаний о функционировании

биологических систем, полученных на лекциях по дисциплине «Моделирование биологических процессов и систем», приобретение навыков работы с программным пакетом MATHCAD.

В 1928 г. голландские ученые Ван-дер-Поль и Ван-дер-Марк разработали динамическую модель сердца в виде трех связанных между собой автогенераторов. С помощью созданного устройства ученые пытались моделировать некоторые болезни сердца, в том числе и аритмию. Повторим такое исследование с помощью компьютерного моделирования трех связанных автогенераторов А-1, А-2, А-З (рис.8).

 

k2

 

k1

A-1

 

 

 

 

 

k5

 

 

k3

 

 

 

 

 

 

 

 

A-3

 

 

 

A-2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k4

k6

 

 

 

 

Рис.8. Структурная схема взаимодействия трех автогенераторов

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

Составим систему дифференциальных уравнений для трех связанных между собой автогенераторов. Такая система как частный случай уравнений (1.9) имеет вид:

d 2u

 

1 (1 a1u1

a2u12 a3u14 )

du

 

1u1 k1u2 k2u3 0,

 

 

 

 

1

 

 

 

 

1

 

 

 

2

 

d

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d 2u

 

2

(1 a1u2

a2u12 a3u24 )

 

du

2

 

2u2 k3u1

k4u3

0,

(1.8)

 

 

 

2

 

 

 

 

 

 

 

 

2

 

 

d

 

d

 

 

 

 

 

 

 

 

 

 

 

 

d 2u

 

 

 

a2u32 a3u34 )

 

du

3

 

 

 

 

 

 

 

 

3

 

3

(1 a1u3

 

 

 

 

 

3u3 k5u1 k6u2 0,

 

 

 

2

 

 

 

d

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

u1 (t), u2 (t), u3 (t)

 

 

 

 

колебания

1-го, 2-го

и

3-го

автогенераторов; k1 k6

—коэффициенты связи между автогенераторами;

1 , 2 , 3

— отношение собственной частоты

i-го

автогенератора в

автономном режиме к частоте автоколебаний

в системе;

t

нормированное время.

 

 

 

 

 

 

 

 

 

 

 

 

 

Система дифференциальных уравнений для i автогенераторов:

d

2

u

 

 

 

du

 

 

 

 

M1

 

 

 

 

1Ô1 (u1 )

1u1

kiui 0,

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

2

 

 

 

 

 

 

d

 

 

 

 

d

 

 

 

 

i 2

 

 

d

2

u

 

 

 

 

 

du

 

 

 

M 2

 

 

 

 

2Ô

 

 

 

2

 

2u2

kiui 0,

 

 

 

 

 

 

1

 

21 (u2 )

 

 

 

 

 

 

 

 

 

 

2

 

d

(1.9)

 

d

 

 

 

 

 

 

 

i 1

 

......................................................................,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

2

u

 

 

 

 

du

 

 

 

 

M n

 

 

 

 

nÔn (un )

n

 

nun

kiui 0,

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

2

 

d

 

 

 

 

d

 

 

 

 

 

 

 

 

 

i 1

 

 

где Ô i (ui ) — нелинейная функция;

k i — коэффициент связи

i-го

автогенератора

с другими;

i — отношение собственной частоты

i-го

автогенератора в автономном режиме работы к частоте автоколебаний в системе.

Программа расчета режима синхронизации автогенераторов согласно (1.9) приведена на рис.9.

ORIGIN 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.2

 

0.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.2

 

 

 

0.1

 

 

0.5

 

 

 

 

 

 

 

0.1

 

 

 

 

 

 

 

 

 

 

 

1.5

 

 

 

 

 

 

 

 

 

 

0.1

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

0.2

 

k

 

 

1

 

 

 

 

 

 

0.1

 

 

 

 

 

 

 

 

1.05

 

 

 

 

0.2

 

 

0.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.01

 

 

 

 

 

 

 

 

0.1

 

 

 

 

 

 

 

 

 

 

 

1.1

 

 

 

 

 

 

 

0.15

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.15

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a1 y1 a2 y1

2

a3 y1

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

y2

y1 1 k1 y3

k2 y5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F(t y)

 

 

 

 

 

 

 

 

 

 

 

 

y

 

2

 

 

 

y

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

1

a

1

y

3

a

2

3

a

3

 

3

y

4

y

3

2 k y

1

k y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

4 5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

2

 

 

 

 

 

5

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

a

 

y

 

a

 

 

y

a

 

 

y

y

 

y

 

 

 

k y

 

k y

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

3

 

3

 

 

 

 

 

 

 

 

1

 

 

5

 

 

2

 

 

 

 

3

 

 

 

6

 

5

 

 

 

5

1

6

 

ORIGIN 0

Z rkfixedy( 0 200 10001 F)

 

 

t:=Z<0>

U1 Z 1

V1 Z 2

 

3

 

5

6

U2 Z

V2 Z 4

U3 Z

V3 Z

Рис.9. Программа расчета колебаний модели сердца (пример 1)

Впрограмме каждое из уравнений (1.9) второго порядка представлено

ввиде двух уравнений 1-го порядка, Ul(t), U2(t), U3(t) — функции колебаний u3(t); u2(t), u3(t); VI, V2, V3 — их производные.

Рис.10. Результаты расчета по программе (пример 1)

ORIGIN 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.2

 

 

0.1

 

0.5

 

 

 

 

 

 

0.1

 

 

 

 

 

 

 

 

 

 

1.5

 

 

 

 

 

 

 

0.2

 

 

 

 

0.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

0.2

 

k

 

0.1

a

1

 

 

 

 

0.1

 

 

 

 

 

 

1.15

 

 

 

 

0.2

 

 

 

0.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.01

 

 

 

 

 

 

0.1

 

 

 

 

 

 

 

 

 

 

1.3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.15

 

 

 

0.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.15

 

 

0.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a1 y1 a2 y1

2

a3 y1

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

y2

y1 1 k1 y3

k2 y5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F(t y)

 

 

 

 

 

 

 

 

 

 

y

 

2

 

 

 

y

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

1

a

1

y

3

a

2

3

a

3

 

3

y

4

y

3

2 k y

1

k y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

4 5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

2

 

 

 

 

 

5

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

a

 

y

 

a

 

 

y

a

 

 

y

y

 

y

 

 

 

k y

 

k y

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

3

 

3

 

 

 

 

 

 

 

1

 

5

 

2

 

 

 

 

3

 

 

 

6

 

5

 

 

 

5

1

6

 

 

Рис.11. Программа расчета колебаний модели сердца (пример 2)

Результаты расчета по программе двух примеров в виде зависимости переходного процесса в автогенераторах U1(t), U2(t), U3(t) и фазовых траекторий VI (Ul), V2(U2), V3(U3) приведены на рис.10 и 11. Исходные данные для первого примера приведены в начале программы рис. 9, для второго — перед графиками на рис. 11.

Из анализа полученных результатов следует, что при относительно небольшом (в пределах 10 %) расхождении собственных частот колебаний автогенераторов в системе устанавливается режим синхронизма, при котором частоты автоколебаний становятся равными, о чем свидетельствуют графики, на которых совмещены функции Ul(t), U2(t), U3(t) (см. рис. 10). При большем расхождении частот автоколебаний в пределах 15...30 % режим синхронизма нарушается и колебания автогенераторов приобретают вид биений (см. рис. 12), напоминающий электрокардиограмму при аритмии сердца.

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

Рис.12. Результаты расчета по программе (пример 2).

Содержание отчета.

Отчет по лабораторной работе должен содержать:

1) краткую теоретическую часть, включающую в себя дифференциальные уравнения для описания данной модели; 2) программу расчета модели;

3)таблицу, содержащую значения решения уравнений в рассчитываемых точках;

4)результат решения модели по программе в виде графиков; 5)выводы по проделанной работе.

Контрольные вопросы

1)Что такое математическое моделирование, цели и задачи?

2)Опишите общий алгоритм математического моделирования.

3)Что такое имитационное моделирование?

4)Расскажите основные принципы имитационного моделирования.

5)Перечислите основные функции математического пакета программ MathCAD, используемые для моделирования биологических систем.

6)Метод Рунге-Кутта.

7)Что такое фазовый портрет.

8)Динамическая модель сердца Ван-дер-Поля и Ван-дер-Марка.

9)Особенности моделирования аритмии сердца.

10)Какие механизмы лежат в основе большинства нарушений ритма сердца?

11)Чем определяется амплитуда и форма электрокардиографических комплексов при различной локализации электродов.

Библиографический список

1)Охорзин, В.А. Прикладная математика в системах MATHCAD: учебное пособие/В.А.Охорзин.-СПБ.:Лань,2009.-с.352.

2)Герман, И. Физика организма человека/ И.Герман.-М.:Интеллект,2011.-

с.992.

3)Дворецкий, С.И. Моделирование систем [ Текст ]/ С.И. Дворецкий, Ю.Л. Муромцев, В.А. Погонин, А.Г. Схиртладзе.-М.: Академия, 2009.-

с.320

4)Советов, Б.Я. Моделирование систем. Практикум / Б.Я. Советов. – М.: Высшая школа, 2009.-с.295.

Математическое моделирование биологических систем

Методические указания к лабораторным работам № 2,3,4,5

по дисциплине «Моделирование биологических процессов и систем»

Составители: Корчагина Вера Анатольевна

Батищева Юлия Николаевна

Редактор Е.А. Федюшина

Подписано в печать . Формат 60x84.1/16. Бумага офсетная. Ризография. Печ.л. 1.8. Тираж 50 экз. Заказ №

Липецкий государственный технический университет. 398600, Липецк, ул. Московская, 30.

Типография ЛГТУ. 398600, Липецк, ул. Московская, 30.

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ

Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования

«ЛИПЕЦКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ»

Кафедра физики и биомедицинской техники

Математическое моделирование биологических систем

Методические указания к лабораторным работам № 2,3,4,5

по дисциплине «Моделирование биологических процессов и систем»

Составители: В.А. Корчагина, Ю.Н. Батищева

Утверждаю к печати

Проректор по учебной работе

Объем

1,8 п.л.

Качановский Ю.П.

Тираж

50 экз.

«

 

»

 

2012г.

Липецк Липецкий государственный технический университет

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]