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

книги из ГПНТБ / Зубов В.И. Математические методы исследования систем автоматического регулирования

.pdf
Скачиваний:
32
Добавлен:
25.10.2023
Размер:
12.11 Mб
Скачать

Выберем строго монотонную функцию L (т), заданную и непре­

рывную

при —схэ <

т < + °°, строго монотонно

убывающую

от + ° °

до 0 при т,

возрастающем от —00 до + °°,

и такую, что

L (t t0) 5s X (t t0). Легко видеть, что построенная таким обра­ зом функция удовлетворяет всем требованиям пункта 1 нашего плана.

Функция W является определенно-отрицательной и допускает

бесконечно малый высший предел, поэтому остается

показать,

что функция

 

 

 

 

 

 

-1-0

 

 

 

 

 

 

 

 

 

 

 

 

 

V(y10»

 

УпОу

U)

Wdx

 

(3.16а)

 

 

 

 

 

 

 

to

 

 

существует и обладает указанными свойствами.

 

 

По числу е' >

0 найдем, как и выше,

такое 6 (е') >> 0, что при

П

 

 

 

 

 

 

 

 

 

Xi у1 < §2 (е')

будет

 

 

 

 

 

 

 

I=1

 

П

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

И У1 (t,

г/10, •

• ., Упо, to) < е'2.

 

 

 

 

1=1

 

 

 

 

 

 

 

Для всех интегральных кривых с такими начальными данными

будет иметь место неравенство

 

 

 

 

 

 

е~~V—to) >

- к

 

 

(3.17)

 

 

е

9

 

которое получается из неравенства (**)

путем применения обрат­

ной функции К к обеим частям.

 

(3.15) получаем

Из неравенства (3.17)

и

равенства

 

 

^

101

• •

мУпо* to)

 

 

 

<

-f00

JП

2

У ю .

 

Упо?.() тe ,~ (x~ U) d x ^

е

' 2 .

 

to

i=1

 

 

 

 

 

 

 

Из этого неравенства следует, что, во-первых, интеграл от правой части (3.15) сходится и в действительности определяет функцию V и, во-вторых, что эта функция допускает бесконечно малый высший предел, так как число е' можно считать произ­ вольно малым.

Покажем теперь, что функция V является положительно-опре­ деленной. Из (3.16а) следует

и+т

п

tfe2

—к ( 2

у])

 

V(y10, . . . , y n0, t ) ^ J

t =l

 

’dx,

(3.18)

to

 

 

 

 

что получается из (3.16a) путем отбрасывания положительного

члена, а именно —

соJ W dx. Вследствие равномерной притягивае-

 

to + T

40

мости нулевого

решения

системы

(1.4) по числу h > О можно

указать такие

Т > О и

а > 0,

что

 

П

 

 

 

{t, Ухо, • • ч Упо, to)

при

1=1

 

 

 

 

 

S у)о: и

7.

1 = 1

 

А тогда подынтегральное выражение в неравенстве (3.18) остается не менее положительного числа у (а), так как функция W опреде­ ленно-отрицательна. Поэтому имеем

П

V (Ухо, ■■■, Упо, t0) ^ y T при £ y lo > h 2,

1=1

что и означает положительную определенность функции V.

Остается показать, что ~ ~ = W. Для этого возьмем точку

У\, . ■■, уп на какой-либо интегральной кривой, соответствующую моменту t0 + At. Тогда, согласно (3.15),

 

 

СО.

 

 

 

V (yu ... ,y 'n ,t 0 + At) = — J Wdx.

(3.19)

Вычитая почленно (3.16а) из (3.19) и деля

обе части на At,

находим

4Т = - Е Г

 

 

(3W)

 

^0I

 

где AV означает приращение функции V вдоль интегральной

кри­

вой. Устремляя At к 0, получим из (3.20)

= W вдоль

инте­

гральной кривой. Этим теорема 8 доказана полностью. Следствие. Пусть нулевое решение системы (1.4) асимптоти­

чески устойчиво. В этой системе всегда можно сделать преобразова­ ние независимого переменного t таким образом, чтобы нулевое решение вновь полученной системы было равномерно притяги­ вающим.

Действительно, в качестве новой переменной времени в си­ стему (1.4) введем длину интегральной кривой в п + 1-мерном пространстве у г. . . ., Уп, t:

ds =

"j/""dt

dyi -(-•••-)- dyn =

 

= d

t ^ 1 +

Д gt {t, УХ, ■■ч Уп)-

(3.21)

41

Тогда система (1.4) перейдет в систему вида

 

dyk _

 

gk (t, УЪ ■; Уп)

------

. 1

ds

 

 

 

 

у

 

 

 

 

 

 

]

