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

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

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

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

а определители Д2 и Д3 получаются круговой заменой верх­ него индекса j в коэффициентах о43) определителя Аг.

Из формулы (8.14) получаем выражение для расстоя­

ния d (т)°) при m = 4 - г 7

 

I AiaW I

(12.38)

d (ц0) = N

/ д [ Г д | + д [

*

Из соотношений (12.37), (12.38) заключаем, что область управляемости Qm (т — 4 ч- 7) представляет собой мно­ жество точек, заключенных между двумя плоскостями

+ Д2т2 + А3з:3 = + N Aja^.

На'этих плоскостях есть точки, принадлежащие области Qm. Рассмотрим теперь случай, когда О, Х3 0. В этом

случае при т = 1,2, 3, уравнения (7.1) имеют вид

Th«i1> = 0,

(12.39)

Ріаз1’ + Ѵ2азѴ( + і7з«з3) = 0.

Система (12.31) вполне управляема, поэтому ранг системы уравнений (12.39) равен двум. Система уравнений (12.36), (12.39) имеет только два решения: rj° и — г]1*; при этом

у (3 )

„(2 )

 

цо = 0, т}« = — -J -, т)®=

,

(12.40)

где

А = "К(а32))2 + (а33))2.

Из выражений (7.19), (12.34), (12.40) получаем формулу для расстояния d (ц0) при т — 1:

d

Щ -=Iа3>2а>з—« 4 W |.

При т = 2 из выражения (7.20) получаем

112

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

[ГЛ. I

При т = 3 из

выражения (7.21) имеем

 

^

I а<»а'я - « S V ’ I [ l - -1- Д ‘

 

при 2К Р > М \

при 2Ік2Р <^М 2.

Область управляемости Qm (т = 1, 2, 3) представляет собой открытое множество точек, заключенных между дву­ мя плоскостями. При т = 1 уравнения этих плоскостей выглядят так:

 

а(з3>х2 — а(32)х3 = + ~

(о43)4 2) — 4 2)а(33)).

 

При т =

2 имеем

 

 

 

4

3)я 2 — с42Уз = +

(о43)4 2) — а 2Ѵ 33)).

 

Уравнения ограничивающих плоскостей при т =

3 выпи­

сывать не

будем.

(8.1)

превращается

в одно

При т = 4 -г- 7 система

уравнение:

 

 

 

 

%<*»5+

+ ^ a h

= °-

(12.41)

Уравнениям (12.36), (12.41) удовлетворяют всевозможные единичные векторы р, ортогональные вектору g =

= (ah ’ с42), а®). Как следует из выражения (12.34), макси­

мум

функции г\0е~АхВ

достигается при т =

0, либо при

т = оо.

Поэтому для

т =

4,

например,

из

формулы

(8.3)

получаем

 

 

 

 

 

 

 

 

N I р?(а)1*+

аЬ) +

р2а22) + Рз<43) |

 

 

о

,

при I рѴі0 К

h l (аГ) +

«21’) + p2f42) +

р"«23) I,

d (Т| )

 

N I рУЬ I

при I pV h I > I р? (ah +

a h ) +

 

 

 

 

 

 

+ P2«h +

p°ah |.

Множество внутренних точек области управляемости Qm при всех т = 4 7 представляет собой внутренность цилиндра, ось которого направлена по вектору g.

3J

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

ИЗ

 

Рассмотрим случай, когда Re А2

0, Re А3 >

0. В этом

случае при т = 1, 2, 3 система (7.1) превращается в одно уравнение

= 0. (12.42)

Величина ф 0, т. е. ранг системы (12.42), равен еди­ нице, поскольку р — п. Следовательно, уравнениям (12.36), (12.42) удовлетворяют всевозможные единичные векторы г]0, у которых

Л°і = 0.

Область управляемости Qm при т — 1, 2, 3 представ­ ляет собой множество внутренних точек цилиндра, осью которого является ось хѵ

При т = 4 -ь- 7 область Qm представляет собой в рас­ сматриваемом случае ограниченное в фазовом простран­ стве Х 3 множество.

