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

Хаотическая динамика и классическая теория устойчивости движения

.pdf
Скачиваний:
0
Добавлен:
15.08.2026
Размер:
1 Мб
Скачать

МАТРИЧНОЕ УРАВНЕНИЕ ЛЯПУНОВА

61

Отсюда следует, что при a > κ

V (t) V (0) exp[2(a − κ)t].

Но тогда по теореме 5.1 получим оценку

|x(t)|2 ε1V (0) exp[2(a − κ)t].

 

Таким образом, здесь получен следующий результат.

Предложение 5.3. Для экспоненциальной устойчивости системы (5.8) с функцией f (t, x), удовлетворяющей условию (5.12), достаточно, чтобы выполнялось неравенство

Λ + κ < 0.

ГЛАВА 6

Критерии устойчивости по первому приближению

Приведем здесь наиболее известные критерии устойчивости по первому приближению для сиcтемы

dx

= A(t)x + f (t, x), t 0,

n

 

(6.1)

dt

x R

.

Здесь A(t) — непрерывная, ограниченная при t 0 (n × n)-матри- ца, f (t, x) — непрерывная вектор-функция, удовлетворяющая в некоторой окрестности Ω(0) точки x = 0 условию

|f (t, x)| κ|x|ν , t 0, x Ω(0).

(6.2)

Здесь κ и ν — некоторые положительные числа, ν 1.

Системой первого приближения будем называть систему

dx

= A(t)x.

(6.3)

dt

 

 

Предположим здесь существование такого числа C > 0 и кусоч- но-непрерывной функции p(t), что для матрицы Коши X(t)X(τ )1 системы (6.3) выполнена оценка

|X(t)X(τ )1 | C exp

t

p(s) ds, t τ 0.

τ

КРИТЕРИИ УСТОЙЧИВОСТИ ПО ПЕРВОМУ ПРИБЛИЖЕНИЮ

63

Теорема 6.1. Если выполнено условие (6.2) с ν = 1 и неравен-

ство

t

1

lim p(s) ds + Cκ < 0,

t→+t

0

то решение x(t) 0 системы (6.1) асимптотически устойчиво по Ляпунову.

Из теоремы 6.1 сразу получаем, что для уравнения первого порядка из отрицательности характеристического показателя системы первого приближения следует асимптотическая устойчивость нулевого решения.

Будем теперь предполагать, что для матрицы Коши X(t)X(τ )1 системы (6.3) имеет место оценка

|X(t)X(τ )1 | C exp[−α(t − τ ) + γτ ], t τ 0, (6.4) где α > 0, γ 0. Об этой оценке уже шла речь в главах 2 и 5.

Теорема 6.2 (Персидский [44, 52]). Если выполнены условия (6.4) с γ = 0, (6.2) с ν = 1 и достаточно малым κ, то решение x(t) 0 системы (6.1) асимптотически устойчиво по Ляпунову.

Теорема 6.2 является следствием теоремы 6.1 при p(t) ≡ −α.

Теорема 6.3 ((Четаев [65], Малкин [9, 44], Массера [45])). Если выполнены условия (6.4), (6.2) и неравенство

(ν − 1)α − γ > 0,

(6.5)

то решение x(t) 0 системы (6.1) асимптотически устойчиво по Ляпунову.

Теорема 6.3 усиливает широко известную теорему Ляпунова об устойчивости по первому приближению для правильных систем [42].

Для доказательства теорем 6.1 и 6.3 нам потребуются хорошо известные леммы Беллмана – Гронуолла и Бихари [12].

64 ГЛАВА 6

Лемма 6.1 (Беллмана – Гронуолла). Если для неотрицательных и непрерывных при t 0 функций u(t) и v(t) и числа C > 0 выполнено неравенство

t

 

 

t 0,

 

u(t) C +

v(τ )u(τ ) dτ,

(6.6)

0

 

 

 

то имеет место оценка

 

 

 

 

t

 

 

u(t) C exp[

 

t 0.

 

v(τ ) ],

(6.7)

 

0

 

 

Лемма 6.2 (Бихари). Если для неотрицательных и непрерывных при t 0 функций u(t) и v(t) и чисел ν > 1, C > 0 выполнены неравенства

u(t) C +

t

 

t 0,

 

v(τ )[u(τ )]ν dτ,

(6.8)

 

 

0

 

 

 

 

 

 

 

(ν

 

1)Cν−1

 

t

 

t

 

0,

(6.9)

v(τ ) dτ < 1,

 

 

 

 

 

 

 

 

 

0

то имеет место оценка

u(t)

 

C[1

 

(ν

 

1)Cν−1

 

t

 

t

 

0. (6.10)

 

v(τ ) ]1/(ν−1) ,

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ДОКАЗАТЕЛЬСТВО ЛЕММ 6.1 И 6.2.

