ISSN 1810-0198. Вестник Тамбовского университета. Серия Естественные и технические науки
Том 23, № 122
2018
DOI: 10.20310/1810-0198-2018-23-122-200-209 УДК 519.71, 517.977.15
ЗАДАЧА ОПТИМАЛЬНОГО УПРАВЛЕНИЯ ДЛЯ ОДНОЙ СИСТЕМЫ ЛИНЕЙНЫХ НАГРУЖЕННЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
<: В. Р. Барсегян1'^, Т. В. Барсегян^
^ Ереванский государственный университет 0025, Армения, г. Ереван, ул. А. Манукяна. 1 E-mail: [email protected] Институт механики HAH Армении 0019, Армения, г. Ереван, пр. М. Баграмяна, 24Б E-mail: [email protected]
Аннотация. Рассмотрена задача управления и оптимального управления одной системы линейных нагруженных дифференциальных уравнений. Сформулировано условие существования программного управления и движения. Приведен явный вид управляющего воздействия для задачи управления и предложен конструктивный подход решения задачи оптимального управления. В качестве приложения построено решение задачи оптимального управления конкретной нагруженной системы.
Ключевые слова: нагруженные дифференциальные уравнения; поэтапно меняющиеся системы; управление; оптимальное управление
Введение
Математическое описание динамических процессов управления, зависящих не только от настоящего, но и от предыстории процесса, осуществляется при помощи обыкновенных дифференциальных уравнений с памятью различных видов, называемых также нагруженными дифференциальными уравнениями. Нагруженными дифференциальными уравнениями в литературе [1-4] принято называть уравнения, содержащие в коэффициентах или в правой части какие-либо функционалы (функции) от решения, в частности, значения решения, в которых фазовое состояние процесса в какой-либо точке и в какой-либо момент может оказывать влияние на динамику процесса в целом.
Наличие в динамике системы нагруженного слагаемого не всегда позволяет непосредственно применять известные методы исследований, развитые при исследованиях обычных (не нагруженных) динамических систем. Это подчеркивает как теоретическую, так и практическую актуальность исследования различных задач управления и оптимального управления для нагруженных дифференциальных уравнений. В последние годы проводится интенсивное исследование нагруженных дифференциальных
уравнений, связанное с различными прикладными задачами механики, биологии, экологии и химии, моделированных с нагруженными уравнениями. Интерес исследователей к задачам управления нагруженных динамических систем связан также с необходимостью использовать наиболее адекватные математические модели рассматриваемых процессов.
Нагруженные обыкновенные дифференциальные уравнения и краевые задачи для таких уравнений рассмотрены в [1-5] и установлены условия их разрешимости различными методами. Значительный вклад в развитие теории нагруженных уравнений внесла работа [1] (и другие работы этого же автора), где даны определения нагруженных дифференциальных, нагруженных интегро-дифференциальных, нагруженных функциональных уравнений и их многочисленные приложения. В монографии [2] нагруженные дифференциальные уравнения интерпретируются как возмущения дифференциальных уравнений. В работе [5] рассмотрена задача управления для одной системы линейных нагруженных дифференциальных уравнений с неразделенными многоточечными промежуточными условиями.
В данной работе рассматривается задача оптимального управления одной системы линейных нагруженных дифференциальных уравнений. Сформулировано условие существования программного управления и движения. Предложен конструктивный подход решения задачи оптимального управления. В качестве приложения построено решение задачи оптимального управления конкретной нагруженной системы.
1. Постановка задач
Рассмотрим управляемый процесс, динамика которого описывается нагруженными линейными дифференциальными уравнениями
X = Aa(t)x + Al (£)a;(ii) + A2(t)x(t2) + A3(t)x(t3) + B(t)u (1.1)
где x(t) / Rn - фазовый вектор системы, Ak(t) (к = 0, 1, 2, 3), B(t) матрицы параметров системы (непрерывны на [t0, Т]), u(t) управляющее воздействие, соответственно с размерностями /U(£) (п О«) , B(t) (п Ог)> u(t) (r Ol)-
В формуле (1.1) слагаемые Ak(t)x(tk) (к = 1,2, 3) как функции влияют на систему, начиная с момента времени t Е t^. Так как значение фазового состояния как
результат измерения, определяется в момент времени t = tk и с этого момента (при i € непрерывно влияет на систему в виде слагаемого
Пусть заданы начальное x(t0) = х0 и конечное х(Т) = хт состояния системы (1.1).
Предполагается, что заданы моменты времени точек нагружения и 0 > i0 < ii < t2 <
t3 < ¿4 = Т. Функция x(i) непрерывна на интервалах [ijt-i, tk) (к = 1,4) ив точках
нагружения tk имеет конечные левосторонние пределы lim x(t) = x(tk).
t—>fjj—о
Рассмотрим следующие задачи.
Задача 1. Требуется найти условия, при которых существует программное управляющее воздействие и = u(t) и программное движение, переводящее движение системы (1.1) из начального состояния x(to) в конечное состояние х(Т) на промежутке времени [£о,Т], а также построить их.
Пусть для отбора оптимальных решений на промежутке времени [i0, Т] задан критерий качества ве[и], который может иметь смысл нормы некоторого нормированного пространства.
Задачу оптимального управления для системы (1.1) с начальной и конечной условиями и критерием качества «[и] можно сформулировать следующим образом.
Задача 2. Требуется найти оптимальное управляющее воздействие w°(i), t / [i0, Т], которое переводит движение системы (1.1) из начального состояния x(t0) в конечное состояние х(Т) и имеет наименьшее возможное значение критерия качества as [и0].
2. О движении нагруженной управляемой линейной системы.
Для построения решений поставленных задач интервал [io,T] разбиваем на части
4
точками нагружения: [fo,Т) = (J [£jt-i, tk). Учитывая последовательность точек на-
к 1
гружения и характер влияния соответствующих нагруженных слагаемых, уравнение (1.1) напишем по отдельности на интервалах (к = 1,4) в следующем виде
iA0(t)x + B{t)u, t / [£0,ii)
A0(t)x + A^xit0 + B(t)u, t / [tut2)
Mt)* + AW^(ii) + A2(t)x(t2) + B(t)u, t/[t2,t3) ' 1 " ;
A0(t)x + Аг(г)х(10 + A2(t)x(t2) + A3(t)x(t3) + B(t)u, t / [£3, T]
которое является поэтапно меняющимися дифференциальными уравнениями [5-7].
Следуя [5], решение системы (1.1) (или (2.1)) для промежутка времени [io,ii) представим следующим образом
t
x(t) = X[t, t0]x{t0) + J H[t, t]u(t)cLt, (2.2)
to
где ii[t,r] = X[i,r]5(r), а через X[t, г] обозначена нормированная фундаментальная матрица решения однородного уравнения х = A0(t)x. Начиная с момента времени движения системы (2.1) представляются в следующем виде
tL t x(t) =X[i,i1]X[i1,t0]s;(i()) + i'[t,i1] I H[tur]u(T)dT + j H[t, т\и(т)<1т (2.3)
при t / [ii,i2),
io
при t /
x(t) = Y[t, t2]X[tu t0]x(t0) + Y[t, t2] J H[tur]u(r)dT+
to
t2 t +X[t,t2]j H[t2,r]u(T)dT+ f H[t,T]u(r)dT ti t2
ti
x(t)=Z[t,t3\X[t1,t0]x(t0) + Z[t,t3] J H[tur]u{T)dr+
to
(2.4)
£ \ ¡2
+ I Х[1, Т]А2(Т)<1Т + Х[1, ¿3]Х[£3, *2] I Н[г2,т]и(т)С1т+ (2.5)
Ч(з / (1
4з £
+Х[Мз] I Н[13,т\и{т)д.т + J Н[г,т\и{т)6.т
при £ / [£3,£4], где
г
ХЫ=Х[1,^] + !Х%т\А^т)дт 0 = 1,2,3), ь
г
Г[М2]=*[М2№Л] + I Х^А^т,
t2
( 4
^[Мз] = I Х[г1т]А1(т)^ + I Х% Т\А2{т)<1ТХ[12, + £3] У[£3, ¿2]. (2.6)
¿3 ¡3
Таким образом, имея начальное состояние ж(£0) системы (1.1) и задавая управляющее воздействие с помощью формулы (2.2)-(2.5) определяется фазовое состояние х(1) системы (1.1) (решение нагруженного уравнения (11)) для промежутков времени (к = 2,...,4), соответственно.
3. Решение задач
Из формулы (2.5) при £ = имеем
т , 4
£о
где
/ евд
и(1)М = х(и) 2[гл,ь]Х[ь,г0]х{г0) Г), (3.1)
С „г I о, ¿о > * < ¿1
1 -* - Г I О, 12>1>Т
Г > £ < ГО £ > £ < £
I 0. £3>£>Г 1 ЩиЛ
Я1[£]=2[£4)£3]Я[£1,£], Я2[£] = у Х[£4,т]А2(Т)^Т + Х[£4,£з]Х[£3,£2] )
ЯзМ = Х[ЬМН[ЬЛ нАщ = Н[ил\.
Функция и{Ь) £ / [¿о,?!) удовлетворяющая интегральному соотношению (3.1), имеет следующий вид [6]
Здесь и(£) - некоторая вектор-функция, удовлетворяющая условию ортогональности
Следующая теорема, аналогичная теореме, приведенной в [6, 9], определяет условия существования программного управления системы (1.1).
Теорема 3.1. Для того, чтобы существовало программное управление и соответствующее ему решение системы (1-1), удовлетворяющее условию (3.1), необходимо и достаточно, чтобы матрица (3-4) была неособой или чтобы ранги матрицы и совпадали между собой.
При заданном критерии качества ее[и] задача оптимального управления с интегральными условиями (3.1) является задачей условного экстремума из вариационного исчисления, где надлежит определить минимум функционала ае[г(] при условии (3.1). Однако, как видно из (3.2), подынтегральные функции в (3.1) являются разрывными, поэтому классические методы вариационного исчисления не применимы для решения этой задачи [8].
Для решения задачи 2 заметим, что левая часть полученного условия (3.1) является линейной операцией, которая порождена функцией и(£) на промежутке времени [¿о,Т]. Условия (3.1) с заданным функционалом ае[и], представляющим собой норму некоторого линейного нормированного пространства, являются проблемой моментов в соответствующем пространстве [6, 8]. Следовательно, решение задачи 2 приводится к решению соответствующей проблемы моментов, которое следует построить с помощью алгоритма решения проблемы моментов [6, 8]. Тогда построенное оптимальное управляющее воздействие и°(1), t / [¿о■ Т], удовлетворяющее условию (3.1) и минимизирующее функционал 8е[и], будет решением задачи 2.
4. Пример
В качестве приложения к вышеизложенному рассмотрим задачу оптимального управления конкретной системы (с одним нагруженным слагаемым), которая состоит из двух этапов и имеет вид
(3.3)
если с^={= 0, где
(3.4)
х = Ао(£)х + В(г)и, I /[io.ii)
(4.1)
X = Аъ(г)х + + В{г)и, I / [£ь Т) (4.2)
где
*-(;;)• *-(")• -с
а - постоянная величина.
Пусть заданы начальное и конечное фазовые состояния
а критерий качества управления имеет вид
1
Т \ 2
ае[и] = |у . (4.4)
Нормированная фундаментальная матрица решения однородной части системы (4.1), (4.2) и их переходные матрицы имеют следующий вид:
*М = (,з V), я[т,4] = (т1().
Из формулы (2.6) для матрицы будем иметь
ХЬЬ] =ХЬЬ]+ IХ\Ь,т\А\{т)(1т = ^ 1 + |Г ^ Т 1 ^ ) .
Согласно обозначениям (3.2) будем иметь
нщ = х[т,г1]й[г1,г]+н[т^] =
где
Й[гъг] =
( /^М] А = Г ¿ь>*<41,
^ Л?'!*!,*] ^ I 0, £1>£>Г,
а = 1 + ^(Г £г)2, /3 = Т 7 = а(Т £,)■
Интегральное соотношение (3.1) запишется в виде
т т
I ЪфифМ =1}Ъ I =т]2, (4.6)
£о ¿о
/ц(£) = оЛ?^,*] + /?/41}М] + М2)[Т,£], = 7М1}М] + Л^М] + 2)[Г,£],
»Та = Ж1(Т) ОХ&о) [«(¿1 £0) + ,5]х2(£0), (4.7)
г]2 = х2(Т) 7ж1(£0) [7(^1 ¿о) + 1]ж2(£0)-
Задача определения функции оптимального управления и°(£), £ / [£о,Т], удовлетворяющего интегральным соотношениям (4.6) и минимизирующего функционал (4.4), рассматриваем как проблему моментов [6, 8].
Решая полученную задачу условного экстремума (4.4) и (4.6) как проблему моментов, для искомого оптимального управляющего воздействия получим следующие выражения:
Г с[(ш*! + 7^2)^1 Ь /[Iо,£0,
\ с(Т £)-Ы2], £/[£ьН
где с = (т^! + т}2<12) + + = ЫЯз Ц2Ч2) , = (т^ т^),
т
41 = / (Ьг(г))2Л = + а2(Ь £0) + + «£(£1 £0)2,
(О
42 = /ь и^т = + (а + + Ж ¿о) +
¡0
Чл = = + ^ ¿о) + (Г £1) + ¿о)2
¡0
Отметим, что при вычислении значения интегралов для функции /^1(£) и /г2(£) (4.7) учтены обозначения (4.5).
Подставляя выражение для и°(£) из (4.8), соответственно, в формулы (2.2) и (2.3), написанные для системы (4.1) и (4.2), получим явные выражения для оптимального движения х°(£) при £ / [£<), Т].
5. Заключение.
Предложен конструктивный подход исследования задач оптимального управления для одной системы линейных нагруженных дифференциальных уравнений. Сформулировано условие существования программного управления и движения. Приведен явный вид управляющего воздействия для задачи управления и предложен конструктивный подход решения задачи оптимального управления. Для иллюстрации указанного способа построено решение задачи оптимального управления конкретной нагруженной системы.
СПИСОК ЛИТЕРАТУРЫ
1. Нахушев А.М. Нагруженные уравнения и их применение. М.: Наука, 2012. 232 с.
2. Дженалиев М.Т., Рамазанов М.И. Нагруженные уравнения как возмущения дифференциальных уравнений. Алматы, 2010. 334 с.
3. Дженалиев М.Т. Оптимальное управление линейными нагруженными параболическими уравнениями // Дифференциальные уравнения. 1989. Т. 25. № 4. С. 641-651.
4. Кожанов А.И. Нелинейные нагруженные уравнения и обратные задачи // Журнал вычислительной математики и математической физики. 2004. Т. 44. № 4. С. 694-716.
5. Барсегян В.Р. Управление составных динамических систем и систем с многоточечными промежуточными условиями. М.: Наука, 2016. 230 с.
6. Барсегян В.Р. Задача управления для одной системы линейных нагруженных дифференциальных уравнений с неразделенными многоточечными промежуточными условиями // Известия Иркутского государственного университета. Серия Математика. 2017. Т. 21. С. 19-32.
7. Barseghyan V.R. Control of stage by stage changing linear dynamic systems // Yugoslav Journal of Operations Resarch. 2012. Vol. 22. № 1. P. 31-39.
8. Красовский Н.Н Теория управления движением. М.: Наука, 1968. 476 с.
9. Зубов В.И. Лекции по теории управления. М.: Наука, 1975. 496 с.
Поступила в редакцию 26 марта 2018 г.
Прошла рецензирование 23 апреля 2018 г.
Принята в печать 5 июня 2018 г.
Конфликт интересов отсутствует.
Барсегян Ваня Рафаелович, Ереванский государственный университет, г. Ереван, Армения, доктор физико-математических наук, профессор кафедры механики; Институт механики НАН Армении, ведущий научный сотрудник, e-mail: [email protected]
Барсегян Тигран Ваняевич, Институт механики НАН Армении, г. Ереван, Армения, кандидат физико-математических наук, научный сотрудник, e-mail: [email protected]
Для цитирования: Барсегян В.Р., Барсегян Т.В. Задача оптимального управления для одной системы линейных нагруженных дифференциальных уравнений // Вестник Тамбовского университета. Серия Естественные и технические науки. Тамбов, 2018. Т. 23. № 122. С. 200-209. БОТ: 10.20310/1810-0198-2018-23-122-200-209
DOI: 10.20310/1810-0198-2018-23-122-200-209
THE PROBLEM OF OPTIMAL CONTROL FOR ONE SYSTEM OF LINEAR LOADED DIFFERENTIAL EQUATIONS
V. R. Barseghyan12), T. V. Barseghyan2)
Yerevan State University 1 Alek Manukyan St., Yerevan 0025, Armenia E-mail: [email protected] 2) Institute of Mechanics of National Academy Science of Armenia 24B Baghramyan Ave., Yerevan 0019, Armenia E-mail: [email protected]
Abstract. The problems of control and of optimal control of a system of linear loaded differential equations are considered. The condition for the existence of program control and motion is formulated. An explicit form of the control action for the control problem is given and a constructive approach to solve the optimal control problem is proposed. As an application, a solution to the problem of optimal control of a concrete loaded system is constructed.
Keywords: loaded differential equations; stage-by-stage changing systems; control; optimal control
REFERENCES
1. Nakhushev A.M. Nagruzhennye uravneniya i ikh primenenie [Loaded Equations and their Applications]. Moscow, Nauka Publ., 2012, 232 p. (In Russian).
2. Dzhenaliev M.T., Ramazanov M.I. Nagruzhennye uravneniya kak vozmushcheniya differen-tsial'nykh uravneniy [Loaded Equations as a Perturbation of Differential Equations]. Almaty, 2010, 334 p. (In Russian).
3. Dzhenaliev M.T. Optimal'noe upravlenie lineynymi nagruzhennymi parabolicheskimi urav-neniyami [Optimal control of linear loaded parabolic equations]. Differentsial'nye uravneniya -Differential Equations, 1989, vol. 25, no. 4, pp. 641-651. (In Russian).
4. Kozhanov A.I. Nelineynye nagruzhennye uravneniya i obratnye zadachi [Nonlinear loaded equations and inverse problems]. Zhurnal vychislitel'noy matematiki i matematicheskoy fiziki -Computational Mathematics and Mathematical Physics, 2004, vol.44, no. 4, pp. 694-716. (In Russian).
5. Barsegyan V.R. Upravlenie sostavnykh dinamicheskikh sistem i sistem s mnogotochechnymi promezhutochnymi usloviyami [Control of Compound Dynamic Systems and of Systems with Multipoint Intermediate Conditions]. Moscow, Nauka Publ., 2016, 230 p. (In Russian).
6. Barsegyan V.R. Zadacha upravleniya dlya odnoy sistemy lineynykh nagruzhennykh diffe-rentsial'nykh uravneniy s nerazdelennymi mnogotochechnymi promezhutochnymi usloviyami [The control problem for a system of linear loaded differential equations with nonseparated multi-point intermediate conditions]. Izvestiya Irkutskogo gosudarstvennogo universiteta. Seriya Matematika -The Bulletin of Irkutsk State University. Series "Mathematics", 2017, vol. 21, pp. 19-32. (In Russian).
7. Barseghyan V.R. Control of stage by stage changing linear dynamic systems. Yugoslav Journal of Operations Resarch, 2012, vol. 22, no. 1, pp. 31-39.
8. Krasovskiy N.N. Teoriya upravleniya dvizheniem [The Theory of Motion Control]. Moscow, Nauka Publ., 1968, 476 p. (In Russian).
9. Zubov V.I. Lektsii po teorii upravleniya [Lectures on the Theory of Control]. Moscow, Nauka Publ., 1975, 496 p. (In Russian).
Received 26 March 2018 Reviewed 23 April 2018 Accepted for press 5 June 2018 There is no conflict of interests.
Barseghyan Vanya Rafayelovich, Yerevan State University, Yerevan, Armenia, Doctor of Physical and Mathematical Sciences, Professor of the Department of Mechanics; Institute of Mechanics of the National Academy of Sciences of Armenia, Leading Researcher, e-mail: [email protected]
Barseghyan Tigran Vanyaevich, Institute of Mechanics of the National Academy of Sciences of Armenia, Yerevan, Armenia, Candidate of Physical and Mathematical Sciences, Researcher, e-mail: [email protected]
For citation: Barseghyan V.R., Barseghyan T.V. Zadacha optimal'nogo upravleniya dlya odnoy sistemy lineynyh nagru-zhennyh differentsial'nyh uravneniy [The problem of optimal control for one system of linear loaded differential equations]. Vestnik Tambovskogo universiteta. Seriya Estestvennye i tekhnicheskie nauki - Tambov University Reports. Series: Natural and Technical Sciences, 2018, vol. 23, no. 122, pp. 200-209. DOI: 10.20310/1810-0198-2018-23-122-200-209 (In Russian, Abstr. in Engl.).