/

1+

i

* ? . ( ' . ■ ■•>Уп)

 

1

 

 

f=i

l

 

 

dt

 

 

 

 

(3.22)

ds

 

 

 

s] (6 yv ■■■•Уп)

 

|

/

1+

£

 

при k =

1, ..., n.

 

.

Если нулевое решение системы (1.4) асимптотически устой­ чиво, то нулевое решение системы (3.22) является равномерно притягивающим, что можно установить тем же способом, как это было сделано в лемме.

Теорема 8 будет давать необходимые и достаточные условия равномерной асимптотической устойчивости для системы (3.22) без дополнительного требования равномерной притягиваемости.

Обратимся к признакам неустойчивости нулевого решения системы (1.4).

Теорема 9

Нулевое решение системы (1.4) неустойчиво, если существуют две функции V (уъ . . ., уп, t) и W (уъ .... ., уп, t), обладающие следующими свойствами:

1) функция V допускает бесконечно малый высший предел и принимает отрицательные значения в сколь угодно малой окрест­ ности полуоси (для любого фиксированного t > Т) у г = у 2 = . . .

■■

=

Уп = 0

и

f

0;

полной производной по t функции V

2) функция

W является

в силу системы

(1.4),

т. е.

W =

 

 

 

 

П

 

a1

W {уъ . . ., уя, t)

 

3)

при £ г/? ^

(/) 0, где сра (/) —

 

 

1=У

 

 

 

 

 

t

 

 

 

 

 

 

 

 

такая

интегрируемая функция,

что

[ фа (т) йт -*■—оо при

t

+°о.

 

 

 

 

 

6

 

 

 

 

 

 

Доказательство. Предположим,

что существуют две функции V

и W, удовлетворяющие условиям теоремы 9. Покажем тогда, что нулевое решение системы (1.4) неустойчиво.

Доказательство будем вести от противного. Предположим, что нулевое решение системы (1.4) устойчиво. Тогда для любого е > 0

 

 

П

можно найти такое число 8 >

0, что при

У20 < 62 будет

п

г=1

 

 

£ y2i(t> Уш • • •,

Упй* h)<&

при t ^ t 0.

1=1

 

 

42

Возьмем положительное число е' <С е. Тогда можно указать

 

 

 

 

П

 

 

 

 

 

 

такое число б',

что при

'£i y \ <

б' будет |

V {уъ . . .,

уп, 0

| < е'

 

 

 

 

1=1

 

 

 

 

 

 

(из-за наличия бесконечно малого высшего предела).

V (у 10, . . .

Выберем точки у 10, . . ., уп0,

t0 так, чтобы

было

 

 

 

П

 

 

 

 

 

 

 

■ ■ •> Упо,

h )

0> з

i=i У ю < С

б .

системы

(1.4)

начинается

при

Если

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

кривая

П

 

 

 

 

 

 

 

 

 

 

X Ую <

б2 и при некоторых t

 

t0 отстоит от

начала

координат

i ~ 1

 

 

 

 

t

будет | V1 (t, у 10, . . .,

ynQ} t0) I <

меньше, чем на б', то для этих

< е'.

 

 

 

 

 

 

 

 

 

 

Вдоль интегральных кривых системы (1.4) имеет место соот-

ношение

dV

W,

откуда

 

 

 

 

 

 

—тг =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

V\{t, Ую,

• • м

Упо,

to) =

V {Ую, • • •> Упо,

to) ~Ь J ^ dx.

(3.23)

 

 

 

 

 

 

 

 

^0

 

 

Положим е' = — V (у10, . . ., уп0, t0). Тогда из соотношения

(3.23) вытекает, что Vx {t, у 10, . . ., уп0, t0) —е' при t t0

и, следовательно, при этих же значениях независимой переменной будет иметь место неравенство

 

У1 (t,

Ую, • ■•> Упо, t0) ^ b .

(3.24)

 

г = 1

 

 

 

 

Из неравенства (3.24) следует, что W (уи . . .,

уп, t) sg фв'

при t

/0 на указанной

интегральной

кривой.

Применив

это

последнее неравенство к (3.23), получим,

что

 

 

 

 

 

t

 

 

Vi V, Ую,

• • ■, Упо, t0) ^ V ( y 10, . . ., уп0, f0) + JФб'(т)*/т. (3.25)

 

 

 

to

 

 

Из (3.25) вытекает, что Vx -»---- 00 при t -*•

а это невоз­

