СЕМИНАР 9
Занятие в компьютерном классе. Колебательные
системы. Локальная модель брюсселятора. Построение
фазовых портретов при разных значениях параметров.
КОЛЕБАТЕЛЬНЫЕ СИСТЕМЫ
Для биологических систем характерно периодическое изменение различных характеристик. С некоторыми
из типов периодических изменений мы уже имели дело
при рассмотрении особых точек типа центр, фокус. Однако часто в живых системах наблюдаются колебания, обладающие особым отличительным свойством: неизменностью во времени периода и амплитуды колебаний, означающей стационарность и устойчивость колебательного режима. В данном случае периодическое изменение величин представляет собой один из типов стационарного поведения системы. Если колебания имеют постоянную период и амплитуду, устанавливаются независимо от начальных условий и поддерживаются свойствами самой системы, без воздействия периодической силы, система называется автоколебательной. На фазовой плоскости такому типу поведения соответствует устойчивый предельный цикл — притягивающее множество. Предельный цикл есть изолированная замкнутая кривая.
Предельные циклы могут быть устойчивыми и неустойчивыми. Предельный цикл устойчив, если существует такая область на фазовой плоскости, содержащая этот предельный цикл — окрестность, — что все фазовые траектории, начинающиеся в окрестности, асмиптотически при t приближаются к предельному циклу.
Если же наоборот, в любой сколь угодно малой окрестности предельного цикла существует по крайней мере одУчебное пособие «Математические модели в биологии»
на фазовая траектория, не приближающаяся к предельному циклу при t, то такой предельный цикл называется неустойчивым. Такие циклы разделяют области влияния разных притягивающих множеств.
Внутри устойчивого предельного цикла обязательно есть неустойчивая точка типа фокус или узел. Необходимым условием рождения устойчивого предельного цикла является переход (бифуркация) типа стационарного состояния от устойчивого фокуса к неустойчивому (действительная часть соответствующих характеристических чисел Re 1,2 через 0 переходит от отрицательных значений к положительным, при этом мнимая часть характеристических чисел Im 1,2 отрицательна). После такой бифуркации возможность существования устойчивого предельного цикла сохраняется до тех пор, пока Re 1,2 сохраняют положительное значение.
ЛОКАЛЬНАЯ МОДЕЛЬ БРЮССЕЛЯТОРА
Брюсселятор — базовая модель, являющаяся классическим примером автоколебательного поведения переменных (концентраций) в системе химических реакций.Брюсселятор представляет собой следующую схему гипотетических реакций:
k k1 k2 k 2X + Y B+ X Y +C, A X, 3X, X R.
k1 k2 k3 k Здесь А, В — исходные вещества, C, R — продукты, X, Y — промежуточные вещества. Пусть конечные продукты C и R немедленно удаляются из реакционного пространства. Это означает, что обратные константы k3 = k4 = 0.
Если субстрат A находится в избытке, то k1 = 0. Предполагают также, что k2 = 0. Значения остальных констант приравнивают 1. Тогда рассматриваемая схема реакций описывается системой уравнений:
Семинар 9. Колебательные системы dx dt = A x( B + 1) + x y, dy = x( B xy ), (9.1) dt A, B > 0, x, y > 0.
Исследуем систему уравнений (9.1). Найдем стационарные значения:
A x ( B + 1) + x 2 y = 0, A x ( B + 1) + x 2 y = 0, x = 0, x = B, x ( B xy ) = 0, y A, B > 0, A, B > 0, x, y > 0, x, y > 0, B B A ( B + 1) + y = 0, B y = A, y y x = B, x = A, A, B > 0, y A, B > 0, x, y > 0.
x, y > 0, Исследуем тип стационарной точки и ее устойчивость. Найдем коэффициенты линеаризации:
Учебное пособие «Математические модели в биологии»
1) Определим знак Px + Q для ответа на вопрос об устойy чивости:
При B > 1 + A2 выполнено неравенство Px + Q > 0, значит, в этой области параметров стационарная точка неустойчива.
При B 1 + A2 выполняется неравенство Px + Q 0, таy ким образом, в этой области параметров стационарная точка устойчива.
ществует.
3) Фокус или узел? Надо определить знак подкоренного выражения в формуле для вычисления характеристических чисел 1,2 :
Это выражение отрицательно при выполнении неравенства (1 A) 2 < B < (1 + A) 2. Таким образом, при выполнении неравенства (1 A) 2 < B < (1 + A) 2 в системе должны наблюдаться колебания.
При B (1 A) 2 и B (1 + A) 2 колебаний в системе не будет.
При Px + Q = 0 происходит смена устойчивости стациоy нарного состояния типа фокус, в этот момент происходит рождение предельного цикла.
Семинар 9. Колебательные системы в виде параметрической диаграммы:
Стационарная точка устойчива, Стационарная точка неустойчива, предельный цикл не существует Учебное пособие «Математические модели в биологии»
ЗАДАЧИ К СЕМИНАРУ
9.1. Определите тип стационарного состояния системы в зависимости от значения параметра a (a > 0). Укажите интервал значений параметра, при котором в системе имеют место автоколебания. Ответ оформите в виде параметрической диаграммы (прямой).