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

книги / Матричное исчисление с приложениями в теории динамических систем

..pdf
Скачиваний:
7
Добавлен:
19.11.2023
Размер:
36.1 Mб
Скачать

составленная в силу системы (18.6.30), является отрицательно­ определенной эрмитовой формой, т.е. ((А*Н + НЛ)х, х) < 0

(11*11 £ 0 ) .

Составим производную по t формы (18.6.31) в силу системы

(18.6.27). Имеем

 

V(t, х) = (Вх, х) + (Их, ф) + (#<р, х),

(18.6.32)

где постоянная матрица В = А*Н 4- НА — отрицательно - опреде­ ленная. Из (18.6.28) следует неравенство ||<p(f, х)|| <е||х|| для до­ статочно малых ||х||:||х|| < /£ < г, где е — произвольно мало. Поэ­ тому из (18.6.32) имеем неравенство

vu . *)< [ Х„„(В) + 2е||Я||]||л||2,

(8.6.33)

откуда следует, что V(t, х) < 0 для всех t > 0, если только выпол­ няются неравенства

0

Ц//11

0 < | | х | | < / в.

 

Из (18.6.32) имеем еще V(t, х)

х=0 ~

Таким образом, для системы (18.6.27) в некоторой окрестности нуля существует положительно-определенная эрмитова форма К(х), явно не зависящая от t и допускающая отрицательно-опреде­ ленную производную по t в силу уравнения невозмущенного дви­ жения (18.6.27). Значит, невозмущенное движение (т.е. тривиаль­ ное решение системы (18.6.27)) устойчиво на [0, «).

Теперь покажем,

что, если ||х(0)|| ^ /е < л, то при выполнении

неравенства И(х(0))

р2 выполняется равенство

 

 

