Математические структуры и моделирование 2002, вып. 10, с. 152-161
УДК 519.6:51-77
АНАЛИТИЧЕСКОЕ ИССЛЕДОВАНИЕ МАТЕМАТИЧЕСКОЙ МОДЕЛИ «ПОЛИТИЧЕСКАЯ ДИФФЕРЕНЦИАЦИЯ - СТЕПЕНЬ АДАПТАЦИИ»
А.А. Лаптев
This article is about analytical qualitative investigation the system of differential equations. This system describes the behaviour of the political and economic society systems. States of the equilibrium are explored. Phase-plane portraits are plotted.
1. Модель «политическая дифференциация — степень адаптации»
В результате последних исследований на основе формализации социологической теории построена математическая модель глобальных социальных процессов [2,3]. В данной работе представлено исследование системы дифференциальных уравнений (1), которые описывают поведение политической и экономической системы (модель «политическая дифференциация-степень адаптации»).
—=kGa(esv-s1-1 )G+(lGE-kGE E 2)-кГж (K+D)G, at
dE
—=kEE (esv-Sl-1 )E-(Ie G-kEoG2)-kEK (K+D)E, dt
G\t=0 = G0, E\t=0 = E0.
(1)
Здесь G и E переменные, характеризующие развитие соответствующих систем: G - политическая дифференциация, E - степень адаптации к окружающей природной среде (экономика). Величины P, K и D - параметры, зафиксированные на некотором уровне; kGG, 8, 81, lG, kGE, kGK, kEE, lE, kEG, kEK - положительные коэффициенты; G0, E0 - начальные условия.
Для упрощения исследования предположим, что kGG = kEE и kGK = kEK. Введем обозначения: A(P) = kGG(eSV-Sl - 1) - kGK(K + D), k1 = kGE, k2 = kEG, li = lG, l2 = lE, x = G, y = E, x = G0, y = E0. Тогда система (1) примет
© 2002 А.А. Лаптев
E-mail: laptev@univer.omsk.su Омский государственный университет
Математические структуры и моделирование. 2002. Вып. 10.
153
следующий вид
dx
dt
Ax + hy — kiy2,
< dy
dt
Ay — l2x + k2x2,
k x|t=o = x> y|t=o = y.
(2)
Таким образом, мы получили систему из двух дифференциальных уравнений с параметром A, который может принимать любые значения. На начальные данные x и у не накладываем никаких ограничений.
2. Качественное исследование модели «политическая дифференциация-степень адаптации»
Рассмотрим систему (2). Предполагаем, что li, ki, l2, k2 — положительные коэффициенты, а A - знакопеременный параметр.
Проследим за изменением качественной картины решения системы (2) в зависимости от параметра A. Для этого изучим возможные состояния равновесия и направления, по которым траектории могут стремиться к ним. Качественное исследование будем проводить, следуя монографии [1].
2.1. Исследование состояний равновесия
Найдем особые точки. В конечной части плоскости - это точки пересечения двух парабол (A = 0)
Ax + liy — kiy2 = 0, Ay — l2x + k2x2 = 0. (3)
Для исследования состояния равновесия (x0,y0) рассмотрим величины P (x,y) = Ax + liy — kiy2, Q(x,y) = Ay — ^x + k2x2,
A = A(xo,yo)
Px(x0,y0) Py (x0,y0)
Qx(x0, y0) Qy (x0,y0)
a = pX (x0,y0) + Qy (x0,y0)-
(4)
(5)
Характер состояния равновесия определяется знаком величин A, a и a2 — 4A. Учитывая, что в точках равновесия P(x0,y0) = 0 и Q(x0,y0) = 0, получаем следующие выражения для A, a2 — 4A и a.
A = A2 + (2kiy0 — li)(2k2x0 — і2), a2 — 4A = —4(2kiy0 — li)(2k2 x0 — ^),
(6)
a = 2A.
(7)
154 А.А. Лаптев. Аналитическое исследование математической модели...
Исключая из (6) 1\ и /2, получаем
y2 х2
А = kik2Xoyo - k\A+ k2Ay0,
(8)
22
a2 - 4А = -4№хоуо + 4kiAgo - 4k2AX0 + 4A2.
Или, исключая в (6) kl и k2,
А = -3A2 - 2/iA00 + 2/2A20 + hh,
l Xo 2 Vo l 2’
a2 - 4А = 16A2 + 8/iA00 - 8/2A00 - 4/1/2.
l xo 2 yo 12
(9)
Рассмотрим точку равновесия (0, 0). В этой точке А = A2 + /1 /2 > 0, a2 - 4А = —4/l/2 < 0, a = 2A. Следовательно, точка (0, 0) - простой (т.к. А = 0) фокус, устойчивый, когда A < 0, и неустойчивый, когда A > 0 [1, с.68-72]. Случай A = 0 разобран далее (см. стр. 156).
Из критерия Бендиксона1 следует, что на плоскости х, у при A = 0 не существует замкнутых контуров, составленных из траекторий, так как PX + Q'y = a = 2A (см. (7)).
A. Исследуем систему (2) при A > 0.
В этом случае имеем два, или три, или четыре состояния равновесия. Их количество зависит от положения парабол (3), которые пересекаются в 2-х, 3-х или 4-х точках (рис.1).
A.I. Пусть имеется 4 состояния равновесия (рис.2, слева). Точка (0,0) -неустойчивый фокус.
В точке Ml: х0 >/2/k2 >0, у0 <0, А<0 (см.(8)). Следовательно, точка Ml -седло [1, с.72].
Для точки M2: x0>/2/2k2>0, y0>/l/kl, А>0, a2 - 4 < 0 (см.(6)). Следовательно, точка М2 - неустойчивый фокус [1, с.72].
В точке M3: 0<х0</2/2k2, у0>/l/kl>0, a2 - 4>0 (см.(6)). Следовательно, точка М3 - седло или узел [1, с.71-72]. х
хКритерий Бендиксона. Если в некоторой односвязной области выражение Pc + QV не меняет знака и не равно нулю тождественно, то в этой области не существует замкнутых контуров, составленных из траекторий [1, с.120].
Математические структуры и моделирование. 2002. Вып. 10.
155
Рис. 2. Параболы (3)
Рис. 3. Параболы (3)
Для определения типа состояния равновесия М3 воспользуемся теоремой Пуанкаре2.
Ниже докажем (стр.158), что на экваторе сферы Пуанкаре лежит один узел. Поэтому N' = 1 и C' = 0. В конечной части плоскости имеем два фокуса и седло (точки O, М2 и Mi). Следовательно, точка М3 может быть только седлом.
A.II. Если в точкеМ2: 0<x0<l2/2k2 (рис.3, слева), то точка (0, 0) - неустойчивый фокус, а точка М1 - седло. Точки М2 и М3 узлы или седла, т.к. о2 — 4Д>0 (см.(6)). Из теоремы Пуанкаре следует, что М2 и М3 неустойчивый узел и седло соответственно.
Итак, если в конечной части фазовой плоскости четыре состояния равновесия, то два из них (М1 и М3) - седла, одно (0(0, 0)) - неустойчивый фокус и одно (М2)- неустойчивый фокус или узел.
A.III. Рассмотрим случай, когда имеются два состояния равновесия (рис.2, справа).
Точка (0, 0) - неустойчивый фокус. В точке М1: х0 >l2/k2 > 0, y0 < 0, Д < 0. Следовательно, точка М1 - седло [1, с.72].
Поэтому, если в конечной части фазовой плоскости два состояния равновесия, то одно из них (М1) - седло, а другое (0(0, 0)) - неустойчивый фокус.
2Теорема Пуанкаре. Если N, Nf и C - соответственно число узлов, фокусов и седел в конечной части фазовой плоскости, а N' и C' - число узлов и седел, лежащих на экваторе сферы Пуанкаре (считая точки, расположенные на концах диаметра, за одну точку (рис. 9.12)), то имеет место соотношение N + Nf + N' = C + C + 1 [1, с.124].
156 А.А. Лаптев. Аналитическое исследование математической модели...
Рис. 4. Две пары параллельных прямых
B. Если A < 0, то исследование системы (2) проводится аналогично (рис.3, справа). При A < 0 положения равновесия (фокусы и узел) устойчивы, так как a = 2A < 0 [1, с.68-72].
C. Исследуем систему (2) при A = 0.
В этом случае параболы (3) вырождаются в две пары параллельных прямых (рис.4)
п l2 п l1 x = °, x = т", У = °, У = у•
k2 ki
Система (2) при A = 0 имеет вид
Г # = hy — kiy2,
\ dt = -І2 X + k2 X2.
(10)
Таким образом, имеем консервативную (гамильтонову) систему (10) с аналитическим интегралом [1, с.136]
дН
ду
hy — hy2
l1 2 k1 3
to| 1 У y
dx dH
dt, dx
h
k2
2 3
l2x — k2x2
2
3
dy
dt
В консервативной системе возможны только простые состояния равновесия типа центра или седла; замкнутые траектории не являются изолированными, а заполняют целые области [1, с.136].
В точках М1 и M3: А =—l1l2 <0 (см.(8)). Следовательно, точки М1 и М3 -седла. А точки М2 и 0(0, 0) - центры (рис.4).
D. Рассмотрим случай, когда в особой точке А=0.
Если А=0, то (см.(6))
A2 + (2kiyo — li )(2k2Xo — I2) = 0. (11)
Точки, удовлетворяющие уравнению (11), располагаются в полуполосах: {0< x0<l2/2k2, y>l1/2k1} и {0<y0<l1/2k1, x>l2/2k2} (рис.5, слева). Таким образом, это может быть только точка М3 (рис.2-3). В этой точке
2k2xo — і2 —A
2k1 Уо — l1
(12)
A
Математические структуры и моделирование. 2002. Вып. 10.
157
Рис. 5. Параболы (3)
Левая и правая часть равенства (12) - угловые коэффициенты касательных к параболам (3). Следовательно, касательные параллельны, а точки касания лежат в интервале: 0 <x0 <l2/2k2 (при A > 0; случай A < 0 рассматривается аналогично). Поэтому (рис.2-3) эти касательные совпадают, т.е. в данном случае параболы касаются (в точке M2), и имеются три точки равновесия в конечной части плоскости (рис.5, справа).
Точка (0, 0) - неустойчивый фокус. В My xo > l2/k2 > 0, y0 < 0, А < 0. Следовательно, точка Mi - седло [1, с.72].
Точка М2 - сложное состояние равновесия (А = 0), которое, по-видимому, является седло-узлом [1, с.90-92].
2.2. Исследование бесконечно удаленных точек
Проведем исследование «бесконечно удаленных» частей плоскости с помощью сферы Пуанкаре [1, с.112-115]. Это сфера единичного радиуса, касающаяся плоскости x,y в начале координат. Каждой точке x плоскости ставится в соответствие две точки сферы, лежащие на прямой, проходящей через центр сферы и точку x. На экватор (большой круг, параллельный плоскости x, у) отображаются бесконечно удаленные точки. При этом интегральные кривые плоскости переходят в соответственные кривые сферы, а седла, узлы и фокусы сохраняют тот же вид.
Преобразование x = 1/z, у = u/z позволяет изучить особые точки на экваторе сферы Пуанкаре, за исключением концов оси у (рис.6).
Исследуемая система (2) содержит в правой части многочлены одинаковой степени: P(x,y) = Ax + l1y — k1y2 и Q(x,y) = Ax — l2x + k2x2. Координаты особых точек на экваторе находим из уравнения
Qn(1,u) — uPn(1,u) = 0,
где Qn и Pn - члены наивысшей степени в Q и P. Получаем: k2 + k1u3 Следовательно, координаты особых точек
u
3
k2 ki'
По критерию Пуанкаре3 находим, что найденная точка - узел.
0.
158 А.А. Лаптев. Аналитическое исследование математической модели...
Рис. 6. Сфера Пуанкаре
Для исследования «концов оси у» полагаем x = v/z, y = 1 /z и рассматриваем точку (0, 0). Получаем
dz -Az2 — k2 v2 z + l2 vz2
dv —ki + li z — k2 v3 + l2 v2 z
Точка (0, 0) не является особой. Следовательно, «концы оси у» не могут быть особыми точками.
Таким образом, мы выяснили, что на экваторе сферы Пуанкаре имеются две диаметрально противоположные точки - два узла.
3. Фазовые портреты и исследование системы на наличие бифуркаций
Проследим за сменой фазовых портретов в зависимости от изменения параметра A системы (2). Коэффициенты li, l2, ki и k2 зафиксированы. Фазовые портреты на рисунках 7-10 построены4 при li = 1, l2 = 1.5, ki = 1.3 и k2 = 1.1. Сплошные линии на рис.7-10 - траектории решений (положительное направление указано стрелками). Короткие штриховые линии показывают векторное
3Критерий Пуанкаре. Пусть Q и P одинаковой степени. Простая особая точка z = 0, и = щ будет седлом, если при изменении и от щ — £ до щ + є выражение
Qn (1, и)
—---:-и
Pn (1,и)
переходит от отрицательных значений к положительным, и узлом, если указанное выражение переходит от положительных значений к отрицательным [1, с.115].
4Для численного решения использовался метод Рунге-Кутта четвертого порядка. Вычисления и построение графиков выполнены с помощью пакета "DynaSys" (http://www.sci.wsu.edu/idea/DSDoc/). Результаты совпадают с решениями, полученными с помощью других пакетов программ. Для проверки результатов мы использовали пакет "DESIR", созданный в Ростовском государственном университете (http://www.math.rsu.ru/mexmat/kvm/MME/), и пакет "ODE", созданный в Московском энергетическом институте (http://www.exponenta.ru/soft/others/ode/ode.asp).
Математические структуры и моделирование. 2002. Вып. 10.
159
Рис. 7. А=-0, 81 и A=-0,13758
Рис. 8. А=—0,1 и A=0
Рис. 10. А=0, 81 и А=-0,13758
160 А.А. Лаптев. Аналитическое исследование математической модели...
Рис. 11. A=-0,13758 и A=-0,1
поле решений.
При больших по модулю отрицательных значениях параметра A в конечной части плоскости два состояния равновесия: т.0 - устойчивый фокус (A< 0) и т.М3 - седло (рис.7, слева). При увеличении A появляется третья особая точка М2 - седло-узел5 (рис.7, справа и увеличенный фрагмент на рис.10, справа), которая распадается на седло М1 и устойчивый (A<0) узел M2. При дальнейшем увеличении A узел становится устойчивым фокусом М2 (рис.8, слева).
При A=0 устойчивые фокусы (т.0 и т.М2) становятся центрами (случай C на стр.156). Вокруг этих точек расположены циклы (рис.8, справа).
При дальнейшем увеличении A фокусы (т.0 и т.М2) становятся неустойчивыми, т.к. A>0 (рис.9, слева). Этот случай соответствует пункту A.I, который разобран на стр.154. Далее фокус М2 переходит в неустойчивый узел (пункт A.II, стр.155), который сливается с седлом М3 (рис.9, справа), превращаясь в седло-узел (пункт D на стр.156). При увеличении A седло-узел распадается, и в конечной части плоскости остаются две точки равновесия (пункт A.III, стр.155): т.0 - неустойчивый фокус и т.Мі - седло (рис.10, слева). При дальнейшем увеличении параметра A характер решения на меняется.
Таким образом, мы исследовали изменение качественного поведения фазовых портретов при изменении параметра вдоль числовой оси и установили, что имеются бифуркационные значения параметра6:
1) при переходе через точку A = 0 происходит смена устойчивости фокусов (т.0 и т.М2) без рождения предельного цикла (бифуркационному значению параметра соответствует консервативная система) [1, с. 198].
2) при появлении седла-узла и его распадении на седло и узел (и наоборот, слияние узла и седла в седло-узел, который исчезает) при прохождении параметра A через те значения, при которых Д=0 (пункт A.II, стр.155) [1, с.193-194].
На рис. 11-12 представлены изображения фазовых портретов на полусфере
5В этой точке Д=0, см. стр.156.
6Это значения параметра, соответствующие негрубой системе, при переходе через которые происходит смена качественной картины фазовых портретов [1, с. 192].
Математические структуры и моделирование. 2002. Вып. 10.
161
Пуанкаре. Внешняя окружность на рисунках - экватор сферы Пуанкаре, на который отображаются бесконечно удаленные точки. Линиями изображены траектории решений. Стрелки показывают движение по траекториям в направлении, соответствующем положительному течению времени. Жирными линиями указаны сепаратрисы седел. Рисунки дают полное представление о характере поведения решений во всей плоскости и в окрестности точек равновесия.
Итак, мы установили количество состояний равновесия исследуемой системы, их тип и расположение на плоскости. Исследовано поведение решений на бесконечности. Выписаны условия на коэффициенты и параметры, которые характеризуют смену качественной картины решения. Построены фазовые портреты.
ЛИТЕРАТУРА
1. Баутин Н.Н., Леонтович Е.А. Методы и приемы качественного исследования динамических систем на плоскости. М.: Наука, 1976. 496 с.
2. Гуц А.К., Коробицын В.В., Лаптев А.А., Паутова Л.А., Фролова Ю.В. Социальные системы: формализация и компьютерное моделирование. Учебное пособие. Омск: Омск. гос. ун-т, 2000. 160 с.
3. Математическое моделирование глобальных социальных процессов: Отчет о НИР(закл.) / Руководитель А.А.Лаптев; №ГР 01.2.00101750; Инв. №02200201162. - Омск, ОмГУ, 2001. 24с.