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

UptoLike

ВУЗ: 

ПГНИУ | Пермь

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

Рубрика: 

ПРИЛОЖЕНИЕ
184
Проверка попадания в начало координат
8x1@TD,x2@TD<
8−0.0000134271, −0.0000148833<
Пример 4.3
Ввод начального и конечного положений фазового вектора
X0 = 8−3, 2, 1<;
XT
= 8−80.7746, −14 7.179, −8 .94415<
8−80.7746, − 147.179, − 8.94415<
Построение фундаментальной матрицы Коши
Resh1 =
DSolve@8x11'@tD 2 ∗ x11@tD+ 2 ∗ x21@tD− 30 ∗ x31@tD,
x21'
@tD 10 ∗ x11@tD− x21@tD− 35 ∗ x31@tD,
x31'
@tD 2 ∗ x11@tD− x21@tD+ x31@tD,x11@0D 1,
x21
@0D 0, x31@0D == 0<, 8x11@tD,x21@tD,x31@tD<,
t
D;
Resh2
=
DSolve@8x12'@tD 2 ∗ x12@tD+ 2 ∗ x22@tD− 30 ∗ x32@tD,
x22'
@tD 10 ∗ x12@tD− x22@tD− 35 ∗ x32@tD,
x32'
@tD 2 ∗ x12@tD− x22@tD+ x32@tD,x12@0D 0,
x22
@0D 1, x32@0D == 0<, 8x12@tD,x22@tD,x32@tD<,
t
D;
Resh3
=
DSolve@8x13'@tD 2 ∗ x13@tD+ 2 ∗ x23@tD− 30 ∗ x33@tD,
x23'
@tD 10 ∗ x13@tD− x23@tD− 35 ∗ x33@tD,
x33'
@tD 2 ∗ x13@tD− x23@tD+ x33@tD,x13@0D 0,
x23
@0D 0, x33@0D 1<, 8x13@t D,x23@tD,x33@tD<,
t
D; 88x11@t, τD<, 8x21@t, τD<, 8x31@t, τD<< =
8Re@x11@tDD ê. Resh1, Re@x21@tDDê. Resh1,
Re
@x31@tDD ê. Resh1<ê. 8t → t −τ<;
88x12@t, τD<, 8x22@t, τD<, 8x32@t, τD<< =
8Re@x12@tDD ê. Resh2, Re@x22@tDDê. Resh2,
Re
@x32@tDD ê. Resh2<ê. 8t → t −τ<;
88x13@t, τD<, 8x23@t, τD<, 8x33@t, τD<< =
8Re@x13@tDD ê. Resh3, Re@x23@tDDê. Resh3,
Re
@x33@tDD ê. Resh3<ê. 8t → t −τ<;
X
@t, τD =
i
k
x11@t, τD x12@t, τD x13@t, τD
x21@t, τD x22@t, τD x23@t, τD
x31@t, τD x32@t, τD x33@t, τD
y
{
;
                                  ПРИЛОЖЕНИЕ

      Проверка попадания в начало координат
8x1@TD, x2@TD<

8−0.0000134271, − 0.0000148833<




Пример 4.3
      Ввод начального и конечного положений фазового вектора
X0 = 8− 3, 2, 1 <;
XT = 8− 80.7746, − 147.179, − 8.94415 <

8 − 80.7746, − 147.179, − 8.94415 <

      Построение фундаментальной матрицы Коши
Resh1 =
 DSolve @8x11 ' @t D 2 ∗ x11 @t D + 2 ∗ x21 @t D − 30 ∗ x31 @t D,
    x21 ' @t D 10 ∗ x11 @t D − x21 @t D − 35 ∗ x31 @t D,
    x31 ' @t D 2 ∗ x11 @t D − x21 @t D + x31 @t D, x11 @0 D 1,
    x21 @0 D 0, x31 @0 D == 0 <, 8x11 @t D, x21 @t D, x31 @t D<,
  t D;

 DSolve @8x12 ' @t D 2 ∗ x12 @t D + 2 ∗ x22 @t D − 30 ∗ x32 @t D,
Resh2 =

    x22 ' @t D 10 ∗ x12 @t D − x22 @t D − 35 ∗ x32 @t D,
    x32 ' @t D 2 ∗ x12 @t D − x22 @t D + x32 @t D, x12 @0 D 0,
    x22 @0 D 1, x32 @0 D == 0 <, 8x12 @t D, x22 @t D, x32 @t D<,
  t D;

 DSolve @8x13 ' @t D 2 ∗ x13 @t D + 2 ∗ x23 @t D − 30 ∗ x33 @t D,
Resh3 =

    x23 ' @t D 10 ∗ x13 @t D − x23 @t D − 35 ∗ x33 @t D,
    x33 ' @t D 2 ∗ x13 @t D − x23 @t D + x33 @t D, x13 @0 D 0,
    x23 @0 D 0, x33 @0 D 1 <, 8x13 @t D, x23 @t D, x33 @t D<,
  t D; 88x11 @t, τD<, 8x21 @t, τD<, 8x31 @t, τD<< =
 8Re @x11 @t DD ê. Resh1, Re @x21 @tDD ê. Resh1,
    Re @x31 @t DD ê. Resh1 < ê. 8t → t − τ<;
88x12 @t, τD<, 8x22 @t, τD<, 8x32 @t, τD<< =
 8Re @x12 @t DD ê. Resh2, Re @x22 @tDD ê. Resh2,
    Re @x32 @t DD ê. Resh2 < ê. 8t → t − τ<;
88x13 @t, τD<, 8x23 @t, τD<, 8x33 @t, τD<< =
 8Re @x13 @t DD ê. Resh3, Re @x23 @tDD ê. Resh3,
    Re @x33 @t DD ê. Resh3 < ê. 8t → t − τ<;
             i x11 @t, τD x12 @t, τD x13 @t, τD y
X @t, τD =     x21 @t, τD x22 @t, τD x23 @t, τD ;
             k x31 @t, τD x32 @t, τD x33 @t, τD {

                                       184

Страницы