Подводя итоги исследования примера (12.31), отметим, что во всех рассмотренных выше трех случаях располо­ жения на комплексной плоскости собственных значений А2 и А3 область управляемости Qm при т = А ~ 1 «мень­

ше», нежели при т — 1,

2, 3. Действительно, в случае,

когда Re А2 < 0, Re А3 <

0,

область Qm при т = 1,

2, 3

совпадает с пространством Х 3, а при т = 4

7 занимает

часть пространства Х 3. В случае, когда Х2

0, А3 <С 0,

область Qm при т =

1, 2,

3 заключена между двумя

пло­

скостями, а при т =

4 -у- 7

заключена внутри цилиндри­

ческой поверхности. Наконец, в случае, когда Re Я2

0,

Re А3 )> 0,

при т =

1, 2,

3

область Qm является неогра­

ниченным

множеством,

а

при

т — 4 -г- 7 — ограни­

ченным.

 

структуры

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

Эти особенности

связаны с наличием в системе (12.31) нулевого собствен­ ного значения.

В главе II (см. § 16), в связи с рассматриваемыми там задачами, будут приведены еще некоторые примеры по­ строения областей Qm при т = 5.

ГЛ Д В А II

НЕКОТОРЫЕ ЗАДАЧИ БЫСТРОДЕЙСТВИЯ

1.ЗАДАЧА БЫСТРОДЕЙСТВИЯ В СИСТЕМАХ

СОГРАНИЧЕННЫМИ ПО ВЕЛИЧИНЕ И ИМПУЛЬСУ УПРАВЛЯЮЩИМИ СИЛАМИ

Вряде реальных систем в качестве органов регулиро­ вания применяются реактивные двигатели. Управляющей силой в таких системах является тяга двигателя, пропор­ циональная секундному расходу топлива. Секундный расход топлива ограничен, поэтому ограничена величина управляющей силы (условие вида (1.2)). Запас топлива двигателя также является величиной конечной, поэтому ограниченным будет импульс управляющей силы (условие вида (1.4)). В настоящем разделе рассматривается задача построения оптимального по быстродействию управле­ ния в случае, когда управление удовлетворяет одновре­ менно двум указанным выше условиям.

При наличии интегральных ограничений на управление задача синтеза оптимального управления обладает опре­ деленной спецификой. В случае ограничения вида (1.4), например, при синтезе управления нужно в качестве до­ полнительной фазовой переменной вводить величину, ха­ рактеризующую «запас топлива на борту объекта» в теку­ щий момент времени.

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

§13. Постановка вопроса

Будем рассматривать управляемую систему, которая описывается матричным уравнением (1.1) (уравнением (2.2) при наличии одной управляющей силы).

В качестве допустимых управлений примем измери­ мые функции us (т) (s = 1, . . ., г), удовлетворяющие од­

1] ОГРАНИЧЕННЫЕ НО ВЕЛИЧИНЕ И ИМПУЛЬСУ СИЛЫ 115

новременно двум неравенствам (1.2), (1.4). В соответствии с терминологией главы I, допустимыми управлениями счи­

таются функции из множества Qm (Q^) при т = 5. Ин­ декс т в настоящем разделе будем опускать, если это не затрудняет понимания, всюду, где т = 5.

Будем рассматривать задачу быстрейшего приведения системы (1.1) (системы (2.2)) в начало координат с по­ мощью управления, удовлетворяющего условиям (1.2), (1.4), в частности, задачу синтеза оптимального управ­ ления.

Из работ, посвященных задаче оптимального быстро­ действия с несколькими ограничениями на управление, отметим статьи [17, 63, 65].

В работах [7, 9, 36а, 47] решена задача оптимального быстродействия в случае, когда на управление наложено только ограничение вида (1.2). Оптимальное управление при этом описывается выражением (4.1). В работе [9] при решении задачи существенно используется, что в этом слу­ чае область достижимости системы (2.2), если она является вполне управляемой (р8 = п), представляет собой строго выпуклое множество. При ограничениях (1.2), (1.4) об­ ласть достижимости системы (2.2), как будет показано ни­ же, не является, вообще говоря, строго выпуклой даже при р„ = п, т. е. граница этой области содержит плоские участки. Это обстоятельство, осложняющее решение за­ дачи быстродействия, не рассматривалось в известных ав­ тору работах. Здесь излагаются результаты, полученные в [55е].

