С - |
неустойчивые точки. |
На |
участке OА F X 0 , |
следовательно |
||||
|
dX |
0 |
и величина X t увеличивается со временем. На участке AB |
|||||
|
dt |
|||||||
|
|
|
|
|
|
|
|
|
|
F X 0 , |
следовательно |
|
dX |
0 и величина |
X t |
уменьшается со |
|
|
|
dt |
||||||
|
|
|
|
|
|
|
|
|
временем.
Рис. DSE_2. Поведение функции |
V X , |
для которой |
F X |
dV |
, а |
|||||||||
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
dx |
|
функция |
F X |
изображена на рис. |
DSE_1. На участке |
OA |
dV |
0 , на |
||||||||
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
dx |
|||
участке |
АB |
dV |
0 , на участке |
BC |
|
dV |
0 . Точка |
|
A – устойчивая, |
|||||
|
|
|
||||||||||||
|
|
dx |
|
|
|
dx |
|
|
|
|
|
|
||
точка, а точки |
B и С - неустойчивые. |
|
|
|
|
|
|
|
|
|
||||
Аналогичные рассуждения |
можно |
сделать и для других точек. |
||||||||||||
Следовательно A − устойчивая точка, а B – неустойчивая |
|
точка, С – также |
||||||||||||
неустойчивая точка. Нетрудно видеть, что тогда устойчивым особым точкам
отвечают |
|
точки |
минимума |
функции |
V X , |
а неустойчивым точкам |
|||||||||||
соответствуют точки максимума функции V X |
или точки перегиба (рис. |
||||||||||||||||
DSE_2). |
|
|
Для |
|
подтверждения |
этого |
факта |
рассмотрим |
функцию |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
2 . |
График |
этой |
функции для |
значений |
|
F |
|
X |
|
|
|
X 1 |
X 2 |
|
X 3 |
||||||||
аргумента |
0.7< X <3.4 |
на рис. |
DSE_2A. |
|
|
||||||||||||
14
Рис. DSE_2A. График функции |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|||||
F |
|
X |
|
|
|
X 1 |
X 2 |
|
|
X 3 |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X |
|
|
|
|
|
|
Первообразную V |
X |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
определим |
|
в виде V |
|
X |
|
|
|
|
F |
|
X dX . Явный |
вид |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
полинома V X X 5 / 5 9X 4 |
/ 4 29X 3 / 3 39X 2 |
/ 2 18X |
приведен |
на |
||||||||||||||||||
рис. DSE_2B.
Рис. DSE_2B. График функции V X для интервала аргумента 0.7< X <3.4.
15
1.2.2. Динамические системы с двумя степенями свободы
Проанализируем качественно возможные сценарии поведения дифференциальных систем уравнений второго порядка. Пусть движение нашей системы описывается уравнениями типа
dX1 t F1 X1; X 2 dt
dX 2 |
t |
F2 |
X1; X 2 |
(DSE.3) |
|
dt |
|||||
|
|
|
|||
Особые точки этой системы, т. е. равновесные состояния определяются из системы уравнений
F1 X1 s ; X 2s 0
F2 X1( s) ; X 2( s) 0 |
(DSE.4) |
и могут быть на плоскости X1; X 2 |
6 типов (см. рис.DSE_3). |
Для анализа особых точек на устойчивость необходимо линеаризовать систему уравнений (DSE.3) относительно малых отклонений от положений равновесия:
x1 t X1 t X1(s) ; x2 t X2 t X2(s) . (DSE.5)
В результате линеаризации системы двух нелинейных уравнений получаем
|
|
dx1 |
a |
|
x a |
|
x |
|
|
||||
|
|
|
|
|
|
|
|
||||||
|
|
dt |
11 |
|
1 |
12 |
2 |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|||
|
dx2 |
a |
|
x |
|
a |
|
x |
|
|
(DSE.6) |
||
|
|
21 |
|
22 |
2 |
|
|||||||
|
dt |
1 |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
||||
Коэффициентами aik в уравнениях (DSE.6) |
представляют собой значения |
||||||||||||
производных от правых частей |
|
(DSE.3), |
взятые при |
X1 X1( s) и |
|||||||||
X 2 X 2( s) : |
|
|
|
|
|
|
|
|
|
||||
a |
F1 X1; X 2 |
|
| ( s ) |
( s ) ; |
|
|
|||
11 |
X1 |
|
X1 |
, X 2 |
|
|
|||
|
|
|
a |
F1 X1; X 2 |
|
| ( s ) |
|
|
|
( s ) |
||
12 |
X 2 |
|
X1 |
, X 2 |
|
|
|||
|
|
|
16
a21 |
F2 X1; X 2 |
|
| |
|
( s ) ; |
a22 |
|
F2 X1; X 2 |
| |
|
|
|
( s ) . |
|
|||
|
X1 |
|
( s ) |
X 2 |
( s ) |
|
|
||||||||||
|
|
|
|
|
|
X1 |
, X 2 |
|
|
|
X1 |
|
, X 2 |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
В уравнении (DSE.6) |
квадратичные слагаемые |
x |
2 , |
|
x2 , их произведения |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
2 |
|
|
|
x1x2 , а также слагаемые более высоких степеней |
|
отбрасываются. |
Ищем |
|||||||||||||
решение |
этой |
|
|
системы |
в |
следующем |
виде: |
xi Ci exp pt , |
|||||||||
|
dxi |
pCi exp pt . |
Сокращая |
все |
слагаемые на |
|
exp pt получаем, что |
||||||||||
|
|
|
|||||||||||||||
|
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
система |
(DSE.6) |
|
является линейной и однородной для |
констант С1 |
и С2 . |
||||||||||||
Из условия существования нетривиального решения однородной системы алгебраических уравнений получаем характеристическое уравнение
|
a11 p; |
a12 |
|
0 . |
(DSE.7) |
||||
|
|
||||||||
|
a21 |
|
; a22 p |
|
|||||
Уравнение для нахождения p |
имеет вид |
|
|
||||||
p2 a11 a22 p a11a22 a12a21 0. |
(DSE.8) |
||||||||
Введем обозначения a11 a22 ; |
a11a22 a12a21 . Корни |
p1 и p2 |
|||||||
определяются из решения квадратного уравнения |
|
||||||||
|
|
|
|
|
|
|
|
||
|
p |
|
|
|
2 4 |
. |
|
||
|
|
|
|
|
|
|
|||
1,2 |
|
|
2 |
|
|
|
|
||
|
|
|
|
|
|
|
|
||
Решение системы (DSE.6) есть линейная комбинация экспонент: x1 t C11 exp p1t C12 exp p2t
x2 t C21 exp p1t C22 exp p2t . (DSE.9)
Свойство решений (DSE.9) зависит от вида корней или коэффициентов
характеристического уравнения. В зависимости от них существует 4 типа |
||
неподвижных точек |
X1( s) ; X 2( s) |
и шесть возможных типов решения, |
которые будут рассмотрены в дальнейшем. Особыми (неподвижными)
точками являются: узлы, фокусы, центры, седла.
Если |
хотя бы один корень имеет |
вещественную положительную часть |
|
Re p1,2 0 |
, то |
стационарное состояние неустойчиво и решения |
|
экспоненциально |
возрастают со |
временем. При отрицательных |
|
|
|
17 |
|
вещественных частях |
Re p1,2 0 стационарное состояние |
асимптотически |
|
устойчиво и решение экспоненциально затухает со временем. |
|||
Если p1 |
и p2 |
являются комплексными числами с отрицательной |
|
вещественной |
частью |
Re p1,2 0 , то решения будут |
носить характер |
колебаний с экспоненциальным затуханием амплитуды. Если вещественная
часть комплексных чисел Re p1,2 0 , |
то фокус становится неустойчивым. |
||
Если |
p1 и p2 |
являются чисто |
мнимыми, но неподвижной точкой |
является устойчивый центр. |
|
||
Если |
p1 и p2 |
вещественные |
числа, имеющие разные знаки, то |
неподвижная особая точка называется седловой. К седловой точке примыкают четыре траектории, которые представляют собой две интегральные кривые. С увеличение времени t по двум траекториям xA t
и xB t изображающая точка приближается к особой точке. В тоже самое время по двум другим траекториям xC t и xD t изображающая точка удаляется от особой седловой точки. Таким образом, седловая точка является существенно неустойчивой точкой. Между указанными траекториями xA t , xB t , xC t и xD t расположены четыре области, содержащие континуум
интегральных кривых гиперболической формы. Эти кривые не достигают особой седловой точки.
Рис. DSE_3. |
Типы особых точек для системы (DSE.6). Величина |
a11 a22 ; |
a11a22 a12a21 . Рисунок приведен в книге Г.Ю.Ризниченко. |
Лекции по математическим моделям в биологии.
18
1 |
|
|
|
|
|
2 |
|
|
|
3 |
|
|
|
|
|
|
|
|
|
||||||||||
Устойчивый узел |
|
(1) |
Неустойчивый узел (2) |
Устойчивый фокус (3) |
||||||||||
|
|
|
|
|
|
|
|
|
||||||
p1 0 |
и |
p2 0 |
p1 0 |
и |
p2 0 |
p1 |
и p2 |
комплексные числа: |
||||||
действительные |
|
|
|
действительные |
|
Re p1 0 и |
Re p2 0 . |
|||||||
отрицательные числа. |
положительные числа. |
|||||||||||||
Периодическое |
затухающее |
|||||||||||||
Апериодические |
|
|
|
|
|
|
||||||||
|
|
|
Апериодическое |
|
колебание, |
приближающееся |
||||||||
затухающие |
колебания. |
|
||||||||||||
самовозбуждающее |
|
к |
положению |
равновесия. |
||||||||||
Область |
|
|
|
|
|
|
||||||||
|
движение. Область |
|
Область |
определяется |
||||||||||
|
|
|
|
|||||||||||
определяется |
условиями: |
определяется |
условиями: |
условиями: |
0 ; |
0; |
||||||||
0; |
|
0; |
||||||||||||
|
0; |
|
0; |
2 4 0. |
|
|
||||||||
2 |
|
|
|
|
|
|
|
|
||||||
4 0. |
|
|
|
|
2 |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
4 0. |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
===========================================================
|
4 |
|
|
|
|
5 |
|
|
|
6 |
|
|
|
|
|
|
|
|
|||||||
Неустойчивый фокус (4) |
Седло |
|
(Неустойчивая |
Центр |
(устойчивая |
|||||||
|
|
|
|
точка) |
(5) |
|
|
|
точка) |
(6) |
|
|
|
|
|
|
|
|
|
||||||
p1 и |
p2 |
комплексные |
p1 и |
p2 |
действительные |
p1 и p2 чисто мнимые |
||||||
числа: |
Re p1 0 |
и |
числа, |
имеющие |
|
различные |
числа: |
|
Re p1 0; |
|||
знаки p1 p2 0 . |
|
|
|
|||||||||
Re p2 0 . |
|
|
При малом |
Re p2 |
0 ; Im p1,2 0 . |
|||||||
|
|
отклонении |
от |
|
положения |
|||||||
Непериодические |
|
|
Незатухающие |
|||||||||
|
равновесия |
|
|
система |
||||||||
самовозбуждаемые |
|
|
|
колебания. |
Область |
|||||||
|
начинает |
от |
него |
|||||||||
колебания с |
нарастающей |
|
|
|
||||||||
удалятьсяв |
неустойчивом |
определяется |
||||||||||
амплитудой. |
Область |
|||||||||||
направлении. Область |
|
|
|
|||||||||
|
|
определяется |
неравенством |
0 ; |
||||||||
|
0 . |
|
|
|||||||||
неравенством |
0 ; |
определяется |
неравенством |
|
|
|||||||
0 . |
|
|
|
|
|
|
|
|||||
0; 2 4 0. |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
19 |
|
|
|
|
|
||
Рис. DSE_3. Шесть типов стационарных состояний в системе в двумя степенями свободы
Существует 6 возможностей, которые соответствуют 6 рисункам DSE_3. 1) p1 и p2 действительные отрицательные числа (устойчивый узел).
Совершенно устойчивый режим узлового типа. Система совершает апериодически затухающее движение, приближаясь к положению равновесия.
2) p1 и p2 − действительные положительные числа (неустойчивый узел).
Неустойчивый режим узлового типа. Система удаляется от положения равновесия, совершая апериодическое самовозбуждающее движение. Оно при определенных условиях может перейти в предельный цикл.
3) p1 и p2 − комплексные числа с отрицательной действительной частью
(устойчивый фокус). Устойчивый фокус. Система совершает периодические затухающие колебания и асимптотически приближается к положению равновесия.
4) p1 и p2 − комплексные числа с положительной действительной частью
(неустойчивый фокус). Неустойчивый фокус. Этот тип движения соответствует периодическим колебаниям с нарастающей амплитудой. В таком случае говорят о режиме самовозбуждающихся колебаний. В реальных системах нарастание амплитуды обычно ограничено физическими процессами, так что возникает устойчивое периодическое движение с постоянной во времени амплитудой. Фазовая траектория при этом накручивается на предельный цикл, соответствующий устойчивым автоколебаниям.
5) p1 и p2 − действительные числа, имеющие разные знаки (неустойчивое седло).
Неустойчивый режим типа седла. При малом отклонении от положения равновесия система начинает от него удалятьсяв неустойчивом направлении. В случае, если решения должны быть конечными, всегда существует два новых стабильных равновесных состояния, к которым приближается система. В свою очередь существование двух устойчивых состояний позволяет создать систему типа тригера. Он может находиться в одном из двух своих устойчивых состояний A или B бесконечно долго. Внешним возбуждением достаточной амплитуды можно переводить его из одного устойчивого состояния в другое. Схематически фазовый портрет тригера представлен на рис. DSE_4.
6) p1 и p2 − чисто мнимые числа (устойчивый центр). Устойчивый режим
центрального типа. Поскольку в этом случае оба корня чисто мнимые, наблюдаются незатухающие колебания и фазовая траектория системы представляет собой эллипс. Типичный пример такой системы − одномерный гармонический осциллятор без трения.
20
Рис. DSE_4: Фазовый портрет триггера с двумя устойчивыми состояниями
A и B.
Для анализа устойчивости необязательно знать явное решение характеристического уравнения. По критерию Гурвица асимптотически устойчивый режим при выполнении неравенств
Re p1 0 и Re p2 0 (DSE.10)
осуществляется, когда все коэффициенты характеристического уравнения (DSE.8) положительны
a11 a22 |
0, |
det |
aij |
a11a22 |
a12a21 0. |
|
|
(DSE.11) |
||
Критерий Гурвица имеет следующее объяснение (рис. DSE_3). Для |
||||||||||
отрицательного |
значения |
0 |
(DSE.8) |
реальная |
часть Re p1,2 0 , |
|||||
|
|
p1,2 |
|
/ 2, |
|
|||||
выражающегося |
соотношением |
2 4 |
также |
|||||||
отрицательна, |
независимо |
от |
знака |
2 4 . |
|
Следовательно, |
||||
линеаризованные решения системы (DSE.6) |
xi Ci exp pt характеризуют |
|||||||||
устойчивый режим.
Если действительные части одного или обоих корней характеристического уравнения равны нулю, то для анализа состояния равновесия необходимо рассматривать члены более высокого порядка малости в разложении в ряд Тейлора правых частей уравнений (DSE.6).
Напомним, что коэффициенты линейной системы aij вычислялись, как соответствующие частные производные правых частей исходных
21
нелинейных уравнений, взятых в стационарном состоянии. Следовательно aij являются комбинациями параметров некой исходной модели. При
изменении какого-либо параметра исходной нелинейной модели будут меняться и коэффициенты соответствующей линейной системы, следовательно, будут меняться величины σ, ∆. При переходе через оси координат σ, ∆ (рис. DSE_3), при которых изменяется знак σ и ∆, характер фазового портрета качественно меняется. Поэтому такие границы называются бифуркационными — по разные стороны от границы система имеет два качественно различных фазовых портрета и, соответственно, два разных типа поведения.
В науке синергетика, которую мы будем изучать в дальнейшем, под термином «бифуркация» мы будем понимать качественную перестройку свойств системы, которая сопровождается изменением ее симметрии.
Нужно отметить, что бифуркационными являются только переходы через оси координат и , поскольку невозможно постепенным непрерывным изменением фазового портрета перейти, например, от узла к седлу или от неустойчивого фокуса к устойчивому. Переходы устойчивый узел – устойчивый фокус (неустойчивый узел – неустойчивый фокус) не являются бифуркационными, поскольку можно постепенным непрерывным изменением фазового портрета перейти от узла к фокусу.
1.2.3. Анализ двух нелинейных дифференциальных уравнений
В качестве примера рассмотрим следующую задачу
dXdt1 1 1 X1 X12 X 2
|
dX 2 |
X |
|
X 2 X |
|
|
(DSE.12) |
||
|
|
1 |
2 |
|
|||||
|
dt |
|
|
1 |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
В нашем случае в |
системе |
уравнений |
(DSE.3) функции F1 X1; X2 и |
||||||
F2 X1; X2 , стоящие в правых частях системы (DSE.12), имеют вид: |
|||||||||
F1 X1; X2 1 1 X1 X12 X2 ; |
F2 X1; X2 X1 X12 X2 |
||||||||
Система (DSE.12) имеет единственный параметр , который может |
|||||||||
принимать любые |
значения. |
|
Этот параметр будет определять характер |
||||||
получающегося решения для зависимостей X1 t и |
X 2 t . |
||||||||
22
Для простоты предположим, что величина 0 . Отметим, что в функции F1 X1; X2 есть два положительных слагаемых 1 X12 X 2 , которые называются слагаемыми «прихода» и одно слагаемое «ухода» 1 X1 , имеющее знак минус. Если бы в первом дифференциальном уравнении были
бы только положительные слагаемые, то в этом случае производная |
dX1 |
0 , |
|||||
|
|||||||
|
|
|
|
|
|
dt |
|
и переменная |
X1 увеличивалась бы со временем. |
Если бы в функции |
|||||
F1 X1; X2 было бы только отрицательное слагаемое |
1 X1 , то это |
||||||
означало, что |
dX1 |
0 и переменная X |
|
увеличивалась бы со временем. |
|||
|
1 |
||||||
|
dt |
|
|
|
|
||
|
|
|
|
|
|
||
Точно такие же рассуждения о положительных слагаемых «прихода» и
отрицательных слагаемых «ухода» можно сделать и для |
второй |
функции |
||||||
F2 X1; X2 . Найдем стационарные решения |
X1( s) ; X 2( s) |
системы (DSE.12), |
||||||
удовлетворяющие условиям F1 X1( s) ; X2( s) 0 ; F2 |
X1( s) ; X 2(s) 0 . Разберем |
|||||||
вначале |
уравнение |
для |
|
|
второй |
|
функции |
|
F2 X1(s) ; X2(s) X1( s) X1( s) X 2(s) |
0 . |
Равенство |
нулю |
функции |
||||
F2 X1( s) ; X 2( s) 0 |
дает два решения: первое решение |
X1( s) 0 |
и второе |
|||||
решение X1(s) X2(s) |
. Если мы первое решение |
X1( s) 0 подставим в |
||||||
первую |
функцию |
F1 X1( s) ; X 2( s) , |
то |
мы |
получим |
F1 X1(s) ; X 2( s) 1. |
||
Следовательно первое решение X1( s) 0 не подходит на роль стационарного
решения. Подставляя второе решение |
X1(s) X2(s) в первую функцию: |
||
F1 X1(s) ; X2(s) 1 1 X1(s) X1(s) X1(s) X2( s) 1 X1( s) 1 1 X1( s) 0 |
|||
Следовательно X1(s) |
1. Учитывая второе решение X1(s) X2(s) находим, |
||
что X 2(s) . Таким образом, для того, |
чтобы одновременно две функции |
||
были бы равны нулю |
F1 X1( s) ; X2( s) 0 |
и F2 X1( s) ; X 2( s) 0 необходимо |
|
чтобы выполнялись два условия |
|
||
X1(s) |
1; |
X 2(s) . |
(DSE.13) |
23
