НАУЧНЫЕ СТАТЬИ SCIENCE ARTICLES
ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ PHYSICAL AND MATHEMATICAL SCIENCES
МАТЕМАТИКА MATHEMATICS
УДК 517.977
В. В. Горячкин1, В. В. Крахотко1, Н. И. Широканова1, В. В. Игнатенко2
1 Белорусский государственный университет 2Белорусский государственный технологический университет
УПРАВЛЕНИЕ АНСАМБЛЕМ ЛИНЕЙНЫХ СИСТЕМ С ЗАПАЗДЫВАНИЕМ ПО УПРАВЛЕНИЮ В КЛАССЕ КУСОЧНО-ПОСТОЯННЫХ ФУНКЦИЙ
Для линейной системы с запаздыванием по управлению с постоянными интервальными коэффициентами (ансамбля систем) получены внешние оценки решений. Исследована задача управляемости ансамблем - перевод пучка траекторий системы из одного множества из Rn в другое за конечное время в классе кусочно-постоянных управлений. В случае интервальной неопределенности нахождение управления сведено к решению задачи линейного программирования, сформулированной по коэффициентам систем ансамбля. Доказаны условия управляемости ансамблем, которые представлены в форме оптимизационной задачи.
Ключевые слова: система дифференциальных уравнений, управляемость, интервальный анализ, линейная система с запаздыванием.
Для цитирования: Горячкин В. В., Крахотко В. В., Широканова Н. И., Игнатенко В. В. Управление ансамблем линейных систем с запаздыванием по управлению в классе кусочно-постоянных функций // Труды БГТУ. Сер. 3, Физико-математические науки и информатика. 2021. № 1 (242). С. 5-8.
V. V. Goryachkin1, V. V. Krakhotko1, N. I. Shirokanova 1, V. V. Ignatenko2
'Belarusian State University 2Belarusian State Technological University
CONTROL OF AN ENSEMBLE OF LINEAR SYSTEMS WITH A DELAY IN CONTROL IN THE CLASS OF PIECEWISE CONSTANT FUNCTIONS
For a linear system with a delay in control with constant interval coefficients (ensemble of systems), external estimates of solutions are obtained. The task of controllability of an ensemble is studied -the translation of a bundle of system trajectories from one Rn set to another in a finite time in the class of piecewise constant controls. In the case of interval uncertainty the control finding is reduced to solving
a linear programming problem formulated due to the ensemble systems coefficients Conditions for the controllability of an ensemble, which are presented in the form of an optimization problem, are proved.
Key words, system of differential equations, controllability, interval analysis, linear system with delay.
For citation: Goryachkin V. V., Krakhotko V. V., Shirokanova N. I., Ignatenko V. V. Control of an ensemble of linear systems with a delay in control in the class of piece-wise constant fnctios. Proseedings of BSTU, issue 3, Physics and Mathematics. Informatics, 2021, no. 1 (242), pp. 5-8 (In Russian).
Введение. Во многих задачах теории оптимального управления, теории систем параметры исследуемых уравнений известны лишь с некоторой точностью, а именно известны лишь интервалы, в которых они принимают значения. В этом плане интервальные модели можно рассматривать как эффективное средство описания широкого круга явлений с неопределенностями. Интервальная модель представляется как континуум моделей с параметрами, принимающими значения из допустимых интервалов. Подход исследования таких интервальных моделей состоит в нахождении одного приемлемого детерминированного решения для всего континуума моделей. Изучение интервальных моделей важно и перспективно для многих прикладных задач, таких как задачи управляемости, стабилизации и т. п. Это приводит к необходимости исследования динамических систем с интервальными коэффициентами [1-6] на управляемость.
В данной статье исследуется управляемость ансамбля систем с запаздыванием по управлению, для которых получены достаточные условия управляемости.
Основная часть. Рассмотрим систему
х(г) = Ах(г) + В1и (г) + В2и (г - к), х (0) = х*, е X*, х(*) = х*е X*, (1) и(-) = {и()г0,ге [-к,0]}.
Здесь хе Я", ие Яг, X*,X* с Я"; Ае ЯпХп, В1, В2 е ЯпХг - неопределенные матрицы со значениями в замкнутых интервалах
А < А < А, В1 < В1 < Д, В2 < В2 < В2; (2)
г* — фиксированный момент времени; к > 0 -число (запаздывание). Матричные и векторные неравенства следует понимать поэлементно. Множества
X*={хе Я": х*< х* < х*},
X*={хе Я": х*< х* <х*}-
это заданные параллепипедные множества (брусы в Я") [2].
Матрицы А, В1, В2, удовлетворяющие (2), будем называть допустимыми.
Ансамблем систем будем называть совокупность систем вида (1), коэффициенты которых -
допустимые матрицы, принимающие значения в интервалах (2) независимо друг от друга. Введем кусочно-постоянные управления
i(t )=,
, Tk_ < t <Tk, k = 1,2,
(3)
с фиксированными моментами разрыва первого рода 0 = т0 <т <... <т = г *, векторы и1, и2,..., и" принадлежат множеству и с Яг, которое определяется конечной системой линейных неравенств.
Обозначим через X(г)е Я" множество фазовых состояний х(г) ансамбля (1) для некоторого кусочно-постоянного управления и (г) в момент времени г, 0 < г < г *, с начальным состоянием х* = х(0)е X* при любых допустимых матрицах А, В1, В2.
Ставится следующая задача: для любого начального состояния из множества X* определить кусочно-постоянное управление (3) (одно и то же для всех систем ансамбля), при котором выполняется включение X(г*) с X*. Если поставленная задача разрешима в этом смысле, то будем говорить об управляемости ансамбля. Очевидно, что поставленная задача имеет решение лишь в исключительных случаях.
Рассмотрим систему (1) с допустимыми матрицами А, В1, В2, условиями х* е X*, х* е X* х* е X * и произвольным кусочно-постоянным управлением и (г), г е |0, г * ^ . По формуле Коши имеем
x (t *) = F (t * ,0)x0 + J F (t * ,т)х
0
x [ Bu (t) + B2u(t_ h)] dT = F(t* ,o)x0
s ( Tk
+Z J [F (t*, t) + F (t *, t + h )B2]dT
+
V Tk _1
где
F (t *, t)
- фундаментальная матрица однородной системы удовлетворяет уравнению
dF (t *, t)
дт
= _F (t *, т)А
с начальными условиями
F ( t*, t *- o) = En, F (t *, t) = 0, T>
t*+ 0.
s
*
В. В. Горячкин, В. В. Крахотко, Н. И. Широканова, В. В. Игнатенко
7
Обозначим
d = F,0)х* ,Ck = f [F(,1)5,
Tk
+
Tk-1 Tk
+
F(,T + h)B2]dт = f F(,t)t(B
Tk-1
Tk _
+ e-AhB2) = f F(,т)TB,k = 1,5.
+
Tk-1
Тогда очевидно, что множество векторов х (*) ансамбля систем (1) в момент времени ^ = I , отвечающее всем допустимым А, В1, В2 и начальному состоянию х*, образует множество X ( *).
Можно показать, используя процедуры оценивания произведения интервальных матриц [5, с. 68], что матрица В = В1 + е В2 принимает значение из замкнутого интервала [В0 -ДВ, В0 +ДВ] (задан в симметричной форме) и интервальные оценки вектора й и матриц Ск имеют вид
й<й<й,Ск <Ск <Ск,к = 1,2,...,5, (4)
где й,й,Ск,Ск определяются по формулам (17)-(20) из работы [4].
Матричные интервалы (4) можно также представить в симметричной форме [2]:
[Ск ] = [Ск, Ск ] = [Сок-ДСк, Сок +ДСк ],
к = 1,2,..., 5,
где Сок = (( + Ск)/2- центр и ДСк= (( - Ск) / 2-
радиус интервальной матрицы [Ск ].
Тогда внешняя интервальная оценка [6, с. 55] X (*) примет ]
множества
' вид
d + £ ( U-ACk\uk\ )< х <
k=1
< d + £ (uk + ACk | U
k=1
(5)
где и понимается как вектор модулей коор-I I к
динат вектора и .
Пусть е> 0,ее Я" -любой вектор. Обозначим X* = {хе Я": х* - е < х < х* + е} замкнутую
е -окрестность бруса X *.
Из полученной оценки (5), следует, что X ( *) с X* имеет место тогда, когда
d + £ {Cokuk-ACk|uk |) > х* -£,
i=1
d + £ (C0kuk + ACk | uk |) < х* + e.
(6)
Составим по условиям (6) задачу нелинейного программирования:
e + e 2
1 2
en ^ min,
d + £ (
k=1
d + £ (C(
-ACk
D> Х-e,
(7)
0kM + ACk
) X
е> 0,ик е и, к = 1,5.
Таким образом, в задаче (7) ищется управление (3) , для которого внешняя оценка множества X ( *) находится в минимальной е-окрестнос-ти бруса X * (под минимальностью е-окрестности бруса понимается окрестность бруса с минималь-
ной нормой вектора
e :1 N1=£ ei.
i=1
Очевидно, задачу (7) можно переписать в виде следующей задачи линейного программирования:
e' - e ^ min,
-£ (CokUk-ACkwk)-e< d -х*,
k=1
_ 5 _
d + £ (C0kuk +ACkwk)-e<d + х*,
k=1
e> 0,-wk <uk < wk,uk e U,wk > 0,
k = 1,5, e = (1,1,...,1),
(8)
где неизвестными будут е и ик, , к = 1,5.
Ясно, что в задаче (8) множество планов не пусто и целевая функция ограничена на нем снизу - значит, решение задачи (8) существует.
Следовательно, откорректированная постановка задачи управляемости ансамбля (1) звучит так: для любого начального состояния из множества X* определить кусочно-постоянное управление (одно и то же для всех систем ансамбля) и минимальный е-брус множества X*, что выполняется включением X (*) с X*.
Поэтому справедлива теорема.
Теорема. Для управляемости ансамбля систем (1) в классе кусочно-постоянных функций достаточно, чтобы разрешимая задача линейного программирования (8) имела оптимальный план (е0,ик0, ^к0, к = 1,5). В этом случае соответствующее управление переводит пучок траекторий системы из X* в минимальную е0 -окрестность бруса X*. Если в решении е0 = 0 , то задача решается точно.
Заключение. Полученные достаточные условия управляемости интервальной системы с запаздыванием по управлению гарантируют попадание пучка траекторий системы (1) в е0 -окрестность бруса X * . Этот результат можно перенести на более сложные объекты, например, интервальные системы с запаздыванием по состоянию.
Список литературы
1. Куржанский А. Б. Управление и наблюдение в условиях неопределенности. М.: Наука, 1977. С. 392.
2. Алефельд Г., Херцбергер Ю. Введение в интервальные вычисления. М.: Мир, 1987. С. 356.
3. Гайшун И. В., Горячкин В. В., Крахотко В. В. Оценка решений двухпараметрической дискретной системы с интервальными коэффициентами // Весщ НАН Беларуси Сер. ф1з.- мат. навук. 2014. № 3. С. 5-8.
4. Гайшун И.В., Горячкин В. В. Робастная и интервальная наблюдаемость двухпараметрических дискретных систем с интервальными коэффициентами // Весщ НАН Беларусь Сер. ф1з.- мат. навук. 2016. № 2. С. 6-9.
5. Ащепков Л. Т. Управляемость интервальной линейной дискретной системы // Известия РАН. Теория и системы управления. 2007. № 3. С. 67-74.
6. Ащепков Л. Т. Внешние оценки и ступенчатая управляемость интервальной линейной системы // Автоматика и телемеханика. 2008. № 4. С. 51-56.
References
1. Kurzhanski A. B. Upravleniye i nablydeniye v usloviyakh neopredelennosti [Control and observation under conditions of uncertainty]. Moscow, Nauka Publ., 1977. 392 p.
2. Alefel'd G., Herzberger Yu. Vvedeniye v interval'nyye vychisleniya [Introduction to interval calculations]. Moscow, Mir Publ., 1987. 356 p.
3. Gaishun I. V., Goryachkin V. V., Krakhotko V. V. Estimation of solutions of a two-parameter discrete system with interval coefficients. Vesti NAN Belarusi [Bulletin of the NAS of Belarus], series Physics and Mathematics, 2014, no. 3, pp. 5-8 (In Russian).
4. Gaishun I. V., Goryachkin V. V. Robust and interval observability of two-parameter discrete systems with interval coefficients. Vesti NAN Belarusi [Bulletin of the NAS of Belarus], series Physics and Mathematics, 2016, no. 2, pp. 6-9 (In Russian).
5. Aschepkov L. T. Controllability of an interval linear discrete system. Izvestiya RAN. Theorii i sistemy upravleniya [News RAS. Theory and control systems], 2007, no. 3, pp. 67-74 (In Russian).
6. Aschepkov L. T. External estimations and step control of an interval linear system. Avtomatika i telemekhanika [Automatics and telemechanics], 2008, no. 4, pp. 51-56 (In Russian).
Информация об авторах
Горячкин Владимир Викторович - кандидат физико-математических наук, доцент, доцент кафедры технологий программирования. Белорусский государственный университет (220030, г. Минск, пр-т Независимости, 4, Республика Белaрусь). E-mail: gorvv@bsu.by
Крахотко Валерий Васильевич - кандидат физико-математических наук, доцент, доцент кафедры методов оптимального управления. Белорусский государственный университет (220030, г. Минск, пр-т Независимости, 4, Республика Белaрусь). E-mail: krakhotko@bsu.by
Широканова Наталья Ивановна - кандидат физико-математических наук, доцент, доцент кафедры общей математики и информатики. Белорусский государственный университет (220030, г. Минск, пр-т Независимости, 4, Республика Белaрусь). E-mail: shirokanova@bsu.by
Игнатенко Василий Васильевич - кандидат физико-математических наук, доцент кафедры высшей математики. Белорусский государственный технологический университет (220006, г. Минск, ул. Свердлова, 13а, Республика Белaрусь). E-mail: ihnatsenko@tut.by
Information about the authors
Goryachkin Vladimir Viktorovich - PhD (Physics and Mathematics), Associate Professor, Assistant Professor, the Department of Programming Technologies. Belarusian State University (4, Nezavisimosti Ave., 220030, Minsk, Republic of Belarus). E-mail: gorvv@bsu.by
Krakhotko Valeriy Vasil'yevich - PhD (Physics and Mathematics), Associate Professor, Assistant Professor, the Department of Optimal Control Methods. Belarusian State University (4, Nezavisimosti Ave., 220030, Minsk, Republic of Belarus). E-mail: krakhotko@bsu.by
Shirokanova Natalia Ivanovna - PhD (Physics and Mathematics), Associate Professor, Assistant Professor, the Department of General Mathematics and Computer Science. Belarusian State University (4, Nezavisimosti Ave., 220030, Minsk, Republic of Belarus). E-mail: shirokanova@bsu.by
Ignatenko Vasiliy Vasil'yevich - PhD (Physics and Mathematics), Associate Professor, Assistant Professor, the Department of Higher Mathematics. Belarusian State Technological University (13a, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: ihnatsenko@tut.by
Поступила после доработки 30.11.2020