−8 +0.845845tCos t +2.Sin t −0.533429tSin t ,
2.Cos t −0.533429tCos t +3.Sin t −0.845845tSin t ,
2.0562Cos t −1.29673Sin t ,−1.29673Cos t −2.0562Sin t
ParametricPlot x10 t ,x20 t , t,0, π ,
Graphics
Исследование поведения решения сопряженной системы
Plot ψ10 t , ψ20 t , t,0, π ,
AxesLabel@8 @ →D8"t","@ D<ψ1,ψ28"<D <
Graphics
Вычисление оптимального значения функционала
3 ∗ x10 π ^2 + 2 ∗ x20 π ^2
0.56252@ D @ D
Проверка постоянства функции Понтрягина
ψ t ∗ik−8x10 t + ψ20 ;
t
Graphics
Вычисление значения функционала на произвольном управлении
Srav =
NDSolve
− @@ ψ10 t 0.1 ^2 ,
y1 0 == −3,y2 0 2 , y1 t ,y2 t , t,0, π ;
y1 t_ ,= y1 t .Srav,y2 t .Srav ;
388∗ y1@@πD^2 + 2 y2@ D^2
0.598124
Пример 2.3.
Интегрирование сопряженной системы с граничными условиями
@ |
, tD,
,
tD ψ1 t_ , t_ , ψ3@ D<<@ D<=
Re.Sopr , Re @ Dê D
@ @ Dê D<решения сопряженной системы
ψ1
Graphics
ψ2
Graphics
700ψ3
Graphics
Определение момента переключения управления
Интегрирование основной системой с полностью определенными оптимальными управлениями и с начальными условиями
Osn =
x2' t 10 x1 t x2 t 35 ∗ x3 tD −
1,
8x1t_ D−, 3.,x2@@0DD 2@,xD3@0D@==<DD1.<,
88x1 @ DDD<@@ DêDê
Re@@x1 t,
8Re@x3 t
Вычислениеоптимального значения функционала
x1 1 + 2 ∗ x2 1 − x3 1
−366.188@ D @ D @ D
Сравнение с неоптимальным управлением
Srav =
tD −
1, y10@DD−3.,y2@@0D@DD<2@@,y@DD3D<<Dê@0D@==D1.D<,
t t ,t ; D , y3 t_ =
Re y1 @y2 t .Srav ,
y1
−365.348@ D @ D @ D
Проверка постоянства функции Понтрягина
Graphics
Пример 2.5.
Построение фундаментальной матрицы Коши
Resh1 =
@ Dt;
x22'@8t x12D@D<tD + 2 ∗Dêx22@@@DD<tD 8 @@DD<<@DDê <,
tx21 t ,
t .Resh2,x21 t .Resh1,
X =.t → 1 − τ
::
ипсилон
XT = Transpose X ;L = l1 ;
;
τ
=
Построениедвух ветвей функции ипсилон
Eps1 l1_; Eps2@@l1_DD = Dob2 − NIntegrate@ @ D 8 <<DD;
Решение задачи математического
u = −1;If Eps1 −1 < Eps2 −1 ,S = Eps2 −1 = u;