можно, так как функция

V допускает бесконечно малый высший

предел.

 

 

 

 

 

Полученное противоречие показывает, что теорема 9 верна, что и требовалось доказать.

Следствие. Если функция V (уи . . ., уп, t) удовлетворяет пер­ вому и второму условиям теоремы 9, а функция W будет опреде­ ленно-отрицательной, то нулевое решение системы (1.4) будет неустойчиво, так как вместо функции фа (t) можно брать отрица­ тельную постоянную. Утверждение, сформулированное в этом следствии, носит название первой теоремы Ляпунова о неустой­ чивости.

43

Теорема 10

Для того чтобы нулевое решение системы (1.4) было неустой­ чивым, необходимо и достаточно, чтобы существовали две функ­ ции V (у 1, . . уп, t) и W (уъ . . ., уп, t), обладающие следующими свойствами:

1) функция V имеет полную производную по t в силу системы (1.4), непрерывную на каждой интегральной кривой, удовлетво­ ряющую соотношению

 

 

 

= W + Г ,

где

Я = Я {t- уи . . ., уп);

(3.26)

 

2)

функция

W (уи

. . .,

уп,

t) — неотрицательна;

 

=

3)

в любой достаточно малой окрестности полуоси у х = у 2 =

. . . = уп =

0; t ^ 0 существуют точки,

в которых функция V

принимает положительные

значения при

любом фиксированном

t >

Т-

 

 

 

 

 

 

 

4)

существует такое число е > 0, что при любом б >■ 0 (6 << е)

неравенство

 

 

 

 

 

 

 

 

 

 

 

 

 

j X(т) dx

 

 

 

\V[yx(t), . . . , y n(t),

t ] \ ^ V ( y w, ...,

ynn,t0)eta

(3.27)

не сохраняется для всех t ^

t0 при любом выборе таких непрерыв­

ных функций у х (t), . .

., уп (t), что

 

 

 

 

 

П

 

 

 

 

 

 

 

 

Е

У21(0 < е2

и у, (t0) =

yi0t

 

где

 

 

i—1

 

 

 

 

 

 

 

П

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ъ у2о<

б2 И V (у10, . . . , упо, Q > 0.

 

 

 

 

»=1

 

 

 

 

 

 

Достаточность. Пусть существуют две функции V и

W, удо­

влетворяющие условиям теоремы 10. Покажем, что нулевое реше­ ние системы (1.4) неустойчиво.

Пусть это не так. Тогда по числу е >

0 можно указать такое

 

П

 

 

 

 

 

 

 

число б > 0, что при 2

У/о

 

б2 будет

 

 

 

 

/=1

 

 

 

 

 

 

 

Е У21

(*. У ю.

• •

•. Упо,

t0) <

е2

( t ^ t0).

(3.28)

i=1

 

 

 

 

 

 

 

 

Выберем y 10,

. . ., #„0,

t0 так,

чтобы

величина

V (у10, . . .

• • •> Упо> t0) была

положительной. Вычислим функцию V на инте­

гральных кривых, удовлетворяющих условию (3.28).

 

Тогда из (3.26) получим линейное дифференциальное уравнение

для функции 1/ х (/, у ш . . .,

упо, g

 

 

 

 

 

^ 1

=

ЯУ1 +

Ц7Ь

 

 

(3.29)

44

где \ГХозначает функцию W, вычисленную на той же интеграль­ ной кривой, что и V. Отбросив в (3.29) неотрицательный член W х

и проинтегрировав полученное неравенство в пределах от

до t,

получим

t

 

 

 

 

 

 

J k(x)dx

 

 

V l W, i/lOi • • •> ^/n0> ^o) ^ (l/lOi • •>

^o)e°

73= ^o))

что противоречит соотношению (3.27), так как в качестве функций yL(t) здесь можно рассматривать интегральную кривую, удовлет­ воряющую условию (3.28).

Таким образом, нулевое решение системы (1.4) неустойчиво. Необходимость. Пусть нулевое решение системы (1.4) неустой­ чиво. Покажем тогда, что существуют две функции V и W, удо­

влетворяющие условиям теоремы 10.

Если нулевое решение системы (1.4) неустойчиво, то существует

такое положительное число г',

что для любого 6 > 0 можно ука-

 

 

 

П

 

 

 

 

зать такую точку у ы,

. . .,

уп0,

tQ, что S

Ую <

Зз 0, и та-

 

П

 

t=l

 

 

 

 

 

 

 

 

 

 

кую, что неравенство

у% (t,

у ы, . . . .

уп0,

t0) >< е'2

