Синергетика_Глава_1_01_Системы_дифференциальных_уравнений
.pdf
Глава 1. Математическое введение
1.1. Системы нелинейных дифференциальных уравнений. Пространство состояний или фазовое пространство. Особые точки и их классификация. Условия устойчивости. Узел, фокус, седло, центр, предельный цикл
Основная идея синергетики состоит в том, что в открытых системах, далеких от термодинамического равновесия, возможно ”спонтанное” нарушение симметрии и образование пространственных структур или периодических колебаний во времени. Связано это математически с тем, что при больших отклонениях от равновесия системы, как правило, ведут себя нелинейно, т. е. описываются нелинейными уравнениями. Поэтому в этой лекции мы приведем некоторые понятия и результаты из теории нелинейных систем.
Для этого рассмотрим автономную систему, состояние которой
характеризуется переменными |
|
X {X1, X2 ,..Xn} |
(VDE.1) |
Множество всех возможных физических пространством состояний, или фазовым переменные удовлетворяют некоторой уравнений первого порядка
состояний системы называется пространством. Пусть эти системе дифференциальных
|
dXi |
Fi X , B |
|
(VDE.2) |
|
dt |
|
||
|
|
|
|
|
где индекс i пробегает |
значения |
i 1,2,...n , а через |
B b1,b2 ,...bl |
|
обозначены, как говорят, управляющие параметры системы. Они могут |
|
описывать как внутренние свойства системы, так и внешние условия. |
|
Автономность системы означает, что функции Fi X , B |
не зависят явно от |
времени t . Если такая зависимость существует, то |
для перехода к |
автономной системе |
необходимо ввести |
дополнительную переменную |
||||
X |
|
t и еще одно уравнение |
dX n 1 |
1. |
Очевидно, что для того, чтобы |
|
n 1 |
|
|||||
|
|
|
dt |
|
||
|
|
|
|
|
||
найти решение системы уравнений (VDE.2), надо еще задать начальные |
||||||
условия, т. е. определить значения функций |
Xi t в некоторый начальный |
|||||
момент времени t0 . |
Обычно таким начальным моментом времени является |
|||||
нулевой момент времени t0 0 . |
|
|||||
|
|
|
4 |
|
|
|
Мы можем рассматривать решение системы уравнений (VDE.2)
|
|
|
|
{X |
1 |
|
|
|
|
2 |
|
|
|
n |
|
|
X |
|
t |
|
|
|
t |
|
, X |
|
t |
|
,..X |
|
t } |
(VDE.3) |
|
(при некотором начальном условии) |
в заданный момент времени |
t как |
||||||||||||||
некоторую точку в n мерном фазовом пространстве системы. С изменением времени t эта точка будет описывать некоторую траекторию рис. VDE_1. Взяв другое начальное условие, мы получим другую траекторию и т. д.
Рис. VDE_1: Траектория системы в трехмерном фазовом пространстве.
Если уравнения в системе (VDE.2) независимы, то при заданном начальном
условии
|
|
|
|
{X |
1 |
|
|
|
|
2 |
|
|
|
|
n |
|
|
X |
|
0 |
|
|
|
0 |
|
, X |
|
|
0 |
|
,..X |
|
0 } |
(VDE.4) |
у этой системы существует единственное решение, удовлетворяющее этим начальным условиям. Среди всех возможных решений системы (VDE.2) важную роль играют так называемые стационарные решения, т. е. решения, не зависящие от времени.
Условие существования стационарных решений системы (VDE.2) имеет вид
|
dXi |
0 ; |
Fi X , B 0 ; |
i 1,2,...n |
(VDE.5) |
||
|
dt |
||||||
|
|
|
|
|
|
|
|
т. е. стационарные решения |
X X (s) |
получаются в результате решения |
|||||
системы |
алгебраических |
уравнений. |
Стационарным |
состояниям |
|||
соответствуют |
фиксированные во времени |
точки фазового |
пространства |
||||
X (s) , которые в теории дифференциальных уравнений называются особыми точками, в то время как все прочие точки носят название обыкновенных точек (регулярных). Если любую из этих особых точек взять в качестве начального состояния, то система будет находиться в них неограниченно долго. Это так называемые состояния равновесия.
5
Бывают, однако, устойчивые и неустойчивые состояния равновесия. Устойчивым называется такое состояние, когда при бесконечно малом отклонении от него система вновь в него возвращается. Неустойчивая особая точка −это соответственно такая точка, при бесконечно малом отклонении от
которой система удаляется от нее навсегда. Если некоторое состояние X (s) является устойчивым, то, по Ляпунову, критерий асимптотической
устойчивости состояния X (s) имеет вид
X t X t X (s) 0 при t (VDE.6)
Вообще характерной особенностью систем нелинейных уравнений (VDE.2)
является то, что они в определенных областях пространства параметров B допускают существование нескольких решений (и необязательно стационарных). Тогда нахождение этих решений не означает еще решение задачи, так как некоторые из решений могут оказаться неустойчивыми. Поэтому в нелинейных задачах надо не только найти решение, но и исследовать его на устойчивость. Если же оказываются устойчивыми несколько решений, то тогда в фазовом пространстве системы надо найти области начальных значений, которые приводят к этим решениям, и определить границы этих областей. Рассмотрим ниже в качестве примера задачу на устойчивость стационарных состояния.
Пусть X (s) есть стационарное решение нашей системы уравнений (VDE.2)
|
|
|
dXi( s) |
Fi X ( s) , B 0 . |
|
(VDE.7) |
|||||
|
|
|
dt |
|
|||||||
|
|
|
|
|
|
|
|
|
|
||
Запишем искомое решение в виде суммы X t X (s) x t где |
x t − |
||||||||||
малое отклонение решения дифференциального уравнения (VDE.2). Разлагая |
|||||||||||
функции Fi X , B в ряд по x t |
и ограничиваясь линейными слагаемыми , |
||||||||||
получим уравнение для возмущения x t |
|
|
|
|
|||||||
|
dx |
n F |
|
|
|
|
|
|
|||
|
i |
|
|
|
i |
| ( s ) x |
|
a x |
|
(VDE.8) |
|
|
|
X |
|
j |
j |
||||||
|
dt |
|
X X |
ij |
|
||||||
|
j |
|
|
|
|
||||||
|
|
|
j 1 |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
||
Вправой части уравнений (VDE.8) использовано правило немого
суммирования Эйнштейна, согласно которому: если в выражении
6
встречается дважды повторяющийся индекс j 1,2,..n , то это означает, что
n
по этому немому индексу проведено суммирование aij x j aij x j .
j 1
По математической терминологии система (VDE.8) представляет собой систему обыкновенных линейных дифференциальных уравнений с
постоянными коэффициентам aij : |
|
|
|
|
a |
Fi |
| |
|
(VDE.9) |
|
( s ) |
|||
ij |
|
|
|
|
X j X X
Теория устойчивости для таких систем линейных дифференциальных уравнений (VDE.8) хорошо развита.
Выберем проект |
решения системы |
уравнений (VDE.8) |
в |
виде |
|||
xi t Ci |
exp pt , где |
Ci − амплитуда, |
а |
p − некоторый постоянный |
|||
параметр. Сокращая все слагаемые на |
множитель |
exp pt , получаем |
|||||
систему |
однородных |
алгебраических |
уравнений |
pCi aijC j |
для |
||
определения амплитуд |
C j . Подставляя |
коэффициент |
Ci C j ij |
в |
виде |
||
суммирования с дельта символом Кронекера ij , можно сделать вывод о том, что нетривиальное решение такой системы алгебраических уравнений для
амплитуд C j |
получается |
в случае, когда детерминант системы |
||||
алгебраических уравнений равен нулю |
|
|
||||
|
|
|
aij p ij |
|
0 . |
(VDE.10) |
|
|
|
|
|||
Таким образом, спектр возможных значений |
p определяется из |
|||||
характеристического алгебраического уравнения (VDE.10). |
||||||
Если среди решений этого алгебраического уравнения найдется хоть |
||||||
одно решение для которого выполняется условие |
Re p 0 , то это означает, |
|||||
что особая |
точка X (s) |
− неустойчивая, |
так |
малое возмущение |
||
xi t Ci exp pt , отсчитываемое от этой особой точки X t X (s) x t ,
экспоненциально возрастает со временем t . Для устойчивости особой точки необходимо, чтобы у всех корней p выполнялось бы условие Re p 0 .
Теперь рассмотрим топологию траекторий в n - мерном пространстве
7
состояний системы: X {X1, X2 ,..Xn}. Предположим, что наша система уравнений (VDE.2) имеет однозначное решение. Это означает, что через каждую обыкновенную (регулярную) точку нашего фазового пространства X t проходит одна и только одна траектория. Следовательно в
обыкновенных точках траектории пересекаться не могут. Для обыкновенных |
||||||||
точек системы дифференциальных уравнений (VSE.2) существует хотя бы |
||||||||
одно значение индекса i m , для которого |
Fm |
X , B 0 . Тогда, исключая |
||||||
время |
t их двух уравнений для индексов i |
и |
m получаем систему из |
|||||
n 1 |
дифференциальных уравнений, |
которые |
|
описывают |
траектории в |
|||
фазовом пространстве X : |
|
|
|
|
|
|
||
|
|
dXi |
|
Fi X , B |
. |
(VDE.11) |
||
|
|
dX m |
Fm X , B |
|||||
|
|
|
|
|
||||
Через обыкновенную точку в фазовом пространстве проходит одна и только одна траектория. Однако, в особых точках однозначность теряется и траектории могут пересекаться. Вследствие однозначности решения
никакая траектория, исходящая из некоторой (обыкновенной) точки X 0 в некоторый момент времени t0 , не может к конечному моменту времени достигнуть какой-нибудь особой точки.Это противоречило бы
единственности решения при заданном начальном условии. Однако это становится возможным при t → ∞. Таким образом, система достигает (или
покидает) стационарные состояния X (s) только в пределе t → ±∞. Не существует движений, ведущих за конечное время к стационарным состояниям. Траектории в фазовом пространстве вблизи устойчивых особых точек могут быть двух типов (или их комбинациями). Если устойчивая особая точка узел, то траектории вблизи нее образуют своеобразный ”ежик” рис. 2. Зависимость координат от времени представляет собой монотонно затухающие кривые.
Рис. VDE_2. Топология траекторий вблизи узла.
8
В том случае, если устойчивая особая точка фокус, то траектории приближаются к ней, закручиваясь вокруг нее по спирали − рис. 3: В этом случае временная зависимость Xi t носит колебательный и затухающий характер.
Рис. VDE_3: Топология траекторий вблизи фокуса.
Таким образом, устойчивые особые точки притягивают к себе траектории. Обратив время вспять, мы получим неустойчивые особые точки, которые отталкивают от себя траектории. В этом случае траектории не сходятся к этим точкам, а наоборот, расходятся от них. Устойчивые особые точки часто называют аттракторами, а неустойчивые − репеллерами.
Существует еще один важный тип особых точек. Это так называемые седловые точки. По одним направлениям они притягивают траектории, а по другим – отталкивают, как изображено на рис. VDE_4.
Рис. VDE_4: Топология траекторий вблизи седловой точки.
Аттракторами могут служить не только точки, но и целые траектории. Классическим примером может служить предельный цикл Пуанкаре C
рис. VDE_5.
9
Рис. VDE_5: Предельный цикл Пуанкаре.
На рис. VDE_5 изображена замкнутая в фазовом пространстве траектория, обладающая тем свойством, что все другие траектории в ее окрестности являются спиралями, закручивающимися в направлении C . Очевидно, что предельный цикл соответствует периодическому движению системы. Предельный цикл может также быть неустойчивым. В этом случае все траектории от него отталкиваются.
В дальнейшем мы будем исследовать еще одно такое состояние в фазовом пространстве, которое называется «центром». Центр представляет собой замкнутые траектории вокруг некоторой точки, которая будет называться центром.
Притягивающая траектория (аттрактор) может также иметь весьма причудливую форму рис. VDE_6. Она незамкнута и выглядит чрезвычайно запутанной, заполняя целиком некоторую область в фазовом пространстве, как правило, с дробной размерностью − фрактал.
Рис. VDE_6: Странный аттрактор.
Движение по такой траектории хаотично и очень чувствительно к начальным условиям. Такой аттрактор поэтому называют странным
аттрактором.
10
Странными аттракторами называются динамические системы, которые характеризуются случайным блужданием своих фазовых траекторий между несколькими неустойчивыми точками. При этом любая траектория спустя некоторое время путешествий в фазовом пространстве притягивается и «садится» на некоторую выделенную область в фазовом пространстве. Эти особые точки в фазовом пространстве и называются странными аттракторами.
Рис. VDE_7. Изображающая точка в фазовом пространстве движется по интересным траекториям, которые имеют области притяжения.
С траекториями странного аттрактора в фазовом пространстве можно ознакомиться на сайте https://nplus1.ru/material/2019/09/06/chaosreigns
11
