§2 Динамические системы

Простейшие свойства автономных (или, как их часто называют, динамических) систем дифференциальных уравнений. В этом очерке мы продолжим изучение специфических свойств автономных систем уравнений:

x′ = f(x), (1)

отличающих такие системы от неавтономных.

Предполагается, что → и удовлетворяет условию Липшица; поэтому решение задачи Коши для уравнения (1) с любыми начальными данными определено на всей оси и единственно. Через мы будем как всегда обозначать оператор сдвига по траекториям уравнения (1), а через — поток, порожденный этим уравнением, т. е. семейство отображений:

{} = {}.

Известно, что:

при всех s ∈ R. В терминах решений уравнения (1) этот факт эквивалентен тому, что сдвиг вдоль оси t любого решения уравнения (1) является его решением.

Напомним, что интегральной кривой уравнения (1) мы называем график Γφ любого решения φ в расширенном фазовом пространстве R×, а траекторией Tφ — проекцию интегральной кривой на фазовое пространство (см. рис. 1). Таким образом, если φ — определенное на D(φ) решение уравнения (1), то по определению:

Γφ = {(t, x) ∈ R× t ∈ D(φ), x = φ(t)},

Tφ = {x ∈ : x = φ(t), t ∈ D(φ)},

Рис. 1.

Задача 1. Докажите, что если две траектории отвечают непродолжимым решениям уравнения (1) и имеют общую точку, то они совпадают.

В дальнейшем мы под траекториями будем понимать только траектории отвечающие максимальным (непродолжимым) решениям системы (1).

Траекторию, соответствующую стационарному (т. е. постоянному) решению, будем называть стационарной точкой (положением равновесия), а траекторию, соответствующую периодическому решению — циклом. Если найдется такое число T > 0, что φ() ≠ φ() при 0 < | – | < T и φ является T-периодическим решением, то о траектории Tφ будем говорить как о цикле без самопересечений (см. рис. 2). Ясно, что если T-периодическое решение φ порождает цикл без самопересечений, то решение φ, рассматриваемое как 2T-периодическое, порождает цикл с самопересечениями; поэтому понятие цикла без самопересечений является характеристикой не только самой траектории, но и порождающего ее решения. Наконец, мы будем говорить, что траектория Tφ не имеет самопересечений, если φ() ≠ φ() при любых неравных и .

Рис. 2.

Имеет место следующая Теорема о типах траекторий автономных систем. Каждая траектория системы (1) является либо стационарной точкой, либо циклом без самопересечений, либо траекторией без самопересечений.

Д о к а з а т е л ь с т в о. Пусть Tφ — отвечающая максимальному решению φ траектория системы (1). Если Tφ имеет самопересечение, т. е. существуют такие τ и θ, что θ > τ и φ(τ) = φ(θ), то решение φ является (θ – τ)-периодическим. Действительно,

φ(t) = [φ(τ)]= [φ(τ)] = [φ(θ)] = φ[t + (θ – τ)]

при всех t ∈ R. Далее в доказательстве нам потребуются следующие утверждения.

Задача 2. Пусть φ — непрерывная функция на R. Тогда множество всех ее периодов:

а) образует аддитивную группу,

б) является замкнутым.

Обозначим через T точную нижнюю грань множества положительных периодов решения φ. Альтернатива такова: либо T > 0, либо T = 0. В первом случае, как нетрудно видеть, траектория образует цикл без самопересечений (периода T). Покажем, что во втором случае Tφ есть стационарная точка. Поскольку T = 0, множество P периодов решения φ содержит сколь угодно малые числа. Поэтому, так как P — аддитивная группа, P всюду плотно на числовой прямой, а так как P замкнуто, P = R. Таким образом, любое вещественное число является периодом решения φ, что очевидно эквивалентно его стационарности.

Задача 4. Докажите, что в одномерном случае траекториями автономной системы могут быть только множества вида: точка, открытый промежуток, открытый луч, вся прямая.

