Добавил:
tg: @Yr66gi4 Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

экзамен / Билеты по синергетике

.pdf
Скачиваний:
0
Добавлен:
08.01.2026
Размер:
13.78 Mб
Скачать

С -

неустойчивые точки.

На

участке 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