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

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

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

Закон управления (21.4.2) скорее всего окажется неоптималь­ ным по интересующему нас критерию, но привлекательна его про­ стота. Теоретически управление u ( t , х) не становится бесконеч­ ным вдоль любой траектории, но при практической реализации за­ кона управления (21.4.2) при I, близких к матрица Ж (/, t L)

будет близка к вырожденной, и управление необходимо прекратить

в момент /, достаточно близкий к

При этом система попадет в

малую окрестность точки Xj.

 

§21.5. Представление критерия управляемости через матрицу управляемости

Использование критерия управляемости, даваемого теоремой 21.3.3, требует знания фундаментальной матрицы Коши. Заманчи­ во получить критерий управляемости, выраженный через исходные матрицы системы A ( t ) и B(t). Оказывается, это возможно в слу­

чае, когда матрица A( t)

(п — 2) раза, а матрица B ( t )

(л — 1) раз

дифференцируемы почти всюду на интервале [/0, *,].

 

 

 

Л е м м а

21.5.1. Если

матрица G(t) размера п х т

дифферен­

цируема п — 1 раз почт и всюду на интервале

[/0,

/ J ,

то для не­

нулевого вектора с равенство

 

 

 

 

 

 

 

 

 

 

 

 

 

 

стС (/) =

0

 

 

 

(21.5.1)

имеет мест о почт и всюду на [/0, /J

тогда и только тогда,

когда

почт и всюду на этом интервале

 

 

 

 

 

 

 

 

 

 

 

 

 

*3

II

0.

 

 

 

(21.5.2)

Здесь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6 ( 0

=

(Г ( Л

dGit)

dn~xG(t)\

 

(21.5.3)

 

 

 

 

 

G W

dt

dtn~ l

/J

 

 

 

 

Д о к а з а т е л ь с т в о .

 

 

 

 

 

 

 

 

 

 

интервале

U0, / J

матрица

 

1.

Пусть

почти

всюду

на

G{t) =

0, с ^

0.

Дифференцируя

равенство

(21.5.1)

no

t

(я — 1)

раз, получаем

 

 

 

 

 

 

 

 

 

 

crG (t) = 0,

 

с

т dG(t) __

 

dn~ lGU)

= 0

 

 

 

 

dt

= 0, . . . , с

п- 1

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

почти

всюду

на

интервале

[f0, /, J. Объединяя последние равенства

в одно матричное и вынося влево вектор-строку ст, имеем

 

 

 

 

 

ст

G (0

dCit)

dn~xG{t)

= 0

 

 

 

 

 

 

 

 

 

 

dt

 

dt л - 1

У

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

crQ ( t ) = 0

почти

всюду на

интервале

[70, /,].

 

 

 

2. Обратное утверждение очевидно, так как равенство (21.5.1) есть часть равенства (21.5.2). ■

Если теперь рассматривать матрицу G(t) в лемме 21.5.1 как матрицу, задаваемую равенством (21.3,11), то можно выразить матрицу Q(t) через матрицы A(t) и

Л е м м а

21.5.2. Если матрица A(t) (п — 2) раза, а матрица

B(t) ( л — 1)

раз дифференцируемы почти всюду на интервале

[/0, 7,], то матрица Q(t) почти всюду на этом интервале может быть представлена в виде

Q(t) = X - l(t,r)Q(t),

(21.5.4)

Q(t) = (B(t)AiBA(t)...

ДГ 'Я (О ),

(21.5.5)

где Дл — дифференциальный оператор

 

 

Дл = - А ( 1 ) + £ .

(21.5.6)

Ранги матриц Q{t) и Q(t) совпадают.

 

 

Д о к а з а т е л ь с т в о .

Вычислим

производную

матрицы

t*)S(l), где 5(0 — матрица размера п х т ,

 

£ (x -'(t.n s (i)) =

5(0+ *-*(».

=

= - 2 Г ‘0 , l‘)A(t)S(t) + X - '( t, f )

= X ~ 1{U f ) b AS(t).

Итак,

 

 

 

£ ( X - ‘(Г, 0 5 ( 0 ) = X - >0. <-)ДлД(0-

(21.5.7)

Полагая в (21.5.7) 5(0 = B(i), получаем

 

Т Г

= x - ‘(t, ОДЛВ(0.

 

Аналогично, полагая S(i) — AAB(t), получаем

 

d2CU)dt2

= х - ‘о , п а 2а*(0-

 

Продолжая этот процесс, имеем

 

 

 

=х->(1, П&"->В(П.

 

Тогда матрица Q(t), определяемая равенством (21.5.3), принимает

ВИД <2(0 = *•'('. OC(Ol

ст ема (21.1.2) вполне управляема на [f0,*,], если хот я бы в одной

т очке эт ого инт ервала матрица Q(t) имеет ранг п.

Д о к а з а т е л ь с т в о . И з условий теоремы и леммы 21.5.2 сле­ дует, что матрица W(tQitv) имеет ранг л, а значит, положительно определена. Тогда по теореме 21,3.3 система (21.3.1) управля­ ема. ■

Теорема 21.5.1 дает достаточный критерий управляемости на интервале

§21.6. Критерий управляемости для линейных стационарных систем

Если система (21.1.2) стационарна, то матрица управляемости Q{t) приобретает простой вид. Условия теоремы 16,5.1 в этом слу­ чае и необходимы. Кроме того, для линейной стационарной систе­ мы понятия вполне управляемой на интервале [f0, / J , управляемой относительно момента tQи абсолютной управляемой системы совпа­ дают.

Л е м м а 21.6.1. Д л я линейной стационарной системы

4Ц- =

А х + Ви,

х G

R n,

u e R m

(21.6.1)

мат рица управляемост и Q(t) постоянна:

 

 

Q(t) ~

Q — ( В А В ...

{ - \ ) п - ' А п~ 1В)

(21.6.2)

и ранги матрицы W { t ^ t x)-u Q совпадают при лю бы х значениях t0,

> t0.

До к а з а т е л ь с т в о .

1.Из леммы 21.5.2 следует, что ранг Q не превышает ранга W .

2.

Пусть матрицы W ( t 0, *,)

имеет ранг

к ^ п . Докажем,

что

ранг

Q не может быть меньше

к, Допустим

противное, что

ранг

Q равен s < п. В этом случае существует п — s линейно независи­

мых векторов ср с2,

с „ _ А таких, что

 

 

 

 

c]Q = (с]Вс]АВ

с]Ап ~ 1В) = 0

(/ =

1,

2, ..., п -

s)

или

 

 

 

 

 

 

 

с]В = с]АВ =

= с]Ап - ' В = 0

( / =

1, 2,

...,

n - s ) .

(21.6.3)

По теореме Гамильтона—Кэли матрица удовлетворяет своему ха­ рактеристическому уравнению

Xй + а1Х'*-1 + + а п- 1Ь + ап = О,

т.е.

An + axAn~l + + ап_ 1А + апЕп —0.

(21.6.4)

Умножая последнее равенство справа на В и слева на с], получаем

с}АпВ + ахс]Ап ~1В + ... + ап _ хс}АВ + апс]В = 0.

(21.6.5)

Из этого равенства и (21.6.3) следует

с}ЛпВ = 0

( / = 1 , 2 , . . . , n - s ) .

(21.6.6)

Если предварительно умножить равенство (21.6.4) на А и по­ вторить процедуру, то будем иметь

с]Ап ~ 1В + ахс)АпВ + + ап _ хс\А2В + апс}АВ = 0,

откуда с учетом (21.6.3) и (21.6.6) с]Ап~{В = 0 и т. д. Таким об­ разом,

с]ЛкВ = 0

