книги / Матричное исчисление с приложениями в теории динамических систем
..pdfЗакон управления (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, /,]. |
|
|
|
ст ема (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) далеко не очевидно, что подпро-