книги из ГПНТБ / Абгарян К.А. Матричные и асимптотические методы в теории линейных систем
.pdf3 9 0 |
Н Е К О Т О Р Ы Е У С Л О В И Я |
У С Т О Й Ч И В О С Т И |
[ГЛ XV |
||
|
Допустим, что К |
(t) |
= (Ki |
К2 ■■■ К п) — невырожденная |
|
и дифференцируемая |
на |
U0, Т] матрица,, столбцы которой |
|||
имеют одинаковую норму, а именно: |
|
||||
и |
II Я/(01 = <*(/)>0 |
(/ = 1,2,...,«) |
|
||
|
|
|
|
|
|
|
dK |
= UK - |
КА + /ѵ, |
|
|
|
dt |
|
|||
где A = diag (^, Ä,2, ..., Kn), a N — некоторая квадратная матрица порядка п.
Замена переменных
je = Ку
приводит уравнение (6.1) к виду
^= A(t)y - M(t)N (t)y М M(t)h(t, Ку).
Полная производная от положительно определенной функции
V(t,x) = (K-'(t)x, K~'(()x) = lyF
по / в силу уравнения возмущенного процесса в данном случае представляется в виде
dV |
= 2 1 У ( 0 II2 Ф (*. У (0 ) + |
2 Re(у*МІі), |
(6 .1 1 ) |
||||
где |
|||||||
|
|
|
|
|
|
||
|
|
|
|
у*Ру |
|
||
4 > ( i , y ( t ) ) = £ R e ^ ™ - + |
|
|
|||||
|
|
0=1 |
11УllT ^ I F F |
|
|||
|
Р = |
— ~ |
(MN + N*M*). |
|
|
||
Интегрируя (6.11) в пределах от t0 до t, |
получаем |
||||||
V (t, х) — V(і0, х0) 11+ |
ехр { 2q>(t', y(t'))dt' — 1 |
+ |
|||||
|
|
|
+ |
R - W |
, у) |
. (6.12) |
|
где |
|
|
|
|
|
|
|
Ф (t, у) - |
|
' |
2J ф(і,у |
(т)) г/т |
|
|
|
|
■ Л » f |
Re (у*М!і) dt'. |
|||||
|
«~<о)ІЫР J |
|
|
|
|||
|
|
Іп |
|
|
|
|
|
s 6] |
ОБ У С Т О Й Ч И В О С Т И |
НА |
К О Н Е Ч Н О М |
П Р О М Е Ж У Т К Е |
391 |
|||
В соответствии с условием |
(6.2) |
|
|
|
||||
|
|
limij)(/, у) = 0. |
|
|
(6.13) |
|||
|
|
(/-►О |
|
|
|
|
|
|
Пусть по-прежнему |
|
|
|
|
|
|
||
|
|
Во (/) = |
таха (ReХа (t)), |
|
|
|||
a vmax (/) — максимальное |
собственное |
значение эрмито |
||||||
вой матрицы Р . |
6.3. Пусть на промежутке П0, Т) |
|
||||||
Теорема |
|
|||||||
и |
|
а |
( 0 |
< |
со(і) |
|
|
|
|
І |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Т ^ Т f [Во (П + |
vmax (Щ |
|
|
(6.14) |
|||
|
1 |
10у |
|
|
|
|
|
|
|
|
*п |
|
|
|
|
|
|
где b — положительное |
число. |
Тогда невозмущенный |
про |
|||||
цесс |
(тривиальное решение |
уравнения |
(6.1)) |
устойчив на |
||||
промежутке U0, |
Т). |
|
|
При условии |
(6.14) сущест |
|||
Доказ ательс тво. |
|
|||||||
вует |
такое 6 > |
О, что |
|
|
|
|
|
|
|
/ |
|
|
|
|
|
|
|
exp j 2ср (?, у (О) dt' — 1< — 26(t — g.
С другой стороны, принимая во внимание (6.13), можно указать такое р0 > 0, что для всех у, удовлетворяющих неравенству |у |< р0, будем иметь |я|з (t, у) |С 26, и тогда (см. (6.12)) V (/, х) С V (і0, х0). Значит, любое решение X = X (t) уравнения (6.1), которое удовлетворяет условию
( К №о)Х Ѵо), К (0) х (^о)) ^ Р2>
где р — произвольное положительное число из промежутка О< р <; ро, в пределах промежутка [t0, Т) будет удовле творять неравенству
(K -'(t)x(t), K ~ [(i)x (t))< р2.
Но тогда будет иметь место и неравенство
G-1(/)х (ОХР2,
где G(0 = К (О —Ш-, так как по условию теоремы а (і) с
< со(0- Условия устойчивости процесса выполняются, по скольку G(0 есть матрица класса К&- Теорема доказана.
392 |
Н Е К О Т О Р Ы Е У С Л О В И Я У С Т О Й Ч И В О С Т И |
[ГЛ . XV |
Следствие. Если на промежутке U„, Т)
а (0 < ш(t)
и
Р-0 (0 + v max (t) |
О, |
то невозмущенный процесс (тривиальное решение уравнения
(6.1)) устойчив на промежутке [/0, Т).
Аналогично теореме 6.2 легко устанавливается и Теорема 6.4. Если на промежутке [(0, (г)
а (О О (О
и
Ро (^о) Ч- v max (^o) < Z. О,
то существует конечный промежуток [/0, Т) cz [/0, /х), на
котором невозмущенный процесс (тривиальное решение урав
нения (6.1)) устойчив.
6.3. Применение алгоритма асимптотического преобра зования уравнения. Наряду с уравнением (6.1) введем в
рассмотрение уравнение
j £ - — U(x )x + h ( t , x ) , |
х = еі, |
(6.15) |
которое при е = 1совпадаете уравнением (6.1). Невырожденным преобразованием
x = |
K im)(T,B)Zy |
(6.16) |
уравнение (6.15) приведем к виду |
|
|
= (z -1A(m>Z — Z~l ~ y -J у — Z~lM{m)N{m)Zy + |
|
|
где M (m) = К {т)~\ |
+ Z~lMim)h(t, K m Zg), |
(6.17) |
|
|
|
N{m) (т, е) == е |
----UKm + К (т)А(т). |
|
Допустим, что U (т) на рассматриваемом промежутке
изменения аргумента является / раз дифференцируемой матрицей. Тогда, используя алгоритм, приведенный в
гл. VIII, можно построить такую матрицу А ' " ' 1 , ч т о матрица
Л(т) будет иметь диагональный или по крайней мере квази диагональный вид, а матрица Л1(т) будет удовлетворять
§ 6] ОБ У С Т О Й Ч И В О С Т И Н А К О Н Е Ч Н О М П Р О М Е Ж У Т К Е |
393 |
условию.
lim /Ѵ(т>(т, в) = 0 (т = О, 1,2.......... / — 1). (6.18)
Мы здесь ограничимся случаем, когда U имеет только простые собственные значения. При этом матрицы К {т) и
Л(ш’ представляются в форме конечных сумм:
К {т)(г, е) = |
т |
|
|
|
Л(ш) (т, е) = 2 |
(т), |
||||
2е'К Ш(т), |
||||||||||
|
|
|
/=0 |
|
|
|
|
|
і=С |
(6.19) |
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
К и] = (К \І]К Р[ |
... |
К п\1 |
л [г] = |
0 |
|
|||||
|
|
|||||||||
|
|
|
|
|
|
|
|
|
іГЧ |
I |
|
|
|
|
|
|
|
|
|
Л-Г7 |
|
причем |
/(„' — столбцовые |
матрицы, |
а К1а 1— скалярные |
|||||||
функции. |
|
|
(k = |
0, |
1, 2, ..., т) таковы, что |
|
||||
Если |
Кік' |
|
|
|
||||||
|
иЛ1*1 |
|
|
|
|
о |
(6.20) |
|||
|
|
|
|
UKm = |
K mAl0\ |
|
||||
UKm = |
K [k]Л[0) + |
/С[0]Л[к1 + |
|
|
(*=1, 2, |
т), |
||||
где |
|
|
|
|
|
|
|
|
|
( 6 . 2 1) |
|
|
|
Л—1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d K } 11 |
|
||
|
|
|
|
|
д - [ А — а ] д [ а ] |
|
||||
ТО |
|
|
|
а=1 |
|
|
|
|
|
|
/Ѵ(ш' = |
еш+' |
/я |
m |
ev-itf'»-‘*+'’]ACa]I |
|
|||||
|
2 |
2 |
|
|||||||
и при условии |
|
|
v=l |
а=ѵ |
/<'[,], |
Atn (/ = 1, 2, |
..., т) |
|||
ограниченности |
||||||||||
требование (6.18) соблюдается.
Равенство (6.20) удовлетворяется, если в качестве диаго
нальных элементов матрицы Лг°] взять собственные значе ния Я2, ..., Хп матрицы U, а в качестве столбцов матрицы
/<[0:і принять отвечающие этим собственным значениям соб
ственные векторы |
К2, ••>Кп матрицы U: |
|
К 12] ^ К о , |
№ = К |
(ст=1,2......... л). |
394 |
Н Е К О Т О Р Ы Е У С Л О В И Я У С Т О Й Ч И В О С Т И |
[ГЛ . X V |
|
Столбцы же матрицы |
и диагональные элементы диа |
||
гональной матрицы А1-*-1из условия (6.21) представляются соотношениями (см. гл. VIII, § 2)
К[ак1= Pcdik~l] + К0Ф 1 W = - л Ы к-Ч
Для наших целей нужно, чтобы столбцы матрицы пре образования (6.16) имели одинаковую норму, определен ную положительной функцией а (/). Это требование мож
но выполнить соответствующим подбором матрицы Z. Примем
Z = diag |
«(О |
а (0 |
’ II |
a (t) |
( 6. 22) |
|
II /СІт) |] |
К т)ІІ ’ |
И |
|
|
Тогда норма столбцов матрицы преобразования (6.16) |
|||||
будет совпадать |
с а (г!), а |
матрица |
Z~‘ A.{m)Z — Z~l |
= |
|
= A — Z — будет иметь диагональный вид.
Возвращаясь к исходному уравнению (6.1), произведем замену переменных
x = |
K im) (t)Z(t) у, |
где К іт) определяется |
суммой (6.19), а Z — формулой |
(6.22) при е = 1. Тем самым уравнение возмущенного дви жения (6.1) приводится к виду
J l - = (л'"0 — Z“1 \ij — Z lMtm)N{m)Zy +
+ Z~lM{m)h(t, K {m)Zy).(6.23)
Согласно теореме 6.3, если
а (t) < со (t),
то для устойчивости невозмущенного процесса (тривиаль ного решения уравнения (6.1)) на заданном промежутке [t0, Т) (Т < оо) достаточно выполнения неравенства
і
- г ~ г |
І К (П + |
vraax(f)]d t ' < ~ b |
( b > 0). |
I 10 |
о |
|
|
|
‘0 |
|
|
В данном случае, как это следует из (6.23), |
|||
р0 (/) = таха |
Re*T (* ) -- !- In |
а ( 0 |
|
ІІА(0"'> (011 |
|||
§ 6] |
ОБ УСТОЙЧИВОСТИ НА КОНЕЧНОМ |
ПРОМ ЕЖ УТКЕ |
395 |
где |
= 2 Л-а , vmax (t) — максимальное собственное зна- |
||
|
/=о |
|
|
чение эрмитовой матрицы |
|
|
|
|
Р<"'> = — - L (M{m)NUn) + |
/Ѵ(т)’М*). |
|
Для устойчивости на [/„, Т) достаточно также выполне |
|||
ния |
неравенства |
|
|
шаха R eC ’ (0 — dt In |
a(t) |
“I- Ттах (t) 0 |
|
IlC’WII |
|||
|
(t G[^o> T)).
Г л а в а XVI
УСТОЙЧИВОСТЬ ПРОЦЕССОВ ОТНОСИТЕЛЬНО
ЗАДАННОЙ ОБЛАСТИ ПРЕДЕЛЬНЫХОТКЛОНЕНИЙ
В определении устойчивости на заданном промежутке времени, введенном в гл. XIV, область предельных откло нений формируется посредством матрицы заданного класса
/С“. При этом матрица G (і), определяющая область пре
дельных отклонений, не рассматривается наперед заданной, так что и сама область предельных отклонений не является наперед заданной: для устойчивости процесса требуется
лишь существование в классе Ад такой матрицы, для ко торой условия (14.3.1) и (14.3.2) соблюдаются.
Наряду с этим представляется целесообразным введе ние и понятия устойчивости на заданном промежутке с априорным заданием конкретной области предельных от клонений. В этом случае, конечно, хотя речь идет не о раз мерах области, а только о ее форме, понятие устойчивости приобретает в большей мере оттенок субъективности, и тем не менее рассмотрение устойчивости по отношению к заданной области представляет и теоретический, и извест ный практический интерес (хотя бы, например, потому, что из устойчивости процесса по отношению к области, опре
деленной через матрицу класса Ад, немедленно следует устойчивость в смысле введенного в гл. XIV понятия устой чивости процесса на заданном промежутке времени).
§ 1. Понятие устойчивости относительно заданной области
О п р е д е л е н и е . |
Если при достаточно малом р > О |
|
любое возмущение х (t) |
процесса, начальное значение |
х0 = |
= X (t0) которого удовлетворяет условию |
|
|
(5 (t0)х0, S (t0)х0) < р2, |
( 1. 1) |
|
5 п |
|
П О Н Я Т И Е У С Т О Й Ч И В О С Т И |
|
|
3 9 7 |
|
на промежутке t0 С ( ■ < Т удовлетворяет условию |
|
|||||
|
|
(S(t)x, S (t) х) < |
р2, |
|
|
(1.2) |
где |
S (і) — заданная ограниченная |
матрица, |
то |
невозму |
||
щенный процесс устойчив на промежутке U0, |
Т). В против |
|||||
ном |
случае — неустойчив. |
|
1, 2, |
.... п) — |
||
Область предельных отклонений xs (s = |
||||||
элементов |
столбцовой матрицы х — задается |
посредством |
||||
функции |
V (/, х) = (5 (/),ѵ, 5 (і)х), определяемой |
матрицей |
||||
5 (і). В зависимости от способа задания 5 (і) |
область пре |
|||||
дельных отклонений приобретает тот или иной вид. На пример, полагая
S(t) = g(t)E ,
где g (I) — некоторая ограниченная скалярная функция, ограниченная снизу положительной константой (g (/) > > а > 0), Е — единичная матрица, получим область пре
дельных отклонений в форме шара с радиусом, равным Р/g(t)-
II * (OIK |
Р |
|
г ( 0 |
Если 5 — постоянная квадратная матрица общего вида, то область предельных отклонений представляет собой ц-мерный эллипсоид с неизменными параметрами. В более общем случае переменной матрицы 5 (t) соотношения (1.1)
и (1.2) представляют эллипсоид с параметрами, изменяю щимися по t.
Таким образом, условия выбора 5 (t) позволяют конст
руировать область предельных отклонений разнообразной формы.
Ниже исследуются условия устойчивости и неустойчи вости невозмущенного процесса на конечном промежутке Н0, Т) (Т < оо) относительно заданной области предельных
отклонений применительно к процессам, которые описы ваются векторно-матричными уравнениями
^ f = U ( t ) x |
(1.3) |
или |
|
U (t)x + h{t, х), |
(1.4) |
где U — квадратная матрица порядка п; х, h — столб цовые матрицы; U и h — непрерывные функции своих
398 |
У С Т О Й Ч И В О С Т Ь П Р О Ц Е С С О В |
[ГЛ . XVI |
аргументов. Элементы матрицы /і — нелинейные функции отклонений xs — таковы, что равномерно по / в пределах промежутка [/Q, Т)
lim h (t, JC) = |
0 |
. |
(1.5) |
лг-И) ~ w |
|
|
§ 2. Устойчивость процесса относительно области, определяемой каноническим преобразованием уравнений
Допустим, что к (і) — ограниченная, невырожденная и дифференцируемая квадратная матрица порядка п и такая,
что замена переменных
х = К (і)9 |
(2 .1) |
приводит линейное уравнение (1.3) (линейную часть урав нения (1.4)) к каноническому виду
|
4 - |
= Л(*)у, |
|
(2.2) |
где Л (0 = diag (A.J, |
...... |
Хп), а А,,-(/= 1, |
2, ..., |
^ — не |
прерывные скалярные функции от t. Полагая 5 (() = |
/С-1 ((), |
|||
область предельных отклонений представим в виде |
|
|||
I/ ((, X)= |
(/С 1(0 Л', К~' (І) х) < |
р2. |
(2.3) |
|
Итак, выясним условия устойчивости невозмущенного процесса относительно области предельных отклонений (2.3), рассматривая К (t) как заданную матрицу. Сначала
рассмотрим вопрос о существовании такого конечного про межутка времени П0, /0 -{-■ А/), на котором условия устой
чивости невозмущенного процесса соблюдаются, а затем займемся условиями устойчивости процесса на заданном промежутке.
2.1. О существовании конечного промежутка устойчи вости .
2.1.1. Линейный процесс. Учитывая диаго нальную структуру матрицы Л, из (2.2) получаем следую щее дифференциальное уравнение относительно нормы век тора у:
d \\ у И |
П |
(2.4) |
|
У] ReК |
|||
dt |
|||
С=1 |
|
||
|
матри |
||
Здесь уа (а = 1, 2...... п) — элементы столбцовой |
|||
цы у. |
|
|
|
§ 2] |
У С Т О Й Ч И В О С Т Ь П Р О Ц Е С С А |
3 9 9 |
Всоответствии с (2.4) полная производная от положитель но определенной функции V (t, х) = Цу |2 в силу уравнения
возмущенного процесса равна
-^-=2 É ИеЯаЫ2- |
(2.5) |
|
Ш |
<J=1 |
|
Теорема 2.1. Если |
|
|
М 0 )< 0 |
(Р-о (О = таха ReХа (t)), |
(2.6) |
то существует конечный промежуток U0, (0 -f АО, на ко
тором линейный процесс (тривиальное решение уравнения
(1.3 ) ) обладает устойчивостью по отношению к области
(2.3) .
Доказ ательство. Из (2.5) следует, что
|
X |
" W < И-о(0 J JIУ°2I= |
Цо ( 0 IУ2II- |
|
|
Отсюда |
видно, |
что если имеет место неравенство |
(2.6), |
||
то в точке |
t = |
ta, |
а по непрерывности и в пределах неко |
||
торого конечного |
промежутка U0, |
t0 + АО er П0, Т } |
пол |
||
ная производная по t от положительно определенной функ ции V (/, X) в силу уравнений возмущенного процесса удов
летворяет условию
что доказывает теорему.
Теорема 2.2. Если |
|
f-lo(0>) 0, |
(2.7) |
то не существует конечного промежутка [f0, t0 + |
АО, на |
котором линейный процесс (тривиальное решение уравне
ния (1.3)) обладал бы устойчивостью по отношению к области
(2.3), т. е. At = 0.
Доказ ательство. Допустим для определеннос
ти, что |
|
Po (0) = |
К (0>)- |
Пусть |
|
Ф(0 у) = |
Е Re |