(* = 0,1,2,...; / = 1, 2, .... п - s).

Тогда для матрицы G(t), определяемой равенством (21.3.11), полу­ пим, что при всех { и / = 1,2, ..., л — s

c}G(t) = с}<ГА<‘- '’>/?= с]\Еп + А(Г - 0 + А (t2~ t}- + ...

+

+ . . . ] B = c ] B + c 1 A B ( f - 0 +

 

+ с}А2В ^ ^ ~ + ...+ с}А кВ ^ р - + ... = 0 .

По лемме 21.6.1 матрица W(tQ,tx) имеет ранг s< t.

Это противоречит предположению, что ранг W равен к. Остает­ ся предположить, что ранг Q также равен к.

3. Так как матрица Q не зависит от моментов t0, tv, то ее ранг

постоянен, а тогда и ранг матрицы W постоянен и не зависит от моментов начала и конца движения. ■

Из леммы 21.6.1 и

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

Т е о р е м а 21.6.1.

Для того чтобы стационарная система

(21.6.1) была вполне управляемой на любом интервале времени

[f0, /J , т.е.

абсолютно управляемой,

необходимо и достаточ­

но,

чтобы

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

Q= (ВАВ...Ап~ 1В) имела

ранг

п.

 

 

Д о к а з а т е л ь с т в о . Оно.следует из теоремы 21.3.3 и леммы 21.6.1, если учесть, что матрица Qc отличается от матрицы Q

(21.6.2) знаками столбцов, что не изменяет ее ранга. ■

Из теоремы 21.6.1 вытекает, что для стационарных систем по­ нятия вполне управляемости на нитервале и абсолютной управля­ емости совпадают.

Далее будет показана эквивалентность этих определений управ­ ляемости определению 21.1.2.

Пр им ер 21.6.1. Линеаризованные уравнения движения летательного аппарата в окрестности круговой траектории в предположении, что возмущенные траектории остаются в плоскости круговой траектории, как показано в примере 20.4.1 на с. 411т имеют вид

d*vr

— — = 2vA Vx+ v^Ar + s,

dt т

Физический смысл переменных ДVr%AVv Дг, Т и постоянной v был также ука­ зан в примере 20.4.1. Проекции 5 и Т силы тяги двигателя Р на радиус и нормаль

крадиусу являются управлениями.

Вматричной записи уравнения движения летательного аппарата будут иметь

вид

d_

(AVry

' 0

2v

v2'

'д к ;

 

<1

=

-v

0

0

д ^

 

 

dt

 

 

V

У

_

О

■v.

Дг

/

 

 

Дг

 

 

 

О

 

 

1. Бели T(t)mо, т е. тяга двигателя всегда направлена по направлению радиуса

г, то матрица управляемости

 

 

 

 

 

 

 

 

 

 

 

 

(\ 0

-v 2)

 

Q = ( B А В А2 В) = 0 V 0

 

 

 

 

 

 

\0 1

0 /

и ее определитель К?£.|= 0 ,

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

Tit) жо. Из этого следует, что летательный аппарат в общем случае нельзя вернуть на круговую траекторию одним лишь радиальным ускорением.

2. Если 5Ц )"0, т.е. тяга двигателя направлена по нормали к радиусу, то мат­

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

 

 

