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

книги из ГПНТБ / Формальский А.М. Управляемость и устойчивость систем с ограниченными ресурсами

.pdf
Скачиваний:
8
Добавлен:
24.10.2023
Размер:
11.01 Mб
Скачать

3l СТРУКТУРА ОБЛАСТЕЙ УПРАВЛЯЕМОСТИ 6t

следовательно, z1 GE RT- Пусть теперь задана точка х1 GE

£= R™ и произвольная

точка

G

тогда х- = -х* +

+ X2 Е Е Q?.

Действительно, во-первых,

х — хх + хг ЕЕ

Е Е XPg, во-вторых, для

всякого

вектора

г )° G E Xjg имеем

d (р0)

I p V I

= I р° (ж1

+ X2) I

= I rfx |.

Таким образом, х G E Q? тогда и только тогда, когда хх GE

е і ? Г (рис. 7.1).

Вследствие того, что неравенства (7.5) являются стро­

гими,

множество 7?™ в пространстве Xjg, а

значит, и об­

ласть

в пространстве Хр , представляет

собой откры­

тое множество.

Если ps = тг (система (2.2) вполне управляема), то Хр =

== Хп. При этом уравнения' (7.1) линейно независимы, по­ скольку система (7.1) являет­ ся частью системы (3.4) ран­ га ps = п. Размерность под­

пространства Xps = Хп в этом

случае равна размерности фундаментальной системы ре­ шений уравнений (7.1). Раз­ мерность фундаментальной системы равна

Г3 г,

п2 Рк ~ 2 Pkt

 

Л=гч+1

(с=1

 

 

т.

е. равна количеству собст-

Рис. 7.1.

 

венных значений матрицы А

частями

с учетом

с

положительными действительными

их

кратностей.

 

 

 

 

Итак, из сказанного вытекает следующая теорема о

структуре областей управляемости.

 

 

 

Теорема 7.1. Облаетъ управляемости Q™ {т =

1, 2, 3) —

цилиндрическое

множество, т. е. Q™ = R™ +

XPg, где

R™ С Xjg — ограниченное открытое множество (сечение цилиндра). При рs = п размерность подпространства

62

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

[ГЛ. I

Хр„ равна числу собственных значений матрицы А с по­ ложительными действительными частями с учетом их кратностей.

Сформулированный результат для случая т = 1 со­ держится, например, в работе 128].

Полученные в настоящем параграфе формулы для рас­ стояний d (rj°) до опорных гиперплоскостей ПГ (т]°) об­ ластей управляемости QГ (т = 1, 2, 3) могут быть исполь­ зованы для построения областей управляемости в кон­ кретных случаях. Это видно в приведенных ниже частных случаях, а также иллюстрируется на примерах, рассмат­ риваемых в конце настоящей главы.

Рассмотрим при условии рs =

п два частных случая.

1. Все собственные значения

матрицы А имеют

не­

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

действительные

части.

Из теоремы

7.1

следует, что (?”

= Х п для т =

1, 2, 3.

Если обратиться

непосредственно к системе уравнений (7.1), то легко уви­ деть, что она совпадает в настоящем случае с системой (3.4), которая при ps = п имеет только нулевое решение. Следовательно, d (р, Т) -> оо при Т —>■ оо для всякого вектора г) ф 0. Для m = 1 настоящий случай при неко­ торых дополнительных предположениях рассматривался,

например, в [9, 26, 47].

2. Матрица А имеет одно положительное действитель­ ное собственное значение Я15 кратность которого p t = 1, все остальные собственные значения имеют неположитель­ ные действительные части. В этом случае система (7.1) состоит из п — 1 уравнений, которые при ps = п линейно независимы. Уравнения (7.1), (7.2) имеют при этом толь­ ко два решения, отличающиеся друг от друга знаком: т)° и —т]°. Для функции I r\°e~Ar:bs1имеет место выражение

 

I

ц0е~АѢ еI

= [ ті°а1)0 bs \ е~^\

Найдем при m =

1, 2, 3 выражения для расстояний d (ц°).

При m =

1 из выражения (7.3)

получаем

d (г|°) -

 

р

I г

м

Мѣ^ I г|0а1)(Д

I т)°а1?0Ь, |. (7,19)

О

3]

СТРУКТУРА ОБЛАСТЕЙ УПРАВЛЯЕМОСТИ

63

При т =

2

