CC BY
39
ISSN 2074-1863 Уфимский математический журнал. Том 3. № 4 (2011). С. 39-42.
УДК 517.938
ПРИБЛИЖЕННОЕ ИССЛЕДОВАНИЕ МНОГОПАРАМЕТРИЧЕСКИХ БИФУРКАЦИЙ В МОДЕЛЯХ ПОПУЛЯЦИОННОЙ ДИНАМИКИ
A.A. ВЫШИНСКИЙ
Аннотация. В работе приводится новая общая схема качественного и приближенного исследования основных сценариев локальных бифуркаций в моделях популяционной динамики. Приводятся основные и достаточные признаки бифуркаций состояний равновесия и Андронова-Хопфа для систем популяционной динамики.
Ключевые слова: прогиб пластины, критические силы, точки бифуркации, асимптотические формулы, состояние равновесия.
1. Введение
Одной из наиболее общих в математической биологии является [1] модель популяционной динамики (или модель Колмогорова), описываемая системой дифференциальных уравнений
(х\ = х1дх(х,^),
•К2 ^'2*72 (•£> АО >
Хп %п9п(,Х, ¡1) ■
В этой системе переменные Х\, Х2, • • • , хп означают численность каждой отдельной биологической популяции, X = {х\, Х2, хп) € Ип, ¡1 € Ик. Функции дг(х,/1) являются гладкими по ж € К+ = {х : ^ 0} и ¡1 € Ик, при этом частные производные ^жДж,/х)
при г ф ] сохраняют знак в октанте К+.
Система (1) зависит от скалярного или векторного параметра ¡л, при изменении которого возможны различные бифуркации. Особый интерес представляют так называемые локальные бифуркации [2], происходящие в окрестностях точек равновесия или циклов системы (1). Одним из основных здесь является вопрос об эволюции бифурцирующих решений при изменении параметра: важно, чтобы эти решения оставались в первом октанте К+. Изучению таких вопросов посвящены работы многих авторов (см. [1] и имеющуюся там библиографию). При этом, как правило, изучаются конкретные модели, исследуются однопараметрические бифуркации, а компьютерное моделирование проводится на основе прямого численного расчета.
В настоящей работе приводится схема качественного и приближенного исследования задач о многопараметрической бифуркации для систем вида (1). Схема основана на операторном методе, изложенном в [3].
A.A. Vyshinskiy, Approximate study of multiparameter bifurcations in models of population
DYNAMICS.
© Вышинский A.A. 2011.
Поступила 15 октября 2011 г.
2. Схема исследования
Предлагается следующая схема исследования задач о многопараметрической бифурка-
Локальные бифуркации системы (1) возможны в окрестностях неподвижных точек (состояний равновесия) и циклов (периодических решений). В этой связи отметим, что система (1) имеет нулевую точку равновесия х = 0, при этом она может иметь неподвижные точки и на координатных плоскостях.
Пусть х* — неподвижная точка системы (1). Вектор х* может иметь несколько нулевых компонент. Пусть для определенности неподвижная точка системы (1) имеет вид х* = (х^,..., х*т) 0,..., 0). Полагая 1г = х — х*, перейдем от (1) к системе
где А(ц) — матрица Якоби правой части системы (1), вычисленная в точке х*, а(/г,,/х) = о(||/г,||) при ||/г,|| —>• 0. Матрица А(ц) здесь имеет вид
Необходимым условием локальной бифуркации в окрестности точки х* является (см., например, [2]) требование, чтобы матрица А(ц) при некотором значении ц = цо имела чисто мнимое собственное значение. Если при этом матрица А(цо) имеет нулевое собственное значение, то значение ц = Цо обычно является точкой бифуркации неподвижных точек системы (1): при близких к ц0 значениях ц у системы (1) в окрестности состояния равновесия х* появляется новое состояние равновесия х**. Если же матрица А(ц0) имеет пару собственных значений вида ±го;о при о;0 > 0, то значение ц = Цо обычно является точкой бифуркации Андронова-Хопфа системы (1): при близких к цо значениях ц у системы (1) в окрестности состояния равновесия х* появляется нестационарное периодичесое решение малой амплитуды.
Для анализа соответствующих бифуркаций в системе (1) воспользуемся операторном методом, изложенным в [3]. Приведем некоторые признаки бифуркаций, при этом ограничимся случаем п = 2.
Теорема 1. Пусть в системе (1) п = 2 и, х* —ее неподвижная точка. Пусть ¡1 = 110 являет,ся точкой бифуркации неподвижных т,очек. Тогда д\(х*, ¡1о) = 0 или д2(х*, ¡1 о) = 0. При этом, если выполнено только одно из этих равенств, то бифуркация имеет, коразмерности 1, если же выполнены оба равенства, то коразмерности бифуркации, равна 2.
Это утверждение дает необходимый признак бифуркации неподвижных точек. Следующая теорема дает достаточный признак бифуркации коразмерности 1.
Теорема 2. Пусть выполнено в точности одно из соот,ношений д\(х*, ¡10) = 0 или, д2(х*, ¡1о) = 0. Пусть д1Х1^(х*, ¡1о) ф 0 или, д^^х*, ¡1о) ф 0 соответственно. Тогда ¡1 = 110 являет,ся точкой бифуркации, неподвижных точек системы (1).
Аналогичное утверждение может быть получено и для бифуркации коразмерности 2. Справедливость этих утверждений и других приводимых в работе теорем следует из более общих результатов, полученных в [3].
Бифуркация Андронова-Хопфа в системах вида (1) возможна лишь в окрестности неподвижной точки с двумя или большим числом ненулевых координат. В частности,
ции для системы (1). Ограничимся рассмотрением случаев х Е В2 и х Е Д3.
к' = А(ц)їі + а(/і, ¡і)
х*ід'іХі{х\ ¡1) Хід'2хі(х*, ¡1,)
хід'пх^х*, и)
9т{% і АО
0
о
9п(х*,ц)
она невозможна для двумерной системы в окрестности точки вида х* = (0, £2) или х* = (ж1,0). Приведем достаточный признак бифуркации Андронова-Хопфа для системы (1) при п = 2.
Теорема 3. Пусть х* = (а^,^) — неподвижная точка системы (1) с положительными компонентами. Значение являет,ся точкой бифуркации, Андронова-Хопфа для этой системы, если выполнены условия:
91x192x2^ > А^о) 0, 92x1^ ) ^0)91x2^ > А^о) ^ 0.
Аналогичное утверждение может быть получено и для систем более высокого порядка.
Приближенное исследование бифуркационных процессов в системах вида (1) предлагается проводить по схеме, описанной в [3]. А именно, на первом этапе необходимо перейти к операторному уравнению, содержащему основные параметры задачи. При этом для исследования сценария бифуркации необходимо выбрать число параметров, соответствующих коразмерности бифуркации. На втором этапе к полученному операторному уравнению применяется метод функционализации параметра. Наконец, для исследования функцио-нализированного уравнения используется метод Ньютона-Канторовича.
3. Система “хищник — две жертвы”
В качестве иллюстрации рассмотрим модель “хищник-две жертвы” [1], описываемую уравнениями
{и[ = щ{а\ — щ — 6 и2 — 4г>),
и2 = и2{о12 — Щ — щ — Юг»), (3)
V' = —v{l — 0, 25«! — 4м2 + у).
Здесь П\ И «2 — численность популяций каждой из жертв, V — численность популяции хищников.
Система (3) содержит два параметра а\ и а^. Рассмотрим бифуркацию Андронова-Хопфа в окрестности ненулевой неподвижной точке с координатами х* = (—11,2 + 8, 2а\ — 4,4ск2; 1,2 — 0, 7а\ + 0,4а2; 1 — 0,75а\ + 0,5о!2). Коразмерность этой бифуркации равна одному, поэтому в качестве бифуркационного рассмотрим один параметр, а именно, а2-Пусть при этом а\ = 3, 7. Циклы в этой системе возникают при изменении параметра (У.2• При 012 ~ 3,874 матрица Якоби правой части системы (3) имеет два чисто мнимых собственных значения.
На рис. 1 продемонстрированы результаты численного исследования бифуркации в системе (3), полученные по разработанной автором программе. На этом рис. изображены проекции фазовых траекторий системы (3) на плоскость (щ, 112) при различных значениях параметра си 2 • Символ * используется для обозначения начальных точек соответствующих периодических решений.
Рис. 1. Семейство периодических траекторий системы (3)
СПИСОК ЛИТЕРАТУРЫ
1. Базыкин А.Д. Нелинейная динамика взаимодействующих популяций. М.-Ижевск: Институт компьютерных исследований, 2003. 368 с.
2. Арнольд В.И. Геометрические методы в теории обыкновенных дифференциальных уравнений. М.: НИЦ «Регулярная и хаотическая динамика», 2000. 400 с.
3. Вышинский A.A., Ибрагимова Л.С., Муртазина С.А., Юмагулов М.Г Операторный мет,од приближенного исследования правильной бифуркации в многопщхш,етрических динамических системах // Уфимский математический журнал, 2010. Т.2. №4. С. 3-26.
Александр Алексеевич Вышинский,
Сибайский институт (филиал) БашГУ, ул. Белова, 21,
453837, г. Сибай, Россия E-mail: sanek3484@gmail. com
