Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
TUD_Nogin / TUD / part2.pdf
Скачиваний:
122
Добавлен:
16.04.2015
Размер:
808.23 Кб
Скачать

ГЛАВА 2. ВТОРОЙ МЕТОД ЛЯПУНОВА_______________________________________________________102

.

 

 

 

dV(ξ(t; x))

 

 

 

 

d

+

 

 

2

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

V(x)

 

 

=

 

 

 

 

 

 

=

 

 

 

 

||ξ(τ; x)|| dτ

 

 

 

 

 

=|| x|| .

 

 

 

 

 

 

 

 

 

 

 

= −|| ξ(τ; x)||

t = 0

 

 

(2)

 

 

dt

 

 

 

t = 0

dt

 

 

 

t

 

 

 

 

 

t = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теперь вычислим производную функции V(x) по t

в точке x в силу

системы (1):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

+ V, ϕ(t, x) =

 

 

 

 

 

 

V(x)

= V, A x + V, ϕ(t, x) = V(x)

 

 

 

 

 

 

 

(***)

(1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= − x 2 + V, ϕ(t, x) .

 

 

 

 

 

 

Напомним, что x – вектор-столбец. Благодаря

 

 

 

 

 

 

 

 

 

 

 

 

 

V(x) = S x, x = Sjkx jxk

 

 

 

(Sjk =Skj) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j,k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

имеем

 

 

 

 

 

 

 

 

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S1kxk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 2S x .

 

 

 

 

 

 

 

 

 

 

 

 

 

V(x) = M

 

 

= 2

 

M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

 

 

 

 

 

S

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

nk

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

||ϕ(t, x) ||

 

Кроме того, из равномерной по

 

t

0

сходимости

 

 

x

 

0

при

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x 0 получаем, что для

любого положительного

ε найдётся

 

 

hε > 0, при ко-

тором

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

||ϕ(t, x)|| < ε||x||

 

 

 

t 0, x : ||x|| hε <

 

H .

 

 

 

 

 

 

Следовательно, из (***)

 

следует

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

≤ −|| x||2 + |

V,

ϕ(t, x)

|≤ −|| x||2 +2 || S|| || x|| ε || x|| =

 

 

 

 

V(x)

 

 

 

 

 

 

 

(1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= −|| x||2 (12ε || S||)< −

1

|| x||2 < 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

если только

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 < ε <

 

, 0 <|| x|| < hε,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 || S||

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ГЛАВА 2. ВТОРОЙ МЕТОД ЛЯПУНОВА_______________________________________________________103

.

= 0.

причём V(0)

(1)

В итоге установлено, что для системы (1) в некоторой окрестности начала координат существует положительно определенная функция V(x) , не завися-

щая от t и допускающая отрицательную производную в силу системы. Согласно теореме Ляпунова об асимптотической устойчивости нулевое решение системы (1) асимптотически устойчиво ■

З а м е ч а н и е 1. Очевидно, для выполнения условия равномерной сходимости (2) достаточно, чтобы при некоторых положительных a и α имело место неравенство

||ϕ(t, x) ||a

 

 

 

x

 

 

 

1+α

t 0, x :

 

 

 

x

 

 

 

< H .

 

 

 

 

 

 

 

 

Из доказанной теоремы вместе с теоремой 1 предыдущего параграфа получаем следующее утверждение.

З а м е ч а н и е 2. В случае, когда вещественные части всех собственных значений матрицы A неположительны и среди них есть собственные числа с нулевой вещественной частью, по линейному приближению на основе теоремы 1 невозможно сделать вывод об устойчивости нулевого решения системы (1). Такие случаи именуют критическими.

Следствие 1. Рассмотрим нелинейную систему

 

 

dx

 

= f(t, x) ,

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

в которой f(t,0) 0

и векторная функция f(t, x)

дважды непрерывно диф-

ференцируема по x

в области t 0,

 

x

 

 

 

H , причём производные

 

 

 

 

2f(t, x)

,

 

 

 

 

 

 

i, j =1,2,...n ,

 

 

 

xix j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ограничены в указанной области, а матрица f(t,0)

x

постоянна. Тогда сис-

 

 

 

 

 

 

 

 

A = f(t,0)

 

тема (3) представима в виде (1), где

 

 

, и для неё справедлива

теорема 1.

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следствие 2. Рассмотрим нелинейную автономную систему

ddty = f(y) , f(y) C2 ({y | yH}) ,

ГЛАВА 2. ВТОРОЙ МЕТОД ЛЯПУНОВА_______________________________________________________104

у которой y = y0 – состояние равновесия, т.е.

 

f(y0 ) = 0 ( ||y0 ||< H) . Если все

собственные значения матрицы Якоби

 

 

 

A = f(y0 ) = (fjk(y0 ))

n x n

 

 

имеют

отрицательные вещественные части,

то положение равновесия

y = y0

данной системы асимптотически устойчиво по Ляпунову.

□ Положим y = y0 +x . Тогда

 

 

ddty = f(y) = f(y0 ) +f(y0 ) x +o(|| x||) = A x +o(|| x||) = ddtx

и можно применить доказанную выше теорему 1 ■

20. Теорема Ляпунова о неустойчивости квазилинейной системы.

Теорема 2. Пусть дана нелинейная система

 

 

 

 

dx

= A x +ϕ(t, x) ,

 

 

(1)

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

где А – постоянная матрица и ϕ(t, x) C {(t, x)| t 0, ||x|| < H}

)

, причём

 

 

 

(

 

 

||ϕ(

t,

x) ||

0 равномерно относительно t 0 при x 0 .

 

x

 

 

 

 

 

 

 

 

Если хотя бы одно из собственных значений λ1,...,λn

матрицы

А имеет по-

ложительную вещественную часть, то нулевое решение системы (1) неустойчиво.

□ Пусть

 

Re λ1,...,Re λm > 0 ; Re λm+1,...,Re λn 0

(1 m n) .

В курсе линейной алгебры доказывается, что матрица А неособенным линейным преобразованием может быть приведена к почти диагональному виду, т.е. существует неособая матрица S (det (S) 0) :

S1 A S = Λ+B ,

где Λ = diag (λ1, ... ,λn ) , B = (bjk )n x n , bjk = 0 для всех j k и ||B|| ε, причём число ε > 0 может быть выбрано сколь угодно малым. Заметим, что матрицы S и B, вообще говоря, комплексные.

Выполним замену

 

x = eαt S y

(0 < α < min Re λj ),

 

1jm

ГЛАВА 2. ВТОРОЙ МЕТОД ЛЯПУНОВА_______________________________________________________105

где вектор y, вообще говоря, комплексный. Имеем

 

 

 

 

eαt S

 

dy

+

eαt α S y = eαt A S y +ϕ(t, eαt S y) ,

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= (ΛαE) y +B y +ψ(t, y) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(*)

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

ψ(t, y) = eαt S1 ϕ(t,eαt S y)

и

матрица

 

ΛαE

не имеет собственных

значений с нулевой вещественной частью.

 

 

 

 

||ϕ(

t,

 

x) ||

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из условия равномерной сходимости

 

 

 

 

 

0

 

 

 

 

относительно t 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при

x 0, следует, что для любого

ε > 0

 

найдётся такое положительное

hε,

что выполнено неравенство

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

||ϕ(t, x)|| ε||x||

 

 

 

 

 

 

t 0,

x :

 

 

 

x

 

 

 

 

 

 

 

< hε

< H .

 

 

 

 

 

 

 

 

 

 

 

 

 

Поэтому

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ(t, y)

 

 

 

 

 

eαt

 

 

 

S1

 

 

 

ε eαt

 

 

 

S

 

 

 

 

 

 

 

y

 

 

 

ε

 

 

 

S1

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

y

 

 

 

 

(**)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

для всех x, удовлетворяющих

 

 

||x|| eαt ||S || || y|| hε, что равносильно неравен-

ству

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

||y|| eαt || S||-1 hε

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(***)

для всех достаточно малых

ε > 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Положим

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y =

 

1

 

 

 

 

 

 

 

μs = λs α,

s

 

=1,2,...,n ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и перепишем систему (*) следующим образом

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy1

 

=μ y

+ b

y

+... + b

y

 

+ψ (t, y),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

1 1

12

2

 

 

 

 

 

 

 

 

 

 

 

 

1n

 

n

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy2

=

 

 

 

 

 

 

 

 

 

μ2y2 +... + b2n yn +ψ2 (t, y),

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

LLLLLLLLLLLLLLLLL

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dyn

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

μn yn +ψn (t, y),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Соседние файлы в папке TUD