Введем функцию Φ(u) = u для леммы 6.1 и Φ(u) = uν для леммы 6.2. В силу возрастания Φ из неравенств (6.6) и (6.8) имеем

КРИТЕРИИ УСТОЙЧИВОСТИ ПО ПЕРВОМУ ПРИБЛИЖЕНИЮ

65

оценку

t

Φ(u(t)) Φ(C + v(τ )Φ(u(τ )) ).

0

Отсюда получим неравенство

v(t)Φ(u(t))

t

v(t).

Φ(C + v(τ )Φ(u(τ )) )

0

t

Обозначая w(t) = C+ v(τ )Φ(u(τ )) и интегрируя последнее нера-

0

венство в пределах от 0 до t, получим оценку

 

 

 

t

w˙ (s)

 

 

 

t

 

 

 

 

 

 

ds

 

 

 

 

 

 

 

 

 

 

 

 

v(s) ds.

 

 

 

 

 

 

 

 

 

 

 

 

0

Φ(w(s))

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Последнее неравенство можно записать в виде

 

 

 

 

 

w(t)

 

 

 

 

t

 

 

 

 

 

 

 

 

 

exp

v(s) ds

 

 

 

 

 

 

 

 

 

 

 

 

 

w(0)

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при Φ(u) = u и в виде

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

1

 

 

t

 

 

 

 

 

 

 

 

v(s) ds

 

 

 

 

 

 

 

 

1 − ν

 

w(t)ν−1

w(0)ν−1

0

при Φ(u) = uν .

Учитывая условие u(t) w(t), w(0) = C и (6.9), отсюда полу-

чим утверждения лемм 6.1 и 6.2.

 

 

ДОКАЗАТЕЛЬСТВО ТЕОРЕМЫ 6.1.

 

 

Представим решение системы (6.1) в следующем виде:

 

x(t) = X(t)(x(0) +

t

 

X(τ )1 f (τ, x(τ )) ).

(6.11)

0

66 ГЛАВА 6

Из соотношения (6.11) и условий теоремы вытекает оценка

t

 

t

t

 

 

 

 

 

 

 

|x(t)| C exp p(s) ds|x(0)| + C

exp(

p(s) ds)κ|x(τ )| dτ.

 

0

0

 

τ

 

 

Эту оценку можно переписать следующим образом:

 

t

 

t

τ

 

 

 

 

 

 

p(s) ds)|x(τ )|) dτ.

(exp(− p(s) ds)|x(t)|) C|x(0)|+

exp(

0

 

0

0

 

 

Применяя лемму 6.1, получим оценку

 

 

 

 

 

t

 

 

 

 

 

 

 

t 0.

 

|x(t)| C|x(0)| exp[

p(s) ds + Cκt],

 

0

 

 

 

 

 

Из этого неравенства следует утвеждение теоремы.

 

ДОКАЗАТЕЛЬСТВО ТЕОРЕМЫ 6.3.

 

 

 

 

 

Из соотношений (6.2), (6.4) и (6.11) вытекает оценка

 

t

 

 

 

 

 

|x(t)|) Ce−αt|x(0)| + C

exp[−α(t − τ ) + γτ ]κ|x(τ )|ν dτ.

 

0

 

 

 

 

 

Эту оценку можно переписать следующим образом:

t

(eαt|x(t)|) C|x(0)| + exp[((1 − ν)α + γ)τ ](eατ |x(τ )|)ν dτ.

0

Применяя лемму 6.2, получим оценку

|x(t)| C|x(0)| exp[−αt][1 (ν − 1)(C|x(0)|)(ν−1) ()×

 

t

 

 

 

 

 

 

 

exp[((1

ν)α + γ)τ ] ]1/(ν−1)

,

 

t

 

0.

×

 

 

 

 

 

0

КРИТЕРИИ УСТОЙЧИВОСТИ ПО ПЕРВОМУ ПРИБЛИЖЕНИЮ

67

Из этого неравенства и из условия (6.5) при достаточно малых |x(0)|

получим утверждение теоремы.

 

 

 

 

Рассмотрим систему

 

 

 

 

dx

t 0, x R

n

 

 

dt = F (x, t),

 

,

(6.12)

где F (x, t) — дважды непрерывно дифференцируемая вектор-функ- ция.

Будем здесь предполагать, что для решений системы (6.12) с начальными данными y = x(0, y) из некоторой области Ω выполнено

следующее условие.

 

 

 

 

Максимальное сингулярное

число

α1(t, y) фундаментальной

матрицы X(t, y) линейной системы

 

 

dz

= A(t)z

 

(6.13)

