ТУ - лекции Овсянникова / 31.Линейная траспортировка пучков.Программное управление.Задачи синтеза
.pdfЛекция 31. Линейная транспортировка пучков. Программное управление. Задачи синтеза.
Программное управление. Рассмотрим управляемую систему
dx=dt = P (t) x + Q(t) u + f(t): |
(1) |
Будем считать, что элементы матриц P (t) (размерностью n n), Q(t) (размерностью n r) и компоненты вектора f(t) заданы при t 0, вещественны
инепрерывны, u r-мерный вектор управления. Пусть заданы векторы x0
иx1 начальное и конечное положение системы (1).
Управление u = u(t) будем называть программным, если решение системы (1), начинающееся при t = 0 в точке x = x0, попадает при t = T в точку x = x1 и при этом интеграл R0T u ( )u( ) d ограничен.
143
Обозначим через Y (t) фундаментальную матрицу решений уравнения
dx=dt = P (t) x
при начальном условии Y (0) = E, где E единичная матрица. Представим решение системы (1) с конечным условием x(t) = x1 по
формуле Коши:
x(t) = Y (t) |
Y 1(T )x1 + |
Tt |
Y 1fQ( )u( ) + f( )gd |
; |
отсюда |
|
Z |
|
|
1 x0 = Z0 T B( )u( )d ; |
|
|||
|
(2) |
где 1 = Y 1(T )x1 R0T Y 1( )f( ); B( ) = Y 1( )Q( ).
144
Следуя В.И.Зубову, решение уравнения (2) будем искать в виде
u = B c + v(t); |
(3) |
где c постоянный вектор, подлежащий определению, а v(t) произволь-
ная функция, суммируемая вместе со своим квадратом на [0; T ] |
и такая, |
|||||||
что |
0T B( )v( )d = 0. |
|
|
|
||||
|
R |
|
|
|
|
|
|
|
|
Подставляя выражение (3) в равенство (2), находим |
|
|
|
||||
|
|
|
|
|
1 x0 = A(T )c; |
|
|
(4) |
где |
A(T ) = |
|
T BB dt |
. Уравнение (4) разрешимо, если |
det(A(T )) = 0 |
либо |
||
|
|
|
0 |
|
6 |
|||
ранг |
матрицы A(T ) совпадает с рангом расширенной матрицы |
|
|
|||||
|
R |
|
|
|
|
|
||
A = |
fA(T ); 1 x0g. Таким образом, имеет место следующая теорема |
В.И.Зубова:
145
Теорема 1 Для того, чтобы существовало программное управление, переводящее систему (1) из любого положения в любое другое положение x0
R T
за время T , необходимо и достаточно, чтобы матрица A = 0
где B = Y 1Q, была неособой. При этом все множество программных управлений дается формулой
u = B c + v;
где v(t) суммируемая с квадратом на [0; T ] функция и
0 T |
Bvdt = 0; |
c = A 1(T ) Y 1(T )x1 |
0 T |
Y 1( )f( )d x0 : |
Z |
|
|
Z |
|
Динамика эллипсоидов. Пусть задан в пространстве Rn эллипсоид S0:
(x0 |
|
0) G0(x0 |
|
0) = 1: |
(5) |
x |
x |
Здесь x0 центр эллипсоида, а G0 симметричная положительноопределенная матрица размерностью n n, задающая эллипсоид.
146
Пусть x(t) = x(t; x0; u) решение системы (1) с начальным условием x(0) = x0 при управлении u = u(t). Очевидно, что любое другое решение x(t) = x(t; x0; u) системы (1) представимо в виде
x(t) = Y (t)(x0 |
|
0) + x(t): |
(6) |
x |
Сдвиг эллипсоида (5) вдоль решений (6) есть также эллипсоид с центром в точке x(t):
(x |
|
x(t)) G(t)(x |
|
x(t)) = 1; |
(7) |
|
|
|
где G(t) = (Y 1(t)) G0Y 1(t). При этом
|
dG |
= P G GP: |
|
(8) |
|||
|
|
|
|
||||
|
dt |
|
|||||
Отметим также, что |
|
|
|
|
|
||
G 1(t) = Y (t)G 1Y (t); |
dG 1 |
+ G 1P : |
(9) |
||||
|
= P G 1 |
||||||
0 |
|
|
|
dt |
|
||
|
|
|
|
|
|
147
Из матричного уравнения (8) следует в частности, что программное управление u = u(t) в системе (1) не влияет на матрицу G(t), определяющую в каждый момент t эллипсоид (7). Программным управлением можно изменять только положение центра эллипсоида, выбирая нужный закон изменения траектории x(t).
Синтез управления. Пусть заданы множества M0 и M1. Требуется построить управление
u = M(t)x + N(t) |
(10) |
так, чтобы оно переводило за время T все траектории, начинающиеся в множестве M0, в множество M1. Относительно такого управления будем говорить, что оно осуществляет транспортировку пучка траекторий из множества M0 в множество M1.
В работе В.И. Зубова дано построение таких управлений. Оно основано на решении задачи транспортировки n-мерного симплекса, включающего множество M0 в другой n-мерный симплекс, содержащий множество M1.
148
Следуя указанной работе, осуществим построение управления вида (10), обеспечивающего транспортировку пучка траекторий из множества M0 в множество M1. Предположим, что эллипсоид (5) охватывает множество M0. Рассмотрим эллипсоид S1:
(x |
|
1) G1(x |
|
1) = 1 |
(11) |
x |
x |
с центром в точке x1, задаваемый симметричной положительно-определенной матрицей G1 и лежащий в множестве M1. Если мы решим задачу транспортировки эллипсоида (5) в эллипсоид (11), то очевидно, что будет решена и исходная задача транспортировки множества M0 в множество M1.
Пусть S неособая матрица, определяющая преобразование эллипсоида (5) в эллипсоид (11): x = x1 + S(x0 x0). Столбцы матрицы S обозначим через S(j), j = 1; n. Построим программное управление u0(t), переводящее точку x0 в точку x1. Обозначим соответствующую этому управлению траекторию x0(t), при этом x0(0) = x0, x0(T ) = x1.
149
Построим также программные управления u(j) = u(j)(t), которым соответствуют траектории x(j)(t), удовлетворяющие условиям
x(j)(0) = e |
+ |
|
|
; |
x(j)(T ) = Sj + |
|
|
; |
x |
x |
|||||||
j |
0 |
|
1 |
|
где ej = (0; : : : ; 0; 1; 0; : : : ; 0) единичный орт.
Указанные программные управления вида u(j) = B (t)c(j) + v существу-
ют тогда и только тогда, когда совместны системы |
|
|||||
|
|
1(j) x0(j) = A(T )cj; |
j = 0; 1; : : : ; n; |
(12) |
||
|
x |
|||||
где |
|
|
|
|
||
|
A(T ) = Z0 |
T |
x0(j) = x(j)(0); |
|
||
|
BB d ; |
|
||||
|
|
x1(j) = Y 1(T )x(j)(T ) Z0 |
T |
|
||
|
|
Y 1fd : |
|
150
Построим матрицы
F (t) = x(1)(t) x(0)(t); x(2)(t) x(0)(t); : : : ; x(n)(t) x(0)(t) ;
U(T ) = u(1)(t) u(0)(t); u(2)(t) u(0)(t); : : : ; u(n)(t) u(0)(t) :
Из построения матриц F и U видно, что
dF (t) = P (t)F (t) dt
+ Q(t)U(t): |
(13) |
Предположим, что F (t) неособая матрица при t 2 [0; T ]. Тогда, если в управлении (10) положим M = UF 1, а N = u(0) Mx(0), то эллипс (5) в силу системы (1) при данном управлении перейдет в момент T в эллипс (11). Исследуем вопрос о существовании обратной матрицы F 1(t), t 2 [0; T ]. Умножая (13) на Y 1 и интегрируя в пределах от 0 до t, получаем равенство
Z0 t Y 1F_ d = Z0 t Y 1P F d + Z0 t Y 1QUd : |
(14) |
151
Интегрируя слева по частям и пользуясь соотношением dY 1=dt = Y 1P , из равенства (14) найдем
Z t
Y 1(t)F (t) = F (0) + Y 1QUd : (15)
0
Учитывая, что u(j) = B c(j) + v, выпишем j-й столбец подынтегральной матрицы в равенстве (15):
fY 1QUgj = Y 1QB (c(j) c(0)); j = 1; n;
отсюда
Z0 t |
Y 1QUd |
j |
= A(t)(C(j) C(0)); j = 1; n; |
||
n |
|
o |
|
|
|
Z t
A(t) = BB d :
0
152