CC BY
46
НАУЧНОЕ ИЗДАНИЕ МГТУ ИМ. Н. Э. БАУМАНА
НАУКА и ОБРАЗОВАНИЕ
Эл № ФС77 - 48211. Государственная регистрация №0421200025. ISSN 1994-0408
электронный научно-технический журнал
Локализация инвариантных компактов двумерных
непрерывных динамических систем
# 07, июль 2013
Б01:10.7463/0713.0583104
Канатников А. Н., Федорова Ю. П.
УДК 517.925.5
Россия, МГТУ им. Н.Э. Баумана [email protected]
Введение
Один из методов качественного анализа динамической системы состоит в оценке положения ее инвариантных компактных множеств. Под инвариантным множеством понимают такое множество в фазовом пространстве системы, для которого любая траектория системы либо целиком принадлежит множеству, либо не пересекается с ним. Инвариантные множества с дополнительным условием компактности тесно связаны с ограниченными траекториями системы (положениями равновесия, предельными циклами, сепаратрисами) и такими их конгломератами как аттракторы и репеллеры.
Для оценки положения инвариантных компактных множеств можно использовать локализирующие множества, т.е. такие множества в фазовом пространстве системы, которые содержат все ее инвариантные компакты [1].
Один из методов построения локализирующих множеств основан на использовании функций, определенных на фазовом пространстве системы и называемых локализирующими [5]. Очевидно, что, рассматривая различные локализирующие функции, можно получить целое семейство локализирующих множеств, пересечение которых дает в итоге более сильный результат.
Использование локализирующих функций позволило исследовать ряд известных систем с хаотическим поведением [4, 6, 7, 8, 9,10,11], в частности, знаменитую систему Лоренца. Отметим, что все эти системы имели 3-й порядок или выше. Исследование систем 2-го порядка имеет особенности. С одной стороны, качественное поведение систем 2-го порядка заметно проще, чем поведение систем более высокого порядка, в частности, в таких системах явления динамического хаоса не наблюдаются. С другой стороны, попытки решения задач локализации для систем 2-го порядка приводили к практическим трудностям.
В настоящей статье исследуются две непрерывные динамические системы второго порядка, описывающие поведение некоторых биологических систем. Для каждой из рассмотренных систем строится семейство локализирующих множеств и вычисляется пересечение полученного семейства локализирующих множеств. При этом для первой системы решение удалось получить чисто аналитически, а для второй системы предложена численная процедура построения локализирующих множеств. Результаты исследования двух систем показаны на рисунках.
Статья организована следующим образом. В разд. 1 описан метод построения локализирующих множеств [1]. В разд. 2 рассмотрена динамическая система ICAISC, представляющая собой модель взаимодействия адаптивной иммунной системы с раковой опухолью [2]. В разд. 3 рассмотрена динамическая система VDFM, которая описывает двухкомпонентную экосистему [3]. В заключении подводятся итоги исследования.
1. Функциональный метод локализации
Основная идея функционального метода локализации в применении к автономной непрерывной динамической системе
ж = f (ж), ж е (1)
состоит в следующем.
Пусть ф: Rn ^ R — гладкая (хотя бы класса C*) функция, определенная на фазовом пространстве системы (эта функция далее называется локализирующей). Рассмотрим множество
= {ж е Rn: ф(ж) = 0} ,
где ф обозначает производную функции ф в силу системы (1). Множество Sv называется универсальным сечением [1]. Все инвариантные компактные множества системы (1), содержащиеся в множестве Q С Rn, содержатся в множестве
QV(Q) = {ж е Rn: (Q) < ф(ж) < Фзир^)} ,
где
ф-mf (Q) = inf ф(ж), фзир = SUp ф(ж).
x£SvnQ xGS^nQ
В качестве локализирующих можно использовать любые гладкие функции, определенные на фазовом пространстве системы, ограничений здесь нет. Однако неудачный выбор локализирующей функции может привести к тому, что локализирующее множество Q^ = совпадет с фазовым пространством, т.е. решение задачи локализации окажется тривиальным. Другая проблема, возникающая при решении задач локализации, — решение двух экстремальных задач, связанных с вычислением значений ф^ и ф8ир. В общем случае требуется предварительный анализ экстремальных задач. Решение этих задач упрощается, если
универсальное сечение является замкнутой кривой, тогда для вычисления значений ^¡^ и ^ир можно использовать метод множителей Лагранжа или численные методы.
Рассматривая несколько локализирующих функций, мы получаем несколько локализирующих множеств. В качестве наилучшей оценки положения инвариантных компактов следует выбрать пересечение найденных локализирующих множеств. Проведенные исследования ряда динамических систем [4,6,7,8,9,10,11] показали, что эффективным является использование параметрических семейств локализирующих функций с последующим построением пересечения найденного параметрического семейства локализирующих множеств.
2. Локализация инвариантных компактов системы 1СЛ18С
Система
х = ах — Ьху — сх2; , ч
о (2)
у = dy - еху - fу2
описывает динамику взаимодействия адаптивной иммунной системы с раковой опухолью [2]. Здесь х — концентрация раковых клеток; у — концентрация клеток иммунной системы; а, Ь, с, d, е, f — параметры системы, представляющие собой положительные числа.
Содержательный смысл фазовые переменные х и у имеют только в первой четверти фазовой плоскости (включая границу). Поэтому в качестве фазового пространства следует рассматривать множество
= {(х, у): х > 0, у > 0}.
Отметим, что это множество инвариантное, поскольку инвариантными являются прямые х = 0 и у = 0, ограничивающие первую четверть плоскости. Однако исследование системы (2) удобно проводить на всей фазовой плоскости, сужая затем полученные локализирующие множества на первую четверть.
Система (2) имеет четыре положения равновесия:
= (»■ 0); Р2 = (0,|); Р = (а, 0); Р. = (^'
При этом [2]:
Р — неустойчивый узел;
„ „ а Ь а Ь
Р2 — устойчивый узел при - < — и седло при - > —;
й / а /
„ „ с а с а
Р3 — устойчивый узел при - < - и седло при - > -.
ей ей
Точка Р.4 может в зависимости от параметров находиться в любой четверти и иметь различный тип: седло, устойчивый или неустойчивый узел или фокус. Если точка Р. попадает в первую четверть плоскости, то либо Р2 и Р3 являются седлами, а Р. устойчивым фокусом, либо, наоборот, точки Р2 и Р3 являются устойчивыми фокусами, а Р. седлом.
Наличие в системе положений равновесия типа седла предполагает существование сепаратрис. Фазовый портрет системы (2) при значениях параметров системы а = 1/16, Ь = 1/8,
с = 1/8, й = 1/32, е = 11/128, / = 1/8 представлен на рис. 1 (красным цветом отмечены положения равновесия, зеленым — сепаратрисы).
Рис. 1. Фазовый портрет системы 1СА18С при значениях параметров
а = 1/16, Ь = 1/8, с = 1/8, й = 1/32, е = 11/128, / = 1/8
В качестве локализирующей выберем линейную функцию ^ = Ах + Ву. Универсальное сечение Б^ для этой функции задается уравнением
Асх2 + (АЬ + Ве)ху + В/у2 - Аах - Вйу = 0
и представляет собой кривую 2-го порядка. Линейная функция на кривой 2-го порядка имеет конечные экстремальные значения лишь в случае, когда эта кривая является эллипсом или параболой. Ограничимся первым случаем. Множество Б^ является эллипсом тогда и только тогда, когда дискриминант соответствующей квадратичной формы положителен, т.е. при условии
(АЬ + Ве)2 - 4АВс/ < 0. (3)
Очевидно, что это условие не выполняется в случаях А = 0, В = 0 и АВ < 0. Поэтому можно считать, что АВ > 0. Поскольку умножение локализирующей функции ^ на число не изменяет множества ^, можем считать, что А = 1 и В > 0. При этом следует ограничиться случаем, когда выполнено условие (3). В случае А = 1 это условие сводится к неравенству
В2е2 - 2(2/с - Ье)В + Ь2 < 0
которое имеет решения относительно В, если
/с - Ье > 0.
Далее ограничимся этим случаем. Тогда возможные значения параметра В локализирующей функции определяются неравенствами
(У/с-у/с-ъе)2 (у/с+у/с-&ё)2
-^- < В < ---
(4)
е
е
Итак, для функции <^(х, у) = х + Ву, где В удовлетворяет (4), получаем следующую задачу на экстремум:
х + Ву ^ ехй",
сх2 + (Ь + Ве)ху + Bfy2 — ах — Bdy = 0.
Поскольку эллипс — ограниченная гладкая кривая, то значения ^¡^ и ^8ир могут быть найдены методом множителей Лагранжа. Другой вариант — приведение квадратичной функции к каноническому виду методом Лагранжа. Выделим в левой части ограничения полный квадрат по переменной х:
( Ь + Ве а)2 А 2 к а2 4х + — 2с) +4Су +2Су — 4с = 0'
где
А = 4В^ — (Ь + Ве)2, к = (Ь + Ве)а — 2Bcd. Затем выделим полный квадрат по переменной у:
Ь + Ве а )2 А ( к )2 _ а2 к2
щх +
2с у 2с/ +4Ду + АУ 4с + 4сА' Выполнив замену переменных
Ь + Ве а к
X = х +--у--, У = у + —,
2с 2с у А'
сведем задачу (5) к следующей:
аХ + вУ + 7 ^ ех№, рХ2 + дУ2 = г2,
(6)
где
Л (2с —е)В —Ь B(2аf — Ьd — аЬ) + В 2(2cd — de — ае)
а = 1, в =-2с-, 7 =-А-,
А 2 к2 + а2 А Р = с, д = —, г =—-—-—. 1 ' 4 4с 4сА
Решения задачи (6) легко записываются аналитически:
/а2 в2 /а2 в2 (В) = 7 — г*/--1--, ^вир(В) = 7 + г*/--1--.
V р д V р д
В результате получаем семейство П^(В) локализирующих множеств, описываемых неравенствами
^(В) < х + Ву<^р(В). (7)
Для построения пересечения П^ = Р| П^ (В) перепишем неравенства (7) в виде
в
—Ву + (В) < х < — Ву + ^зир(В).
рис.2
Тогда множество ^ можно представить неравенством
шах(-Ву + (В)) < х < шт(-Ву + ^8ир(В)). (8)
в в
Из множества следует выделить часть , попадающую в первую четверть. Для этого в (8) достаточно добавить еще два неравенства:
шах |шах(-Ву + (В)), о| < х < шт(—Ву + (В)), у > 0. В качестве примера рассмотрим систему 1СЛ18С с параметрами
а = - Ь =1 с =1 й =— е = — а 16' 8' С 8' 32' е 128'
При этих значениях параметров неравенство (4) принимает вид
/ =
(9)
(10)
336 - 128^5 336 + 128^5 -121-<В<-121-•
Максимум и минимум в неравенствах (9) можно вычислять по конечному набору значений параметра В. Это несколько расширит локализирующее множество, но сохранит главное его свойство — оно будет локализирующим для инвариантных компактов системы 1СЛ18С. Результаты вычислений приведены на рис. 2.
0.3
0.25
0.2-
0.15
0.1
0.05
-0.05
с N
\
0.1
0.2
0.3
0.4
0.5
0.6
Рис. 2. Локализирующее множество для системы 1СА18С при значениях параметров а = 1/16, Ь = 1/8, с = 1/8, d = 1/32, е = 11/128, / = 1/8
Замечание 1. Локализирующие множества для инвариантных компактов системы 1СА18С получены лишь в частном случае /с — Ье > 0. Итоговый результат дает компактное локализирующее множество. Из этого, в частности, следует, что система имеет максимальный инвариантный компакт. Действительно, объединение всех инвариантных компактов системы есть инвариантное множество. Замыкание этого объединения также является инвариантным множеством. Существование компактного локализирующего множества означает,
1
8
рис.3
что указанное замыкание представляет собой инвариантный компакт, который содержит в себе все другие инвариантные компакты.
Для системы с параметрами (10), как следует из анализа фазового портрета, максимальный инвариантный компакт представляет собой замкнутую область, ограниченную сепара-трисным контуром, натянутым на точки покоя Pb P2, P3.
Замечание 2. Анализ фазового портрета при различных значениях параметров показывает, что максимальный инвариантный компакт (если ограничиться первой четвертью фазовой плоскости) существует и в случае f c — be < 0, хотя локализация с помощью линейных функций приводит к тривиальным результатам. Однако, если при вычислении экстремальных значений ^¡nf и ^sup учесть ограничения x > 0, y > 0, то нетривиальное локализирующее множество дает любая линейная функция <^(x, y) = Ax + By с неотрицательными коэффициентами A и B. В частности, при <^(x, y) = x получаем множество Q = {(x, y): 0 < x < a/c}, а при <^(x, y) = y — множество Q2 = {(x, y): 0 < y < d/f}. Таким образом, все инвариантные компакты первой четверти содержатся в прямоугольнике 0 < x < a/c, 0 < y < b/f. Однако аналитическое решение экстремальных задач при дополнительных ограничениях типа неравенства существенно труднее, чем решение при отсутствии таких ограничений.
3. Локализация инвариантных компактов системы VDFM
Система VDFM
x = x(1 — x) — sx — axy; , ч
• ( . ) y; (11) y = sx — by,
описывает процесс взаимодействия зеленой x и мертвой y биомасс внутри экосистемы [3]. В системе s — коэффициент старения зеленой биомассы (0 < s < 1); a > 0 — коэффициент подавления роста зеленой биомассы, вызванного ростом мертвой биомассы; b > 0 — коэффициент разложения мертвой биомассы.
Содержательный смысл система (11) имеет лишь в первой четверти фазовой плоскости. Однако для этой системы, в отличие от системы ICAISC, первая четверть не является инвариантным множеством, поскольку на оси абсцисс векторное поле системы не касается этой оси. Можно лишь утверждать, что первая четверть для системы VDFM является положительно инвариантным множеством: траектории, при t = 0 начинающиеся в R+, при t > 0 остаются в R+.
Система (11) имеет два положения равновесия:
' b(1 — s) s(1 — s)
Pi = (0, 0); P2
b + as ' b + as
При этом точка Р\ является седлом, а точка Р2 — устойчивым узлом. Фазовый портрет системы представлен на рис. 3 (красным цветом обозначены положения равновесия, зеленым — сепаратриса, идущая из Р1 в Р2).
Рис. 3. Фазовый портрет УБРМ-системы при значениях параметров а = 3, Ь = 3, £ = 0,35
Линейная локализирующая функция. Использование линейных функций на всей фазовой плоскости приводит к тривиальным результатам. Однако ситтуация меняется, если ограничиться первой четвертью фазовой плоскости. Действительно, выбрав функцию ^о(х, у) = х, получим универсальное сечение
Б^0: х(1 — х — в — ау) = 0,
состоящее из пары прямых. Множество П М+ ограничено, а функция (х, у) достигает на нем наименьшего и наибольшего значений ^0тт(М+) = 0, ^0тах(М+) = 1 — в. Это позволяет в задачах локализации ограничиться множеством
Q = {(х, у): 0 < х < 1 — в, у > 0} .
В качестве локализирующей рассмотрим функцию вида ^(х, у) = Ах + у. Тогда универсальное сечение будет описываться уравнением
Ах2 + Ааху — (А(1 — в) + в)х + Ьу = 0.
Возникают экстремальные задачи
Ах + у — 1п£,
Ах2 + Ааху — (А(1 — в) + в)х + Ьу = 0; Ах + у — 8Ир,
Ах2 + Ааху — (А(1 — в) + в)х + Ьу = 0,
(12) (13)
решение которых будем искать в полосе 0 < х < 1 — в.
Для решения экстремальных задач (12), (13) из уравнения связи выразим переменную у
через х:
у=
(А(1 — в) + в)х — Ах2 Аах + Ь
С помощью этого соотношения исключим из целевой функции переменную у. В результате получим следующие, эквивалентные (12) и (13), задачи:
(A(1 - s) + s)x - Ax2
Ax + —-{---► inf,
Aax + b
(A(1 - s) + s)x - Ax2
Ax + —-{---► sup,
Aax + b
0 < x < 1 - s. (14)
Решения задач (14) ^¡^ (А) и ^8ир(А) зависят от соотношений между параметрами. При А > 1/а целевая функция
(А(1 - в) + в)х - Ах2
—(х) = Ах + ^-7—Т1-
Аах + Ь
возрастает на отрезке [0' 1 - в]. Значит, в этом случае ^¡^(А) = — (0) = 0, ^8ир(А) = -0(1 - в). При 0 < А < 1/а функция — (х) выпукла вверх при х > 0 и достигает максимума в точке
_ b 1 /b(Aa(1 - s) + as + b)
xm _ - Aa + Aav 1-Aa .
Следовательно, минимальное значение на [0, 1 - s] функция достигает в одном из концов отрезка, т.е. ^inf(A) _ min {0, -0(1 - s)}. Максимальное значение достигается либо в точке xm, либо в точке 1 - s, а именно в той, что на числовой оси находится левее. Таким образом, ^sup(A) _ -(min {xm, 1 - s}).
При A _ 0 функция — (x) линейная, так что в этом случае <^inf(A) _ 0, ^sup(A) _ in \ s(1 - s) _ —(1 - s) _ —b—.
При —b—у < A < 0 функция — (x) выпукла вниз. Ее максимальное значение на
[0, 1 - s] достигается в одном из концов отрезка, т.е. ^sup(A) _ max {0, —((1 - s)}. Минимальное значение достигается в точке локального минимума xm, если эта точка попадает на отрезок [0, 1 - s], иначе в одном из концов отрезка. Таким образом, в этом случае
, „ч I — (xm), 0 < xm < 1 - s;
^inf(A)_^ • rn in u
min {0, —(1 - s)}, иначе.
Если же A < ————-, то особая точка -функции —(x) попадает на отрезок
a(l — s) Aa
[0, 1 - s]. В результате в этом случае получаем тривиальный результат <^inf(A) _ -то,
^sup(A) _ то.
Итак, для инвариантных компактов множества R+ получено семейство локализирующих множеств
b
<Anf(A) < Ax + y < ^sup(A), a > —---.
a(1 - s)
Пересечение этого семейства можно записать в следующем виде:
sup (-Ax + Pinf(A)) < y < inf (-Ax + ^sup(A)). (15)
b '' (1-s)
A>--Щ—1
рис.4
0.08 0.07 0.06 0.05 0.04 0.03 0.02 0.01 0
-0.01
0 0.1 0.2 0.3 0.4 0.5
Рис. 4. Локализирующее множество для инвариантных компактов системы УОБМ при значениях параметров а = 3, Ь = 3, £ = 0,35, полученное с помощью линейных функций
Неравенства (15) удобны для компьютерного расчета локализирующего множества. Эти расчеты для случая а =1, Ь =1, в = 0,35 представлены на рис. 4.
Квадратичная локализирующая функция. Если локализирующая функция является квадратичной, то универсальное сечение для системы УЭБМ в общем случае будет кривой 3-го порядка. Выделим случай, когда универсальное сечение есть кривая 2-го порядка, причем это эллипс.
Теорема 1. Для локализирующей функции
<р(х, у) = Ах2 + 2Вху + Су2 + 2(Вх + Еу) (16)
производная в силу системы УЭБМ является квадратичной функцией тогда и только тогда, когда
А = В = 0, С = 0. (17)
При выполнении этих условий универсальное сечение является эллипсом тогда и только тогда, когда дискриминант соответствующей квадратичной формы будет положителен, т.е.
в2С2 — 2(ав + 2Ь)СВ + а2В2 < 0. (18)
Доказательство. Найдем производную в силу системы для функции (16):
^ = —2Ах3 + (2А — 2Ав — 2В + 2Вв)х2 + (2В — 2Вв + 2Ев)х — 2СЬу2 — 2ЕЬу +
+ (—2Аа — 2В )х2 у — 2Ваху2 + (2В — 2Вв — 2Ва — 2ВЬ + 2Св)ху.
Это выражение задает квадратичную функцию, если коэффициенты при степенях выше двух равны нулю. В результате приходим к равенствам А = В = 0. При выполнении этих равенств локализирующая функция имеет вид
^у) = Су2 + 2(Ох + Еу) а соответствующее универсальное сечение Б^ описывается уравнением
Ох2 + (Да - Св)ху + СЬу2 - (.0(1 - в) + Ев)х + ЕЬу = 0. (19)
Уравнение (19) является кривой 2-го порядка в том случае, если один из коэффициентов С и О отличен от нуля. При этом уравнение (19) задает эллипс тогда и только тогда, когда дискриминант квадратичной формы отрицателен (при этом учитываем, что случаи мнимого и вырожденного эллипса невозможны). Отрицательность дискриминанта эквивалентна неравенству
(Оа - Св)2 - 4С.Ь < 0' (20)
которое, в свою очередь, равносильно неравенству (18). Теорема доказана.
Очевидно, что случаи С = 0, О = 0, СО < 0 не удовлетворяют условию (18). Поэтому будем рассматривать только случай СО > 0.
Поскольку эллипс — ограниченная гладкая кривая, то значения ^¡^ и р8ир могут быть вычислены с помощью метода множителей Лагранжа. Однако уравнение на критические точки функции Лагранжа имеет 4-ю степень, так что его решения, хотя и могут быть представлены в аналитическом виде, оказываются очень громоздки и для анализа малопригодны. Поэтому для получения результатов остановимся на разработке соответствующей вычислительной процедуры.
Рассмотрим квадратичную функцию, полагая С = 1 :
р = у2 + 2Ох + 2Еу. При этом универсальное сечение будет эллипсом, если
ав + 2Ь - 2^Ь(аз + Ь) < О < ав + 2Ь + 2^Ь(аз + Ь)'
или
(УЬ + ав-Уб )2 <О< (УЬ + ав + Уб )2. (21)
Для решения задачи
у2 + 2Ох + 2Еу ^ ех^ Ох2 + (Оа - в)ху + Ьу2 - рх + ЕЬу = 0'
где р = О(1 - в) + Ев, выразим из ограничения переменную х и подставим в целевую функцию. Ограничение представляет собой квадратное относительно х уравнение. Ясно, что
при поиске максимума целевой функции необходимо выбирать больший корень квадратного уравнения, а при поиске минимума — меньший. В результате придем к задачам
y2 + 2D£i(y) + 2Ey ^ min, y2 + 2D&(y) + 2Ey ^ max, (22)
где
р - (Ра - в)у - У(р - (Ра - в)у)2 - 4Р(Ьу2 + ЭД
Ыу) =-2р-'
р - (Ра - в)у + У(р - (Ра - в)у)2 - 4Р(Ьу2 + ЁЬУ)
С2 (у) =-2Р-•
Область изменения переменной у находим из условия положительности дискриминанта квадратного уравнения:
(р - (Ра - з)у)2 - 4Р(Ьу2 + ЕЬу) > 0,
или
[4РЬ - (Ра - в)2] у2 + 2 [(Ра - в)р + 2РЕЬ]у - р2 - 0. (23)
Из уравнения (23), полагая
5 = [(Ра - в)р + 2РЕЬ]2 + р2 [4РЬ - (Ра - в)2]
и учитывая неравенство (20) при С = 1, находим
-(Ра - в)р - 2РЕЬ - -(Ра - в)р - 2РЕЬ +
4РЬ - (Ра - в) - у - 4РЬ - (Ра - в) •
На указанном интервале изменения переменного у решения задач (22) можно найти любым стандартным численным методом.
Решив задачи (22), получим двухпараметрическое семейство П(Р, Е) локализирующих множеств
Ум(Р, Е) - у2 + 2Рх + 2Еу - ^вир(Р, Е),
или
Еу + ^(Р, Е) ^ ЕЁ + ^вир(Р,Е) (24)
2Р Р + 2Р < Х < 2Р Р + 2Р • ( )
Неравенства (24) позволяют записать формулу для пересечения Р| П(Р, Е) локализирующих множеств:
( У! Еу + ^(Р,Е)^ ^ ( Еу + ^(Р,Е)
^-2Р-Р + 2Р Х < ЩЕ^ -2Р-Р + 2Р
где диапазон изменения Р определяется неравенством (21), а диапазон изменения Е — вся числовая ось. Область изменения переменной у включает те ее значения, при которых левая часть двойного неравенства не превышает правой части. Для учета ограничения динамической системы на первую четверть необходимо добавить условие у > 0 в область
изменения y и ввести в двойное неравенство ограничение x > 0:
in ( У2 Ey + yinf(D,E)\\ .. ( y2 Ey + ySUp(D, E)^
max< 0, sup--—--— +--—- > < x < mi--—--— +---- , (25)
\ d,E\ 2 D D 2D yj - - D,E\ 2 D D 2D )
В качестве примера рассмотрим систему VDFM (11) при значениях параметров а = 3, b = 3, s = 0,35 (см. рис. 3). Локализирующее множество для ее инвариантных компактов, построенное в соответствии с неравенством (25), представлено на рис. 5.
Рис. 5. Локализирующее множество для инвариантных компактов системы УББМ при значениях параметров а = 3, Ь = 3, £ = 0,35
Замечание 3. Получение компактного локализирующего множества для системы УЭРМ, как и в случае системы 1СА18С, указывает на существование у системы максимального инвариантного компакта. Анализ фазового портрета системы УЭРМ при указанных значениях параметров приводит к выводу, что максимальным инвариантным компактом этой счистемы является объединение двух положений равновесия и соединяющей их сепаратрисы. Поэтому полученные локализирующие множество (см. рис. 4, 5) не являются точными. Тем не менее отметим, что нижняя граница локализирующего множества, полученного с помощью квадратичных функций, оказывается очень близкой к сепаратрисе.
Заключение
В статье рассмотрена задача локализации инвариантных компактных множеств для непрерывных динамических систем второго порядка. Эти системы отличаются относительно простым поведением, однако решение задач локализации для таких систем наталкивается на определенные трудности.
Исследование проведено для двух динамических систем, представляющих собой простейшие биологические модели. В каждом случае построено параметрическое семейство
локализирующих множеств и найдено пересечение этого семейства. В обоих случаях решение задачи локализации строилось на использовании функционального метода локализации с привлечением как аналитических вычислений, так и численных расчетов.
Работа выполнена при финансовой поддержке РФФИ (гранты 11-01-00733 и 12-07-00267) и Программы Президента РФ по государственной поддержке ведущих научных школ (грант НШ-3659.2012.1).
Список литературы
1. Канатников А.Н., Крищенко А.П. Инвариантные компакты динамических систем. М.: Изд-во МГТУ им. Н.Э. Баумана, 2011. 231 с.
2. Chrobak J.H., Herrero H. A mathematical model of induced cancer-adaptive immune system competition // Journal of Biological Systems. 2011. Vol. 19, no. 3. P. 521-532.
3. Wang K., Zhang N., Niu D. Periodic oscillations in a spatially explicit model with delay effect for vegetation dynamics in freshwater marshes // Journal of Biological Systems. 2011. Vol. 19, no. 2. P. 131-147.
4. Канатников А.Н. Локализация инвариантных компактов ПРТ-системы // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2007. № 1. С. 3-18.
5. Крищенко А.П. Локализация инвариантных компактов динамических систем // Дифференциальные уравнения. 2005. Т. 41, № 12. С. 1597-1604.
6. Krishchenko A.P., Starkov K.E. Localization of compact invariant sets of the Lorenz system // Phys. Lett. A. 2006. Vol. 353, no. 5. P. 383-388.
7. Krishchenko A.P., Starkov K.E. Localization of compact invariant sets of nonlinear systems with application to the Lanford systems // Int. J. of Bifurcation and Chaos. 2006. Vol. 16, no. 11. P. 3249-3256.
8. Krishchenko A.P., Starkov K.E. Estimation of the domain containing all compact invariant sets of a system modeling the amplitude of a plasma instability // Phys. Lett. A. 2007. Vol. 367, no. 1. P. 65-72.
9. Starkov K.E. Bounds for the domain containing all compact invariant sets of the system modeling dynamics of acoustic gravity waves // Int. J. of Bifurcations and Chaos. 2009. Vol. 19, no. 10. P. 3425-3432.
10. Starkov K.E. Bounding a domain which contains all compact invariant sets of the Bloch system // Int. J. of Bifurcations and Chaos. 2009. Vol. 19, no. 3. P. 1037-1042.
11. Coria L.N., Starkov K.E. Bounding a domain containing all compact invariant sets of the permanent-magnet motor systems // Comm. Nonlin. Sci. and Numer. Simul. 2009. Vol. 14. P. 3879-3888.
SCIENTIFIC PERIODICAL OF THE BAUMAN MSTU
SCIENCE and EDUCATION
EL № FS77 - 48211. №0421200025. ISSN 1994-0408
electronic scientific and technical journal
Localization of invariant compact sets of two-dimensional
continuous dynamical systems
# 07, July 2013
DOI: 10.7463/0713.0583104
Kanatnikov A. N., Phedorova Yu.P.
Bauman Moscow State Technical University 105005, Moscow, Russian Federation [email protected]
One method of the qualitative analysis of a dynamical system is to estimate the position of its compact invariant sets closely associated with bounded trajectories of the system. As a solution to such a problem, one can use a localizing set, i.e. a set in the phase space containing all invariant compact sets of the system. In this article two continuous two-dimensional dynamical systems describing behavior of some biological systems are explored. For each of these systems a family of localizing sets is constructed, and then the intersection of the family is calculated. For the first system the solution was obtained analytically and for the second one the numerical procedure of constructing localizing sets was proposed. The investigation results are shown in figures.
References
1. Kanatnikov A.N., Krishchenko A.P. Invariantnye kompakty dinamicheskih system [Invariant compact sets of dynamical systems]. Moscow, Bauman MSTU Publ., 2011. 231 p.
2. Chrobak J.H., Herrero H. A mathematical model of induced cancer-adaptive immune system competition. Journal of Biological Systems, 2011, vol. 19, no. 3, pp. 521-532.
3. Wang K., Zhang N., Niu D. Periodic oscillations in a spatially explicit model with delay effect for vegetation dynamics in freshwater marshes. Journal of Biological Systems, 2011, vol. 19, no. 2, pp. 131-147.
4. Kanatnikov A.N. Lokalizatsija invariantnyh kompaktov PRT-sistemy [Localization of the invariant compact sets of the PRT system]. Vestnik MGTU im. Baumana. Ser. Estestvennye nauki [Bulletin of the Bauman MSTU. Ser. Natural Sciences], 2007, no. 1, pp. 3-18.
5. Krishchenko A.P. Lokalizatsija invariantnyh kompaktov dinamicheskih system [Localization of the invariant compact sets of the dynamical systems]. Differentsial'nye uravnenija, 2005, vol. 41, no. 12, pp. 1597-1604. (English version: Krishchenko A.P. Localization of In-
variant Compact Sets of Dynamical Systems. Differential Equations, 2005, vol. 41, no. 12, pp. 1669-1676).
6. Krishchenko A.P., Starkov K.E. Localization of compact invariant sets of the Lorenz system. Phys. Lett. A, 2006, vol. 353, no. 5, pp. 383-388.
7. Krishchenko A.P., Starkov K.E. Localization of compact invariant sets of nonlinear systems with application to the Lanford systems. Int. J. of Bifurcation and Chaos, 2006, vol. 16, no. 11, pp. 3249-3256.
8. Krishchenko A.P., Starkov K.E. Estimation of the domain containing all compact invariant sets of a system modeling the amplitude of a plasma instability. Phys. Lett. A, 2007, vol. 367, no. 1, pp. 65-72.
9. Starkov K.E. Bounds for the domain containing all compact invariant sets of the system modeling dynamics of acoustic gravity waves. Int. J. of Bifurcations and Chaos, 2009, vol. 19, no. 10, pp. 3425-3432.
10. Starkov K.E. Bounding a domain which contains all compact invariant sets of the Bloch system. Int. J. of Bifurcations and Chaos, 2009, vol. 19, no. 3, pp. 1037-1042.
11. Coria L.N., Starkov K.E. Bounding a domain containing all compact invariant sets of the permanent-magnet motor systems. Comm. Nonlin. Sci. and Numer. Simul, 2009, vol. 14, pp. 3879-3888.