из выражения (7.7)

получаем формулу

(Л°) =

 

(40ahA )2e - ^ d x

= ] /

| ifciiA

|- (7-20)

 

 

о

 

 

 

 

 

Пусть

теперь

т = 3.

Для

определения расстояния

d (р0) по

формуле

(7.9)

нужно

из уравнения (4.9) при

Т = оо найти величину о0.

 

 

 

График

функции ф (т) =

| p°ai,ofrs

I

показан

на рис. 7.2.

Легко видеть (см. рис. 7.2), что если искомая

величина а0 удовлетворяет

условию a0Ms <^

| p°ai,o&3 l>

то

 

 

 

 

 

 

 

 

 

 

Es ( oo, (To)

 

To),

Gs ( oo,

n0)

(Tq,

oo),

 

 

T0 =

----r1,ln

„----r-

> 0 .

 

 

 

 

 

 

h

 

I ^ ai,0bs

 

 

 

При этом уравнение (4.9)

записывается в виде

 

 

aaM .

 

 

W

 

 

 

 

 

 

 

i_ l*

 

 

 

= Ря

- І Г М‘ 1п

I B°«i,ob.

 

„2 ^

(ті0«і,о^)2 e~aW t

После

вычисления

интеграла

получаем

 

 

 

 

1 Л/Го .

 

soМ

_1

M \

 

 

 

Ms ln

 

 

L — p

 

 

p°ai,obJ ^

2^i

r *‘

 

Отсюда

имеем

 

 

 

 

ö0 = I П°*іA I

 

 

 

Tn = — —

2A.1

 

 

 

 

,\p

 

 

 

M.

 

 

64 ОБЛАСТИ УПРАВЛЯЕМОСТИ [ГЛ. I

Эти равенства имеют смысл при 2ХгР а 7> М\. Из выраже­

ния (7.9)

при атом получаем

 

 

 

т0

оо

ПУП.

d (т|°) = Ms $ | 4°ah0bs I e-^dx +

^

(л0«!A ) 2 e~2^xdx =

 

О

-c0

 

 

Л°аіА I (1 — Vo) + 2x7 M s I

 

 

 

8.\

 

Xi Ms I if«! A

M2s ]•

Если

2куРа ^ Ml, то уравнение (4.9) записывается

в виде

СЮ

 

 

 

 

 

-V $ (лп<*іA )2e~2^xdx = Ps-

ao 0

Пользуясь этим соотношением, из (7.9) получаем

d(T|°) =

№ ь і , А У е 2ХіТ^ = У А - 111°«!А |.

Таким образом, при т = 3 имеем

*>р,

ІИ2в при 2ХіPs)>M 2i

при 2A,iPs Mf. (7.21)

Нижняя строка в (7.21) совпадает с выражением (7.20).

Область (?Г (иг = 1, 2, 3), как следует из теоремы 7.1, представляет собой в рассматриваемом частном случае множество точек х ЕЕ Х п, заключенных между двумя ги­

перплоскостями ПГ(л°) и ПГ (—г)°), ортогональными век­ тору т]0 и находящимися на расстоянии d (л°) от начала координат. Этот результат для т = 1 содержится, на­ пример, в работах [26, 30].

з]

СТРУКТУРА ОБЛАСТЕЙ УПРАВЛЯЕМОСТИ

65

В случае рs = п можно выписать систему уравнений, эквивалентную системе (7.1) и не зависящую от вектора bs. Система (3.4) г гоит из rs групп уравнений

r\aklbs = О

(I — 0, 1,. .

p h — 1),

(7.22)

соответствующих корням Хк уравнения (3.1). В группе (7.22) содержится p h уравнений. Если ранг системы (3.4)

Ps = п, то ранг системы (7.22) равен p k, поскольку

п

2

Ph — л,

следовательно, столбцы

aklbs

{1 =

0,

1, . . .

fe=i

 

р к — 1)

линейно независимы.

 

Эти

столбцы

пред­

. . .,

 

ставляют собой линейные

комбинации столбцов матриц

а к1

{1 = 0,

1,

. . .,

p h — 1).

Среди столбцов

матриц а кІ

{1 =

0 ,1 , .

. .,

p h — 1) содержится ровно р к линейно не­

зависимых

[21]. Отсюда следует,

что система р к уравне­

ний

(7.22)

эквивалентна системе

прк уравнений

 

 

 

Н «кі ~ О

