Правительство Российской Федерации
Федеральное государственное бюджетное образовательное учреждение
высшего профессионального образования
«Санкт-Петербургский государственный университет»
Кафедра прикладной кибернетики
Лежнина Юлия Геннадьевна
А ТТРАКТОРЫ
И КЛАССИЧЕСКИЕ ОПРЕДЕЛЕНИЯ
НЕУСТОЙЧИВОСТИ
Дипломная работа
Допущена к защите
Заведующий кафедрой:
д. ф.-м. н., профессор Г. А. Леонов
Научный руководитель:
д. ф.-м. н., профессор Г. А. Леонов
Рецензент:
д. ф.-м. н., профессор А. Л. Тулупьев
Санкт-Петербург
2016
SAINT PETERSBURG STATE UNIVERSITY
Department of Applied Cybernetics
Lezhnina Yulia
A TTRACTORS
AND CLASSICAL STABILITY THEORY
Graduation Thesis
Admitted for defence
Head of Department:
Doctor of Physics and Mathematics,
Professor Leonov G. A.
Scientific Supervisor:
Doctor of Physics and Mathematics,
Professor Leonov G. A.
Reviewer:
Doctor of Physics and Mathematics,
Professor Tulupyev A. L.
Saint Petersburg
2016
T3
Содержание
Введение .........................................................................................4
1. Предварительные сведения ...................................................6
1.1. Определения аттракторов .................................................................6
1.2. Основные определения в классической теории
устойчивости движения ....................................................................9
2. Взаимоотношения между основными понятиями
устойчивости ...........................................................................13
3. Чувствительность траекторий к начальным данным
и основные понятия неустойчивости .................................19
4. Сведение задачи к исследованию нулевого решения .....26
Заключение ...................................................................................27
Список литературы ......................................................................28
T4
Введение
Основоположникам геометрической (качественной)
теории дифференциальных уравнений и теории устойчивости
движения А. Пуанкаре [8, 9], Н. Е. Жуковскому [5],
А. М. Ляпунову [1, 6], Дж. Биркгофу[7] принадлежат методы
исследования и фундаментальные результаты, положившие
начало дальнейшему развитию в России и за рубежом.
В то время как теория устойчивости Ляпунова достигла
высокого уровня развития, теория прочности Жуковского
не получила должного развития и многие вопросы в ней
оказались открытыми.
Кроме вопроса о структуре аттракторов, интенсивному
изучению был подвергнут вопрос о зависимости аттракторов
от возмущения системы. Считается, что определяющим
свойством странного аттрактора является чувствительность
его траекторий по отношению к начальным данным. В этой
работе эти проблемы рассматриваются с точки зрения
классической теории устойчивости движения.
Многие известные исследователи предполагают, что
такая чувствительность по отношению к начальным данным
адекватна неустойчивости по Ляпунову [2, 3]. Однако это
справедливо только для дискретных динамических систем.
Для непрерывных систем необходимо вспомнить понятие
неустойчивости по Жуковскому. Один из основателей
современной аэродинамики, выдающийся русский учёный
Николай Егорович Жуковский ввёл своё понятие устойчивости движения в 1882 году [4, 5], на десять лет раньше
публикации исследований А. М. Ляпунова (1892 год [1]).
Понятие неустойчивости по Жуковскому адекватно
для чувствительности траекторий по отношению к начальным данным для непрерывных динамических систем.
T5
В этой работе были рассмотрены взаимосвязь основных
понятий теории устойчивости, а также примеры динамических систем, неустойчивых в одном смысле, но имеющие
устойчивость в другом.
Текст данной работы организован следующим образом.
В первом разделе вводятся такие понятия, как динамическая
система и аттрактор, даётся классификация аттракторов, а
также описаны необходимые определения классической
теории устойчивости, такие как устойчивость по Ляпунову,
по Красовскому, асимптотическая и экспоненциальная
устойчивость, устойчивость траекторий по Жуковскому и
по Пуанкаре. В следующих двух разделах представлены
обязательные теоремы, демонстрирующие связь основных
понятий, описанных в предыдущем разделе, и приведены
примеры, имеющие различные характеристики по основным
понятиям классической теории устойчивости. В последнем
разделе описан основной метод, используемый при исследовании устойчивости (или неустойчивости) траекторий
системы.
T6
1. Предварительные сведения
В этом разделе мы в основном следуем книгам и обзорам
[1, 11, 17-20].
1.1. Определения аттракторов
Аттрактор динамической системы — это притягивающее
замкнутое инвариантное множество в её фазовом
пространстве. Рассмотрим динамические системы, порождённые дифференциальными
dx =f(x), t∈R 1, x∈R n, (1.1)
(1.1)
x(t+1)=f x(t) , x∈Z, x∈R n (1.2)
(1.2)
и разностными
уравнениями. Здесь R n — евклидово пространство, Z — множество целых чисел, f(x) — вектор-функция: R n → R n.
Определение 1.1. Будем говорить, что уравнение (1.1)
или (1.2) порождает динамическую систему, если по любому
начальному состоянию x 0 ∈ R n однозначно определена траектория x(t, x 0 ) при t ∈ [0, +∞).
Здесь x(0, x 0) = x 0.
Хорошо известно, что в рассматриваемом случае для решений системы (1.1) справедливо следующее полугрупповое
свойство
x(t + s, x 0) = x t, x(s, x 0) (1.3)
(1.3)
при всех t ≥ 0, s ≥ 0.
Динамическую систему, порожденную уравнением (1.1),
будем называть непрерывной, а уравнением (1.2) — дискретной динамической системой.
T7
Введем определения аттракторов, следуя в основном
работе [10].
Определение 1.2. Будем говорить, что K инвариантно,
если x(t,K) = K, ∀t ≥ 0.
Здесь
x(t,K) = {x(t,x 0)|x 0 ∈ K}.
Определение 1.3. Будем говорить, что инвариантное
множество K является локально притягивающим, если
для некоторой ε-окрестности этого множества K(ε) выполнено
соотношение
lim ρ K, x(t, x 0) = 0, ∀ x 0 ∈ K(ε).
t→+∞
Здесь ρ(K, x) — расстояние от точки x до множества K,
которое определяется по формуле
ρ(K, x) = inf |z − x|.
z∈K
| · | — евклидова норма в R n, K(ε) — множество точек x,
|·|
для которых ρ(K, x) < ε.
Определение 1.4. Будем говорить, что инвариантное
множество K является глобально притягивающим, если
lim ρ K,x(t,x 0 ) = 0, ∀x 0 ∈ R n .
t→+∞
Определение 1.5. Будем говорить, что инвариантное
множество K является равномерно локально притягивающим,
если для некоторой его ε-окрестности K(ε) и для любых числа
δ > 0 и ограниченного множества B существует число t(δ,B) > 0
такое, что
x t,B∩K(ε) ⊂K(δ), ∀t≥t(δ,B)
Здесь
x t,B∩K(ε) ={x(t,x 0)|x 0 ∈(B∩K(ε))}.
T8
Определение 1.6. Будем говорить, что инвариантное
множество K является равномерно глобально притягивающим,
если для любых числа δ > 0 и ограниченного множества B ∈ R n
существует число t(δ, B) > 0 такое, что
x(t, B) ⊂ K(δ), ∀t ≥ t(δ, B).
Определение 1.7. Будем говорить, что инвариантное
множество K устойчиво по Ляпунову, если для любого числа ε > 0
существует число δ > 0 такое, что
x(t, K(δ)) ⊂ K(ε), ∀t ≥ 0.
Заметим, что, если K состоит из единственной траектории, то последнее определение совпадает с классическим
определением устойчивости решения по Ляпунову. Если
такое K является также локально притягивающим, то имеет
место асимптотическая устойчивость в смысле Ляпунова.
Определение 1.8. Будем говорить, что K — аттрактор,
если K является инвариантным, замкнутым, локально притягивающим множеством.
Определение 1.9. Будем говорить, что K — глобальный
аттрактор, если K является инвариантным, замкнутым,
глобально притягивающим множеством.
Определение 1.10. Будем говорить, что K — B-аттрактор,
если K является инвариантным, замкнутым, равномерно
локально притягивающим множеством.
Определение 1.11. Будем говорить, что K — глобальный
B -аттрактор, если K является инвариантным, замкнутым,
равномерно глобально притягивающим множеством.
T9
1.2. Основные определения в классической теории
устойчивости движения
Рассмотрим систему
(1.4)
T
где
– непрерывная вектор-функция и систему
(1.5)
T
Обозначим через
решение системы (1.4) или
системы (1.5) с начальными данными
T
Определение 1.12. Решение
называется устойчивым по Ляпунову, если для любых
и
существует
число δ(ε, t 0 ) такое, что
1) все решения x(t, t 0 , y 0 ) , удовлетворяющие условию
|x 0 − y 0 | ≤ δ,
определены при t ≥ t 0 ,
2) для этих решений справедливо неравенство
J
Если δ(ε, t 0 ) не зависит от t 0 , то устойчивость по Ляпунову
называется равномерной.
Определение 1.13. Решение x(t,t 0 ,x 0 ) называется асимптотически устойчивым по Ляпунову, если оно устойчиво по Ляпунову и для любого t 0 ≥ 0 существует число ∆(t 0 ) > 0 такое, что
решение x(t, t 0 , y 0 ), удовлетворяющее условию |x 0 − y 0 | ≤ ∆
обладает свойством
T
T10
Определение 1.14. Решение x(t, t 0 , x 0 ) называется устойчивым по Красовскому, если существуют такие положительные
числа δ(t 0 ) и R(t 0 ), что для любого y 0 , удовлетворяющего условию
|x 0 − y 0 | ≤ δ(t 0 ) решение x(t, t 0 , y 0 ) определено при всех
t ≥ t 0 и выполнено неравенство
J
Если δ и R не зависят от t 0 , то устойчивость по Красовскому
называется равномерной.
Определение 1.15. Решение x(t,t 0 ,x 0 )
называется экспоненциально устойчивым, если существуют такие положительные
числа δ(t 0 ) , R(t 0 )и α(t0), что для любого y 0 , удовлетворяющего
условию
|x 0 − y 0 | ≤ δ(t 0 ) решение x(t, t 0 , y 0 ) определено при
всех t ≥ t 0 и выполнено неравенство
J
Если δ, R и α не зависят от t 0 , то экспоненциальную
устойчивость называют равномерной.
Перейдём теперь к рассмотрению динамических систем
(1.1) и (1.2).
Введём следующее обозначение
T
Определение 1.16. Траектория x(t,x 0) динамической системы
называется устойчивой по Пуанкаре (или орбитально устойчивой), если для любого числа ε > 0 существует число δ(ε) > 0 такое,
что при всех y 0 , удовлетворяющих неравенству |x 0 − y 0 | ≤ δ(ε),
выполнено соотношение
ρ(L + (x 0 ), x(t, y 0 )) ≤ ε, ∀t ≥ 0
T11
Если, кроме того, для некоторого числа δ 0 и всех y 0 , удовлетворяющих неравенству |x 0 − y 0 | ≤ δ 0 , выполнено соотношение
lim ρ(L + (x 0 ),x(t,y 0 ))=0, t→+∞
то говорят, что траектория x(t,x 0 ) асимптотически устойчива
по Пуанкаре (или асимптотически орбитально устойчива).
Заметим, что приведённом выше определении для непрерывных динамических систем t ∈ R, а для дискретных динамических
систем t ∈ Z.
Введем теперь определение устойчивости по Жуковскому
для непрерывных динамических систем. Для этого потребуется рассмотрение следующего множества гомеоморфизмов:
Hom{τ(·)|τ : [0,∞) → [0,+∞),τ(0) = 0}.
Функции τ(t) из множества Hom будут играть роль перепараметризации времени для траектории системы (1.1).
Определение 1.17. Траектория x(t,x 0 ) системы (1.1) называется устойчивой по Жуковскому, если для любого числа ε > 0
существует число δ(ε) > 0 такое, что для любого вектора y 0 ,
удовлетворяющего неравенству |x 0 − y 0 | ≤ δ(ε),найдется функция
τ(·) ∈ Hom, при которой выполнено неравенство
|x(t,x 0 )−x(τ(t),y 0 )|≤ε, t≥0.
Если, кроме того, для некоторого числа δ 0 > 0 и любого y 0 из
шара{y||x 0 −y 0 | ≤ δ 0 }найдется функция τ(·) ∈ Hom,при которой
выполнено соотношение
lim |x(t,x 0 )−x(τ(t),y 0 )|=0,
t→+∞
T12
то будем говорить, что траектория x(t,x 0 ) асимптотически
устойчива по Жуковскому.
Другими словами, устойчивость по Жуковскому — это
устойчивость по Ляпунову при подходящей перепараметризации каждой из возмущенных траекторий.
По определению неустойчивостью по Ляпунову (по Красовскому, по Пуанкаре, по Жуковскому) считается отрицание
соответствующего вида устойчивости.
T13
2. Взаимоотношения между основными
понятиями устойчивости
Сформулируем следующие очевидные утверждения.
Предложение 2.1. Для непрерывных динамических систем
из устойчивости по Ляпунову следует устойчивость по Жуковскому, а из устойчивости по Жуковскому следует устойчивость
по Пуанкаре.
Предложение 2.2. Для дискретных динамических систем
из устойчивости по Ляпунову следует устойчивость по Пуанкаре.
Предложение 2.3. Для состояний равновесия все введенные
выше определения Ляпунова, Жуковского и Пуанкаре эквивалентны.
Предложение 2.4. Для периодических траекторий дискретных систем эквивалентны определения устойчивости по Ляпунову
и по Пуанкаре.
Доказательство. Пусть периодическая траектория дискретной системы устойчива по Пуанкаре. Выберем ε-окрестность U(ε, L+(x0)) периодической траектории L+(x0) так,
чтобы
U(ε, u) ∩ U(ε, z) = ∅.
для любых точек u∈L +(x 0), z∈L +(x 0), u =
̸ z.
множество.
(2.1)
Здесь – пустое
Рассмотрим y 0 из δ-окрестности точки x 0 и соединим x 0 и
y 0 отрезком x 0+γ(y 0−x 0), γ∈[0,1].
T14
Из непрерывности f (x), устойчивости по Пуанкаре и
соотношения (2.1) следует, что
x(t,x 0+γ(y 0−x 0))∈U ε(x(t,x 0)), ∀t≥0, ∀γ∈[0,1].
(2.2)
Включение (2.2) при γ = 1 означает устойчивость по Ляпунову траектории x(t, x 0). Отсюда следует эквивалентность
определений устойчивости по Пуанкаре и по Ляпунову
периодических траекторий дискретных систем.
Аналогично доказывается эквивалентность асимптотической устойчивости по Пуанкаре и по Ляпунову периодических траекторий дискретных систем.
Предложение 2.5. Д л я п е р и о д и ч е с к и х т р а е к т о р и й
непрерывных систем эквивалентны определения устойчивости
по Пуанкаре и по Жуковскому.
Доказательство. Пусть периодическая траектория
непрерывной системы устойчива по Пуанкаре. Выберем
ε-окрестность U(ε,L +(x 0)) периодической траектории L +(x 0)
так, чтобы
+
+
U(ε,L (x 0 ),t 1 )∩U(ε,L (x 0 ),t 1 )=∅,
(2.3)
∀t 1 =
̸ t 2 , t 1 ∈[0,T), t 2 ∈[0,T)
где
U(ε,L +(x 0),t)={y||y−x(t,x 0)|≤ε, (y−x(t,x 0)) ∗f(x(t,x 0))=0}.
По ε выберем δ так, чтобы из неравенства |x 0 −y 0| ≤ δ
следовало соотношение
x(t, y 0) ∈ U(ε, L +(x 0)), ∀t ≥ 0.
Не умаляя общности можно считать, что y 0 ∈ U(ε,L +(x 0),0).
Тогда из соотношения (2.3) следует, что можно определить
T15
перепараметризацию τ(t) траектории x(t, y 0 ) следующим
образом
x(τ(t),y 0) ∈ U(ε,L +(x 0),t).
Ясно, что τ(·) ∈ Hom и выполнено неравенство
|x(τ(t),y 0)−x(t,x 0)|≤ε, t≥t 0.
Последнее означает устойчивость по Жуковскому.
Аналогично доказывается эквивалентность асимптотической устойчивости по Жуковскому и по Пуанкаре. В этом
случае хорошо известна [11] специальная перепараметризация τ (t) = t + c называемая асимптотической фазой.
Хорошо известны примеры периодических траекторий
непрерывных систем, которые не устойчивы по Ляпунову, но
устойчивы по Пуанкаре.
Приведём один из таких примеров.
Пример 2.1. Рассмотрим уравнение Дуффинга
T
или эквивалентную ему двумерную систему
T
(2.4)
Легко видеть, что все решения системы (2.4) периодичны
и удовлетворяют соотношению
T
(2.5)
Здесь y_0=y(0), \ x_0=x(0).В дальнейшем, не умаляя
общности, будем считать, что y_0=y(0). Последнее возможно,
поскольку любая периодическая траектория системы (2.4)
T16
проходит в некоторый момент времени через точку
с координатой y=0 (рис. 2.1).
T
Рис. 2.1
На рисунке 2.1 траектории системы (2.4) — это графики
функции y(x)=±√{c-x^2-(x^4)/2}, где параметр с принимает
положительные значения.
Обозначим через T(x0) период решения, проходящего
через точку y=0, x=x0.
Из соотношений (2.4) и (2.5) получим, что в полуплоскости {y≥0} выполнено равенство
dt/dx=(x0^2+(x0^4)/2-x^2-(x^4)/2)^(-1/2)
Отсюда и из равенства x(T(x0)/2)=-x0
получим формулу для периода T(x0):
(см. рис. 2.1)
T(x_0)=2\int^{x_0}_{-x_0} (x_0^2+1/2x04-x2-1/2x4)^{-1/2}dx
(2.6)
T17
В отличие от линейной системы
x=y,y=-x,(2.7)
где период решения с начальными данными
не зависит от x0 и равен 2pi:
(2.7)
y=0, x=x0
период T(x0) решения системы (2.4), вычесленный по формуле (2.6) зависит от x0.
График T(x_0) изображен на рисунке 2.2.
T
Рис. 2.2
Из соотношения (2.5) следует, что при малых изменениях
начальных данных траектории системы (2.4) целиком
остаются в малых окрестностях друг друга и, следовательно,
являются устойчивыми по Пуанкаре. Поскольку все рассматриваемые траектории периодичны, то они устойчивы
по Жуковскому (Предл. 2.5).
T18
Однако здесь нет устойчивости по Ляпунова в отличие
от траекторий системы (2.7). Это следует из следующих
соображений. Если рассмотреть решения системы (2.4)
с начальными данными x=x_0, y=0 и x=x_1, y=0, где
x_0-x_1 =d и d мало, то можно подобрать натуральное число
N так, чтобы |N(T(x_0)-T(x_1))| было близко к T(x_1)/2.
Отсюда следует, что x(NT(x_0),x_1) близко к значению -x1.
Из приведённых здесь рассуждений следует, что решение
x(t,x_0) не является устойчивым по Ляпунову: для малых
x_0-x_1 разность |x(NT(x_0),x_0)-x(NT(x_0),x1) близка к 2|x0.
T19
3. Чувствительность траекторий
к начальным данным и основные
понятия неустойчивости
Перейдем теперь к сравнению приведенных выше определений и эффекта чувствительности траекторий по отношению
к начальным данным на странных аттракторах.
Как было только что показано в примере 2.1 для периодических траекторий непрерывных динамических систем
возможна устойчивость по Пуанкаре и по Жуковскому и
в тоже время неустойчивость по Ляпунову. Таким образом,
неустойчивость по Ляпунову не может характеризовать
"отталкивание" друг от друга непрерывных траекторий
при малых вариациях начальных данных.
Оказывается, что и неустойчивость по Пуанкаре не может
характеризовать "отталкивание" траекторий друг от друга, но
уже по другой причине. Здесь возмущённое решение может
выходить из ε-окрестности некоторого куска невозмущённой
траектории (эффект отталкивания), но одновременно это
решение может входить в ε-окрестность другого куска невозмущённой траектории (свойство устойчивости по Пуанкаре).
Таким образом, отталкивающиеся друг от друга траектории
могут быть устойчивыми по Пуанкаре.
Рассмотрим более подробно эти эффекты.
В компьютерных экспериментах часто оказывается, что
траектории, находящиеся на неустойчивом многообразии
седловой особой точки всюду плотно заполняют B-аттрактор
(или его часть, состоящую из ограниченных траекторий). Это
наблюдается на B-аттракторе системы Лоренца [12]
x ̇= − σ ( x − y )
(3.1)
y ̇ = rx − y − xz
z ̇ = −bz + xy,
где σ = 10, r = 28, b = 8/3 (рисунки 3.1, 3.2).
T20
T
Рис. 3.1. Неустойчивое многообразие седла системы Лоренца.
Первые пятьдесят оборотов.
T
Рис. 3.2. Неустойчивое многообразие седла системы Лоренца.
Следующие пятьдесят оборотов.
Аналогичная ситуация имеет место для системы Хенона [13]
x(t + 1) = a + by(t) − x(t) 2
y(t + 1) = x(t),
где a = 1.4, b = 0.3 (рис. 3.3).
(3.2)
T21
T
Рис. 3.3. Неустойчивое многообразие седла системы Хенона.
Затемненная часть квадрата — начальны данные, при
которых траектории стремятся к бесконечности.
Таким образом для этих систем
ρ(L +(x 0), x(t, y 0)) ≡ 0,
(3.3)
∀ x 0 ∈ M, y 0 ∈ K, ∀ t>0. Здесь M — неустойчивое многообразие седла, K– упомянутый выше B-аттрактор (или его
часть, состоящая из ограниченных траекторий).
Из соотношения (2.6) следует, что на аттракторе K
траектории x(t, x 0) асимптотически устойчивы по Пуанкаре.
Последнее означает, что неустойчивость по Пуанкаре
не может характеризовать чувствительность по отношению
к начальным данным на странных аттракторах.
Приведём теперь аналогичные примеры, носящие
аналитический характер.
Пример 3.1. Рассмотрим линеаризованные уравнения
двух несвязных маятников
x1’=y1, y1’=-w1^2x1,
x2’=y2, y2’=-w2^2x2.
(3.4)
T22
Решениями этих уравнения являются функции
x1(t)=Asin(w1t+𝜑1(0)),
y1(t)=Aw1cos(w1t+𝜑1(0)),
x2(t)=Bsin(w2t+𝜑2(0)),
y2(t)=Bw2cos(w2t+𝜑2(0)).
При фиксированных A и B траектории системы (3.4)
расположены на двумерных торах
w1^2x1^2+y1^2=A^2, w1^2x2^2+y2^2=B^2.
(3.5)
В том случае, когда число w1/w2 является иррациональным, при любых начальных условиях 𝜑1(0) и 𝜑2(0) траектории системы (3.4) всюду плотно расположены на торах (3.5).
Таким образом, при иррациональном w1/w2 для любых
точек x0∈R4, y0∈R4, расположенных на торе (3.5), справедливо равенство (3.3). Отсюда следует асимптотическая
устойчивость по Пуанкаре траекторий динамической системы
с фазовым пространством (3.5). Однако движение точек x(t,x0)
и x(t,y0) вдоль траекторий происходит так, что они
не стремятся друг к другу при t->\infty. Нет "прижимания"
друг к другу и траекторий. Следовательно, наше интуитивное
представление об асимптотической устойчивости, как о процессе сходимости объектов друг к другу, противоречит
формальному определению Пуанкаре. Легко видеть, что
аналогичного эффекта нет для понятия устойчивости
по Жуковскому: в рассматриваемом случае асимптотическая
устойчивость по Жуковскому не имеет места.
Пример 3.2. Вновь рассмотрим динамическую систему
(3.4) с иррациональным w1/w2. Изменим здесь поток траекторий на торе следующим образом. Произведём разрез
T23
поверхности тора вдоль некоторого куска фиксированной
траектории от точки z1 до точки z2 (рис. 3.4).
T
Рис. 3.4
Далее диффеоморфно растянем рассматриваемую
поверхность вдоль тора так, чтобы разрез преобразовался
в окружность с фиксированными точками z1 и z2 (рис. 3.5).
Рис. 3.5
Внутренность окружности обозначим через H.
Изменим нашу динамическую систему так, чтобы z1 и z2
стали седловыми стационарными точками, а полуокружности,
соединяющие z1 и z2, — гетероклиническими траекториями,
стремящимися при t->\infty и t->-\infty соответственно к z2 и
z1 (рис. 3.5).
T24
Вне "дырки" H оставляем прежним расположение траекторий на торе, которое получилось после диффеоморфного
растягивания.
Рассмотрим теперь поведение траекторий полученной
здесь системы с точки зрения Пуанкаре и Жуковского.
Вне дырки H траектории всюду плотны на торе и
поэтому, как было показано ранее, они асимптотически
устойчивы по Пуанкаре.
Рассмотрим теперь некоторую 𝛿-окрестность точки z0,
расположенной на торе и вне множества H. Траектория,
выпущенная из точки z0, либо всюду плотна, либо совпадает
с сепаратрисой S седла z1, стремящейся к z1 при t->\infty
(рис. 3.5). Поэтому существует момент времени t, когда часть
траекторий, выпущенных из 𝛿-окрестности z0 будут
находиться в малой окрестности точки z1 и справа от сепаратрисы S. Другая же часть траекторий, выпущенных из этой
окрестности z0 будут в момент времени t находиться слева
от сепаратрисы S. Ясно, что после этого траектории,
расположенные справа и слева от сепаратрисы S, будут
огибать "дырку" H соответственно справа и слева. Но тогда
ясно, что при таком огибали траектории будут отталкиваться
друг от друга и, следовательно, траектория, выпущенная
из точки z0 будет неустойчива по Жуковскому.
Таким образом, траектория может быть асимптотически
устойчива по Пуанкаре и неустойчива по Жуковскому.
В этом примере траектории чувствительны по отношению
к начальным данным и могут расходиться друг от друга
в некоторые моменты времени на значительные расстояния.
Понятие неустойчивости по Жуковскому адекватно такой
чувствительности.
T25
Заметим, что множество таких чувствительных траекторий расположено на гладком многообразии "тор минус
дырка H". Таким образом, ограниченное, инвариантное
множество траекторий, чувствительных по отношению
к начальным данным не всегда имеет нецелую хаусдорфову
размерность или структуру канторовского множества.
Из сказанного выше можно сделать вывод о том, что
среди классических понятий неустойчивости для дальнейшего
изучения странных аттракторов наиболее адекватны
неустойчивость по Жуковскому (в непрерывном случае) и
неустойчивость по Ляпунову (в дискретном случае).
Элементы теории устойчивости по первому приближению
в смысле Жуковского развивались в [14, 15].
Недавно был получен следующий замечательный
результат [16]. ω-предельное множество асимптотически
устойчивой по Жуковскому траектории состоит из периодической траектории.
Таким образом, предельные множества асимптотически
устойчивых по Жуковскому траекторий устроены достаточно
просто. C другой стороны, как было показано выше,
неустойчивость по Жуковскому может быть одной из характеристик странных аттракторов.
T26
4. Сведение задачи к исследованию
нулевого решения
Первая процедура, которая применяется чаще всего [17]
при исследовании устойчивости (или неустойчивости)
решения x(t,t0,x0) системы (1.4) является следующее
преобразование x = y + x(t,t0,x0).После такой замены получим следующее уравнение
уравнение 4.1,
(4.1)
которое часто называют дифференциальным уравнением
возмущенного движения. Очевидно, что задача об
исследовании устойчивости решения x(t,t0,x0) преобразуется
в задачу об исследовании устойчивости тривиального
решения y(t)==0 уравнения (4.1).
При этом считают, что правая часть (4.1) известна,
поскольку известны F(x,t) и решение x(t,t0,x0). В настоящее
время трудности вычисления x(t,t0,x0) часто преодолеваются
с помощью численных методов и проведения компьютерных
экспериментов.
При исследовании устойчивости по Красовскому и
экспоненциальной устойчивости поступают аналогичным
образом. Устойчивость по Жуковскому есть ничто иное как
устойчивость по Ляпунову перепараметризованных траекторий.
Поэтому для исследования устойчивости по Жуковскому есть
возможность применить методы и приёмы, разработанные для
изучения ляпуновской устойчивости.
T27
Заключение
В этой работе были приведены основные определения
устойчивости, и рассмотрена взаимосвязь основных понятий
теории устойчивости и неустойчивости, а также показана
чувствительность траекторий к начальным данным. Для неё
наиболее адекватны неустойчивость по Жуковскому
(в непрерывном случае) и неустойчивость по Ляпунову
(для дискретных динамических систем).
Список литературы
1. Ляпунов А.М. Общая задача об устойчивости движения.
Харьков, 1892, 250 с.
2. Анищенко В.С. Сложные колебания в простых системах.
М., Наука, 1990, 311 c.
3. Fradkov A.L., Pogromsky A.Y. Introduction to Control
of Oscillation and Chaos. World Scientific. Singapore, 1998,
390 p.
4. Жуковский Н.Е. О прочности движения. Собрание
сочинений. Т.1. М., Гостехиздат, 1948, c. 67–160.
5. Жуковский Н.Е. О прочности движения // Ученые
записки Московского университета. Отд. физ.-мат. 1882,
вып. 4, с. 1 104.
6. Ляпунов A.M. Об устойчивости движения в одном
частном случае задачи о трех телах // Сообщения
Харьковского Матем. Об-ва, 2 серия, 1889.
7. Биркгоф Дж. Динамические системы. М.: Гостехиздат,
1941.
8. Пуанкаре А. О кривых, определяемых
дифференциальными уравнениями. M.-JL: ГТТИ, 1947.
9. Пуанкаре А. Избранные труды. Т. 1,2.- М.: Наука, 1971,
1972.
10. Ладыженская О.А. О нахождении минимальных
глобальных аттракторов для уравнений Навье-Стокса и
других уравнений с частными производными. УМН, 1987,
т. 42, вып. 6, с. 25–60.
T29
11. Демидович Б.П. Лекции по математической теории
устойчивости. М., Наука, 1967, 472 с.
12. Странные аттракторы. Сборник статей. Под редакцией
Я.Г. Синая и Л.П. Шильникова. Математика. Новое
в зарубежной науке N 22. М., Мир, 1981, 253 с.
13. Леонов Г.А., Полтинникова М.С. О ляпуновской
размерности аттрактора диссипативного отображения
Чирикова. Труды Санкт-Петербургского математического
общества, 2002, т. 10, с. 186–198.
14. Леонов Г.А. Об устойчивости по первому приближению.
ПММ, 1998, т.69, вып.4, с. 548– 555.
15. Дунаева О.В. Признаки асимптотической прочности и
непрочности движения динамической системы. Докл. АН,
1997, т. 315, N4, с. 476–478.
16. Yang X. Lyapunov Asymptotically Stability and Zhukovskij
Asymptotically Stability. Chaos, Solutions and Fractals, 2000,
N 11, p. 1995–1999.
17. Малкин И.Г. Теория устойчивости движения. М., Наука,
1966, 525 с.
18. Четаев Н.Г. Устойчивость движения. М., Гостехиздат,
1955, 207 c.
19. Далецкий Ю.Л., Крейн М.Г. Устойчивость решений
дифференциальных уравнений в банаховом пространстве.
М., Наука, 1970, 534 с.
20. Леонов Г.А. Странные аттракторы и классическая теория
устойчивости движения. Успехи механики. 2002, т.1, N 3,
с. 3–42.
Отзывы:
Авторизуйтесь, чтобы оставить отзыв