|
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), напоминающий электрокардиограмму при аритмии сердца.
Таким образом, с помощью компьютерного анализа подтверждается идея Ван-дер-Поля и Ван-дер-Марка о возможности моделирования работы сердца с помощью трех связанных между собой автогенераторов.