КАЧЕСТВЕННОЕ ИССЛЕДОВАНИЕ СИСТЕМЫДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙНА ИНВАРИАНТНОМ ТРЕУГОЛЬНИКЕ ЧАСТОТ Текст научной статьи по специальности «Математика»

ФУНКЦИОНАЛЬНЫЙ АНАЛИЗ И ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ

Научная статья УДК 517.925.5

Б01: 10.18101/2304-5728-2023-3-3-13

КАЧЕСТВЕННОЕ ИССЛЕДОВАНИЕ СИСТЕМЫ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ НА ИНВАРИАНТНОМ ТРЕУГОЛЬНИКЕ ЧАСТОТ

© Бортковская Мария Романовна

кандидат физико-математических наук, доцент, Санкт-Петербургский политехнический университет Петра Великого Россия, 195251, г. Санкт-Петербург, ул. Политехническая, 29 mbort@mail.ru

Аннотация. В статье рассматривается 6-параметрическое семейство автономных систем обыкновенных дифференциальных уравнений, применяемых в классических моделях конкурентного взаимодействия трех групп представителей одной популяции. Такие системы могут использоваться для моделирования взаимодействия трех конкурирующих групп в задачах разной тематики. Системы рассматриваются на инвариантном треугольнике частот; исследование проводится для произвольных значений параметров, за исключением случаев их специальных соотношений. Дано определение множества приближения к особой точке (удаления от нее) относительно семейства отрезков. Исследуются границы этих множеств для особых точек, расположенных в вершинах треугольника частот.

Ключевые слова: автономная система, обыкновенные дифференциальные уравнения, параметрическое семейство систем, трехмерная система, динамическая модель, особые точки, фазовый портрет системы, инвариантное множество, численный пример.

Для цитирования

Бортковская М. Р. Качественное исследование системы дифференциальных уравнений на инвариантном треугольнике частот // Вестник Бурятского государственного университета. Математика, информатика. 2023. № 3. С. 3-13.

Введение

Будем рассматривать систему трех дифференциальных уравнений вида

х7 = X; (Ах), - хТАх), 7 = 1,2,3, (1)

(0 а Ь

где А =

у 0 8 е й 0

матрица с нулевыми диагональными элементами,

х(Г) =

( х^) ^

Х2^ )

V х3({) 0

— искомая вектор-функция. Система вида (1) с конкретны-

ми числовыми значениями элементов матрицы А приведена в книге [1] как пример динамической модели конкуренции представителей трех групп в пределах популяции животных одного вида. Также система такого вида применялась автором данной статьи для моделирования конкуренции трех групп личных имен в русском языке (см. [2]). Система (1) рассматривается на инвариантном треугольнике частот

О = Х' = 1'Х' ~ 01 (хг(^) — относительная частота в популяции представителей 1 -й конкурирующей группы в момент времени t). Ненулевые элементы матрицы А естественно рассматривать как параметры, а множество систем вида (1) с произвольными значениями параметров как параметрическое семейство систем. В книге [1] разбирается численный пример модели конкуренции в популяции и объясняется, почему в такой модели всегда можно, не умаляя общности, считать диагональные элементы матрицы А нулевыми, а также утверждается инвариантность О относительно системы (1), которую легко проверить непосредственно (см. [2]).

1 Постановка задачи

В настоящей статье рассмотрим особые точки системы (1), которые могут в зависимости от значений параметров а,р,у,8,е,^ возникать внутри и на сторонах треугольника частот О (его вершины — особые точки системы (1) при любых значениях параметров)' исследуем (предварительно дав необходимые определения) множества приближения к особым точкам и удаления от них для особых точек системы (1), находящихся в вершинах треугольника частот; разработаем способ нахождения уравнений кривых - границ множеств приближения (удаления) и исследуем геометрическое расположение найденных кривых по отношению к особым точкам системы (1). Цель работы: не решая систему численно, изучить качественно ее фазовый портрет на треугольнике частот для произвольных значений параметров системы, за исключением особых случаев (когда параметры связаны специальными соотношениями), выделить такие случаи. После доказательства общих утверждений и описания способа нахождения границ множеств приближения к особым точкам (удаления от них) рассмотрим пример с конкретными числовыми значениями параметров системы (1), иллюстрирующий результаты исследования.

2 Особые точки системы (1)

Как уже отмечено, вершины инвариантного треугольника частот О являются особыми точками системы (1) при любых значениях параметров, что легко проверяется. Эти особые точки всегда сложные (матрица линеаризованной системы в этих точках имеет по крайней мере одно нулевое собственное число). На стороне треугольника частот {х3 = 0, х + х2 = 1, х1 > 0,1 = 1,2} в случае, когда ау > 0 , имеется особая

а у Л

точка -,-,0 . Она простая (седло либо узел), если

ка + У а + у )

