⎢⎣t0 ⎥⎦ ⎢⎣t0 ⎥⎦
Второй этап процедуры, состоящий в построении «минимального» элемента h0 ()⋅ , сводится к задаче математического программирования следующего ви-
да:
t0i=1 i=1i j, =1 ⎢⎣t0 ⎥⎦
Или по-другому:
n n
i∑, j=1αijlil j → min, ∑i=1 cili =1,
где
T
t0
⎛l10 ⎞
⎜ ⎟
Пусть вектор l 0 = ⎜ ⎟ – решение этой задачи. Тогда «минимальный» элемент
⎜⎝ln0 ⎟⎠
вычисляется по формуле
n h
i=1
Заключительный третий этап построения оптимального управления снова сводится к решению изопериметрической задачи вариационного исчисления
0 ρ
Ее решением является вектор-функция
0 1 0 1
U ( )⋅ = −
где
µ ⎢∫
⎢⎣t0 ⎥⎦
Таким образом, оптимальное управление, решающее задачу теории оптимального управления, имеет вид
U
(ρ0 ) (ρ0 ) i=1 i=1
0 li0
где νi =
(ρ0 )
0 ⎡T 0 0⎤1
Воспользовавшись соотношениями (3), укажем другой способ формирования оптимального управления, который быстрее и проще приводит к цели. Подставив выражение (3) в равенство (1.3), получим
t0=1j=1 t0 j=1 j
Таким образом, решение задачи об оптимальном управлении можно начинать непосредственно с решения системы линейных алгебраических уравнений
(4), а оптимальное программное управление вычислять по формуле (3).
Итак, оптимальное управление для случая минимизации «энергии» единственно, непрерывно в каждой точке интервала [t0 ,T], пропорционально «минимальному» элементу h0 ()⋅ и зависит линейно от краевых условий. Пример 3*. Рассмотрим линейный управляемый динамический объект
x1 = 2x1 + 2x2 −30x3 +u1, x2 =10x1 − x2 −35x3 +u2, x3 = 2x1 − x2 + x3 +u3, t0 = 0,T =1;
x10 = −3, x20 = 2, x30 =1; xT1 = −80.7746, xT 2 = −147,179, xT3 = −8.94415;
I U⎣⎡ ( )⎦⎤ u ( ) u ( ) u ( ) . (5)
⎣ 0 ⎦
В рассматриваемом примере дифференциальные уравнения движения динамического объекта и начальное положение фазового вектора взяты из примера 2.3. Конечное положение фазового вектора совпадает с той точкой фазового пространства, в которую переводит в конечный момент времени фазовый вектор оптимальное управление из примера 2.3.
Последовательно вычисляем
⎛ x11 [t,τ] x12 [t,τ] x13 [t,τ]⎞
X t[ ,τ] = ⎜⎜ x21 [t,τ] x22 [t,τ] x23 [t,τ]⎟⎟, H t[ ,τ] = X t[ ,τ]B = X t[ ,τ],
⎜⎝ x31 [t,τ] x32 [t,τ] x33 [t,τ]⎟⎠
⎛ x11 [1,τ]⎞ ⎛ x21 [1,τ]⎞ ⎛ x31 [1,τ]⎞
h[ ]1 ( )τ = ⎜⎜ x12 [1,τ]⎟⎟ , h[ ]2 ( )τ = ⎜⎜ x22 [1,τ]⎟⎟ , h[ ]3 ( )τ = ⎜⎜ x32 [1,τ τ]⎟⎟, ∈[0,1],
⎜⎝ x13 [1,τ]⎠⎟ ⎝⎜ x23 [1,τ]⎟⎠ ⎝⎜ x33 [1,τ]⎠⎟
1 ⎛α α α11 12 13 ⎞ ⎛4770.34 8405.45 469.812⎞
0 ⎜⎝α α α31 32 33 ⎟⎠ ⎜⎝469.812 832.786 48.3062⎟⎠
⎛c1 ⎞ ⎛ xT1 ⎞ ⎛ x10 ⎞ ⎛−67.4743⎞ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟
⎜c2 ⎟ = ⎜ xT 2 ⎟− X [1,0]⎜ x20 ⎟ = −⎜ 115.885⎟.
⎝⎜c3 ⎠⎟ ⎝⎜ xT3 ⎟⎠ ⎜⎝ x30 ⎟⎠ ⎜⎝−5.34546⎠⎟
Запишем систему линейных алгебраических уравнений (4)
α11 1ν αν αν+ 12 2 + 13 3 = c1, αν αν αν21 1 + 22 2 + 23 3 = c2, . αν αν αν31 1 + 32 2 + 33 3 = c3
Ее решением будут числа
ν10 = 0.0682083, ν20 = −0.0954224, ν30 = 0.871024 .
Тогда оптимальное управление определяется по следующей формуле:
U
Вычислим функционал (5) на оптимальном управлении
0 ⎡ 1 0 0⎤
⎣ 0 ⎦
Проверим проведенные вычисления. Покажем, что полученное управление переводит фазовый вектор из положения x0 в положение xT за время [0,1]. Действительно проинтегрируем систему дифференциальных уравнений