не

выпол-

i=i

 

 

 

 

 

 

няется при всех t t0.

 

 

момент после t0, когда

Пусть t (у10, . . ., уп0, t0) есть первый

нарушается указанное

выше неравенство.

Разобьем

множество

 

П

 

б'2; t ^ O на два множества, отнеся

точек в полуцилиндре

^

г/2<

 

i=i

 

 

интегральные

кривые,

к первому те точки, из которых исходят

остающиеся при неограниченном возрастании времени в этом

полуцилиндре.

На этом множестве точек положим V (уъ . . .

уп , t) =

0 . К другому множеству точек отнесем все те, из

которых выходят интегральные кривые, пересекающие поверх­

ность

полуцилиндра

при возрастании

времени.

Если у 10, . . .

■■•>

Упо, такая точка,

то положим

 

 

 

 

 

V(y10, • •

.,

упо,

t0) = й*0~* («ю. ■

- W

o

) .

Ясно, что в первом случае

Ех (t, у 10,

. . .,

уп0,

t0) тождественно

равна нулю при t ^

t0.

 

 

 

 

 

 

Во втором случае

 

 

 

 

 

 

 

 

VAt,

Ую ,

Упо,

t0) =

 

 

 

ч

откуда следует,

что

 

= V, т. е. здесь X = 1,

a

W = 0. Функ-

 

 

 

 

 

 

t

 

 

 

ция V неотрицательна и ограниченна, a J X (т)

dx

 

расходится при

t t0 -► °°. так как

 

 

^0

 

 

 

X (т) = 1.

 

 

 

 

45

Следовательно, условие 3 теоремы 10 выполнено, что и требо­ валось доказать.

Следствие 1. Нулевое решение системы (1.4) неустойчиво, если:

1)

существуют две функции

V и W, удовлетворяющие первым

двум

условиям

теоремы 10;

 

2)

функция

V ограничена;

 

 

t

 

 

 

3)

J к dr ->

-j- 00

при

 

 

tо

 

зависит от

переменных у г, у 2, . . ., уп, то

Если функция к

интеграл вычисляется вдоль ограниченных решений. Действительно, все условия теоремы 10, за исключением усло­

вия 3, выполнены автоматически; но условие 3 также выполняется благодаря тому, что его левая часть ограничена при ограничен­ ности функции V, а правая часть не ограничена. Впрочем, можно потребовать лишь, чтобы

t

1im J X (т) dx = + оо при t —+оо. to

Следствие 2. Если выполнены условия 1 и 2 теоремы 10, функ­ ция V ограничена и функция к является положительной постоян­ ной, то нулевое решение системы (1.4) будет неустойчиво.

Доказательство этого утверждения вытекает из следствия 1. Утверждение, сформулированное в следствии 2, носит название второй теоремы Ляпунова о неустойчивости.

Следствие 3. Нулевое решение системы (1.4) неустойчиво, если существуют две функции V и W, удовлетворяющие следующим условиям:

1)обе функции удовлетворяют условию 3 теоремы 10;

2)функция V непрерывно дифференцируема вдоль интеграль-

ных кривых системы и dV = W;

3)

по числу е >■ 0 можно указать такое число б >■ 0, что при

V >■ е

будет W >• б.

Действительно, положим к = - у \ тогда

— = kV dt '

При Е0 > 0 будет ] / > V0, а следовательно,

