Динамический хаос (ИПИС, ФКС)
.pdfУстойчивость и неустойчивость траекторий
Как говорилось ранее, наличие хаотической динамики тесно связано с неустойчивостью, присущей фазовым траекториям системы. В качестве примера, на рис. 107 представлены временные зависимости одной из динамических переменных для модели Лоренца, задаваемой системой уравнений (110)
x |
|
x y |
|
|
|
|
|
y |
r x y xz . |
(110) |
|
|
|
|
|
z |
bz x y |
|
|
|
|
|
|
Рис. 107. Чувствительность решений уравнения Лоренца к
изменениям начальных условий в хаотическом режиме: r 28, |
||||
10 и b 8 3 . |
Начальные |
условия: |
x0 |
10, 485 , |
y0 10, 485 , z0 |
27 , где |
, и |
|
– случайные |
добавки порядка 0,01. |
|
|
|
Все эти зависимости z(t) получены решением одной и той же системы (110), демонстрирующей хаотическую динамику, при слегка отличающихся друг от друга начальных условиях. Из рис. 107 видно, что все эти зависимости практически неотличимы на начальном участке, но с течением времени различие между ними увеличивается и в итоге картина движения размазывается. Такое расхождение траекторий спустя длительное время обуславливается различиями в начальных условиях, несмотря на то, что они были малы. В этом случае говорят о неустойчивости фазовых траекторий системы.
Таким образом, как видно из рис. 107, движение системы Лоренца можно с хорошей точностью прогнозировать лишь на небольшом начальном интервале времени (в нашем случае до t 30 ). При t 30 описывать движение системы обычным способом, т.е. решая систему уравнений (110), не имеет смысла и следует использовать уже другие методы описания, например, статистические.
Связь между хаосом и неустойчивостью фазовых траекторий можно отследить и для других модельных систем, изученных нами ранее. Именно присутствие неустойчивости даёт возможность соединить, казалось бы, несоединимое – динамическую природу системы (предсказуемость) и хаос (непредсказуемость).
Устойчивость по Ляпунову
Рассмотрим динамическую систему, описываемую автономной системой дифференциальных уравнений
dx |
F(x) , |
(241) |
dt |
|
|
где x – N-мерный вектор. Рассмотрим одну из траекторий x(t) системы с началом в точке x0 . Выберем ещё одну траекторию y(t) с началом в произвольной точке y0 , близкой к x0 , так, чтобы в начальный момент времени траектории были близки. Если обе эти траектории остаются близкими и в любой последующий момент времени, то траектория x(t) называется устойчивой по Ляпунову.
Другими словами, траектория |
x(t) устойчива, если для любого |
||||||
малого |
|
0 существует такое |
0 , |
что для любой точки |
|||
старта из |
|
-окрестности x0 |
, т.е. при |
y0 x0 , имеем |
|||
для всех t 0 : y(t) x(t) |
|
, т.е. (см. рис. 108) |
|||||
0, 0, y0 , y0 |
x0 |
, t 0 : |
y(t) x(t) , (242) |
||||
где |
|
|
|
|
|
|
|
|
|
x |
x2 |
x2 |
x2 |
||
|
|
|
1 |
|
2 |
|
N |
– длина вектора x , а
y x ( y x )2 |
( y |
2 |
x )2 |
( y |
N |
x |
N |
)2 |
1 1 |
|
2 |
|
|
|
– расстояние между точками с координатами x и y .
Рис. 108. Качественная картина устойчивости по Ляпунову
– две близкие в начальный момент времени траектории остаются близкими в любой последующий момент времени.
Соответственно, если условие (242) не выполняется, то траектория x(t) называется неустойчивой по Ляпунову.
Траектория x(t) может обладать и более сильным свойством – асимптотической устойчивостью. Траектория x(t) называется
асимптотически устойчивой, если она устойчива по Ляпунову и |
|||
если существует |
такое 0 , |
что для любой точки старта |
|
из -окрестности |
x0 , т.е. при |
y0 x0 |
, имеем: |
lim y(t) x(t) |
0 . |
(243) |
t |
|
|
Теорема Ляпунова. Ляпуновские показатели
Рассмотрим подробнее, как проводится анализ траекторий системы на устойчивость по Ляпунову. Пусть нам необходимо исследовать на устойчивость траекторию x(t) . Для этого выберем близкую к ней траекторию
y(t) x(t) x(t) , |
(244) |
соответствующую немного изменённому начальному условию. Здесь x(t) – некоторая малая добавка, определяющая различие между траекториями.
Т.к. y(t) является также траекторией системы, то она должна удовлетворять дифференциальному уравнению (241), т.е.
|
d y |
F(y) . |
|
|
dt |
|
|
Подставляя сюда выражение (244), получим |
|
||
d x |
d x |
F(x x) . |
(245) |
dt |
dt |
|
|
Учитывая малость отклонения x , разложим функцию F в ряд по x .