
ТУ - лекции Овсянникова / 32.Управление пучком траекторий в линейном случае
.pdfЛекция 32. Управление пучком траекторий в линейном случае.
Рассмотрим функционал
Z T Z Z
I(u) = '(t; xt; u(t; xt)) (t; xt)dxtdt + g(xT ) (T; xT )dxT ; (1)
0 Mt;u Mt;u
заданный на сечениях пучка траекторий системы, исходящих из множества M0 Rn, и учитывающий плотность распределения траекторий. Считаем, что начальная плотность (0; x) = 0(x) задана. Далее полагаем
g= x 0x + x 0;
'= x A(t)x + u B (t)x + x B(t)u + u C(t)u + x A1(t) + u C1(t);
где A; B; C; 0 матрицы, A1; C1; 0 векторы соответствующих размерностей. Предполагаем также, что 0 постоянная матрица, 0 постоянный вектор, а компоненты остальных упомянутых матриц и векторов
156
определены и непрерывны на [0; T ]. Считаем также, что матрицы A и C симметричны и квадратичная форма u Cu является положительно определенной, т.е. существует постоянная c > 0, такая, что cu u u C(t)u, t 2 [0; T ]. Требуется построить управление вида u = Mx + N(t), минимизирующее функционал (1).
Функционал (1) можно представить в виде
I(u) = ZM0 |
J(u; x0) 0(x0)dx0; |
(2) |
где |
|
|
J(u; x0) = Z0 T '(t; x(t; x0); u(t; x(t; x0)))dt + g(x(T; x0)): |
(3) |
Управление вида u = Mx + N(t), где матрица M(t) и вектор N(t) имеют вещественные и непрерывные компоненты, определенные на [0; T ], называют оптимальным по отношению к функционалу (3), если оно доставляет этому функционалу наименьшее возможное значение на любом движении
157
x = x(t; x0; u), x = x0 при t = 0.
Введем матрицу (t) и вектор (t), удовлетворяющие уравнениям
_ = (P QC 1B ) (P QC 1B ) + + QC 1Q A + BC 1B ;
_ = QC 1( Q + C1 ) 2 f P + (Q + B)C 1Q + +BC 1( Q + C1 ) ( P + B)C 1( Q + C1 )
A1 + ( Q + B)C 1C1
при конечных условиях
(4)
(5)
= 0; |
t = T; |
(6) |
= 0; |
t = T; |
(7) |
Имеет место следующая теорема В.И.Зубова:
158
Теорема 1 Для того, чтобы существовало оптимальное управление в системе dx=dt = P (t) x + Q(t) u + f(t) в смысле функционала (3), необходимо и достаточно, чтобы уравнение (4) с начальным условием (6) имело решение, непрерывное при t 2 [0; T ], при этом оптимальное управление представимо в форме
u0 = M0x + N0; M0 = C 1(Q + B ); |
1 |
C 1( Q + C1 ) ; (8) |
|
N0 = |
|
||
2 |
|||
где вектор (t) удовлетворяет уравнению |
(5) с начальным условием (7). |
Из теоремы 1 и представления (2) для функционала (1) следует, что управление, оптимальное по отношению к функционалу (3), является оптимальным и по отношению к функционалу (1).
Дадим еще одно представление функционала (1), позволяющее обойтись без интегрирования по сечениям пучка траекторий в функционале (1) или без интегрирования по множеству M0 в функционале (2) и сводящее бесконечномерную задачу к конечномерной.
159

Введем следующие характеристики пучка: |
|
|||
|
|
= ZMt;u (t; xt)dxt = const; |
|
|
|
|
(9) |
||
x(t) = ZMt;u xt (t; xt)dxt; |
||||
1 |
|
|
Z
D(t) = (xt x(t)) (xt x(t)) (t; xt)dxt:
Mt;u
Здесь масса (заряд) частиц; x(t) центр тяжести пучка;
D(t) матрица вторых моментов. Изменение x(t) и D(t) в силу системы dx=dt = P (t) x + Q(t) u + f(t) при управлении u = Mx + N описывается уравнениями
_ |
|
|
|
(10) |
|
|
|
||
x = (P + QM)x + QN + f; |
||||
D_ = (P + QM)D + D(P + QM) |
(11) |
160

с начальными условиями x(0) = x0, D(0) = D0. Функционал (1) можно записать теперь в виде
I(u) = Z0 T Spf(A + M B + BM + M CM)Dgdt+ |
|
|||||||||||||||||||
+ |
|
Z0 T ( |
|
Ax |
+ u B |
|
+ |
|
Bu + u Cu + |
|
A1 + u C1)dt+ |
(12) |
||||||||
|
x |
x |
x |
x |
|
|||||||||||||||
+ Sp( 0D(T )) + |
|
( |
|
(T ) 0 |
|
(T ) + |
|
(T ) 0(T )); |
|
|||||||||||
|
x |
x |
x |
|
где u = Mx(t) + N(t), x(t) и D(t) определяются уравнениями (10), (11). Таким образом, управление, оптимальное по отношению к функционалу
(3), является оптимальным и по отношению к функционалу (12). Особый интерес здесь представляет случай, когда x(t) 0.
161

Замечание 1 Пусть sii(t) диагональные элементы матрицы G 1 влетворяющей уравнению
|
dG 1 |
|
|
G 1(t) = Y (t)G0 1Y (t); |
|
= P G 1 |
+ G 1P : |
|
|||
|
dt |
|
. Нетрудно показать, что sii(t) равны квадратам максимальных нений xi xi(t) на эллипсоиде
(x x(t)) G(t)(x x(t)) = 1;
где G(t) = (Y 1(t)) G0Y 1(t). , т.е.
p
sii(t) = max zi; z = (z1; z2; : : : ; zn) :
z Gz=1
(t), удо-
(13)
откло-
(14)
(15)
162