Линейные задачи оптимизации. Ч.2. Оптимальное управление линейными динамическими объектами. Лутманов С.В. - 133 стр.

UptoLike

Составители: 

Рубрика: 

4. ЛИНЕЙНАЯ ЗАДАЧА ТЕОРИИ ОПТИМАЛЬНОГО УПРАВЛЕНИЯ
КАК ФУНКЦИОНАЛЬНАЯ ПРОБЛЕМА МОМЕНТОВ
133
[]
[
]
[
]
[
]
[
]
[] [] [] []
[] [] [] []
[] [] [] []
[] []
[
]
[
]
[] []
[] []
[] []
11 12 13 14 13 14
21 22 23 24 23 24
31 32 33 34 33 34
41 42 43 44 43 44
,,,, ,,
00
,,,, ,,
00
,,,,
,,,, ,,
10
,,,, ,,
01
xt xt xt xt xt xt
xt xt xt xt xt xt
Xt Ht Xt
xt xt xt xt xt xt
xt xt xt xt xt xt
τ
τττ ττ
τ
τττ ττ
τττ
τ
τττ ττ
τ
τττ ττ
⎛⎞
⎛⎞
⎜⎟
⎜⎟
⎜⎟
⎜⎟
===
⎜⎟
⎜⎟
⎜⎟
⎜⎟
⎜⎟
⎜⎟
⎝⎠
⎝⎠
⎜⎟
⎜⎟
⎜⎟
[]
()
[
]
[]
[]
()
[
]
[]
[]
()
[
]
[]
[]
()
[
]
[]
[]
12 34
13 23 33 43
14 24 34 44
1, 1, 1, 1,
,,,,0,1
1, 1, 1, 1,
xxxx
hh hh
xxxx
ττττ
τττττ
ττττ
⎛⎞ ⎛⎞ ⎛⎞
====
⎜⎟ ⎜⎟ ⎜⎟
⎝⎠ ⎝⎠ ⎝⎠
,
[]
()
[]
()
1
1
0
, , , 1,2,3,4
i
ij
hhdij
ατττ
==
11 12 13 14
21 22 23 24
31 32 33 34
41 42 43 44
0.766436 0.303602 1.51071 0.779061
0.303602 0.557218 0.764328 1.12012
1.51071 0.764328 3.3364 1.89574
0.779061 1.12012 1.89574 2.6037
αααα
αααα
αααα
αααα
⎛⎞
⎛⎞
⎜⎟
⎜⎟
⎜⎟
⎜⎟
=
⎜⎟
⎜⎟
⎜⎟
⎜⎟
⎜⎟
⎜⎟
⎝⎠
⎝⎠
,
[]
1011
2022
3033
40
44
0.640532
0.491302
1, 0
1.61672
1.31002
T
T
T
T
xcx
xcx
X
xcx
x
cx
⎛⎞
⎛⎞
⎛⎞
⎜⎟
⎜⎟
⎜⎟
⎜⎟
⎜⎟
⎜⎟
=− =
⎜⎟
⎜⎟
⎜⎟
⎜⎟
⎜⎟
⎜⎟
⎜⎟
⎜⎟
⎜⎟
⎝⎠
⎝⎠
⎝⎠
.
Запишем систему линейных алгебраических уравнений (4)
11 1 12 2 13 3 14 4 1
21 1 22 2 23 3 24 4 2
31 1 32 2 33 3 34 4 3
41 1 42 2 43 3 44 4 4
,
,
,
.
c
c
c
c
α
ναναναν
αν αν αν αν
αν αν αν αν
αν αν αν αν
+
++=
+++=
+++=
+++=
Ее решением будут числа
0000
1234
0.982227, 0.833006, 0.790821, 0.579604
ννν
=− =− = = .
Тогда оптимальное управление определяется по следующей формуле:
(
)
[
]
()
[
]
(
)
[
]
(
)
[
]
(
)
[
]
1234
00 0 0 0
1234
,0,1Ut h h h h t
ντντντντ
=+++ .
Вычислим функционал (5) на оптимальном управлении
() () ()
1
1
2
000
0
, 0.999712IU U U d
τττ
⎡⎤
⎡⎤
⋅= =
⎢⎥
⎣⎦
⎣⎦
.
Заметим, что для оптимального управления из примера 2.7 функционал
(6) принимает значение
1.0 0.999712> . Такой результат является естествен-
                         4. ЛИНЕЙНАЯ ЗАДАЧА ТЕОРИИ ОПТИМАЛЬНОГО УПРАВЛЕНИЯ
                                КАК ФУНКЦИОНАЛЬНАЯ ПРОБЛЕМА МОМЕНТОВ
            ⎛ x11 [t ,τ ]      x12 [t ,τ ]     x13 [t ,τ ]       x14 [t ,τ ] ⎞                            ⎛0        0 ⎞ ⎛ x13 [t ,τ ]     x14 [t ,τ ] ⎞
            ⎜                                                                ⎟                            ⎜            ⎟ ⎜                             ⎟
               x [ t ,τ ]      x22 [t ,τ ]     x23 [t ,τ ]       x24 [t ,τ ] ⎟                               0      0 ⎟ ⎜ x23 [t ,τ ]     x24 [t ,τ ] ⎟
