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

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

&

x

(t 0).

 

 

 

 

x = − t +1

 

Его решение в форме Коши имеет вид

x(t) = x0 (t0 +1)

. Следовательно, нулевое

 

 

 

t +1

 

решение ξ = 0 асимптотически устойчиво. Однако оно не является экспонен-

циально устойчивым, поскольку не существует положительных N и α, для которых неравенство

| x(t) |=| x0 (tt +0 1+1) | N | x0 | eα(tt0 )

справедливо при x0 0 и всех t 0.

Теорема 1. Предположим, что существует положительно определенная квадратичная форма V(x) = xT P x , производная которой в силу системы (1) при некотором h > 0 удовлетворяет неравенству

&

t 0 , x : x < h < H ,

V(x(t)) W(x(t))

где W(x) = −xTQ x – отрицательно определенная квадратичная форма и P, Q постоянные симметричные матрицы. Тогда нулевое решение системы (1) экспоненциально устойчиво.

□ Обозначим через λj (P) , λj (Q) ,

 

j =1,2,...,n , собственные значения

матриц P и Q соответственно.

 

 

Введём числа

 

 

a = min λj (P),

a1

= max λj (P) ,

j

 

j

b = min λj (Q),

b1

= max λj (Q) .

j

 

j

Очевидно,

0 < a a1, 0 < b b1 .

Из курса линейной алгебры известно, что

a x 2 V(x) = xT P x a1 x 2 , b x 2 ≤ −W(x) = xT Q x b1 x 2 .

Далее

&

T

 

 

2

 

b

 

 

 

 

V(x) W(x) = −x

 

Q x ≤ −b

x

 

≤ −

a1

V(x) .

 

 

 

 

 

 

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

Пусть здесь x = x(t) − произвольное ненулевое решение системы (1) с начальными данными в области Z0. Интегрируя последнее неравенство в пределах от t0 0 до t, находим

 

 

 

 

 

 

 

V(x(t)) V(x(t0 )) e

2α(t t0 )

t t0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

α =

 

b

 

> 0 . Но

V(x) a

 

 

 

x

 

 

 

2 ,

V(x) a1

 

 

 

x

 

 

 

2 . Поэтому

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2a1

|| x(t)||2

1

V(x(t))

a1

|| x(t0 ) ||2 e2α(t t0 )

 

 

 

t t0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А значит,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|| x(t)|| N || x(t0 ) || e

α(t t0 )

t

t0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

N =

 

a1

и норма || x(t0 ) ||

 

 

 

мала настолько,

что

 

 

 

x(t0 )

 

 

 

< h и

 

 

 

x(t)

 

 

 

< H

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

для всех t ≥ t0

20. Экспоненциальная устойчивость линейных нестационарных сис-

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

dx

= A(t) x

(t 0)

(5)

dt

 

 

 

с непрерывной на промежутке [0, +∞) матрицей A(t). Для таких систем нередко вводят несколько более «жёсткое» понятие экспоненциальной устойчивости, чем это было сделано выше.

О п р е д е л е н и е 2. Линейную систему (5) называют экспоненциально устойчивой, если существуют такие положительные константы α1 2 12 ,

что для всех начальных данных t0 0, x0 Rn выполняются неравенства

β1 x0 eα1 (tt0 ) x(t; t0 , x0 )β2 x0 eα2 (tt0 ) t t0 . (6)

Нетрудно понять, что из экспоненциальной устойчивости системы (5) вытекает её асимптотическая устойчивость (но не обратно). В случае, когда матрица A(t) = A постоянна, оба указанных понятия эквивалентны.

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

x& = −2t x , обладающее решением x(t) = x0 et2 , для которого, как легко проверить, не будет выполняться левое из неравенств (6).

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

В следующей теореме представлен критерий экспоненциальной устойчивости системы (5).

Теорема 2. Для того чтобы система (5) была экспоненциально устойчивой, необходимо и достаточно, чтобы существовали две квадратичные формы

V(t,x) = xT P(t) x , W(t,x) = xT Q(t) x , удовлетворяющие следующим двум условиям

1) для некоторых положительных констант

a1,a2 , b1, b2 при всех t 0,

x Rn имеют место неравенства

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a1

 

 

 