уС + ае — ау Ф 0. Если ау = 0 , а2 + у2 Ф 0, то эта особая точка сложная и совпадает с вершиной треугольника частот. Если ау < 0 или а2 + у2 = 0, то на соответствующей стороне треугольника О также нет особых точек системы (1), за исключением вершин треугольника. Анало-

( 8 С ^

гично рассматриваются особые точки е Л

-;0;-

0;-

'8 + £'8+£

и

, которые при соответствующих значениях параметров

Ь+е ь

имеются на двух других сторонах треугольника частот (подробно об этом в статье [2]). Стороны треугольника О в свою очередь являются инвариантными множествами системы (1), и на них легко из вида правых частей системы (1) определить направления движения. Исследование возможных сложных особых точек здесь не проводится, но из дальнейших рассуждений будет видно, что в ряде случаев оно может и не понадобиться для получения качественного фазового портрета системы (1) на О .

Теорема 1. Для того чтобы точка N (х1, х2, х3), такая, что х1 > 0, х2 > 0, х3 > 0, была особой точкой системы (1), расположенной внутри треугольника частот О , необходимо и достаточно, чтобы ее координаты х1, х2, х3 удовлетворяли системе линейных алгебраических уравнений

ах2 + /Зх3 =ух1 +8х3,

ех1 + Сх2 = ух1 +8х3, (2)

I х2 I х3 — 1.

Доказательство. Система (2) получается приравниванием нулю правых частей системы (1) при условиях х1 • х2 • х3 Ф 0 и х1 + х2 + х3 — 1. Отсюда очевидна необходимость заявленного в теореме условия неподвижности точки N. Докажем его достаточность. Пусть х1,х2, х3 - решение системы (2). Тогда имеем

3

ух1 + 8х3 — (ух1 + 8х3) • ^ х1 — ухе12 + 8х1 х3 + ух1 х2 +8х2х3 + ух1х3 + 8х32 .

I—1

Это равенство, поскольку в силу системы (2) ух1 — ах2 + Ьх3 — 8х3 и

8х3 — еёс1 +Сх2 — ух1, означает, что

ух1 + 8х3 — х1(ах2 + ¡Зх3 — 8х3) + 8х1 х3 + +ух1 х2 + 8х2 х3 + ух1 х3 + х3 (ех1 +С,хг — ух1).

Отсюда уХ1 + 5Х3 = Х1(аХ2 + /ЗХ3) + Х2(уХ1 +5Х3) + Х3(еХ1 + ^Х2) = ХТАХ, а поскольку в силу системы (2) верны равенства УХ +5Х3 =аХ2 + /ЗХ3 =еХ1 + ^Х2, то будут верны и равенства аХ2 + ЬХ3 = ХТАХ и еХ1 + ^Х2 = ХТАХ, то есть в точке N все правые части системы (1) обратятся в нуль. Теорема доказана.

Замечание 1. Вычислив определитель системы (2), получаем, что при условии (у+5)(-5-£) - (у -е + £)(-а + Ь -5) Ф 0 система (1) может иметь внутри О не больше одной особой точки.

Замечание 2. Легко убедиться, что, если система (2) имеет такое решение, что две его компоненты (координаты) положительны, а третья равна нулю, то это решение совпадает с координатами одной из особых точек системы (1), лежащих на сторонах (но не в вершинах) треугольника О, найденных ранее.

3 Границы множеств приближения к особым точкам и удаления от особых точек системы (1), расположенных в вершинах треугольника частот

Пусть {1С }, 0 < С < 1, — непрерывное однопараметрическое семейство

отрезков параллельных прямых; все отрезки лежат в треугольнике частот О и через каждую внутреннюю точку О проходит ровно один отрезок. Пусть Р - особая точка системы (1), расположенная в вершине треугольника О, и имеет место следующее: каждый луч, исходящий из точки Р,

3

лежащий в плоскости ^ Х1 = 1 и направленный внутрь О, пересекает ка-

1=1

ждый отрезок 1С ровно в одной точке, причем значения С , соответствующие пересекаемым лучом прямым, меняются вдоль луча строго монотонно.

Определение 1. Такое семейство отрезков {1С} будем называть правильным по отношению к особой точке Р .

Пусть точка Р и семейство отрезков {1С } таковы, как в определении 1. Пусть для определённости при движении вдоль каждого луча, исходящего из точки Р, значения С убывают. Рассмотрим точку Т внутри О, она лежит на некотором (единственном) отрезке &{1С} . Пусть траектория

системы (1), проходящая через точку Т, имеет участок ё , содержащий внутри точку Т и такой, что при движении по этому участку траектории в соответствии с системой (1) значения С при пересечении ё с отрезками {1С } строго возрастают.

Определение 2. Такую точку Т назовем точкой приближения к особой точке Р относительно правильного семейства отрезков {1С } . (В этом определении, очевидно, можно поменять местами убывание и возрастание

параметра C вдоль произвольного луча, исходящего из P, и вдоль участка траектории, проходящей через T , соответственно).

Аналогично определим точку удаления от особой точки P относительно семейства отрезков {lC } .

Определение 3. Плоское множество G с D , состоящую из точек приближения к особой точке P (удаления от особой точки P) относительно правильного семейства отрезков {lC }, назовем, соответственно, множеством приближения (удаления) точки P относительно этого семейства отрезков.

Ниже, если оговорено, о каком правильном семействе отрезков идет речь, будем для краткости называть множества приближения к особой точке (удаления от нее) относительно этого семейства просто множествами приближения (удаления) особой точки.

Теорема 2. Пусть g — пересечение внутренности int D треугольника

3

частот D с кривой g, лежащей в плоскости ^ xt = 1, и такой, что ее про-

¿=1

екция на плоскость x1Ox2 задается уравнением второго порядка (ß + e) x12 + (S + Z) x22 + (-a-g+ß+S+e+ Z) x1 x2 + +(-2ß -e)x1 + (a-ß-S-Qx2 + ß = 0.

Тогда g разделяет множество приближения к особой точке S1 (1;0;0) и множество удаления от особой точки S1(1;0;0) системы (1) относительно некоторого правильного семейства отрезков, если только данное уравнение второго порядка не задает прямую, точку или пустое множество.

Если g nint D = 0, то int D — множество приближения или множество удаления особой точки S1 .

Доказательство. Рассмотрим всевозможные плоскости, задаваемые уравнениями x2 + x3 = C, 0 < C < 1. Их пересечения с инвариантным треугольником D — множество отрезков, параллельных стороне треугольника D, противолежащей вершине S1(1;0;0) (особой точке из формулировки теоремы 2). Ясно, что это множество отрезков — правильное семейство отрезков по отношению к S1 . Если траектория системы (1), проходя через точку P , лежащую внутри D на одном из этих отрезков, пересекает с ростом времени t данный отрезок «в направлении» роста значений C , то точка P лежит в множестве удаления от особой точки S1 , а если траектория пересекает этот отрезок «в направлении» убывания значений C , то точка P лежит в множестве приближения к особой точке S1 . Точки внутри D , лежащие на границе множеств приближения и удаления особой точки S1 , могут быть только такими, в которых обращается в нуль производная функции x2 + x3 в силу системы (1). Поскольку мы

исследуем систему (1) на треугольнике D, это то же самое, что найти все точки внутри D, в которых обращается в нуль производная функции x1 = 1 - (x2 + x3) в силу системы (1). Приравняв нулю правую часть первого уравнения системы (1) и подставив в полученное равенство x3 = 1 - x1 - x2, получим уравнение второго порядка из формулировки доказываемой теоремы. При переходе точки плоскости x1Ox2 через соответствующую кривую второго порядка производная функции x1 в силу системы (1) меняет знак, за исключением вырожденных случаев, оговоренных в условии теоремы. Если кривая второго порядка не пересекает треугольник {x1 > 0, x2 > 0, x1 + x2 < 1} , то вся int D — множество приближения или множество удаления особой точки S1. Теорема доказана.

Замечание 3. Аналогично можно получить границу множеств приближения и удаления для двух других особых точек S2(0;1;0); S3(0;0;1) системы (1), расположенных в вершинах треугольника частот (относительно семейства отрезков, параллельных стороне треугольника, противолежащей изучаемой особой точке).

Замечание 4. Где именно по отношению к g будет множество приближения, а где множество удаления, легко понять по знаку левой части полученного в теореме 2 уравнения второго порядка.

Замечание 5. Для построения фазового портрета системы (1) удобно для всех особых точек, расположенных в вершинах D , спроецировать границы множеств приближения и удаления относительно соответствующих семейств прямых на одну и ту же координатную плоскость (неважно, на какую). Ниже в числовом примере для особых точек S2(0;1;0) и S3 (0; 0; 1) обозначим соответственно как g2 и g3 кривые, аналогичные кривой g для точки S1, и выпишем уравнения проекций этих кривых на плоскость x1Ox2. В общем виде их также легко записать.

Замечание 6. Рассмотрение линий уровня функции x2 + x3 и исследование ее монотонности вдоль траекторий системы (1) опирается на идеи второго метода Ляпунова ([3], [4]). Но здесь мы не исследуем как таковую устойчивость особой точки. Рассуждения, подобные примененным в доказательстве теоремы 2, демонстрируются в книге [1] с помощью числового примера, где рассмотрены линии уровня нелинейной функции, там такой подход назван рассуждениями «типа Ляпунова».

Теорема 3. Кривая g,i = 1,2,3, разделяющая в D множества приближения и удаления одной из особых точек Si системы (1), рассмотренные выше, проходит через особые точки системы (1), которые являются внутренними точками сторон треугольника D , смежных с вершиной Si , если такие особые точки существуют.

Доказательство. Рассмотрим для определенности 51(1; 0; 0). На смежных с этой вершиной сторонах треугольника частот могут быть особые

( в 8 Л

точки системы (1): К -,0,- (если только Ве> 0) и

УР+8 р+8)

М

( у 0Л

(если только ау > 0 ). Для того, чтобы убедиться, что

а + у а+у

кривая у1 проходит через эти точки, достаточно убедиться, что ее проекция на координатную плоскость х1Ох2 проходит через проекции точек К и М на эту плоскость, поскольку на О любые две координаты точки од-

3

нозначно определяют точку в силу равенства ^ х7 = 1. Подставим

¿=1

В п

х1 =-, х2 = 0 в уравнение второго порядка из теоремы 2, получим

В +8

В2 (2В + е)В (в е Л

В = -В'--1--+ В = 0, то есть у проходит

В+8 В+8 \Р+8 В+8 0

через точку К .

Аналогичные выкладки показывают, что кривая у1 проходит через точку М. Для вершин 52, 53 доказательство полностью аналогично. Теорема доказана.

Теорема 4. Пусть кривая у{, 7 = 1,2,3, разделяет в О множества приближения и удаления одной из особых точек 5, 7 = 1,2,3, системы (1), рассмотренные выше. Если N(Х1, Х2,Х3) — особая точка, расположенная внутри треугольника частот О, то кривая у7,7 = 1,2,3 проходит через точ -ку N.

Доказательство. Для определенности рассмотрим кривую у1 , разделяющую множества притяжения и отталкивания особой точки 5 (1; 0; 0). Кривая у1 является пересечением поверхности второго порядка (Ах) - хТАх = 0 и треугольника О . Из доказательства теоремы 1 видно, что точка N(х1, х2, х3) принадлежит тому же пересечению. Для кривых у2, у3 доказательство аналогично. Теорема доказана.

Замечание 7. Если кривая у1 не содержит внутренних точек треугольника частот, то и особых точек, расположенных внутри О , система (1) не имеет. Обратное неверно.

Замечание 8. Изучение особых точек параметрического семейства систем можно продолжить с помощью методов теории бифуркаций. В работах такой тематики, где рассматривается динамическое моделирование с помощью систем дифференциальных уравнений, можно встретить

применение термина «неподвижная точка» для особой точки (положения равновесия) системы дифференциальных уравнений ([5]). Это понятие отлично от понятия неподвижных (в отличие от подвижных) особых точек решений дифференциальных уравнений ([6], [7]), также применяемого в задачах динамического моделирования.

4 Числовой пример

Рассмотрим систему (1), где параметры имеют следующие значения: а —11,4, Ь — —1,2, у— 3,4,8 — 5,4, е — — 2,4, С —12,7 . Она получена в работе [2] при построении модели конкуренции трех групп женских личных имен в русском языке. Здесь ау > 0,8С > 0, /е> 0, и на сторонах треугольника частот, помимо его вершин, есть особые точка системы (1): К(0,3; 0; 0,7) (неустойчивый узел), ¿(0; 0,3; 0,7) и М(0,8;0,2;0) (устойчивые узлы); координаты точек даны приближенно (с округлением); тип точек определяется по первому приближению. Внутри треугольника частот одна особая точка N(0,342; 0,306; 0,352) (трехмерное седло), найденная в соответствии с теоремой 1. Уравнения проекций на плоскость х1Ох2 границ множеств приближения и удаления для особых точек системы (1), расположенных в вершинах треугольника, ищем согласно теореме 2 и замечанию 3:

1. —3,6х12 + 18,1х2 — 0,3х1х2 + 4,8х1 — 5,5х2 —1,2 — 0 (у'1, проекция у1) для ^(1;0;0);

2. —3,6х12 + 18,1х2 — 0,3х1 х2 + 1,6х1 — 23,5х2 + 5,4 — 0 (у'2, проекция у2) для ^(0; 1; 0);

3. —3,6х12 + 18,1х2 — 0,3х1 х2 +1,2х1 — 5,4х2 — 0 (у'3, проекция у3) для 5"3(0; 0; 1) . (Все три кривые -гиперболы, см. рис.1).

На рисунке 1 точка К' — проекция точки К на плоскость х1Ох2, Ь' — проекция точки Ь и так далее. На рисунке 2 изображены кривые у, у2, у3 и несколько траекторий системы (1) (стрелками указано направление движения по траектории).

Рис. 1. Сплошная линия — кривая у '1, крупный пунктир — кривая у '2,

на плоскость х1Ох2

мелкий пунктир — кривая у 3

Бз

Рис. 2. Сплошная жирная линия — кривая у1, крупный пунктир — кривая у2, мелкий пунктир — у3 . Сплошные линии со стрелками — траектории системы (1)

Рисунки 1 и 2 иллюстрируют утверждения теорем 2, 3 и 4.На рисунке 3 приведен качественный фазовый портрет системы, обобщающий

полученные результаты.

Х1=1М

Рис. 3. Фазовый портрет системы из примера 11

Заключение

Полученные результаты позволяют анализировать поведение траекторий системы (1) не только качественно, но и численно, за исключением специальных (вырожденных) случаев, когда параметры системы связаны определенными соотношениями, отмеченными в статье. В остальных случаях, зная координаты любой неособой точки треугольника частот, мы можем установить поведение проходящей через эту точку траектории по отношению к любой особой точке системы (1), расположенной в одной из вершин треугольника частот, а именно, — указать, в каком из множеств (приближения или удаления) по отношению к данной особой точке расположена выбранная точка. Этот факт удобно использовать для сравнения эспериментальных данных с результатами моделирования, применяя системы вида (1) для динамического моделирования процессов конкуренции. Теоремы о взаимном расположении множеств приближения (удаления) и особых точек системы (1) помогают построить фазовый портрет системы на треугольнике частот. Перспектива дальнейшего исследования системы (1) состоит в дополнительном рассмотрении вырожденных случаев (возможно с помощью методов теории бифуркаций), а также в применении метода, разработанного в настоящей статье, для нахождения границ множеств приближения и удаления для особых точек, не расположенных в вершинах треугольника частот. Поскольку модели популяционной динамики используются для моделирования процессов конкуренции в различных областях ([8], [9]), исследование системы (1) имеет прикладную актуальность.

Литература

1. Эрроусмит Д., Плейс К. Обыкновенные дифференциальные уравнения. Качественная теория с приложениями: пер. с англ. Москва: Мир, 1986. 248 с.

2. Бортковская М. Р. Пример моделирования конкуренции в языке с помощью системы дифференциальных уравнений // Анализ, моделирование, управление, развитие социально-экономических систем (Амур-2022): сборник научных трудов XVI Всероссийской с международным участием школы-симпозиума. Симферополь, 2022. С. 74-79.

3. Бибиков Ю. Н. Курс обыкновенных дифференциальных уравнений: учебное пособие для ун-тов. Москва: Высшая школа, 1991. 303 с.

4. Демидович Б. П. Лекции по математической теории устойчивости: учебное пособие для вузов. Санкт-Петербург: Лань, 2023. 480 с.

5. Лукин Д. С., Гайдук А. Р. Исследование нелинейной системы с бифуркацией // Математические методы в технологиях и технике. 2021. № 2. С. 7-10. Б01: 10.52348/2712-8873_ММТТ_2021_2_7.

6. Арнольд В. И., Ильяшенко Ю. С. Обыкновенные дифференциальные уравнения // Итоги науки и техники. Сер. Современные проблемы математики. Фундаментальные направления. 1985. Т. 1. С. 7-140.

7. Орлов В. Н., Гасанов М. В. Влияние возмущения подвижной особой точки на структуру аналитического приближенного решения одного класса нелинейных

дифференциальных уравнений третьего порядка в комплексной области // Вестник МГТУ им. Н. Э. Баумана. Сер. Естественные науки. 2022. № 6. С. 60-73. DOI: 10.18698/1812-3368-2022-6-60-76.

8. Макаренко С. И. Динамическая модель двунаправленного информационного конфликта с учетом возможностей сторон по наблюдению, захвату и блокировке ресурса // Системы управления, связи и безопасности. 2017. № 1. С. 60-89. Doi: 10.24411/2410-9916-2017-10106.

9. Альрефаи В. А. Моделирование и численный анализ динамических систем с конкурентным взаимодействием // Технологический аудит и резервы производства. 2015. № 2/5(22). С. 36-41. Doi: 10.15587/2312-8372.2015.41156.

Статья поступила в редакцию 02.09.2023; одобрена после рецензирования 25.09.2023; принята к публикации 27.09.2023.

PHASE PORTRAIT OF A DIFFERENTIAL EQUATIONS SYSTEM

MODELING COMPETITIVE INTERACTION

Maria R. Bortkovskaya

Cand. Sci. (Phys. and Math.), A/Prof,

St. Petersburg State Polytechnic University

29 Polytekhnicheskaya St., Saint-Petersburg 195251, Russia

mbort@mail.ru

Abstract. The 6-parametric family of 3-dimentional autonomous ordinary differential equations systems is considered. Systems of that type are classic in the dynamic modeling of competitive interaction in populations, when there are three interacting groups. Such systems can be used in other spheres where we model the dynamics of competitive processes for three groups. We consider the systems on the invariant triangle of frequencies and investigate them for arbitrary parameter values except their special relations. Definitions of approaching and retiring areas of a singular point of system relatively to a family of line segments are formulated. Next, we investigate these areas bounds for singular points located in the triangle vertexes.

Keywords: аutonomous system, ordinary differential equations, parametric family of systems, 3-dimentional system, dynamic model, singular points, phase portrait, invariant set, numeric example.

For citation

Bortkovskaya M. R. Phase Portrait of a Differential Equations System Modeling Competitive Interaction // Bulletin of Buryat State University. Mathematics, Informatics. 2023. N. 3. P. 3-13.

The article was submitted 02.09.2023; approved after reviewing 25.09.2023; accepted for publication 27.09.2023.