/0 2v

O '

Qc

1 0

—2V2

\

1 0

2v t

и ее определитель \Qc\ = — 4у2^ 0 . Система является вполне управляемой, т.е. ле-

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

§ 21.7. Структура пространства состояний линейной управляемой системы

Изучим сначала структуру пространства состояний для преобра­ зованной системы (21.2.3).

Т е о р е м а 21.7.1. Если матрица W(tQ) tx) имеет ранг к, то

ортогональным преобразованием у = А*Гу система (21.2.3) может

быть приведена к виду

 

 

dA

= о

(21.7.1)

dt

и’

 

dh

G2(t)u,

(21.7.2)

dt

 

 

(уЛ

 

вектор у2

где у — - , вектор у{ имеет размерность п — к,

V 7

размерность к, а система (21.7.2)вполне управляема на интерва­ ле [/0, /,]. Очевидно, что система (21.7.1) вообще неуправляема.

До к а з а т е л ь с т в о .

1.При к = 0 теорема верна, так как в этом случае равенства (21.3.8) возможны только, если <7 = 0, и тогда исходная система

имеет вид (21.7.1), а система (21.7.2) отсутствует.

2.

Будем считать, что к > 0. Совершенно ясно, что векторы

с1>с2»

c« - t в лемме 21.3.3 могут быть выбраны ортонормиро-

ванными и порождают (п — к)-мерное подпространство Ry~к про­ странства Rn. Образуем матрицу М — (Л/, М2), где столбцы мат­ рицы М, составлены из векторов ср с2, ..., сп_ а столбцы матри­

цы Мг представляют собой орты подпространства Rk,

дополняющего Rny~k до Rn. Тогда можно положить у — Af*y или

 

r)V

'Щу\

 

 

 

А )