(Z =

0 ,

1 ,

 

. . . ,

р к — 1 ) .

 

 

Значит, при рs =

п система

(7.1) эквивалентна системе

 

 

 

 

 

г\ак1 = 0

 

 

 

 

 

 

(7.23)

 

 

{I = 0, 1, . . ., р к — 1;

 

к = гг + 1,. . ., г3).

 

Из рассмотрения уравнений (7.23) видно, что векторы

ц0 при рs =

п определяются только матрицей И и не зави­

сят

 

от столбца bs.

 

 

 

 

 

 

 

 

 

 

 

Перейдем к выяснению структуры области управля­

емости Qm {т = 1, 2, 3).

 

 

 

 

 

 

 

 

 

 

 

Система

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r\aklbs =

0

 

 

 

 

 

 

(7.24)

{I = 0, 1,. . .,

Рй1\к =

гг +

1,. . .,

r3;s

=

1,

. . ., г),

полученная из системы (7.1) при s =

 

1, . . . ,

г,

содержит

 

Т*з

 

 

 

 

 

 

 

 

 

 

 

 

 

г

2

Рк уравнений.

 

 

 

 

 

 

 

 

 

 

fc=r,+i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Система (3.15) состоит из г3 групп уравнений

 

r\aklbs = 0

{ 1 = 0 ,

1, . . ., р к — 1;

s = 1, . . ., г),

(7.25)

соответствующих корням Хк уравнения (3.1). Среди

столбцов матриц a kl {I = 0, 1, . . ., р к — 1) содержится

3 А. М. Формальский

66

 

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

 

[ГЛ. f

p k линейно

независимых,

следовательно,

среди столбцов

а * А

(/ = 0,

1, . .

р к — 1;

s =

1,

и среди урав­

нений (7.25)

содержится

не

более р к линейно незави­

симых. Значит, в системе

(7.24)

содержится не

более

г,

Pu линейно

 

 

 

 

 

 

2

независимых уравнений.

 

fc=r,+l

 

 

 

 

 

п, то в системе (7.25) —

Если ранг системы (3.15) р =

ровно p h линейно независимых уравнений.

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

 

 

 

 

 

 

 

п

 

при р = п в системе (7.24) содержится ровно 2

Рк ли-

 

 

 

 

 

 

 

/С=Г,+1

 

нейно независимых уравнений, и система (7.24) эквива­ лентна системе (7.23).

Обозначим векторы ц, которые являются решениями уравнений (7.24), (7.2), через ц°.

Если вектор ц° удовлетворяет системе уравнений (3.15), то область управляемости Qm удовлетворяет условиям (3.19). Если вектор т)° не удовлетворяет системе (3.15), то область Qm удовлетворяет условиям вида (7.5), в ко­ торых d (г)°) есть расстояние до опорной гиперплоскости ГГ (rj°) области Qm. Это расстояние представляет собой

сумму расстояний до опорных плоскостей П™ (г|°) областей

Q? (s = 1,. . ., г).

В § 3 через Хр обозначено подпространство размерно­ сти р пространства Х п векторов ж, удовлетворяющих ра­

венствам (3.19). Через X J обозначим подпространство, на­ тянутое на векторы rj°, принадлежащие пространству Х р.

Через Хр обозначим ортогональное дополнение подпро­ странства Хр до пространства Хр, т. е. Х \ _]_ Хр, Xj +

+ Хр = Хр. Если р = п (система (1.1) вполне управля­ ема), то Хр = Х„. При этом размерность подпростран­

ства Хр равна размерности фундаментальной системы решений уравнений (7.24),

Га

г ,

П - 2

Рк = 2 Рк-

/с=г,+і

/с=і

Обозначим через R m множество точек х еЕ Xj, удо­ влетворяющих неравенствам вида (7.5). Тогда имеет место следующая теорема.

3]

СТРУКТУРА ОБЛАСТЕЙ

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

67

 

Теорема 7.2. Область управляемости Qm (т =

1,2, 3) —

цилиндрическое множество, т.

е. Qm = R m +

Х 2Р, где