dt

 

 

 

 

удовлетворяет неравенству

 

 

 

 

α1(t, y) α(t),

 

t 0,

y Ω.

(6.14)

Здесь матрица

∂F (x, t)

A(t) = ∂x |x=x(t,y)

является матрицей Якоби вектор-функции F (x, t) на решении x(t, y),

X(0, y) = I.

Теорема 6.4. [30]. Пусть функция α(t) ограничена на интервале (0, +). Тогда решение x(t, y), y Ω, устойчиво по Ляпунову.

Если, кроме того,

lim α(t) = 0,

t→+

то решение x(t, y), y Ω, асимптотически устойчиво по Ляпунову.

∂x(t, w)
∂w

68

ГЛАВА 6

ДОКАЗАТЕЛЬСТВО.

Хорошо известно, что

∂x(t, y)

= X(t, y), t 0.

∂y

Также хорошо известно [18], что для любых векторов y, z и числа t 0 существует вектор w, для которого выполнены соотношения

|w − y| |y − z|,

|x(t, y) − x(t, z)| | ||y − z|.

Поэтому для любого вектора z из шара с центром в y и целиком расположенного в Ω справедлива оценка

|x(t, y)−x(t, z)| |y−z| sup α1(t, w)| α(t)|y−z|, t 0. (6.15)

Здесь супремум берется по всем w из шара {w| |w − y| |y − z|}. Из оценки (6.15) сразу следуют утверждения теоремы 6.4.

Обсудим теперь условия теоремы 6.4. Теорема 6.4 устанавливает асимптотическую устойчивость решений с начальными данными из Ω, если соответствующие уравнения (6.13) имеют отрицательные ляпуновские экспоненты (или отрицательные характеристические показатели). Здесь требование равномерной по Ω отрицательности ляпуновских экспонент заменяет требование малости коэффициента неправильности в теореме 6.3.

Таким образом, перроновские эффекты, описанные в главе 4, возможны только на границах устойчивого по первому приближению потока.

В самом деле, рассмотрим поток решений системы (4.1) с начальными данными в некоторой окрестности точки x1 = x2 = 0:

x1(0, x10, x20 ) = x10, x2(0, x10 , x20) = x20.

Легко видеть, что

x1(t, x10 , x20) = exp[−at]x10.

КРИТЕРИИ УСТОЙЧИВОСТИ ПО ПЕРВОМУ ПРИБЛИЖЕНИЮ

69

Поэтому матрица A(t) системы (6.13) в этом случае примет вид

 

A(t) =

−a

0

.

 

2 exp[−at]x10

[sin(ln(t + 1)) + cos(ln(t + 1)) 2a]

 

Решениями системы (6.13) будут функции

z1(t) = exp[−at]z1(0)

z2(t) = exp[(t + 1) sin(ln(t + 1)) 2at](z2(0) + 2z1(0)x10 ×

t

×exp[(τ + 1) sin(ln(τ + 1))] ).

0

Как было показано в главе 4, при выполнении соотношений (4.2) и при

z1(0)x10 = 0

характеристический показатель z2(t) является положительным. Следовательно, в сколь угодно малой окрестности тривиально-

го решения x1(t) ≡ x2(t) 0 существуют начальные условия x10, x20 такие, что система первого приближения для x1(t, x10, x20), x2(t, x10, x20) имеет положительный старший характеристический показатель (и ляпуновскую экспоненту μ1).

Таким образом, эффект Перрона для системы (4.1) имеет место в силу того, что здесь нет равномерной по Ω оценки (6.14). Здесь Ω — некоторая область, содержащая точку x1 = x2 = 0.

ГЛАВА 7

Критерии неустойчивости

Рассмотрим систему

 

 

 

dx

= A(t)x + f (t, x), t 0,

n

 

(7.1)

dt

x R

,

где A(t) — непрерывная, ограниченная на [0, ∞) (n ×n)-матрица. Будем предполагать, что вектор-функция f (t, x) непрерывна и в некоторой окрестности Ω(0) точки x = 0 выполнено неравенство

|f (t, x)| κ|x|ν , t 0, x Ω(0).

(7.2)

Здесь числа κ и ν таковы,что κ > 0, ν > 1.

Введем в рассмотрение нормальную фундаментальную матрицу

Z(t) = (z1(t), . . . , zn(t)),

(7.3)

состоящую из линейно независимых решений zj (t) линейной системы первого приближения

dz

= A(t)z.

(7.4)

dt

 

 

Для анализа неустойчивости одним из основных подходов является приведение линейной части системы (7.1) к треугольному виду. Здесь наиболее эффективным оказывается метод Перрона – Винограда [9, 12] триангуляции линейной системы. Опишем здесь этот метод.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]