Теперь мы переходим к исследованию поведения решений при t → ∞. Пусть φ — решение уравнения (1), определенное на всей оси. Частичные пределы φ при t → +∞ называются ω-предельными точками траектории Tφ (и решения φ). Множество всех ω-предельных точек траектории Tφ называется ω-предельным множеством этой траектории и обозначается через Ω(Tφ). Другими словами, точка x лежит в Ω(Tφ) в том и только том случае, когда существует такая последовательность → +∞, что φ() = x. Аналогичная конструкция при t → –∞ приводит к понятию α-предельной точки и α-предельного множества A() траектории Tφ (см. рис. 3).

Рис. 3.

Задача 5. Покажите, что:

где черта означает замыкание множества в .

Структуру ω-предельных множеств частично выясняет следующая

Теорема о структуре ω-предельных множеств. Для любого решения φ системы (1) множество Ω(Tφ) замкнуто. Если, кроме того, решение φ ограничено вправо, т. е. множество φ([0, +∞)) ограничено, то Ω(Tφ) не пусто, компактно и связно.

Д о к а з а т е л ь с т в о. Замкнутость Ω(Tφ) очевидна в силу задачи 5. Не пустота вытекает из относительной компактности ограниченных подмножеств конечномерного пространства, а компактность Ω(Tφ) — из только что упомянутого факта и замкнутости ω-предельного множества.

Для доказательства связности Ω(Tφ) предположим противное. Тогда, по определению, Ω(Tφ) представимо в виде объединения двух замкнутых непересекающихся подмножеств и . Расстояние:

ρ = inf{ : x ∈ , y ∈ }

в силу замкнутости и ограниченности и положительно. Пусть ∈ (i = 1, 2) — произвольные точки.

По определению Ω(Tφ) найдутся такие , → +∞, что:

φ()→ ∈ и φ()→ ∈ , при k → ∞.

Не теряя общности, можно считать, что при всех k ∈

и

inf{ : x ∈ } < ρ/2,

inf{ : x ∈ } < ρ/2.

В силу непрерывности φ и замкнутости и ограниченности и найдутся такие ∈ (,), что

inf{ : x ∈ } < ρ/2,

inf{ : x ∈ } < ρ/2

(2)

(см. рис. 4). Поскольку последовательность φ() ограничена (и следовательно, относительно компактна), некоторая ее подпоследовательность φ() сходится, скажем, к точке z (ω-предельной точке траектории Tφ). Заменяя в (2) k на и переходя к переделу при m → ∞, получим противоречие: z ∉ ∪ = Ω(Tφ).

Аналогичная теорема имеет место и для α-предельных множеств.

Рис. 4.

Множество M ⊆ называется положительно (соответственно, отрицательно) инвариантным множеством уравнения (1) (и порожденного им потока ), если (M) ⊆ M при всех t > 0 (соответственно, t < 0). Одновременно положительно и отрицательно инвариантное множество называется инвариантным.

Задача 6. Докажите, что объединение любого семейства траекторий является инвариантным множеством. Наоборот, любое инвариантное множество есть объединение некоторого семейства траекторий.

Теорема об инвариантности предельных множеств. ω- и α- предельные множества любой траектории инвариантны.

Д о к а з а т е л ь с т в о. Пусть x ∈ Ω(Tφ), т. е. существует такая последовательность → ∞, что φ()= x.

Тогда при любом t ∈ R в силу непрерывности оператора сдвига:

= Ω(Tφ),

что и требовалось доказать.

Таким образом, предельные множества любой траектории состоят из траекторий системы (1).

Продолжая изучение поведения решений при t → ∞, перейдем к исследованию устойчивости решений и траекторий динамических систем. Начнем с одного наводящего примера. Рассмотрим автономную систему, фазовый портрет которой изображен на рис. 5, причем, пусть время обращения фазовой точки по круговым траекториям зависит от радиуса этой траектории. Похожий фазовый портрет имеет, например,