X [t ,τ ] = ⎜ 21                                                                , H [ t ,τ ] = X [ t ,τ ] ⎜             =
            ⎜ x31 [t ,τ ]      x32 [t ,τ ]     x33 [t ,τ ]       x34 [t ,τ ] ⎟                            ⎜1        0 ⎟ ⎜ x33 [t ,τ ]     x34 [t ,τ ] ⎟
            ⎜⎜                                                               ⎟⎟                           ⎜⎜           ⎟ ⎜                            ⎟
             ⎝ x41 [t ,τ ]     x42 [t ,τ ]     x43 [t ,τ ]       x44 [t ,τ ] ⎠                             ⎝0       1 ⎟⎠ ⎝⎜ x43 [t ,τ ]   x44 [t ,τ ] ⎟⎠


              ⎛ x [1,τ ] ⎞              ⎛ x23 [1,τ ] ⎞            ⎛ x33 [1,τ ] ⎞            ⎛ x43 [1,τ ] ⎞
  h[ ] (τ ) = ⎜ 13                                                                                       ⎟ , τ ∈ [ 0,1] ,
                                [ 2]                      [3]                       [ 4]
                           ⎟ , h (τ ) = ⎜            ⎟ , h (τ ) = ⎜            ⎟ , h (τ ) = ⎜
    1

              ⎝ x14 [1,τ ] ⎠            ⎝ x24 [1,τ ] ⎠            ⎝ x34 [1,τ ] ⎠            ⎝ x44 [1,τ ] ⎠
                                                    1
                                          α ij = ∫ h[1] (τ ) , h[i] (τ ) dτ , i, j = 1, 2,3, 4 ⇒
                                                    0


                   ⎛ α11       α12     α13      α14 ⎞ ⎛ 0.766436 0.303602 1.51071 0.779061⎞
                   ⎜
                   ⎜ α 21      α 22    α 23     α 24 ⎟⎟ ⎜⎜ 0.303602 0.557218 0.764328 1.12012 ⎟⎟
                                                       =                                         ,
                   ⎜ α 31      α 32    α 33     α 34 ⎟ ⎜ 1.51071 0.764328 3.3364      1.89574 ⎟
                   ⎜⎜                                 ⎟ ⎜                                      ⎟
                    ⎝ α 41     α 42    α 43     α 44 ⎟⎠ ⎜⎝ 0.779061 1.12012 1.89574    2.6037 ⎟⎠

                                              ⎛ c1 ⎞ ⎛ xT 1 ⎞               ⎛ x10 ⎞ ⎛ 0.640532 ⎞
                                              ⎜ ⎟ ⎜            ⎟            ⎜ ⎟ ⎜                ⎟
                                              ⎜ c2 ⎟ = ⎜ xT 2 ⎟ − X [1, 0] ⎜ x20 ⎟ = ⎜ 0.491302 ⎟ .
                                             ⎜ c3 ⎟ ⎜ xT 3 ⎟               ⎜ x30 ⎟ ⎜ 1.61672 ⎟
                                             ⎜⎜ ⎟⎟ ⎜⎜         ⎟⎟           ⎜⎜ ⎟⎟ ⎜⎜             ⎟⎟
                                              ⎝ c4 ⎠ ⎝ xT 4 ⎠               ⎝ x40 ⎠ ⎝ 1.31002 ⎠

       Запишем систему линейных алгебраических уравнений (4)
                                                    α11ν 1 + α12ν 2 + α13ν 3 + α14ν 4 = c1 ,
                                                    α 21ν 1 + α 22ν 2 + α 23ν 3 + α 24ν 4 = c2 ,
                                                    α 31ν 1 + α 32ν 2 + α 33ν 3 + α 34ν 4 = c3 ,
                                                    α 41ν 1 + α 42ν 2 + α 43ν 3 + α 44ν 4 = c4 .

Ее решением будут числа
                       ν 10 = −0.982227, ν 20 = −0.833006, ν 30 = 0.790821, ν 40 = 0.579604 .

Тогда оптимальное управление определяется по следующей формуле:
                             U 0 ( t ) = ν 10 h[ ] (τ ) +ν 20 h[                                               , t ∈ [ 0,1] .
                                                                   2]
                                                1
                                                                        (τ ) +ν 30 h[3] (τ ) +ν 40 h[4] (τ )
Вычислим функционал (5) на оптимальном управлении
                                                                                                1
                                                        ⎡    1
                                                                                             ⎤2
                                      I ⎡⎣U 0 ( ⋅) ⎤⎦ = ⎢    ∫   U 0 (τ ) , U 0 (τ )      dτ ⎥ = 0.999712 .
                                                        ⎣    0                               ⎦
         Заметим, что для оптимального управления из примера 2.7 функционал
(6) принимает значение 1.0 > 0.999712 . Такой результат является естествен-



                                                                             133