При наличии только ограничения (1.2) оптимальное по времени (минимальное время в этом случае обозначим через Ѳ = Ѳ (х)) управление является релейным. Задача синтеза такого управления заключается в «разбиении»

пространства Х п, состоящего из фазовых координат

хи ...

...,

хп, поверхностями переключения на области,

в кото­

рых

управления

и$ (т) принимают значения M s и —Ms

(s =

1, . . ., г).

Если такое разбиение произведено, то

известно оптимальное управление как функция фазовых координат и — и' (х).

При наличии ограничений (1.2), (1.4) постановка за­ дачи синтеза видоизменяется. В начальный момент вре­ мени оптимальное управление зависит не только от вектора X (0) начального состояния, но также и от вектора N (0) =

116

 

ЗАДАЧИ

БЫСТРОДЕЙСТВИЯ

 

 

 

ІГЛ. II

N =

I Ns I , T.

e.

и = и (x (0),

jV (0)).

В текущий

момент

времени t

имеем

и =

и (x(t),

N

(t)), где

х (t)

состояние

системы

(1.1)

в момент

t,

а

вектор

N (t)~

= I N s (t) I

характеризует

ограничение на величину им­

пульса управляющей силы в этот момент,

 

 

 

 

ОО

 

 

 

 

 

 

 

 

 

 

 

^ \ u s{ x ) \ d x ^ N s(t)

(s =

l,...,r).

(13.1)

 

 

t

 

 

 

 

 

 

 

 

 

Введем координаты xn+s (t), определяемые дифферен­

циальными

уравнениями

 

 

 

 

 

 

 

 

 

dx

 

 

 

(s = l,...,r).

 

(13.2)

 

 

—^ i = - K ( * ) l

 

Решение уравнения (13.2) имеет вид

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

tfn+s (t) =

ж„+8 (0) — §1щ (т) I dx.

 

(13.3)

 

 

 

 

 

 

о

 

 

 

 

 

Пусть

xn+s (0) = Ns

= Ns (0)

(s =

1,

. . .,

г), тогда из

соотношений (1.4),

(13.3)

имеем

 

 

 

 

 

 

 

Оо

 

 

 

t

 

 

 

 

 

 

 

\ I us (т) I d t <

xn+s (t) -f § I us (t) I dx.

 

 

 

0

 

 

 

0

 

 

 

 

 

Отсюда получаем

 

 

 

 

 

 

 

 

 

 

 

оо

 

 

 

 

 

 

 

 

 

 

 

$ 1u$(т) I dx < xn+a(f)

(s =

1,...,

r).

(13.4)

 

 

t

 

 

 

 

 

 

 

 

 

Из

соотношений

(13.1),(13.3),(13.4) следует,

что

xn+s (t)

= Ns (t). Величина xn+s (t) характеризует ограни­

чение на импульс управляющей силы us (т) в текущий мо­ мент времени t.

Таким образом, задача синтеза оптимального управ­ ления при ограничениях (1.2), (1.4) состоит в построении

функции и = и (жц . . .,

хп, хп+1, . . ., хп+г).

Обозначим

через Хп+Г пространство,

состоящее из фазовых координат

хх, . . .,

хп системы (1.1)

и фазовых координат

жп+1, . . .

, . ., хп+г

системы (13.2).

 

 

1І ОГРАНИЧЕННЫЕ ПО ВЕЛИЧИНЕ И ИМПУЛЬСУ СИЛЫ Ц 7

В области пространства Х„+г, обусловленной неравен­

ствами

1

 

 

 

 

%п+s

(^ — 1, ..., г),

(13.5)

ѳ (^i? •••) %п) ~Tj

 

yl s

 

 

 

