Решение задач оптимального управления с использованием математической системы MATLAB и пакета имитационного моделирования SIMULINK. Сивохин А.В - 129 стр.

UptoLike

129
()
= 1exp
2
max1
T
tt
uktx
&
. (5.40)
Теперь мы имеем две неизвестные
1
t и
2
t . Для их определения
применим метод стыковывания уравнений. В точке
1
tt
=
()
tx
1
, вычисленные
по формулам (5.16) и (5.18), должны совпадать. В этой точке значения
(
)
tx
1
&
,
вычисленные по формулам (5.17) и (5.19), также должны совпадать. Имеем
следующую систему уравнений
=
++=
+
.1expexp
exp1exp
12
max
1
maxmax
12
12max
1
max1max
T
tt
uk
T
t
ukuk
T
tt
TTttukx
T
t
Tuktuk
n
, (5.41)
В первом уравнении раскроем скобки и приведем подобные члены, а
правую и левую часть второго уравнения разделим на
max
uk , получим
=
+
+=
+
.02expexp
,exp2exp2
112
12
max2max
1
max1max
T
t
T
tt
T
tt
Tuktukx
T
t
Tuktuk
n
(5.42)
Перепишем первое уравнение в виде:
2max
112
max1max
2expexp2 tukx
T
t
T
tt
Tuktuk
n
+=
+
+
. (5.43)
Так как выражение в скобках равно 0, то из последнего уравнения
находим
max
12
2
uk
x
tt
n
=
. (5.44)
Второе уравнение сначала умножим на
T
t
1
exp , а вместо
2
t подставим
его найденное значение, получим
01
2
expexp2
max
1
1
=
T
uk
x
t
T
t
n
. (5.45)
Обозначим
y
T
t
=
1
exp ,
(
)
1>y . Тогда следует
01exp2
max
2
=
Tuk
x
yy
n
. (5.46)
                                                 ⎛    ⎛t −t ⎞ ⎞
                         x&1 (t ) = k ⋅ u max ⋅ ⎜⎜ exp⎜ 2 ⎟ − 1⎟⎟ .                                                  (5.40)
                                                 ⎝    ⎝ T ⎠ ⎠

     Теперь мы имеем две неизвестные t1 и t 2 . Для их определения
применим метод стыковывания уравнений. В точке t = t1 x1 (t ) , вычисленные
по формулам (5.16) и (5.18), должны совпадать. В этой точке значения x&1 (t ) ,
вычисленные по формулам (5.17) и (5.19), также должны совпадать. Имеем
следующую систему уравнений
⎧                                  ⎛     ⎛ t1 ⎞ ⎞                      ⎛                       ⎛ t 2 − t1 ⎞ ⎞
⎪k ⋅ u max ⋅ t1 + k ⋅ u max ⋅ T ⋅ ⎜⎜ exp⎜ − ⎟ − 1⎟⎟ = x n + k ⋅u max ⋅⎜⎜ t 2 − t1 + T − T ⋅ exp⎜          ⎟ ⎟⎟
⎪                                  ⎝     ⎝ T  ⎠   ⎠                    ⎝                       ⎝     T    ⎠⎠
⎨                                                                                                              , (5.41)
                               ⎛     t ⎞
⎪k ⋅ u − k ⋅ u ⋅ exp − 1 = k ⋅u ⋅⎜ exp 2 1 − 1⎟.  ⎛    ⎛ t  − t  ⎞      ⎞
⎪      max         max         ⎜       ⎟      max ⎜    ⎜         ⎟ ⎟
⎩                              ⎝     T ⎠          ⎝    ⎝   T     ⎠ ⎠

     В первом уравнении раскроем скобки и приведем подобные члены, а
правую и левую часть второго уравнения разделим на k ⋅ u max , получим
⎧                                      ⎛    ⎛ t1 ⎞ ⎞                                                 ⎛ t 2 − t1 ⎞
⎪2 ⋅ k ⋅ u max ⋅ t1 + k ⋅ u max ⋅ T ⋅ ⎜⎜ exp⎜ − ⎟ − 2 ⎟⎟ = xn + k ⋅ u max ⋅ t 2 − k ⋅ u max ⋅ T ⋅ exp⎜          ⎟,
⎪                                      ⎝    ⎝ T⎠ ⎠                                                   ⎝ T ⎠
⎨                                                                                                                    (5.42)
⎪exp⎛ t 2 − t1 ⎞ + exp⎛ − t1 ⎞ − 2 = 0.
⎪⎩ ⎜⎝ T ⎟⎠                 ⎜
                           ⎝ T⎠
                                   ⎟

        Перепишем первое уравнение в виде:
                                               ⎛    ⎛t −t ⎞      ⎛ t ⎞ ⎞
         2 ⋅ k ⋅ u max ⋅ t1 + k ⋅ u max ⋅ T ⋅ ⎜⎜ exp⎜ 2 1 ⎟ + exp⎜ − 1 ⎟ − 2 ⎟⎟ = xn + k ⋅ u max ⋅ t 2 .             (5.43)
                                               ⎝    ⎝ T ⎠        ⎝ T⎠ ⎠

     Так как выражение в скобках равно 0, то из последнего уравнения
находим
                                                              xn
                                      t 2 = 2 t1 −                  .                                                (5.44)
                                                          k ⋅ u max

        Второе уравнение сначала умножим на exp⎛⎜ 1 ⎞⎟ , а вместо t 2 подставим
                                                 t
                                                                                    ⎝T ⎠
его найденное значение, получим
                                                     ⎛          xn          ⎞
                                                     ⎜ 2t1 −                ⎟
                                          ⎛t ⎞               k ⋅ u max
                                     2 exp⎜ 1 ⎟ − exp⎜                      ⎟ −1 = 0 .                               (5.45)
                                          ⎝T ⎠       ⎜       T              ⎟
                                                     ⎜                      ⎟
                                                     ⎝                      ⎠

        Обозначим exp⎛⎜ 1 ⎞⎟ = y , ( y > 1) . Тогда следует
                       t
                                  ⎝T ⎠
                                                      ⎛        xn       ⎞
                                      2 y − y 2 ⋅ exp⎜⎜ −               ⎟⎟ − 1 = 0 .                                 (5.46)
                                                      ⎝   k ⋅ u max ⋅ T  ⎠


                                                              129