Rm f— ограниченное открытое множество (сечение

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

Формулировка теоремы 7.2 получается из формули­ ровки теоремы 7.1, если в последней всюду опустить индекс s.

Результаты, полученные выше для области Q™ при условии ps = п в двух частных случаях, могут быть легко получены в этих же частных случаях для области Qm при условии р = п.

Сформулированные в теоремах 7.1 и 7.2 результаты,

связанные со структурой областей Q™ и Qm (т = 1, 2, 3), получены без использования возможности приведения систем (2.2) и (1.1) к каноническому виду. В § 9 эти ре­ зультаты будут сформулированы в терминах, связанных

сканоническими переменными.

§8. Структура областей управляемости при т = 4 —7

Полученные в предыдущем параграфе выводы о струк­

туре

областей управляемости Q™ и Qm показывают, что

при

т = 1,

2, 3 собственные значения %k матрицы А

с нулевыми

действительными частями (eft = 0) влияют

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

частями.

При т — 4

7

собственные значения Kk,

для

которых Eft = 0, играют,

как будет видно ниже, другую

роль, нежели при т — 1, 2,

3.

 

 

Структура областей управляемости QГ,

Qm при т =

= 4-^ -7

отличается

от

структуры этих

областей

при

т = 1,

2, 3.

 

 

 

 

Q™

Для

того чтобы

выяснить

структуру

области";

(т — 4 -ч- 7), будем, так же как в § 7, искать векторы ц,

для

которых расстояние с? (г|, Т) остается ограниченным

при

Т

оо, и векторы ц, для которых d (ц, Т)

оо при

Т

ОО.

 

 

3*

68

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

[ГЛ. I

Корни kh уравнения (3.1) пронумерованы так, что содержащиеся в выражении це-АхЬв (см. формулу (3.3))

функции f-XftTT!, в которых

1 = 0 ,

1, • • •,Pk

1 при

ft = 1, . . ., /р и 1 =

0 при к = 7-J +

1, . . г2,

являются

ограниченными при 0

т

со. Остальные функции вида

е~х*т X1, содержащиеся в выражении (3.3), не являются ограниченными на полуоси [0, оо).

Приравняем в выражении (3.3) нулю коэффициенты при неограниченных функциях е"г^х1:

 

 

цак1Ь$ =

0

 

 

(8.1)

/ I =

1,

. .. , р к — 1

при ft =

щ +

1,. . . , г2, \

 

V I =

0,

1, . . . ,р к — 1

при

ft =

г2 +

1, .. . , г3 /

'

Будем рассматривать соотношения (8.1) как уравнения относительно компонент вектора г|. Количество уравнений, входящих в систему (8.1), определяется числом

Г 2 г 3 г 3

2

(Р/с — 1 )+

2 Р к =

2

Л - ( г2 — о).

(8.2)

й'= г1+ і

 

к = г г+1

7г=г,+1

 

В системе (8.1),

являющейся

частью системы

(7.1),

на (г2 — Tj)

уравнений меньше,

чем в системе (7.1). Ве­

личина г2 — гг равна числу корней уравнения (3.1) с ну­ левыми действительными частями без учета их кратностей.

Векторы р, которые являются решениями уравнений

(8.1), (7.2), обозначим через

р0. Тогда функция | xfe~Mbs |

остается

ограниченной при

т -> о о .

 

Функция d (т), Т) не убывает с ростом Т, поэтому

функция

d (т)°, Т) при пѣ = 4 (см. формулу (5.8)) при

Т — оо

стремится

к конечному пределу

 

 

d (р°, оо) =

d (р°) =

N s max | р°е~ЛтЬ8 ].

(8.3)

 

 

 

 

т е [ о , оо)

 

Если

вектор

р°

такой,

что р°е_Лт&5 = 0,

то d (р°,

Т) == 0,

d (р°) =

0.

Координаты точек х GE Qi удовлетво­

ряют условиям (3.10), т. е. область Q\ принадлежит под­ пространству Х р .

Если yfe~Mbs ф 0, то d (р°) Ф 0, область управляе­ мости Qi заключена при этом между плоскостями

3]

СТРУКТУРА ОБЛАСТЕЙ УПРАВЛЯЕМОСТИ

69

nf (г)°) (7.4). Для тех векторов ц0, для которых функция d (ц0, Т) строго монотонно возрастает с ростом Т, мно­

жество <2s (Т) «достигает» плоскостей

(т]°) лишь

при

Т —>

оо, координаты точек

х е= (І* удовлетворяют

при

этом

строгим неравенствам

(7.5).

 

 

Функция d (ц0, Т) не при всяком векторе т)° строго монотонно возрастает. Если, например, вектор ц0 такой,