Уравнение математического маятника:

y′′ + ω2sin y = 0

в окрестности нижнего положения равновесия. Выберем произвольную траекторию (отличную от стационарной точки) и порождающее ее решение φ. Заметим теперь, что это решение не является устойчивым по Ляпунову. Действительно, как бы близко к φ(0) мы не взяли начальные данные решения ψ (отличного от φ), при некоторых t > 0 в силу различия периодов решений φ и ψ точки φ(t) и ψ(t) будут почти диаметрально противоположными и, как следствие, расстояние между φ(t) и ψ(t) не может быть сделано сколь угодно малым. Положение же равновесия этой динамической системы очевидно устойчиво по Ляпунову.

Рис. 5.

С асимптотической устойчивостью дело обстоит еще хуже, а именно, всякое периодическое нестационарное решение φ автономной системы никогда(!) не бывает асимптотически устойчивым. В самом деле, начальные данные решений φ(t) и ψ(t) = φ(t + τ) при малых τ сколь угодно близки, а функция не может стремиться к нулю при t → ∞, т. к. α является периодической ненулевой функцией.

В то же время, представляется физически естественным считать колебания математического маятника вблизи нижнего состояния равновесия устойчивыми: малые возмущения начальных данных приводят к малым изменениям амплитуды, фазы и частоты колебаний. В этом направлении мы и модифицируем понятия устойчивости и асимптотической устойчивости.

Полутраекторией , порожденной решением φ, будем называть часть траектории Tφ, отвечающую положительным значениям t. Решение φ системы (1) и соответствующую ему траекторию Tφ называют орбитально устойчивыми, если для каждого ε > 0 найдется такое положительное δ, что из неравенства ≤ δ вытекает неравенство:

Inf { : y ∈ }≤ ε (3)

при всех t ≥ 0. Соотношение (3) означает, что при малых возмущениях начальных данных траектория решения остается в ε-окрестности полутраектории (см. рис. 6).

Таким образом, требование:

в определении устойчивости по Ляпунову заменяется более слабым требованием (3).

Рис. 6.

Орбитально устойчивое решение φ и соответствующая траектория Tφ называются орбитально асимптотически устойчивыми, если найдется такое Δ > 0, что Inf { : y ∈ } при k → ∞ как только Inf { : y ∈ } Δ. Другими словами, расстояние от до полутраектории стремится к нулю при t → ∞, если решение начинается достаточно близко к (см. рис. 7).

Рис. 7.

Предельным циклом называют цикл Tφ, для которого либо Tφ = Ω(Tφ), либо Tφ = A(Tφ).

Задача 7. Является ли цикл {(, ) ∈ : + = 1} системы:

= + ( ),

= + ( )

предельным? Орбитально асимптотически устойчивым?

При исследовании устойчивости предельных циклов весьма полезным оказывается рассмотрение так называемой функции последования. Строится она следующим образом. Пусть Tφ — предельный цикл системы (1). Тогда f(x) ≠ 0 при всех x ∈ Tφ (почему?).

Зафиксируем точку ∈ Tφ и обозначим через E (n – 1)-мерную гиперплоскость, проходящую через точку x0 и ортогональную вектору f() (см. рис. 8). Будем считать, что E наделена аффинной структурой, индуцированной из , и кроме того, будем считать, что начало координат в E фиксировано в точке . Наконец, будем считать E нормированным пространством с нормой опять же индуцированной из . Пусть U — малая окрестность нуля в E. Для точки x ∈ U обозначим через G(x) первую точку пересечения кривой t → (x) (t > 0) с гиперплоскостью E, т. е. такую точку y ∈ E, что y = (x) при некотором > 0, причем таково, что (x) ∉ E при t ∈ (0, ) (см. рис. 8). Через D(G) обозначим множество тех x ∈ U, для которых G(x) определена.