[к (т) dx >оо при t -j- °о.

Отсюда из следствия 1 вытекает справедливость утверждения, которое носит название теоремы Четаева о неустойчивости.

46

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

правые части системы предполагаются непрерывными справа.

З а м е ч а н и е 2. В этом параграфе все теоремы относились к тому слу­

чаю,

когда правые части системы дифференциальных уравнений заданы в об-

 

п

ласти

0 и ^ (/? <1 г1. Однако в реальных системах процессы и соответствую-

1=1

щие им системы дифференциальных уравнений задаются лишь при t^-T. Тогда путем введения новой переменной времени т = t Т можно привести рассмо­

трение к системе вида (1.4).

З а м е ч а н и е 3. Рассмотренную здесь теорию можно также применять к исследованию устойчивости движения в том случае, когда правые части си­ стемы дифференциальных уравнений заданы лишь на конечном интервале вре­ мени длительности Т.

Это исследование можно проводить путем замены конечного интервала бес­ конечным в другом отсчете времени.

Действительно, пусть правые части системы (1.4) заданы лишь при 0

t<C_ Т.

Положим t =

<р (т), где ф (т) — такая непрерывно дифференцируемая

функция,

что при изменении аргумента т от 0 до + о о ф (т) монотонно изменяется

от 0 до Г.

 

Тогда система (1.4)

примет вид

 

= ф '

(т )& [ф (т ), Уи ■■-,Уп] (s = 1, • • -,«)•

Отметим, что для подробного ознакомления с теорией второго метода Ляпу­ нова необходимо использовать работы [42, 46].

4. Исследование вопроса об устойчивости нулевого решения по отношению к отдельным заданным координатам

При решении технических вопросов весьма часто ставится за­ дача исследования переходных процессов (неустановившихся движений) по некоторым заданным координатам.

Пример, приведенный в [43], показывает, что весьма важно знать общие методы исследования переходных процессов по от­ дельным заданным координатам. Для более определенной поста­ новки задачи введем ряд определений.

 

Определение

11

 

Нулевое решение системы (1.4) назовем устойчивым по коор­

динатам у ъ . . ., ук\

k <

п если для

каждого е > 0

можно ука­

зать такие два числа

8Х>

0 и 62 >

0 (б1< е ) , что при

£ 4 & > < б ? и S

0? о < б 2

(4.1)

;=1

 

(=*+1

 

47

будет

k

 

 

 

 

(4.2)

S

 

У! (^» УMi • • •>

<o) ^

^ 3^ ^0

1=1

 

 

 

 

для любого t0 ^

0.

 

 

 

Определение 12

Нулевое решение системы (1.4) называется асимптотически устойчивым по координатам г/1; . . ук, если по любому е > 0 можно указать такие числа 6Х> 0 и 62 > 0, что при выполнении (4.1) имеет место (4.2) и, кроме того,

к

 

Е У1 (t, ут, ■Упо, to) —0 при t~* -f оо.

(4.3)

i=i

 

Если система (1.4) имеет свойство, противоположное свойству устойчивости по координатам у ъ . . ., yk, то она называется не­ устойчивой по этим координатам.

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

dx

 

=*(*>*;

dy

<4 -4 )

где g (t), f (t) и h (t) — функции, заданные при 1 ^ 0 и кусочно­

непрерывные. t

Если J g (т) dx ограничен сверху при t ^ t0, то нулевое ре-

шение системы (4.4) будет устойчивым по координате х. Дей­ ствительно, система (4.4) имеет решение:

 

t

 

 

J g (T) dx

 

x — x0e

;

t

t

(4.5)

J

j h (x)dx

Ig (в)—A(0)]rf0

У — e io

/ (t) e*»

x0dx

Из вида этого решения следует наше утверждение, ибо \х \ <С < k |х„|, где

t

J g (Т)dx

48

i

Если при этом в системе (4.4) j h (т) dr не ограничен сверху,

^0

то нулевое решение системы (4.4) не будет устойчивым по обеим

координатам х и у

(в смысле Ляпунова), так как система имеет

 

 

 

t

 

 

 

 

 

 

\ h (t) dx

 

 

 

решение х

=

0, у —

eto

у0, которое принимает значения, ле­

жащие сколь

угодно далеко

от

положения

равновесия.

 

t

 

 

 

 

 

Если j

g (т) dr —>----со при

t ->- -f °°,

то нулевое решение

tо

системы (4.4) будет асимптотически устойчивым по координате х, в то время как устойчивость по обеим координатам в смысле Ляпу­ нова может и не существовать.

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

Факт неустойчивости по заданным интересующим нас коорди­ натам также не может служить критерием непригодности, ибо малым отклонениям в случае неустойчивости на конечном проме­ жутке Т могут соответствовать малые отклонения переходного процесса от исследуемого невозмущенного движения. При этом считается, что переходные процессы, возникающие в системе, инте­ ресуют нас при t г? Т.

Этот предварительный анализ определений устойчивости и неустойчивости по заданным координатам приводит снова к не­ обходимости введения понятия технической устойчивости. Это понятие было сформулировано в п. 1 и благодаря его общности может быть использовано при изучении устойчивости по заданным координатам.

Приведем условия устойчивости и асимптотической устойчи­ вости по заданным координатам у и у^ нулевого решения системы (1.4). Пусть из каких-либо физических или механических соображений известны оценки (может быть, даже сколь угодно грубые) поведения переходных процессов в системе (1.4) по по­ следним п— k координатам

\уi

(^> У1о» • ■•> Уп01

^о)| ^ Ф (^> Уш • • •> Уп о> Q

{t),

 

i =

k -f—1» • • •)

(4.6)

4 В.

и. Зубов

 

49

Соседние файлы в папке книги из ГПНТБ