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

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

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

Г л а в а ///

ОЦЕНКА ПОВЕДЕНИЯ ПЕРЕХОДНЫХ ПРОЦЕССОВ В НЕЛИНЕЙНЫХ И УПРАВЛЯЕМЫХ СИСТЕМАХ

12. Метод оценок. Техническая устойчивость

Пусть заданы интервал изменения времени U0, t0 + Т] и два множества Gt h G( в пространстве п измерений yv . . ., уп. Если

существует такое множество F в пространстве допустимых зна­ чений параметров, что интегральные кривые системы (1.4), на­ чинаясь при t = 10 в множестве G{ , остаются в множестве G

при / £ П0, t0 -f Т], то нулевое решение системы (1.4) называется технически устойчивым.

Как уже отмечалось в п. 1, множество G( определяется, как

правило, техническими условиями работы реальной системы. Мно­ жество Gt формируется из определенных требований, предъявляе­ мых к расчетной системе. В большинстве случаев можно считать, что множества Gt и G( являются сферами с радиусами г( и г(.

Тогда задача технической устойчивости сводится к нахождению оценок поведения решений системы (1.4) по всем координатам; при этом координаты считаются равноправными в отношении предъявляемых к ним требований.

В п. 4 были приведены примеры, в которых было важно про­ вести детальное исследование поведения решений по некоторым заданным координатам (устойчивости по заданным координатам). В большинстве случаев такого рода в качестве множеств Gf и G(

можно задать «-мерные эллипсоиды.

Оказывается возможным свести решение вопроса о технической устойчивости к нахождению таких оценок на некоторой функции R

(yv . . .,

уп, t),

что из /?, (t, у10, . . ., уПо,

t0) ^ f ( t ,

y1Q, . . ., уПо, /0)

следует,

что

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

кривая ys

= ys (t, у10, . . ., уп0,

tо)

(s = 1 , . . . , « )

системы (1.4), начинающаяся во

множестве

Gt0,

попадает во множество Gt, t 6

[t0, i0 +

Т]. Через

(t, y 10, . . .

• ■•. Уno уtо) обозначено значение функции R (уъ . . ., уп, t) на этой интегральной кривой.

150

Ясно, что функция R (уъ . . ., уп, t) определяется в некотором смысле геометрическими свойствами множеств Gt и Gt.

и

Предположим, что существуют две функции

V (уъ

. . . ,

уп, t)

W (уI, . .

уп,

t),

удовлетворяющие

 

следующим

условиям:

 

 

 

 

V

и W заданы

 

 

t £ [tQ, t0 + Т]

 

П

у)

 

 

1)

функции

 

при

и £

< а 2;

где

а

и Т — положительные

постоянные;

 

 

 

 

 

 

 

2)

функция

V

{уъ

. . .,

уп, t) удовлетворяет неравенствам,

 

Фх (О R 1' (*> Уъ

 

Уп) <

V «S ср2 (0

 

tf'2 (f, г/ь

. .

уп),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( 12. 1)

где

 

(t) и ср2 (0 — кусочно-непрерывные положительные функ­

ции, заданные в промежутке

[t0,

t0 +

T J ;

функция R

(t, y x, . . .

 

 

 

 

 

 

 

 

 

 

 

П

 

 

 

 

 

 

 

 

. . ., yn) задана при t £

\tQ, tQ+

ТЛ,

£

y]

О 2,

a

R

 

( t, yv . . .

 

 

 

 

ПРИ У i

L

 

 

 

J

1=1

 

и

 

 

V

 

поло­

• • •, Уп) =

0 —

=

• =

 

Уп

=

0 ;

h

/ 2 —

 

такие

жительные

числа,

что

 

^

/2;

 

 

 

 

 

 

 

 

 

 

 

 

 

3)

функция

W (у1,

. . ., уп, i) удовлетворяет

неравенствам

 

 

 

 

-------

y n ) ^ W ^ — %(t)Rk2(t, У1..........Уп),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( 12.2)

где Фх (0 и ф2 (t) — кусочно-непрерывные функции в промежутке

[^о, t0 + T 1]-, &i и к 2— такие положительные числа,

что

функции Е

 

4)

функция

W является

полной

 

производной

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

(1.4):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т

=

 

 

 

 

 

 

 

 

02.3)

Теорема 32

Если существуют две функции Е и W, удовлетворяющие предыдущим четырем условиям, функции ф,- (t) положительны

и

ч

> 1 (i

= 1, 2),

то

любая интегральная кривая си-

 

 

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

стемы (1.4)

 

 

 

 

 

I—А-1

 

I

 

 

 

ф2 1(0 Ео

1 +

(Ki — 1) Ео1” 1 J ф! (т) фТх' (т)£Йг

 

 

 

 

 

: Rl (t,

Ую>

• •> УПО) А>)

при

 

t £ [tQ,

10 +

7 \],

 

(12.4)

 

 

 

где

 

Eq = V

( ух о,

. , -| у по*

to)>

 

151

Доказательство. Возьмем интегральную кривую системы (1.4) и подставим ее в равенство (12.3); тогда оно превратится в тож­

дество = W. Умножая обе части этого тождества на функ­

цию V~%, где X >

1, и интегрируя в пределах

от t0 до t, получаем

 

t

 

t

 

 

 

j

dV =

j WV~Kdt,

 

 

 

to

 

to

 

 

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

 

 

 

 

 

 

 

 

 

(12.5)

Производя тождественные преобразования в равенстве (12.5),

находим

 

 

 

 

 

У =

,

 

t

(12.6)

 

я-1 г

 

v

'

 

J/

1+ (1

j WV~x dt

 

 

 

 

tо

 

 

Оценим подынтегральное выражение в формуле (12.6). По­ ложим X = Х2.

Функция W удовлетворяет неравенствам

— Фг (t) R kl (t, у ъ

. . .,

у „ )<

<

— ф2 (0 R k* (t,

у ъ

. . ., уп).

Функция W~%2

удовлетворяет

неравенствам

4>Tu (t)R - ^h (t, уи . . . , Уп) ^ У ~ Хг^:

^<pTX2(t)R-x‘l'(t, уи . . ., Уп).

Из этих неравенств следует, что

 

WV~U ^

ф2 (0 ц>ТХг (t),

(12.7)

так как

Х2 = ■

 

 

Из

неравенства (12.7) и

соотношения (12.6)

получим при

X = Х2, что

 

 

V*

( 12.8)

 

V

(t) dt

152

Положив в (12.6) к = k t и

используя

неравенства

(12.1) и

(12.2), найдем

 

 

 

WV~K <

<pi (0

(/)•

(12.9)

Из неравенства (12.9) и соотношения (12.6) при к = кх по­

лучим

 

■- - ■

(12.Ю)

У И - ( ^ - 1 ) V'S1-1 j b V ) v r X l V)dt

 

to

Из неравенств (12.7) и (12.9) следует, что функция R (t, у и . . .

■■Уп) удовлетворяет неравенствам

^ ФГ1

yv . . ., Уп) ^ У

(t) V.

(12.11)

Подставив в правую часть неравенства (12.11) вместо функ­ ции V правую часть (12.8) и в левую часть (12.11) вместо функ­ ции V правую часть (12.10), получим равенство (12.4).

З а м е ч а н и е . Оценка (12.4) имеет место лишь вдоль тех интегральных кривых, которые лежат в области задания функций V и W.

Это замечание не носит характер ограничения на область применимости оценки (12.4), а лишь предполагает, что функции V и W продолжены по пере­ менным г/1; . . . . уп на все фазовое пространство.

Теорема 33

Если существуют две функции V и W, удовлетворяющие усло­ виям 1—4, приведенным выше, и kL > 1 (г = 1, 2), то при =

— /2 = I любая интегральная кривая удовлетворяет неравенствам вида

1-?.,

ф2 1(0 Ко ] /

1 +

-

1) Ко1” 1 } и (т)

:^ l( ^ >

J/lOl

• • •!

Упй1 *<>) =

j

Фт {t) ф1Я' (0.

если фх (t) ;== 0

1 1

Фь (0 ф2~Л‘ (*),

если

(t) < 0

 

|

Фг (0 Ф'Г*"2(0,

если ф2 (t) 5- 0

~

1

ф2 (t) срГЯг (0,

если

ф2 (0 < 0

153

Доказательство. Возьмем интегральную кривую системы (1.4) и подставим ее в равенство (12.3); тогда теми же преобразованиями, какие были сделаны при доказательстве предыдущей теоремы, получим при Я > 1 соотношение вида (12.6).

Произведем оценку величины — WV~X, стоящей под интегра­ лом в (12.6).

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

если ф2 (t) > 0 , то имеет место неравенство (12.7)

при Я2 = Я.

Если ф2 (О <С 0, то величина — WV~X удовлетворяет неравенству

W \ T X*S& ЬЧ>ТХг (t).

(12.12)

Таким образом, из неравенств (12.12) и (12.7) следует, что имеет место неравенство

 

 

WV~X2 ^

f 2 (t),

(12.13)

 

.

( ФзфГ*2 при 1р2 ^

0

/

'

\ г|)2ф1

при ф2 <

0

Из неравенства (12.13) и соотношения (12.6)

вытекает, что

 

V =

И,

 

(12.14)

 

 

 

 

] / "

l + ^ - l )

^ - 1 J / 2 (x) dx

Произведя аналогичные оценки сверху для величины — WV~% при Я = Я1( получим

V :

1 1 +

 

(12.15)

] /

q1(Я1)- 1Ef i"J (т)

 

 

 

dx

где

 

 

 

=

 

при -ф! ^

0

1

i ф1фГЯ‘

при ф1<

0

В предыдущей теореме было отмечено, что функция R удовле­ творяет неравенству (12.11). Используя это неравенство и нера­ венства (12.14) и (12.15), получим требуемые оценки.

З а м е ч а н и е . Оценки, приводимые здесь, справедливы для тех кривых, которые лежат в области задания функций V и W, и лишь при тех значениях ta< ПРИ которых подкоренные выражения, встречающиеся в формулах,

неотрицательны.

Теорема 33 является более сильной, чем предыдущая, так как функции фг могут принимать значения разных знаков. Это замечание было бы неполным, если бы мы не отметили также, что теорема 32 относится к более широкому классу функций V. Однако вид этих функций каждый раз зависит от способа, по кото­ рому мы их выбираем. От нашего выбора зависит также функция R, поэтому условие li = /2 не является слишком ограничительным.

154

Теорема 34

Если существуют две функции V и W, удовлетворяющие усло­ виям 1—4 и = Л2 = 1; ф,- ^ 0, то любое решение системы (1.4) будет удовлетворять неравенствам

—J t iT f 1 dx

 

 

Фа'ЧОУое

 

yw, ■■ Упо, to) ^

j

Ф2(|:12_1 dx

ф ! 1 ( 0 Кое

to

(12.16)

Доказательство. Пусть существуют две функции К и W, удо­ влетворяющие условиям теоремы 34. Возьмем интегральную кри­ вую системы (1.4) и подставим ее в соотношение (12.3). Умножим полученное тождество на функцию К-1 и проинтегрируем в пре­ делах от tо до t\ тогда получим

t t

J k'MK = J WV~4t,

to

to

 

 

откуда

t

 

 

 

 

 

V =

J WV~‘ dt

 

 

K0e *°

.

(12.17)

Отметим, что интеграл,

состоящий

в правой

части равен­

ства (12.17), можно оценить тем же способом, как это было сде­ лано при доказательстве теоремы 33. Тогда будем иметь

 

 

—ф1фГг^

IKK"”1^ —фгфГ1-

(12.18)

Из неравенств (12.18) и (12.11) получим неравенство (12.16).

З а м е ч а н и е .

Если ф1 и ф2

меняют знак,

то неравенства (12.18) имеют

место лишь при тех

t, при которых ф г (t) > 0, ф2 (t) >

0. Если же фг <С 0 и

ф2 <С 0, то функция

WV_1

будет

удовлетворять

неравенствам:

для ф2 <

0

 

 

 

 

 

 

 

ДЛЯ фх <

О

 

 

'tv'- 1^

зфг1;

 

 

 

 

 

 

 

 

 

 

 

 

 

w v - 1^

—Ф ^г1.

 

 

Из этих

неравенств

и

неравенств (12.11) получим

 

f

 

 

 

—J Ф!фГ1 dx

 

 

ФУ1(0

Т0е и

 

ic Rj (i!, у10,-

уп0, <0) ^

 

 

 

 

 

 

 

 

_

 

 

 

 

 

 

 

-

—(

dx

 

 

 

 

 

Я ъ 1 {t)V& >

.

(12.19)

155

Здесь

 

<Pi =

<*i<Pi (0 +

о2ф2 (О;

 

?2 =

ОГзф! (О +

04ф2 (О,

где

 

 

 

 

 

0 1 =

I + sign .

 

 

— sign .

2

 

 

 

2

 

 

 

 

 

I — sign \|J3 .

о,=

+ sign ф2

 

 

 

 

Теорема 35

Если существуют две функции V и W, удовлетворяющие усло­ виям I—4 и = I, А.2 > I и lx — Z2 = I, то любое решение си­ стемы уравнений (1.4) будет удовлетворять неравенствам

где

J Ф2Ф2 при фг О

I фгфГ^2 при фг < О

Доказательство. Левая часть неравенства (12.20) может быть найдена из замечания к теореме 34.

Правая часть этого неравенства может быть получена тем способом, который применялся при доказательстве теоремы 33.

З а м е ч а н и е . Аналогичная теорема при

> 1, К2 = 1, = /2 = /

также имеет место. В этом случае любое решение системы (1.4) будет удовлетво­

рять неравенствам

 

i Г

 

 

Г

I

 

1 /

Vj-'WVo

у

1+ (^1 -

1) Vj»-1 J /, (т) dx;

 

 

 

 

 

-J 1|>2Ф2 dx

; R\ (t,

У10. • ■ • . УпО>

h) « = \

4 4 1 ( 0 Е 0е

( 12.21)

где

 

 

 

 

 

 

|

1М 0 ф1 Х‘ (0 при % ( 0 ^ 0

 

1 Ф 1(0 фГЯ‘ (0

при Ф г(0 <

0

156

Отметим, что неравенства, о которых говорится во всех теоремах, сформули­ рованных в этом параграфе, имеют место лишь при тех значениях t0, пока

интегральные кривые системы (1.4) лежат в области, в которой существуют функ­ ции V и W, удовлетворяющие условиям (12.1)—(12.3). Кроме того, эти неравен­

ства имеют смысл лишь до тех пор, пока подкоренные выражения, встречающиеся в них, неотрицательны.

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

Теорема 36

Пусть даны: число Г > 0 и два множества Gto и Gt.

Если существуют две функции V и W, удовлетворяющие усло­ виям 1—4 и такие, что одно из неравенств (12.4), (12.4а), (12.16)

(12.20), (12.21) обладает тем свойством, что если

(у10, . . ., уп0) £ х

X Gtо, то (ух, . . , уп) 6 Gt при tG U„, tо +

Т], тогда нулевое

решение системы (1.4) является технически устойчивым.

Эта теорема будет давать, вообще говоря, лишь достаточные условия технической устойчивости; при этом неудачный выбор функций R и V может повлечь к жестким достаточным условиям. Для окончательного выяснения этих условий необходимо перебрать несколько экземпляров функций R и V и после этого продол­ жить рассмотрение полученных достаточных условий по возмож­ ности путем применения электронной модели системы (1.4).

Ниже приводятся примеры использования приведенного здесь метода.

13.Устойчивость на конечном интервале

Вработах Г. В. Каменкова [17] и А. А. Лебедева [23] была раз­ вита теория устойчивости нулевого решения на конечном интер­ вале времени. Эта теория имеет весьма тесную связь с понятием технической устойчивости. В настоящем параграфе мы изложим основные результаты, относящиеся к устойчивости на конечном интервале.

Впредыдущих параграфах отмечалось, что значения некоторой

функции V на интегральных кривых системы дифференциальных уравнений (1.4) позволяют иногда судить о поведении этих кривых. Эта же идея лежит в основе теории устойчивости на конечном ин­ тервале.

 

 

 

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

 

 

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

данном 10

по отношению к положительно-определенной функ­

ции V (ух,

. . .,

уп, t)

на интервале т, если

из V (у10, ■

</ло>

tо) = а следует,

что

{t, y 1Q, ■. уп0, t0) < а „при t0 <

* «=s

«s; t0 + т

и любом выборе достаточно малой

положительной

ве­

личины а.

157

Здесь, как и выше, Vx (t, у 10, . . ., уп0, /„) означает величину функции V (уи . . ., уп, t) на интегральной кривой системы (1.4), проходящей при t — t0 через точку (у10, . . ., уп0).

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

Нулевое решение системы (1.4) называется устойчивым на ко­ нечном интервале при данном t0, если существует положительное число т > 0 и такая положительно-определенная функция V (Ух, . . ., уп, t), что по отношению к ней нулевое решение си­ стемы (1.4) устойчиво на конечном интервале т.

Рассмотрим систему дифференциальных уравнений

%

= S Psi (0 Vi + /s (t, уъ .. ., yn){s = 1, .. ., n).

(13.1)

ul

(=i

 

Здесь psi (t) — вещественные функции, заданные при t ^ 0 и непрерывные там; функции fs (t, у и . . уп) заданы при t ^ 0

и ^

У21

“2 (а — положительная

постоянная, а>»0), непре­

рывны и

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

 

 

 

1 М < * ( 0

( £ * ) ' .

(13.2)

где

с

ф5(^) — непрерывные

положительные

функции, за­

данные при t ^ 0.

Систему (13.1) можно рассматривать как один из возможных вариантов системы (1.4).

Теорема 37

Для того чтобы нулевое решение системы (13.1) при данном t 0 было устойчиво на конечном интервале по отношению к некоторой фиксированной положительно-определенной квадратичной форме V при любом выборе функций fs, удовлетворяющих неравен­ ствам (13.2), необходимо и достаточно, чтобы корни уравнения

\ P (t0) — kE\ = 0

(13.3)

имели отрицательные действительные части, где Р (t0) — матрица, элементами которой являются числа psi (t0).

Прежде чем перейти к доказательству этой теоремы, заметим, что нулевое решение линейной системы

= S P s i ( t 0) y i {s = \ , . . . , n )

(13.4)

спостоянными коэффициентами асимптотически устойчиво тогда

итолько тогда, когда все корни уравнения (13.3) имеют отрица­ тельные вещественные части.

458

Известно, что в случаё асимптотической устойчивости любое решение удовлетворяет неравенству

i

y U c x £ J & e-',

/=i

i=i

где clt c2— положительные постоянные. А тогда, как было по­ казано в п. 9, существуют две квадратичные формы V и W, удовле­ творяющие следующим неравенствам:

П

П

ai 2 y l ^ V ^ c i i 2 У)\

i=i

;=i

~ b xf 1 y l ^ W ^ - b ^ у],

t=i

i=i

где at, bt (i = 1,2) — положительные

постоянные.

Отметим, что эти формы можно выбрать таким образом, чтобы

их коэффициенты были постоянными

числами.

Действительно, пусть нулевое решение системы (13.4) асимпто­ тически устойчиво. Возьмем произвольную определенно-отрица­

тельную квадратичную форму W.

Следуя тому

способу,

который

был изложен

при доказательстве

теоремы 27,

положим

 

 

 

 

 

+

CD

 

 

 

 

V =

j

W dt.

 

(13.5)

 

 

 

 

 

о

 

 

 

Легко видеть, что функция V является квадратичной формой

переменных у 10, . . ., уп0,

положительно-определенной

и не за-

,

при этом

dV

=

TV7

 

dV

 

 

висящей от г;

 

IV,

где-^- есть полная производная

формы V в силу системы (13.4). Отсюда следует, что это соотно­

шение имеет вид

 

 

 

 

 

 

 

 

T i W t H Plk{to)yk = W-

 

(13'6)

 

1=1

k=1

 

 

 

 

Форма V, определяемая (13.5), является единственным реше­ нием уравнения (13.6) в классе квадратичных форм. Таким обра­ зом, получено следующее утверждение.

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

Это утверждение впервые сформулировано А. М. Ляпуновым и доказано им при помощи теории производных определителей. Отметим, что весьма полезна алгебраическая форма этой теоремы Ляпунова.

159

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