Отображение G: D(G) → E и называется функцией последования или оператором последования Пуанкаре.

Рис. 8.

Теорема об орбитальной асимптотической устойчивости. Пусть Tφ — предельный цикл уравнения (1) и G — какая-либо функция последования, отвечающая этому циклу.

Допустим, что при некотором натуральном k в окрестности U нуля в E определена k-ая степень функции последования, причем, существует q < 1 такое, что при всех x ∈ U:

≤ q (4)

Тогда цикл Tφ орбитально асимптотически устойчив.

Один достаточный признак асимптотической орбитальной устойчивости предельного цикла дает формулируемая ниже теорема Андронова — Витта. Ее можно рассматривать как аналог теоремы Ляпунова об устойчивости по первому приближению.

Пусть f непрерывно дифференцируема и пусть φ — нестационарное T-периодическое решение системы (1). Наряду с уравнением (1) рассмотрим линеаризованное вдоль решения φ уравнение:

x′ = A(t)x, (5)

где .

Очевидно матрица A является T-периодической и поэтому, в частности, система (5) является правильной по Ляпунову. К сожалению, в случае нестационарного решения φ условия теоремы Ляпунова об устойчивости по первому приближению никогда не бывают выполненными.

Действительно, как показывает дифференцирование тождества

φ′(t) ≡ f[φ(t)] по t, (ненулевая) периодическая функция ψ(t) = φ′(t) является решением линеаризованного уравнения (5). Поэтому нулевое решение этого уравнения не является асимптотически устойчивым и, следовательно, спектр системы (5) содержит нулевой характеристический показатель.

ТЕОРИЯ ФЛОКЕ

- теория о строении пространства решений и о свойствах самих решений линейной системы дифференциальных уравнений с периодическими коэффициентами:

(1)

матрица A(t) периодическая по с периодом и суммируемая на каждом компактном интервале из .

1) Любая фундаментальная матрица системы (1) имеет представление:

(2)

называется представлением Флокe (см. [1]), где F(t) - некоторая -периодическая матрица, К - некоторая постоянная матрица. Существует базис , ..., пространства решений системы (1) такой, что в этом базисе матрица. имеет жорданову форму; этот базис можно представить в виде:

где - многочлены относительно t с - периодическими коэффициентами, - характеристические показатели системы (1). Любая компонента решения системы (1) является линейной комбинацией функций вида (решений Флоке) В случае когда все характеристич. показатели различны (или среди них есть кратные, но им отвечают простые элементарные делители), функции суть просто - периодические функции. В представлении (2) матрицы и К, вообще говоря, комплекснозначные. Если ограничиться только действительным случаем, то F(t) может не быть - периодической, но обязательно будет - периодической.

2) Систему (1) можно привести к дифференциальному уравнению с постоянной матрицей у'=Ку с помощью преобразования Ляпунова:

(3)

где F(t) и K -матрицы из представления Флоке (2) (см. [2]). Представление (2) вместе с подстановкой (3) часто называют теоремой Флоке - Ляпунова.

3) Пусть - спектр матрицы К. Для каждого такого, что в силу представления (2) пространство распадается в прямую сумму двух подпространств и :

таких, что

здесь V(t) - нормированная в нуле фундаментальная матрица системы (1). Отсюда следует экспоненциальная дихотомия системы (1), если

Re

ни для какого j= 1,...,i.

Теорема Андронова — Витта. Пусть нулевой характеристический показатель линеаризованной системы (5) однократен, а остальные характеристические показатели этой системы отрицательны. Тогда решение φ орбитально асимптотически устойчиво.

Более того, для любого решения ψ, начинающегося достаточно близко к траектории Tφ, найдется число τ (называемое асимптотической фазой решения ψ) такое, что → 0 при t → ∞.

Заметим, что наличие асимптотической фазы означает, что при t → +∞ решения "наматываются" на предельный цикл Tφ как спирали.