111
нутой кривой, называемой предельным циклом. Таким образом, в системе
устанавливается автоколебательный процесс.
При
c
BB <
положение равновесия устойчиво, и система ведёт себя со-
вершенно иначе. Концентрации X и Y совершают затухающие колебания,
приближаясь к своим равновесным значениям: при
AX
,
ABY /
. Изображающая точка на фазовой плоскости в этом случае асим-
птотически стремится к положению равновесия, т.е. к точке с координа-
тами (A;B/A). При этом возможны два случая. Если
22
,
точка равновесия (A;B/A) является устойчивым фокусом, а концентрации
X и Y совершают затухающие колебания около своих равновесных значе-
ний. Если же
2
, то точка равновесия является устойчивым уз-
лом, и затухающие колебания концентраций X и Y носят апериодический
характер. Сравним теперь результаты решения системы (1) при различных
значениях B, полагая
. В этом случае
5
=
c
B
. Решение проводилось с
помощью пакета Mathcad методом Рунге – Кутта четвёртого порядка.
На рис. 1 представлен результат решения системы (1) при
;
6
B
с начальными условиями
8,1)0(
X
;
3,3)0(
Y
. Поскольку в данном случае
c
BB >
, то, как следует из вышесказанного, образуется предельный цикл.
Это видно из верхнего графика рис. 1. Этот график представляет собой
фазовую траекторию. Светлой точкой на фазовой плоскости изображено
положение равновесия, т.е. точка с координатами (2; 3). Кроме того, по-
скольку выполнено условие
2
, точка равновесия (2; 3) пред-
ставляет собой неустойчивый фокус. Действительно, на фазовой траекто-
рии видна раскручивающаяся спираль. На втором и третьем графиках рис.
1 представлены соответственно зависимости концентраций X и Y от вре-
мени. Видно, что амплитуды колебаний этих концентраций возрастают,
приближаясь к своим предельным значениям. На рис. 2 представлен ре-
зультат решения системы (1) при тех же значениях A и B, что и на рис. 1
(
;
6
B
), но с новыми начальными условиями:
3,1)0(
X
;
9,1)0(
Y
.
Параметры A и B в данном случае имеют такие же значения, что и на рис.
1, поэтому образуется такой же предельный цикл. Это иллюстрирует тот
факт, что образующийся предельный цикл не зависит от начальных усло-
вий, а определяется лишь параметрами системы. Анализ двух последую-
щих графиков рис. 1 и 2 показывает, что и периоды установившихся ко-
лебаний концентраций X и Y на обоих графиках также одинаковы. На рис.
3 представлен результат решения системы (1) при
;
10
B
с началь-
ными условиями
8,1)0(
X
;
2,5)0(
Y
. Поскольку в данном случае
c
BB >
и
2
, то образуется предельный цикл, и точка равновесия (2; 5)
является уже неустойчивым узлом. На рис. 4 представлен результат ре-
шения системы (1) при тех же значениях A и B, что и на рис. 3 (
;
10
B
), но с новыми начальными условиями:
8,1)0(
X
;
3)0(
Y
. Посколь-
ку в этом случае параметры A и B имеют такие же значения, что и на рис.
3, то образуется такой же предельный цикл.