оптимальное

управление

не зависит от координат

xn+s

(s = 1,

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

оптимальное управление,

найденное при начальных условиях, удовлетворяющих неравенствам (13.5), и при ограничении (1.2), удовлетворя­ ет, очевидно, ограничению (1.4). Следовательно, в области (13.5) синтез осуществляется с помощью функции и —

= и' (X).

При решении задачи быстродействия будем использо­ вать области достижимости Qs (Т), Q (Т) и их свойства, рассмотренные в предыдущей главе. Кроме того, докажем некоторые новые свойства этих областей, необходимые для решения задачи. В следующем параграфе как раз иссле­ дуются новые свойства областей достижимости.

§ 14. Плоские участки и угловые точки границ областей достижимости

Если

Т <; xn+s/Ms, то

$ (Т) = Ql (Т). При

этом

множество Ql (Т), так же как и Q\

(Т), будет строго вы­

пуклым в подпространстве

X Pg. В настоящем параграфе

будем предполагать

выполненным

неравенство

(5.10),

т. е. неравенство Т

xn+s/M s-

максимизирующего •

При

подстановке

управления,

функционал (2.10), в выражение (2.1) получается вектор (Т) (см. рис. 2.1) координат точки, которая является общей для опорной гиперплоскости П8 (ц, Т) и множе­ ства Qs (Т) (точки касания). Если вектор ц такой, что Y\e~Azbs ф const, то, как показано в § 5, функция us (т, %0) (5.12) будет единственным управлением, которое макси­ мизирует функционал (2.10). Гиперплоскость П« (ц, Т), ортогональная такому вектору ц, имеет только одну общую

с областью Qs (Т) течку.

Нумерация собственных значений Xk введена в § 3 та­ ким образом, что Re Kk = 0 при к = гг -f- 1, . . ., г2. Ус­ ловимся, кроме того, о том, что при к = гг 4-1 собственное значение = 0. Тогда из соотношения (3.3) следует, что

118

ЗАДАЧИ БЫСТРОДЕЙСТВИЯ

[ГЛ. ir

т|e~Mbs =

<ss = const при тех и только

тех векторах rj,

которые удовлетворяют системе линейных алгебраических

уравнений

ЩкА = 0

 

(14.1)

 

 

=

— і

при

А = г х + 1 \

 

\l = 0, 1, . . рк — 1

при

к ф Г і + і )

 

В системе (14.1) содержится п — 1 уравнений.

Обозна­

чим решения уравнений (7.2),

(14.1)

через ц°. Могут пред­

ставиться два

случая:

 

 

 

 

т]°е Axbs =

<3S =

0,

(14.2)

 

■x\°e~Ai:bs =

crs ф

0.

(14.3)

В обоих этих случаях управление, максимизирующее функционал (2.10), не является единственным. В случае (14.2) множество Qs (Т) целиком принадлежит гиперпло­

скости г)°х — 0.

Подобная ситуация имеет место для об­

ластей Q™ (Т) при

всех

значениях

тп

1,

. . ., 7. При

условии

рs — п

может

представиться

только

случай

(14.3).

 

п и, кроме того, матрица А не имеет нулевого

Если ps =

собственного

значения,

 

то цe~Axbs ф const

ни

при ка­

ких векторах

ц Ф 0.

Поэтому

имеет

место

следующая

теорема.

 

 

Если,

 

система

(2.2) вполне управляема

Теорема 14.1.

 

. (ps = п),

матрица

А

не

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

собственного

значения,

то множество

Qs (Т)

является

строго вы­

пуклым.

 

п,

то систему (2.2) можно представить в ви­

Если ps <

де (3.11),

(3.12).

При этом из теоремы 14.1 следует, что

если матрица Ап не имеет нулевого собственного значения, то множество Qs (Т) является строго выпуклым в подпро­ странстве Х р (подпространство совпадает с фазовым

пространством системы (3.11)).