что ц°e~Mbs

0

при т ->

сю (r\°e~Ar:bs ф 0),

то, очевидно,

существует такое значение Т , что d (ц0, Т) =

d (ц°, Т') =

= d (ц0) при всех значениях Т

Т’. Подобная ситуация

имеет

место

и

в

случае, когда

yfe~Azbs =

const Ф 0;

в этом

случае d (ц°,

Т) =

d (г)°)

=

N s | т)°e~Mbs | при всех

значениях Т >

0.

Для

таких

векторов

ц0

множество

Ql (Т)

с ростом

Т

«достигает»

плоскостей

(ц0) при

конечном значении Т и больше «не расширяется в направ­ лении т]°». При этом на плоскостях П* (т)°) существуют

точки, принадлежащие множеству (Д- Из сказанного выше следует существование таких

векторов т]°, что координаты точек х ЕЕ Qt удовлетворяют

нестрогим неравенствам

 

 

 

I фх I

d (ц0).

(8.4)

Для расстояний d (ц,

Т) при значениях

индекса

m = 5, 6, 7 имеет место оценка

 

d(p, T)*ElNs

шах \ r\ë~Axbs \.

(8.5)

 

те[о, т]

 

Если r\e~Axbs =

const, то неравенство (8.5) обращается

в равенство (см.

формулы

(5.18) и (6.1)).

 

Из неравенства (8.5) следует, что условиям (3.10),

(7.5), (8.4)

удовлетворяет не только область Qi, но также

и области

Q™ (т — 5, 6, 7). Формулы

для расстояний

d (ц0) получаются из выражений (5.17), (6.14), (6.16).

Пусть

теперь

т)

ц°. Докажем, что

при этом d (ц,

Т) —>■ оо,

если Т —*■ оо. Для случая m =

4

этот факт не

вызывает

сомнений,

поскольку

функция |

т}e~A^bs | при

т —*■ оо не остается ограниченной, если ц ф

ц0, и, следо­

вательно,

max

I т\ё~м ЪѣI —> оо

при Т

 

 

 

оо.

 

 

te[o,T]

 

 

 

 

 

70

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

[ГЛ. I

Рассмотрим случай т = 5 (при условии (6.2) такое рассмотрение охватывает также и случай т — 7).

Из соотношений (5.17), (5.13), (5.16) вытекает нера­ венство

d (т),

Т) > М Л о (7 > Я . (Г, х0(Г))

= N sXo (Л-

(8.6)

Покажем, что %о (Т) —ѵ оо при Г

оо. Если ц ^ ц°,

то T]e_ATbs

const. При этом функция \iEs (Т, %)

строго

монотонно возрастает с ростом Т для каждого фиксиро­ ванного значения %. При каждом фиксированном значе­

нии Т функция \iEs (Т, %) строго

монотонно убывает

с ростом величины %(если только %

%'). Отсюда сле­

дует, что функция Хо (Л> являющаяся решением урав­ нения (5.16), строго монотонно возрастает с ростом ве­

личины

Т.

 

 

 

 

 

Предположим, что при ц Ф ц0 функция хо (Т) огра­

ничена

некоторой

константой;

 

 

 

Ъ

<

С

при

0 <( Т < оо.

(8.7)

Тогда

при всех

0

Т

оо

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

 

цЯ5

(Т,

%0 (Т)) >

[iE, (Т, С).

(8.8)

Выражение для функции |r\e~Axbs | можно записать в виде (7.15). При этом, в отличие от случая m = 3, в рас­ сматриваемом”*' здесь случае пг = 5 имеем: если ц уь ц0,

то всегда е~гк'~ т1' -а- оо при т ->■ оо. В §’7 доказано, что если е~е]с'~ г 1' —>- оо при т —>■оо, то

\xEs (Г, С) оо при Г —V оо.

Из неравенства (8.8) получаем

(Л Хо (Л ) 00 при Г -> оо,

а это противоречит тому, что Хо (Л — решение уравнения (5.16). Предположение (8.7) приводит к противоречию, зна­ чит, Хо (Л при Т -*■ оо. Из неравенства (8.6) полу­ чаем, что d (г), Т) -а- оо при Т ->• оо.

Рассмотрим при условии (6.3) случай m = 7.

Пусть при некотором векторе ц Ф ц0 для всякого сколь угодно большого значения Т' найдется такое значение

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