x

 

 

 

2 V(t, x) a2

 

 

 

x

 

 

 

2 ,

b1

 

 

 

x

 

 

 

2

W(t,x) b2

 

 

 

x

 

 

 

2 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2)матрица P(t) непрерывно дифференцируема на промежутке [0, +∞) и справедливо равенство

dV

 

= −W(t,x) .

(7)

dt

 

(5)

 

 

 

□ Необходимость. Любое решение системы (5) с начальными данными

(0, x0) можно представить в виде

 

 

 

x(t) = X(t) X1 (0) x0

(t 0),

где X(t) – фундаментальная матрица данной системы. Введем квадратичную форму W(t,x) := xT x и функцию

+∞

+∞

 

V(t,x) := W(τ,x)dτ = xT xdτ.

(8)

t

t

 

Очевидно, квадратичная форма

W(t, x) удовлетворяет условию 1) доказывае-

мой теоремы при b1 = b2 =1 и, кроме того, выполнено равенство (7).

Остается убедиться в том, что функция V(t,x),

будучи квадратичной

формой, удовлетворяет условию 1). В самом деле, используя (8), при произвольно выбранном решении x = x(t) = x(t; t0 , x0 ) системы (5) получаем

+∞

+∞

V(t,x) = x(τ)T x(τ)dτ = x0T (X1(0))T XT (τ) X(τ) X1(0)dτ x0T =

t0

t0

+∞

= xT0 (X1(0))1 XT (t) ((X1(t))T XT (τ) X(τ) X1(t)dτ

t0

+

X(t) X1(0) xT0 = xT BT (τ, t) B(τ, t)dτ xT ,

t0

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

где B(τ, t) = X( ) X1 (t) . Изτполученного представления для V(t,x) благодаря произвольности выбора начальных данных (t0 , x0 ) (а значит, и произвольности

x) следует,

что функция V(t,x) является квадратичной формой с матрицей

+∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P(t) = BT (τ, t) B(τ, t)dτ.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из (6) для всех x0 Rn вытекает

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(β )2

 

 

 

x

0

 

 

 

2 e1 t

 

 

 

x(t;0,x

0

)

 

 

 

2

2

)2

 

 

 

 

 

 

 

x

0

 

 

 

2 e2 t

 

 

 

 

t 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Интегрируя эти неравенства в пределах от 0

 

 

до +∞, получим

 

 

 

 

 

 

 

 

 

2

+∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+∞

 

 

 

 

 

2 2 )2

 

 

 

 

 

 

 

2

+∞

1 )2

 

 

 

x0

 

 

 

e1 t0 ) V(0, x0 ) =

 

 

 

x(τ)

 

 

 

 

 

 

x0

 

 

 

e2 t0 )dτ.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

Поскольку

 

 

 

 

 

 

 

 

 

 

 

 

 

+∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ei τ

=

 

 

 

 

,

 

 

 

 

 

 

 

 

 

i = 1,2,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

полученное неравенство в силу произвольности

x0 Rn означает, что квадра-

тичная форма V(t,x) удовлетворяет условию 1).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Достаточность. Пусть существуют квадратичные формы V и W, удовле-

творяющие условиям 1) и 2) теоремы. Тогда применимо следствие 1 из § 3, в соответствии с которым линейная система (5) является асимптотически устой-

чивой.

Остается показать, что при условии t0

 

0,

 

x0 Rn

любая интегральная

кривая

x = x(t; t0 , x0 )

системы (5)

 

удовлетворяет неравенствам (6).

 

Интегрируя равенство (7), предварительно делённое на V, в пределах от

t0 до t, получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

W(τ,x(τ))

 

 

V (t;t

 

,x

 

) := V(t,x(t;t

 

,x

 

)) = V(t

 

,x

 

 

 

t0

V(τ,x(τ) )

(9)

 

0

0

0

0

0

0

) e

 

,

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где функция

V1, записанная в начале строки, является значением функции

V

на выбранной интегральной кривой x(t; t0 , x0 )

 

системы (5).

 

Рассмотрим подынтегральное выражение из (9):

 

 

 

 

 

 

 

 

 

 

 

W

=

xT Q(t) x

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

 

xT P(t) x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Используя условие 1), для него можно записать

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