Ига ||х(0Н = 0 ,

(18.6.34)

т.е. невозмущенное движение (тривиальное решение системы (18.6.27) ) асимптотически устойчиво. В самом деле, в силу нера­ венства (18.6.33) функция V(x(t)) — монотонно убывающая, и по­ этому она при t-+ со будет монотонно стремится к некоторому пределу а, оставаясь все время больше этого предела, так что для всех /> 0

К(х(/)) > а.

(18.6.35)

Докажем, что а = 0. Пусть а =/= 0, следовательно, а > 0. Так как Г(х(/)) есть функция непрерывная, то из (18.6.35) вытекает, чго нетривиальное решение x(f) уравнения невозмущенного движения (18.6.27) удовлетворяет неравенству

В силу (18.6.36) получим V(t, x(t)) < —b < 0 для всех t > 0. Следо­ вательно, при всех / > 0 будет выполняться неравенство V(x(t)) ^ < К(х(0)) —Ы, что, очевидно, невозможно.

Таким образом,

мы приходим к заключению, что

lim K(x(f)) = 0, откуда,

вследствие знакоопределенности F(x(/)),

вытекает (18.6.34), что и доказывает теорему. ■ Т е о р е м а 18.6.2.2. Если хотя бы одно собственное значение

матрицы А обладает положительной вещественной частью, то невозмущенное движение (т.е. тривиальное решение уравнения (18.6.27)) неустойчиво.

Д о к а з а т е л ь с т в о . Система (18.6.27) не может быть устой­ чивой, так как при ограниченном ш(/) из устойчивости в смысле определения 15.1.1 следовала бы устойчивость по Ляпунову, а при условии теоремы — невозмущенный процесс (т.е. тривиальное ре­ шение уравнения (18.6.27) неустойчиво по Ляпунову). Теорема до­ казана. ■

18.6.3. Об устойчивости нелинейного процесса по первому приближению. Ниже устанавливаются некоторые условия устой­ чивости движения, представленного тривиальным решением х = 0 векторно-матричного уравнения

§ = A ( t ) x + h(t,x),

(18.6.37)

где A(t) — л х «-матрица на [а, »), непрерывная и ограниченная (11-4(01!) < и0 < оо (а0> 0 — некоторое число) и такая, что каждое решение системы

(18.6.38)

ft "

ограничено на [а, »); Л(г, х) — л х 1-матрица, элементы которой (нелинейные функции отклонений xs) таковы, что равномерно по i на [а, оо) выполняется соотношение

 

(18.6.39)

Т е о р е м а 18.6.3.1. Если на [/0, оо)

 

t

 

|ехр £ Sp A{x)dx\ 2* mt > 0,

(18.6.40)

t

где jij(£),..., (<) — все собственные значения эрмитово-симмет-

ризованной матрицы AN(t) = ^ (Л*(*) + Л(01» т0 невозмущенное

движение (тривиальное решение уравнения (18.6.37)) устойчиво.

Д о к а з а т е л ь с т в о . Пусть преобразование x = K(t)y приво­ дит уравнение (18.6.38) к диагональному виду

4 £ = Л (/)у.

(18.6.42)

Тогда

 

K(t) = X(t)CZ(t)

(18.6.43)

(см. гл. 12), причем можно, используя имеющуюся произвольность

выбора, сделать так, что столбцы Kg(t) (cr= 1, 2, .... п),

матрицы

K(t) удовлетворяли бы условию

 

||/Со(0Н = о на [f0, оо) (а = 1, 2, ..., л),

(18.6.44)

а матрица А(0 имела бы вид Л(t) =diag {X,, ..., Хп}, причем

К = jjin \\X c J .

(18.6.45)

Будем предполагать, что в (18.6.43) матрицы X{i)y С, Z(t) по­ строены таким же образом, как и при доказательстве теоремы 18.6.1.5. Из (18.6,44) следует, что K(t) € Кд. Следовательно, для

доказательства теоремы достаточно показать, что все решения x(t) Ф 0 уравнения (18.6.37), удовлетворяющие условию

(*~Ч*о)*(*о)> *~Ч*0)*(*о)) ^Р 2

( л < / < “ ),

(18.6.46)

для всех t > t0 удовлетворяют условию

(K~l(t)x(t), K~l(t)x(t)) К р\

(18.6.47)

Преобразование х = K(t)y приводит уравнение (18.6.37) к виду

% = A(t)y + M(t)h(t,Ky)

(М = К~1). (18.6.48)

Производная по t положительно-определенной эрмитовой формы

V(U х) = (K-'(t)x, K-'(l)x) = 1Ы12,

вычисленная в силу уравнения (18.6.48), представляется в виде

^7 = 2 2 *о1>и2 + 21*е (у*МЛ),

(18.6.49)

О■»1

 

где уа (cr = 1, 2, п) — элементы столбцовой матрицы y(t)> Интег­ рируя (18,6.49) вдоль решения возмущенного движения, получим

V(ty х) = К(*0,дс0) [ 1 + 2 (exp $ 2 Xad x -

1) 1уо>|

+

 

0=1

 

 

Ну(/п>11:

 

 

 

 

 

 

 

 

+

(< -<0Ж<. У)]>

(18.6.50)

где

 

/

t

 

 

 

 

 

 

4>0, У) = (t- tQ)\\y(tQf

J Re [у* exp

J 2A(x)dx

Mh]dt'.

i

 

 

 

 

Покажем, что на [?0, oo)

равномерно no t

 

 

Iim ip(f, у) =

0.

 

 

(18.6.51)

y-*Q

 

 

 

 

Действительно,

t

П

 

llyll2

m

V(t> у) < j —г $

2 exP J 4 V *

| | м | | а д \ Г > " '

*0 \a=l

 

 

 

 

С другой стороны

 

t n

 

 

1Ы12

 

 

 

exp

$ 2 2

+ $ 2|M*|| $ j| Л -

ll^oll2

 

<0o=l

 

 

 

 

 

 

В силу (18.6.44) имеем на [/0, «»)

 

 

 

 

||АГ(ОН='/ла).

 

Для матрицы dK(t)/dt имеем

 

 

 

II ^

II = Л о ( М || + ИЛИ).

 

Из (18.6.45) имеем на [tQ, »)

 

 

M'minO) < K W < ^шах(0

= 1, 2, ..., rt),

(18.6.52)

(18.6.53)

(18.6.54)

где fimin(0 и цшах0 ) — соответственно минимальное и максималь­ ное собственные значения матрицы AH(t). Но, так как на [f0, <»)

ц , ( 0 < 1 М я ( 0 1 К 1 М ( 0 И

0 = 1 , 2 , . . . ,

Л),

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

||Л 0)|| «£ Vna0.

Поэтому в силу (18.6.54) имеем на [f0, »)

 

dK(t)

п)ша0 < ».

(18.6.55)

II dt 1 (Уп +

 

На основании (18.6.40) имеем на [а, »): |det Х(*)| > /и, > 0. Для постоянной невырожденной матрицы С имеем: | det С| 5* т 2 > 0, а так как каждое решение уравнения первого приближения (18.6.38)

ограничено на [а, <»), то | det Z(t) \ ^ т ъ > 0 на ,

оо). Значит, на

[а, оо)

 

| del K(t) | > m > 0.

(18.6.56)

Из (18.6.53), (18.6.55) и (18.6.56) следует, что в (18.6.43) К есть матрица Ляпунова. Но так как M(t) = К~1также есть матри­ ца Ляпунова, значит, она ограничена по норме на [а, »). Кроме того,

Ит Л(/, К » /||у || = 0

у-0

в силу (18.6.39) равномерно по t на [а, оо). В этих условиях ||у||/||у0|| — ограниченная величина, и так как ограничены и все другие множители подынтегрального выражения в (18.6.52), то

 

$

2

exp J 4 \ adx

ID

 

 

а=1

где г — некоторая положительная постоянная.

Итак,

|-ф(/, у)| < г sup ||Л||/||у||,

откуда следует (18.6.51). На

основании

(18.6.41) и (18.6.54) следует, что на [а, °°)

 

 

 

г

 

 

- Ц -

и „ ( т ) Л « - / < 0

(0 = 1 , 2 ...„ в )

 

г —Гп

J

 

и существует такое 6 > 0, что

 

*

'

\yef

2

exp J 2 Xadx — 1

 

О= 1V

/

В Д 5 "

С другой стороны, учитывая (18.6.51), можно указать такое р> 0,

что для

всех

y(t), удовлетворяющих неравенству V{t, х(0) <

V(t0i x(t0)),

будем

иметь

|г|>(*, y(i) \ < 26. Тогда (см. [73])

V(t, x(t)) < V(tQ, х ( / 0) )

на [f0 ,

оо), а это означает, что любое реше­

ние x(t)

уравнения возмущенного движения (18.6.37), удовлетворя­

ющее условию (18.6.46), где 0 < р < р0 для всех / > t0, удовлетворяет неравенству (18.6.47). Теорема доказана. ■

^«.УСТОЙЧИВОСТЬ НА КОНЕЧНОМ ИНТЕРВАЛЕ ВРЕМЕНИ И ТЕХНИЧЕСКАЯ УСТОЙЧИВОСТЬ

В настоящей главе будет показано, что теория устойчивости на конечном интервале времени по Каменкову (и в других постанов­ ках), а также теория технической устойчивости могут быть получе­ ны в рамках теории /^-устойчивости.

§ 19.1. Ад-устойчивость на конечном интервале времени

Общее определение понятия Ад-устойчивости, приведенное в гл. 15, представляется следующим образом:

О п р е д е л е н и е 19.1.1. Если в заданном классе Ад существует такая матрица. G(t), совпадающая в момент t = t0 с заданной по­ стоянной матрицей G0 класса А£, что при достаточно малом р > О любое возмущение x(t), начальное значение х = x(tQ) которого удовлетворяет условию

(GQ1X0,GQ1X0) ^ P2,

(19.1.1)

на интервале Д = (t0, Т)

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

 

(G

-l(t)x,G -l( t ) x ) ^ p 2,

(19.1.2)

то невозмущенное движение устойчиво; в противном случае — не­ устойчиво. ■

Класс

Ад вполне

определяется заданным интервалом

А = [/0, Т)

и заданной

положительной функцией со(<). Рассмот­

рим класс Ад-постоянных матриц следующим образом.

Положим

сo(t) = со = const,

А = [*0, Т), приняв Т = tQ+ x (т> 0

— неко­

торое конечное

положительное число).

При таком выборе со и

Т матрица класса Ад — постоянная

матрица, и определение

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

[*0, 10 + т)

при­

нимает следующий

вид:

 

 

О п р е д е л е н и е

19.1.2. Если в заданном

классе

где

Л = [<0, /0 + т), существует такая матрица G, что при достаточно ма­ лом р > 0 любое возмущение x(t), начальное значение х0 = x(t0) ко­ торого удовлетворяет условию (C7_I JC0, С7_1л0) ^ р2, на конечном ин­ тервале [/0, Т) удовлетворяет условию

(<Г1х(0,СГ1х ( 0 ) « | *

то невозмущенное движение устойчиво на конечном промежутке времени [/0, tQ+ т). В противном случае — неустойчиво. ■

Сравнивая последнее определение понятия устойчивости с опре­ делением понятия устойчивости по Каменкову (см. § 14.6), легко видеть, что они эквивалентны и отличаются лишь обозначениями. Роль постоянной матрицы А в определении Каменкова здесь играет

постоянная матрица G~l.

Продолжим наши рассуждения. Применим общую теорему о Ад-устойчивости (см. теорему 16.2.1) к понятию устойчивости со­

гласно определению

19.1.2 в условиях теоремы Каменкова об ус­

тойчивости (теорема

14.6.1).

об у с т о й ч и в о с т и 19.1.1. Если

Т е о р е м а К а м е н к о в а

матрица Р0 = P(t0),

не имея

кратных собственных значений,

имеет только отрицательные собственные значения или комп­ лексные с отрицательными вещественными частями, то невоз­ мущенное движение обладает устойчивостью на конечном интер­ вале времени [/0, t0 + т).

Аналогичная теорема о Ад-устойчивости формулируется так:

Т е о р е м а

о Ад-устойчивости

19.1.2.

Если существует

положительно-определенная эрмитова

форма

V(t, х) = x*A(t)x

такая, что

 

 

 

Г.

A(i0) =

(G j')-c ;1;

 

 

Т .

-п Sp /Г 1 < ш2;

 

 

3е. ^ <0 при V /е [t0, *0 + т),

то система

% = Р 0х + ДP(t)x + h(l, х)

устойчива.

Пусть А — матрица, преобразующая матрицу системы Р0 == = P(t0) к каноническому виду Л. Предположим, что столбцы мат­

рицы К пронормированы так, что ||Яу.|| = <о (/ = 1, 2

, л); будем

иметь К Е Яд, так что G = K . Примем V(t, JC) =

(K~lx, K~lx).

Матрица этой формы

 

А = ( К ~ 1УК-1= (C7-l)*G_1.

 

Так как в данном случаем G = К = const и согласно определению понятия устойчивости имеем G0 = G, условие 1° теоремы о Яд-ус- тойчивости выполняется:

A(t0) = (К -'УК “ = (G0-')-G -'.

Далее согласно теореме 16.1.1 о скелетном разложении положи­ тельно-определенная эрмитова матрица А может быть представлена

в виде А~1= Я Я \ где Я — квадратная матрица, все столбцы ко­

торой имеют одну и ту же эрмитову норму or =

Sp A~l .

Из выражения матрицы А имеем

 

А = (G~l)*G~l = (G‘)~xG~i - GG\

так что, очевидно, Я = G = А, и значит, в данном случае ст = со, поэтому условие 2° теоремы также выполняется.

Остается проверить условие 3°. Очевидно, имеем К(г, х) = (К~'х, JT'x) = ||уН2.

Надо из уравнения возмущенного движения определить d\\y\\/dt, и тогда при условиях теоремы Каменкова будет ясно, что

% =Чу\\^ « о .

Тем самым будет показана эквивалентность теоремы Каменкова об устойчивости движения на конечном интервале времени [<0, t0 + T) и общей теоремы Яд-устойчивости в условиях теоремы

Каменкова. Аналогично можно установить эквивалентность и дру­ гих теорем Каменкова и общей теоремы о Яд-устойчивости.

§19.2. Ад-устойчивость и техническая устойчивость

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

«г Ч О * . с - ‘ ( / ) * ) « р ? ,

§ 19.2. ^-УСТОЙЧИВОСТЬ И ТЕХНИЧЕСКАЯ УСТОЙЧИВОСТЬ

399

посредством матрицы G{t) класса К%. Полагая

% = {•V (G"'(!o)*0G->(t0)x0) « р2},

ш(() = {x:(G-'(r)x. G-‘(r)x) Sp2}

определение понятия /^-устойчивости можно представить так.

О п р е д е л е н и е 19.2.1. Невозмущенное движение (тривиаль­ ное решение уравнений возмущенного движения) называется К%- устойчивым на 1 — [/0, /0 + т) относительно заданных областей

о)0 и со в том и только в том случае, если всякое возмущенное дви­ жение уравнений, определенное начальным условием

x(t0) Е u)0 = {*0:(G-i(f0)x0. G~l(lo)xo) < Р2>>

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

х(0 Е со(И) = {x:(G~l(t)x(t),G~l(/)x(t)) « S p 2} при V/ Е /. ■

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

ГЛАВА 20

УПРАВЛЯЕМЫЕ СИСТЕМЫ

§ 20.1. Понятие об управляемых системах

Управление — это такой процесс, который охватывает чрезвы­ чайно обширные области деятельности человека. По сути дела, прак­ тически любое сознательное действие человека означает некоторое воздействие с его стороны на ход событий с той или иной целью. В данном параграфе мы коснемся вопросов управления лишь в очень небольшой части этой обширной области, а именно — будем рассмат­ ривать управление такими объектами, которые при допустимых иде­ ализациях можно описать системами линейных обыкновенных диф­ ференциальных или интегро-дифференциальных уравнений.

Состояние объекта в каждый момент времени обычно можно охарактеризовать посредством некоторого числа (конечного или бесконечного) параметров. Так, поступательное движение твердого тела по прямой можно задать некоторой линейной координатой точки на прямой, по которой движется тело, и скоростью его дви­ жения. В некоторых задачах этими двумя параметрами можно ог­ раничится, например, при рассмотрении движения автомобиля (или другого подобного технического объекта) по прямой. Для опи­ сания движения материальной точки на плоскости требуется уже четыре параметра, в трехмерном пространстве — шесть.

Параметры, которые характеризуют состояние объекта, при­ нято называть фазовыми координатами. В приведенных приме­

рах

роль фазовых координат выполняют линейная координата

s и

скорость v = s. Фазовые координаты могут иметь самую

разную природу (геометрические, кинематические, динамиче­ ские, термодинамические, химические и т.д.).

Объект является управляемым, если имеется возможность путем вмешательства «извне» влиять на изменение фазовых

координат х1 (/'= 1, 2 ,..., л). Например, на величину пройден­ ного пути и скорость автомобиля можно воздействовать путем изменения положения педали газа и педали тормоза. Его дви­ жением по дороге можно управлять еще путем поворота руле­ вого колеса. Те величины, посредством которых осуществляется

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