АР2у

 

 

 

 

к j

 

 

Действительно, из (21.3.8) следует

 

 

ТГ = ^

ш 0.

 

 

О"-

3. Если бы система (21.7.2), т.е. система

 

 

 

7 Г - * * “ %■ =Щ О(О »,

 

не была вполне управляема на

интервале

[*0, <,],

то по лемме

21.3.3 существовал бы ненулевой вектор с

такой, что

стЛ^<7(0 = 0 почти всюду на интервале [/0, /,], или

 

 

с"п- к +1^(0

Cn-Jfc + l

М2с.

(21.7.3)

Вектор сп_ к+, представляет собой линейную комбинацию столбцов матрицы М2 и поэтому ортогонален всем векторам с,, с2, ...» сп_ к.

А тогда из равенства (21.3.8) и (21.7.3) вытекает, что матрица G( t) имеет меньше, чем к линейно независимых строк, и поэтому

по лемме 21.3.3 матрица

W ( t 0, t x) имеет ранг меньше чем к,

что

противоречит предположению теоремы. ■

 

 

Интересно

проследить изменение

ранга матрицы W ( t Q, г,)

при

изменении ее аргументов.

 

 

 

Т е о р е м а

21.7.2. Ранг матрицы

W ( t 0, /,)

не убывает при уве­

личении fj

(уменьш ении

t0) и не возрастает

при уменьш ении /,

(увеличении

t0).

 

 

 

 

 

Д о к а з а т е л ь с т в о .

Оно легко следует из того, что функции,

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

на интервале [/0, t {), линейно зависимы

и на

любом подынтервале

[t'Q,

t\ ] С [г0> /,],

а также

и из следствия лем ­

мы 21.3.3. ■ Так как ранг — целочисленная переменная, то, например, при

увеличении tx ранг W (t0, tx) в некоторые моменты времени может

измениться скачком, так что зависимость ранга W от 10 представ­

ляет собой ступенчатую функцию .

Из доказанной выше теоремы 21.7.1 видно, что пространство R n

векторов у расслаивается на подпространство А" ” * и А*. Подсисте­

ма, леж ащ ая в А”, вполне управляема, а леж ащ ая в А " - * неуправ­

ляема. Если сделать обратное преобразование переменных

у = Х ~ *(<, Г)х, то подпространство А” ~ к перейдет такж е в (п — к )-

мерное подпространство R xn ~ k(t) С А”, но зависящее от t. Аналогич­

но подпространство А* перейдет в подпространство А^(/).

Так как

у х = Щ у = M \ X ~ l(t, Г ) х ,

уг = Щ у = Щ Х ~ \ и

0 * 1

то многообразие L k точек ( t , х ) G R n + [, которые описывает под­ пространство R * ~ k(t) при изменении времени /, задается равенст­ вом

Щ Х - ' Ц , О * = 0,

(21.7.4)

а многообразие L k точек (/, х) £ Ал + 1, описываемых подпростран­ ством А *(г), равенством

M \ X ~ l ( t , O x = 0.

(21.7.5)

Следствием теоремы 21.7.1 для исходной системы (21.1.1) явля­ ется следующая теорема.

Т е о р е м а 21.7.3. Пусть матрица W{tQ, /,) имеет ранг к< п.

Тогда существует единственное многообразие L C Rn + l размер­ ности к, такое, что:

а) многообразие Lk инвариантное многообразие для системы (21.1.1), т.е. каждая траектория, начинающаяся в любой момент t G [/0, /, J в этом многообразии, не выходит из него, и каждая тра­

ектория, начинающаяся в любой момент т G [/0, /х] вне Lk, не мо­

жет попасть в него;

только в этом мно­

б) если рассматривать систему (21.1.1)

гообразии, то она управляема на интервале

т.е. из любой

точки дс0 такой, что (/0, л'0) G Lk, можно перевести систему за время — /0 любую точку x v (t, V) G ZA

Справедливость теоремы следует их теоремы 21.3.3 и взаимной однозначности преобразования у= МТХ~1(1, 1*)х.

Заметим, что на меньшем интервале [xg, xj G [*0, i{] система

может оказаться неуправляемой даже в многообразии Lk, так как ранг матрицы 1К(т0, х,) может оказаться меньше к. Если же матри­

ца W{т0, х,) имеет ранг к для любых х0, xt G [/0, /,](х0 < т,), то си­

стема (21.1.1) будет абсолютно управляемой в многообразии ZA

§21.8. Структура пространства состояния стационарной системы

Весьма интересно, что для стационарных линейных систем мно­ гообразия Lk и Ln~k являются цилиндрическими по t, т.е. подпро­ странства Rk(t), не зависят от t. Для стационарной линей­ ной системы

a t = Ах + Ви, х G R '1,

и G R m,

(21.8.1)

 

подпространство R*~k задается равенством

 

д^те-А(< - Ох = 0,

 

(21.8.2)

а подпространство R k — равенством

A/re-A (/-ox —о,

(21,8 3)

получающимися из (21.7.4), (21.7.5) подстановкой выражения для фундаментальной матрицы линейной системы X{t, t*) — еЛ^~^К Из вида равенств (21.8.2), (21.8.3) далеко не очевидно, что подпро-

странства Я” R* не зависят от t, так как матричный коэфф ици­

ент при векторе JC является функцией времени. Однако имеет место следующая теорема, являющаяся конкретизацией теоремы 21.7.1 для стационарных систем.

Т е о р е м а 21.8.1! Пусть матрица Qc = ( B A R . . А п ~ хВ ) име­

ет ранг к К п. Тогда существует единственное линейное подпро­ ст ранст во R k С R ’1 такое, что:

а) подпространство R * порождается системой векторов, со­

ст авленной из к линейно независимых столбцов мат рицы Qc;

б) R k

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

подпространство, т.е.

т раект ория

системы

(21.8.1), начинаю щ аяся

в

точках,

принадлеж ащих

R k,

целиком лежит

в R k, и никакая траектория, начинаю щ аяся

вне

Л*, не может привести в /?*;

 

 

 

 

 

 

 

в) если систему (21.8.1)

рассматривать т олько в

R k, т о она

абсолю т но управляема.

При к = п теорема

 

 

 

 

 

Д о к а з а т е л ь с т в о .

очевидно

верна,

поэ­

тому будем считать, что к < п. Так

как по условию теоремы

ранг

Qc равен

к < п,

то по лемме

21.6.1

ранг матрицы

W (tQ, ^ ) также

равен к.

Образуем матрицу

М =

( М { М 2) точно так

же,

как

при

доказательстве

теоремы

21.7.1.

Покажем,

что

п — к

векторов

Pk + v

Рк. + г>

образующих строки матрицы М \е ^ ~ г ^ в равен­

стве

(21.8.3), определяющем подпространство

R k,

в любой момент

времени t ортогональны системе к линейно независимых постоян­ ных векторов qv q2, ..., qk, представляющих собой к линейно неза­

висимых столбцов матрицы Q.

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

 

 

 

M]G(t) = 0,

(21.8.4)

а по лемме 21.5.1, применяемой к каждой строке

матрицы М \ в

(21.8.4),

такж е и

 

 

 

 

M \ Q ( t ) = 0.

(21.8.5)

М атрица

Q связана с матрицей Q, которая для стационарной сис­

темы имеет вид Q — ( В — А В ... ( \ ) п ~ хА п ~ хВ),

соотношением

(см. лемму 21.5.2)

Q ( t ) = X ~ l(t, t') Q или, учитывая, что для ста­

ционарной системы

X ~ l( t , /*) = е~л(-1~ 1\ — соотношением

(21.8.6)

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