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

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

 

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

 

2

(1 a1u2

a2u12 a3u24 )

 

du

2

 

2u2 k3u1 k4u3 0,

(1.8)

 

 

 

 

 

 

 

 

 

 

 

 

 

d

d 2

 

 

 

 

 

 

 

 

 

d 2u

3

 

3

(1 a1u3

a2u32 a3u34 )

du

3

 

3u3 k5u1 k6u2 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

d 2

 

 

 

 

 

 

 

 

 

 

 

 

где

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

 

 

 

колебания 1-го,

2-го и 3-го

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

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

между

автогенераторами;

1 , 2 , 3

— отношение собственной частоты i-го автогенератора в автономном

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

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

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

d

2

u

 

 

 

du

 

 

 

 

M1

 

 

 

 

 

1Ô1

 

1u1

kiui 0,

 

 

 

 

 

1

 

(u1 )

 

1

 

 

 

 

 

2

 

 

 

 

 

d

 

 

 

 

d

 

 

 

 

i 2

 

 

 

d

2

u

 

 

 

 

 

du

 

 

 

M 2

 

 

 

 

 

2Ô21 (u2 )

 

2

 

2u2

kiui

0,

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

(1.9)

d

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

0.2

 

k

 

a

1

 

 

 

 

 

 

 

0.1

 

 

 

 

 

 

 

1.05

 

 

 

 

0.2

 

 

 

0.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.01

 

 

 

 

 

 

 

 

 

 

0.1

 

 

 

 

 

 

 

 

 

 

 

1.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.15

 

 

0.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.15

 

0.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y1

2

a3 y1

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

a1 y1 a2

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

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

25

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

ORIGIN

 

 

 

 

 

 

 

 

 

F)

 

 

U3V1

 

 

 

Z

 

 

 

 

 

 

V3

 

Z

 

 

 

 

 

 

 

 

 

 

Z rkfixedy(

 

00U1200 10001Z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t:=Z<0>

U2 Z 3

4

 

V2 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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

1

 

 

 

 

 

 

 

0.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

0.2

 

k

 

0.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1.15

 

 

 

 

 

 

 

 

0.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.2

 

 

 

 

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), напоминающий электрокардиограмму при аритмии сердца.

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

Источник: https://studfile.net/preview/16724881/