CC BY
17
ISSN 1810-0198. Вестник Тамбовского университета. Серия Естественные и технические науки
Том 23, № 122
2018
DOI: 10.20310/1810-0198-2018-23-122-158-167 УДК 517.977.5
ВЫРОЖДЕННАЯ ЛИНЕЙНО-КВАДРАТИЧНАЯ ЗАДАЧА
С ЗАПАЗДЫВАНИЕМ
< И. Ю. Андреева, А. С. Шляхов
ФГАОУ ВО «Уральский федеральный университет им. первого Президента России Б.Н. Ельцина» 620002. Российская Федерация, г. Екатеринбург, ул. Мира, 19 El-mail: ira. iyandreeva. [email protected], [email protected]
Аннотация. Рассматривается вырожденная линейно-квадратичная задача оптимизации с постоянным запаздыванием в системе и в функционале. В ранее полученных результатах авторов удалось сократить число уравнений, описывающих параметры оптимального управления. Благодаря этому впервые для таких задач просчитаны модельные примеры.
Ключевые слова: вырожденная линейно-квадратичная задача; оптимальное управление; импульсное управление; последействие
Введение
Исследование вырожденных линейно-квадратичных задач оптимизации с запаздыванием является актуальным, так как для практических задач факт вырожденности функционала встречается достаточно часто. Данная задача рассматривалась в работах [1], [2]. Невырожденные линейно-квадратичные задачи для систем с последействием рассматривались в [3], [4].
1. Основные понятия
Рассмотрим задачу минимизации функционала
~Т ~0
„Ти\,ъ т„/\л I JT/
= X1 (T)Sx(T) + X1 (i)*b(t)x(i) + X1 (t) <L\(t,6)x(t + 9)d0+
to
о о
„Tu i /)\ат/1 mjû„f+\ , т/
+ X1 (t + 0)Ф{ (i, 9)<Юx(t) + X1 (t + з)Ф2(£, s)x(t + s)ds+
Работа выполнена при поддержке Российского фонда фундаментальных исследований (проект № 16-01-505а).
о о
т
+ х1(г+в)Ф3(^в,р)х(г + р)Шр+х1(г т)ФА{€)х{г т) dt, (1)
—т —т
где Ф^^, Ф2(хЯ, Ф4(^ — симметричные (за исключением Ф3(^X^и
Ф^Х^) непрерывные по совокупности переменных матрицы-функции размерности п О71? $ — неотрицательно определенная симметричная матрица с постоянными элементами вдоль траекторий системы дифференциальных уравнений
~о
x(t) = A(t)x(t) + AT(t)x(t т) + G(i, 6)x{t + 9)d0 + B(t)v(t) (2)
— T
с начальным условием x(t) = <p(t)y при ¿о т € i e io, которую назовем задачей 1. Здесь x(t), <p(t) — вектор-функции размерности п, v(t) — вектор-функция размерности m, A(t), AT(t) —непрерывные матрицы-функции размерности п Qn, G(t, в) — непрерывная по совокупности переменных пQn, матрица^функция, B(t) —непрерывно дифференцируемая матрица-функция размерности п Q)m.
В работах |1], [2] данная задача была решена, и результаты этого решения отражены в следующей теореме.
Теорема 1. Пусть для t / [¿о, Т] выполнены условия
1) det)BT(t)<bo(t)B(t) + P6(t,0){ =0,
2) матрица ^
1Ф0(£) Ф^в) ФФ2{Ь,9)
неотрицательно определена при в / [ т. 0],
3) матрица Ф3(t,9,p) имеет структуру
4) матрица Ф4 (i) неотрицательно определена,
а матрицы P(t), Pd(t,s), Pe(f,r), Pi(t,<?i), PB(t,p), Q(t,0), P2(t,Mi), P{£, 9, p) = RT(t, p, в), , #2) = Pj(t,02,d 1) являются решением системы уравнений
dP{t) + Q(t, 0) + QT(t, 0) + PT(t)A(t) + AT(t)P(t) + P4(t, 0) + Ф0(г) =
dt
= Fj(t)H~1(t)F0(t),
r\ r\
) dt w (g(i) в) + m в) + P(i)G(i'e) + + & в) =
= FT(t)H~i(t)F1(t,e)7
) dt W Щ p) + GT(i'p) + QT{t>9)G{t'p) + Фз{£'p) =
= F?(t,e)H~1(t)F1(t,p),
О г\
) dt Щ (Pl(i'9l) + 9l) + AT(i)Pi(i> öl) + 9+
Ö!)G(i, e2)B(t + 02) + BT(t + вг)Ф3(1, вг, e2)B(t+02) = = Fj(t191)H~\t)F2(t192),
О о л о
)öt + öp(Pb(t,p) + $2(t,p)B(t+p)= 0,
) it îr (Pe(i'r) + ßr(i + r)$2(i'r)jB(i + r) = 0 (3)
тгрн граничных условиях
P(T) = ЛГ, Q(T,0) = Я(ТДя) = Pi(T, öi) =
= Р2(Т,0А) = Рз(Т, 9ъ92) = Pi(T. s) = Р5(Т,р) = P6(T,r) = 0, P(t)AT(t)B(t r) = P1(t, r), BT(t r)ATr{t)Q{t19) = pT{t,e, т),
ATT(t)P{t) = QT{t1 t), ÂÏ(t)Pl(t,e1) = P2(t, rA), 5r(i T^WAMO + ^MOArfW r) = Pj(£,fl1, r) + P3(i, $4(t) = P4(i, г), Ф4(£)Б(£ r)=JP5(i, r), ßr{i г)Ф4(£)В(£ r)=P6(t, r), для t < T, те M !, 02, r, s € 0, где
Po(i) = AT(i, 0) + Bf(t)P(i) + p5r(i; 0) + 5r(i)5>o(t),
F1(t,9) = B'[(t)Q(t,e) + P^t,e,0) + BTm1(t,e), F2M) = BÏ(t)p1(t,el) + pT(t,el,o) + BT(t)<ï>l(t,e1)B(t+e1), H(t) = BT(t)<b0(t)B(t) + pe(t,0), B^t) = A(t)B(t) B(t). Тогда на [io,T] оптимальное управление примет вид
v{t) = A»(to, S(t t0) + + Z2{t)x{t т)+
~o
+ Z3{t, 0)x(i + в) дв + Дг>(Т, 0)) S(t Г),
где величины скачков функции v(t), определяются формулами
Av{t0: = W0(t0)<p(t0) + W^to, 0)<p(to + 9) d9,
и
Av(T,x(T 0))= ^Bt(T)SB(T)(^ Bt{T)Sx{T 0) + +) E ) Bt{T)SB(T) ^BT(T)SB(T) |p
wiM) = н (t) + + ,
=
ctlJl i / (i0,P)-
dt
+ Wo(t)A(t) + W1(t,0) , Z2(t)= \\V0(t)AT(t) W^t, T) ,
2. Упрощение системы для нахождения коэффициентов оптимального управления
Основным результатом по упрощению системы уравнений (3) является ее приведение к системе с меньшим количеством уравнений. Рассмотрим второе и четвертое уравнение этой системы. Умножим второе уравнение на B(t + 9) справа, а затем из полученного уравнения вычтем четвертое уравнение, для которого в качестве второго параметра возьмем 9. В результате выполнения указанных действий получим следующее уравнение:
it ^(QM^it + e) ^(Pi(tA) + R(tA9)B(t + 9)
P2(t, 0, вг) + AT(t)Q(t19) B(t + 9) AT(t)P1(t,9) = = Fj(t)H~1(t)F1(t,e)B(t + 0) F^{t)H-l(t)F2{t.,9). Подставим в это уравнение:
Pi(i, 9) = В?(t)Q(t, 9) + P[(t,9; 0) + ВТ(ф1(г,в), F2(t,e1) = Bf(t)P1(t,e1) + Pl(t,e1,0) + BT(^1(t,91)B(t + e1).
В результате получим уравнение
P2(t, О, в,) + AT(t)Q(t1 9) B(t + 9) AT(t)P1(t,9) =
= F^(t)H-1(t)[Bj(t)Q(t,e)B(t + 9) Bl(t)P1(t79)+pT(t,910)B(t + 9) JfM,0)].
Очевидно, что данное уравнение станет тождеством, если выполняются следующие условия:
B(t + s) = const, P1(t,0)=Q(t,0)B(t+0), P2(t,0.9) = R(t,O,0)B(t + 0)
PÏ(t,9,0) = PÏ(t,9,0)B(t + 9).
Проведя аналогичные рассуждения с остальными уравнениями системы и граничными условиями, можно прийти к выводу, что при выполнении равенств:
Рг(t, 9) = Q(t, 9) B{t + 9), P2{t, О,9) = R{t, 0,9) B(t + 9), P3r(i, 0,0) = P2r(i, 9,0) B(t + 9),
P5(t, s) = P4(f, s) B(t + s), P6(£, s) = BT(t + s) P4(i, s) B(t + s) система примет вид
+ Q(i) 0) + QTfa 0) + PT(t)A(t) + AT(t)P(t) + P4(i, 0) + Ф0(г) =
at
= Fj(t)H~1(t)F0(t),
) Ж Ъв (Q(i'в) + т е) + ^ + P(i) Gв) + в) =
= F^{t)H~\t)F1{tt0),
) ш ъё щ р) + GT(i'р)+QT{t>9) G(i'р) + Фз(£'р) =
= Ff(t,0)H~1(t)F1(t,p),
с граничными условиями
Р(Т) = N, Q(T,0) = P(T,^) = P4(7>)=O, Ф4(£)=Р4(£, г),
2^(i)Q(i,0)=JRT(i,0, r) + P(t, т,0), ^(i)P(i)=Qr(i, т), при t <Т, те 0,01,^2) р, s е 0, где
P0(i) = Q(t, 9) + P(i) + P4(i, 0) + Ф0(£), Fi(t, 0) = 0) + P(£, 0) + Фг(*,
Я(£) = Фо(£) + Р4(г,0).
Пример 1. Рассмотрим задачу минимизации функционала (1) вдоль траекторий системы дифференциальных уравнений
*М = )? + ) о1 £(»<«
с начальными условиями: <р(Ь) = 1, £ / [ 5;0]. Из постановки задачи имеем, что
I ; п I; = I ^ ; I; я = ) I ( ¿о = 07 т = 5, г=1, с = А 0 0
) 1 2 Q1 0
1 " ' ' ' /00
Для данной задачи будут иметь место следующие равенства:
р4<м)) ] [ = рб(м), (1 1 )p4(t,s) 1 ; f =p6(M), qm ) : f =Pi(t,0),
1
с учетом которых система для нахождения неизвестных функций примет вид:
+СЭМ+дт&0)+Ртю) ; I ЛЛ? I /Р(£) + Р4(£,0)+Фо(£) =
где
Л. ' 7 ^ у ' 7 2^12
= ^(£)Я-1(£)Р0(£),
^ ^ + Ф3(г, в, р) = в)!!1 (1)Рг (£, р),
)т в;(Р*(*>*)+ **(*,*) = о,
Р0(г) = (1 1)<5г(^0) + (2 2)Р(£) + (1 1)Р4г(^о) + (1 1)*ь(*). Р,М) = (2 2)(3(г,9) + (1 1)Рт(£,М) + (1 1)Ф1(£,6>),
Я(£) = (1 1)Ф0(£)) | ^+(1 )
Граничные условия для этой системы будут следующими:
\ пЛ 1 "1 I а1"2 а1а2+а1а2 /• \ п г п г
Р(Т) = М= I +а'2 <Ч+«2 (01+а.2)а / = I ' '
V ' I 2а 102 а^а-г+аш^ /-, а2 \2 °1°2 .7 / 0 5 0 5
д(Т,0) = Р(ТДр) = Р4(Т>) = О, ) ^ ^ [ Р(£) = <?' (£, 5),
) 0! I + ^ ^ = 5) + Р(£, 5,0), Ф4(0=Р4(4, 5)
для ^ < 1, 5 € в, р, .ч е 0. В результате решения данной системы уравнений с указанными граничными условиями были получены неизвестные функции для нахождения оптимального управления. Тогда вид оптимального управления будет следующим:
¿,(£) = ф) + 2х1(г) + х2{г) Х1(г т)
+х2(1 т)Х¥№)
величины скачков функции «(£), соответственно, будут определятся формулами:
Д»(*,) = 1,(0) = 0) )) ^ ) ^ ^ = ^(0)^(0 ), Д ь{Т,х(Т 0))= (Вт(Т)ЗВ(Т))1Вт(Т)Зх(Т 0),
Ро(£)=)в^-> + 1 ФШ) ^ 0
Я(*) = ф1(*)+ «£(*) +«3(4)+«£(*) (г Г)2, Я-1(^)Ро(0= 0(,
16е8Ст"() + 10е5(т"*> 40е4(т"£) + 5е2(т"£) 18е(т"{> ~ 8 + 6еГ-£ + 16е<т"() +
(8е4(г-^) + 4е(т-о 4) (Г ¿)2 + 8 + 8 + + 16е4(т"£)
ФШ)=Ф&г)= е^ е—
т) =
ИоИ =
2 '
8е5(т"{> 32е4<т"() 4е2(т"*> ет"((6 (£ Г)2)
б4е8 (т-0 2е5(Т_Г) 4е2(т"*> + Т)2
8е^т"() 32е4(т"*> 4е2<т-£> + ¿т~1\Т I)2 64е8 (Т-») + 32е5(т-£) + 4е2(Т-г) ^
СПИСОК ЛИТЕРАТУРЫ
1. Андреева И.Ю., Сесекин А.Н. Вырожденная линейно-квадратичная задача оптимизации с запаздыванием по времени // Автоматика и телемеханика. 1997. № 7. С. 43-54.
2. Сесекин А.Н., Фетисова Ю.В. Вырожденная линейно-квадратичная задача оптимизации для систем с последействием // Аналитическая механика, устойчивость и управление движением: тр. 9 Междунар. Четаевской конф., посвящ. 105-летию Н.Г. Четаева. Иркутск, 2007. Т. 3. С. 196-205.
3. Колмановский В.Б., Майзенберг Т.Л. Оптимальные оценки состояния системы и некоторые задачи управления уравнениями с последействием // Прикладная математика и механика. 1977. Т. 41. № 3. С. 446-456.
4. Янушевский Р.Т. Управление объектами с запаздыванием. Москва: Наука, 1978. 416 с.
Поступила в редакцию 19 марта 2018 г. Прошла рецензирование 24 апреля 2018 г. Принята в печать 5 июня 2018 г. Конфликт интересов отсутствует.
Андреева Ирина Юрьевна, Уральский федеральный университет им. первого Президента России Б.Н. Ельцина, г. Екатеринбург, Российская Федерация, кандидат физико-математических наук, доцент кафедры прикладной математики, e-mail: [email protected] Шляхов Александр Сергеевич, Уральский федеральный университет им. первого Президента России Б.Н. Ельцина, г. Екатеринбург, Российская Федерация, аспирант, кафедра прикладной математики, e-mail: [email protected]
Для цитирования: Андреева И.Ю., Шляхов А.С. Вырожденная линейно-квадратичная задача с запаздыванием // Вестник Тамбовского университета. Серия Естественные и технические науки. Тамбов, 2018. Т. 23. № 122. С. 158—167. DOI: 10.20310/1810-0198-2018-23-122-158-167
DOI: 10.20310/1810-0198-2018-23-122-158-167
DEGENERATED LINEARLY-QUADRATIC PROBLEM WITH DELAY
I. Yu. Andreeva, A. S. Shlyakhov
Ural Federal University named after the first President of Russia B.N. Yeltsin 19 Mira St., Yekaterinburg 620002, Russian Federation E-mail: ira. iyandreeva. [email protected], [email protected]
Abstract. We consider a degenerate linear-quadratic problem optimization with a constant delay in the system and in the functional. In the previously obtained results of the authors, it was possible to reduce the number of equations describing the parameters of optimal control. Thanks to this, for the first time model examples are calculated for such problems.
Keywords: degenerate linearly-quadratic problem; optimal control; impulse control; aftereffect
REFERENCES
1. Andreeva I.Y.. Sesekin A.N. Vyrozhdennaya lineyno-kvadratichnaya zadacha optimizatsii s zapazdyvaniem po vremeni [Degenerate linear-quadratic optimization problem with time delay]. Avtomatika i telemekhanika - Automation and Remote Control, 1997, no. 7, pp. 43-54. (In Russian).
2. Sesekin A.N.. Fetisova Y.V. Vyrozhdennaya lineyno-kvadratichnaya zadacha optimizatsii dlya sistem s posledeystviem [Degenerate linear-quadratic optimization problem for systems with aftereffect]. Trudy 9 Mezhdunarodnoy Chetaevskoy konferentsii «A n atitich e shay a mekhanika, ustoy-chivost' i upravlenie dvizheniem», posvyashchennoy 105-letiyu N.G. Chetaeva [Proceedings of the International Chetaev Conference "Analytical Mechanics, Stability and Motion Control" dedicated to the 105th anniversary of N.G. Chetayev]. Irkutsk, 2007, vol. 3, pp. 196-205. (In Russian).
3. Kolmanovskiy V.B., Mayzenberg T.L. Optimal'nye otsenki sostoyaniya sistemy i nekotorye zadachi upravleniya uravneniyami s posledeystviem [Optimum estimates of the state of the system and some problems of control of equations with aftereffect]. Prikladnaya matematika i mekhanika -Journal of Applied Mathematics and Mechanics, 1977, vol. 41, no. 3, pp. 446-456. (In Russian).
4. Yanushevskiy R.T. Upravlenie ob"ektami s zapazdyvaniem [Management of Objects with Delay], Moscow. Nauka Publ.. 1978, 416 p. (In Russian).
Received 19 March 2018
Reviewed 24 April 2018
Accepted for press 5 June 2018
There is no conflict of interests.
The work is partially supported by the Russian Fund for Basic Research (project № 16-01-505a).
Andreeva Irina Yurevna, Ural Federal University named after the first President of Russia B.N. Yeltsin, Yekaterinburg, the Russian Federation, Candidate of Physics and Mathematics, Associate Professor of Applied Mathematics Department, e-mail: [email protected]
Shlyakhov Alexander Sergeevich, Ural Federal University named after the first President of Russia B.N. Yeltsin, Yekaterinburg, the Russian Federation, Post-Graduate Student of the Department of Applied Mathematics, e-mail: [email protected]
For citation: Andreeva I.Yu., Shlyakhov A.S. Vyrozhdennaya lineyno-kvadratichnaya zadacha s zapazdyvaniem [Degenerated linearly-quadratic problem with delay]. Vestnik Tambovskogo universiteta. Seriya Estestvennye i tekhnicheskie nauki -Tambov University Reports. Series: Natural and Technical Sciences, 2018, vol. 23, no. 122, pp. 158-167. DOI: 10.20310/18100198-2018-23-122-158-167 (In Russian, Abstr. in Engl.).