Рассмотрим теперь при условии ps = п случай (14.3). При условии рз = п все п — 1 уравнений, входящих в си­ стему (14.1), будут линейно независимыми. При этом урав­ нения (7.2), (14.1) имеют только два отличающиеся друг от друга знаком решения: т)° и — ц°.

П ОГРАНИЧЕННЫЕ ПО ВЕЛИЧИНЕ И ИМПУЛЬСУ СИЛЫ Цд

Обозначим через соя (Г) множество измеримых функций us (т), удовлетворяющих одновременно неравенству (1.2) и условиям

 

0'US(т)>0

(те [О, Г]),

 

 

 

т

 

 

 

 

 

 

\

щ (т) dx =

xn+s sgn C5S .

 

(14.4)

 

О

 

 

 

 

 

При Т =

хп+$/Ms класс

cos (Т) состоит из одной функции

Ms sgn Os

(остальные отличаются от нее

на

множестве

меры нуль). При

Т )>

 

 

 

 

> xn+s/Ms класс cos (Т)

 

 

 

 

содержит

неединствен­

 

 

 

 

ную функцию.

что

 

 

 

 

Легко

видеть,

 

 

 

 

соs{T) сі Qs (Г)-

Каж­

 

 

 

 

дая функция us (т) е

 

 

 

 

ecos (Г) доставляет мак­

 

 

 

 

симум интегралу

(2.10)

 

 

 

 

при т] = ц°. Этот макси­

 

 

 

 

мум определяется

фор­

 

 

 

 

мулой (5.18).

через

 

 

 

 

Обозначим

 

 

 

 

P s (Т) множество точек,

 

 

 

 

каждая из которых при­

 

 

 

 

надлежит одновременно

 

 

 

 

области Qs (Т) и опор­

 

 

 

 

ной гиперплоскости Пя (т]°, Т) (рис. 14.1).

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

точек касания определяется выражением

 

 

 

Ps (Т)

= {г;8

{Т):

us (т) €Е со* (Г)}.

(14.5)

Через Ps (Т) обозначим совокупность точек, симметрич­ ных точкам множества P s (Т). Гиперплоскость Пя (т)°, Т) остается без изменения при всех Т xnJrSIMs, т. е. Пя (т]°, Т) = Пя (ц°). Уравнение плоскости П8 (ц0) имеет вид

f]°x = d (Т]°) = жп+я I Пя |.

(14.6)

Множество Ps (Т) представляет собой сечение области Qs (Т) плоскостью (14.6), поэтому оно, так же как и Qs (Т), обладает свойствами Iе, 2°, 4е.

Докажем следующую лемму.

120

ЗАДАЧИ

БЫСТРОДЕЙСТВИЯ

 

[ГЛ. II

Лемма 1.

Множество P s (Т)

при Т

xn+s/M s

содер­

жит более одной точки.

 

 

 

 

Будем считать, для определенности, что as 0. Рассмот­

рим однопараметрическое семейство функций (рис.

14.2):

 

Ms

при

т

а, а +

■^n+S

 

и,(т,

а) =

 

 

 

м ,

(14.7)

 

т ф

 

жп+,

 

0

при

j^ct,

 

 

м,

 

 

 

 

 

 

 

где параметр 0 а летворяют условиям

, utcv,aу

1

1

1

1

. 1

Т xn+s/M s. Функции (14.7) удов­

(14.4), поэтому us (т, а) GE cos (Т).

і

1

1

1

1

т *

Ms

Рис. 14.2.

При подстановке управления (14.7) в формулу (2.1) полу­ чаем векторы vs (Т, а) GE Ps (Г),

 

а+

 

vs (Т, а) = Ms

e -^ h ß т.

(14.8)

 

а

 

Для доказательства леммы достаточно показать, что мно­ жество (14.8) содержит более одного вектора.

Продифференцируем соотношение (14.8) по пара­ метру а:

dvs ( Г , а)

 

-А а +

xn-\-z

 

 

Ms

Мв

b, е~АаЪя

(14.9)

da

 

 

 

 

 

 

 

Если множество (14.8) состоит из одного вектора, то этот вектор Vs {Т, а) не зависит от значения а. При этом

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