5. Малыгина В.В. Некоторые признаки устойчивости функционально-дифференциальных уравнений, разрешённых относительно производной // Изв. вузов. Математика. 1992. № 7. С. 46-53.
6. Yoneyama T. The 3/2 stability theorem for one-dimensional delay-differential equations with unbounded delay //J. Math. Anal. and Appl. 1992. № 165. P. 133-143.
7. Азбелев Н.В., Максимов В.П., Рахматуллина Л.Ф. Введение в теорию функционально-дифференциальных уравнений. М.: Наука, 1991.
8. Малыгина В.В. Об устойчивости решений некоторых линейных дифференциальных уравнений с последействием // Изв. вузов. Математика. 1993. № 5. С. 72-85.
БЛАГОДАРНОСТИ: Работа выполнена в рамках госзадания Минобрнауки РФ (задание №2014/152, проект №1890) и при финансовой поддержке Российского Фонда Фундаментальных Исследований (грант 13-01-96050 р_урал_а).
Поступила в редакцию 1 июня 2015 г.
Malygina V.V. ON THE ROBUST STABILITY OF NONAUTONOMOUS FUNCTIONAL-DIFFERENTIAL EQUATIONS
Effective conditions of robust stability are obtained for a linear nonautonomous differential equation with several delays. As an example, the equation with two delays is considered, whose domain of robust stability is constructed in the parameter space of the problem.
Key words: functional-differential equation; robust stability; effective conditions.
Малыгина Вера Владимировна, Пермский национальный исследовательский политехнический университет, г. Пермь, Россия, кандидат физико-математических наук, ведущий научный сотрудник научно-исследовательского центра «Функционально-дифференциальные уравнения», e-mail: [email protected]
Malygina Vera Vladimirovna, Perm State National Research University, Perm, the Russian Federation, Candidate of Physics and Mathematics, Leading Researcher of the Research Center «Functional-Differential Equations», e-mail: [email protected]
УДК 517.977
ЭЛЛИПСОИДАЛЬНОЕ ОЦЕНИВАНИЕ ТРУБОК ТРАЕКТОРИЙ УПРАВЛЯЕМЫХ ИМПУЛЬСНЫХ СИСТЕМ С КОНИЧЕСКИМ ОГРАНИЧЕНИЕМ НА УПРАВЛЕНИЕ
© О.Г. Матвийчук
Ключевые слова: множества достижимости; импульсное управление; эллипсоидальное исчисление.
В работе рассматриваются задача оценивания множеств достижимости управляемых динамических систем, описываемых дифференциальными уравнениями с импульсным управлением и неполной информацией о начальных данных. Предполагается, что импульсные управления в динамической системе находятся в пересечении специального конуса и обобщенного эллипсоида в пространстве функций ограниченной вариации. Целью данной работы является исследование свойств множеств достижимости управляемой системы заданного вида и построение алгоритма внутреннего эллипсоидального оценивания траекторий динамических систем импульсного типа.
Рассматривается задача импульсного управления трубкой траекторий линейной дифференциальной системы
^ж = + ^и(г), х € МП, £ € [¿о, Т] (1)
с неизвестным, но ограниченным начальным состоянием ж(£о — 0) = жо € Хо = Е(г, К) , где Е(г, К) = {ж € Мп| (ж — г)'К-1(ж — г) ^ 1} — эллипсоид в Мп с центром г € Мп и симметричной положительно определенной матрицей К. Здесь А(£) — непрерывная матричная функция размерности п х п, и(-) € у — импульсное управление, где у обозначает пространство п-векторных функций ограниченной вариации на [¿о, Т] [1—2]. В отличие от известных постановок, на управляющие импульсные воздействия и(£) , помимо обычного требования ограниченности вариации, наложено специальное ограничение и(£) € € и = Е * П К * , где Е * — обобщенный «эллипсоид» в пространстве у :
т
Е* = {и(-) € УП| У у(*)'йи(*) < 1, Уу(-) € СП, у(*) € Ео = Е(0,ф-1), ]}
¿0
и К * — обобщенный конус в пространстве у:
т
К * = {и(-) € УП | У у(*)'йи(*) ^ 0, Уу(-) € СП, у(-) € Ко, Т]},
¿0
Ко = {и € МП | и = (и1, ...,иП), и1 ^ 0, и € М, г = 1,...,п}.
Здесь матрица фо является симметричной положительно определенной матрицей, С^ обозначает пространство непрерывных п-векторных функций на [¿о,Т] с нормой ||у(-)||те,д = = шах ШНд (1 < 9 < ^ , 9 = (1 — Р-1)-1).
При ограничении и(£) € и = Е * П К * векторы скачков обобщенных управлений Ди = и(^+1) — и(^) обязаны лежать в пересечении Е(0, фо) П Ко , т. е.
Ди € Е(0,фо) П Ко = {г € МП| ¿'д-1,? < 1, г = (21, ..., 2„), 21 ^ 0, г* € М, 2 = 1, ..., п}.
Предполагается, что заданы фазовые ограничения в моменты ¿у
ж (¿у) € Е (у ), ^ : ¿о <¿1 < ••• <4 = т, ; = 0,..., Л, (2)
где Е(у-, ^) — заданные эллипсоиды с симметричными положительно определенными матрицами ^ и центрами у- € Мп .
Предполагается, что для рассматриваемой системы существует по крайней мере одна траектория ж(£) = ж(£; и(-), жо) , удовлетворяющая условиям жо € Хо , и(-) € и и фазовым ограничениям.
Трубку всех возможных траекторий системы (1) из начального состояния {¿о,Хо} при ограничениях жо € Хо , и(-) € и и фазовых ограничениях (2) далее будем обозначать символом
X(■) = X(■;и,Хо) =и {ж(•; и(-),жо) | жо € Хо, и(-) € и}.
Отметим, что сечение трубки траекторий X(■) в момент времени Т совпадает с множеством достижимости X (Т; и, Хо) системы (1) в момент Т , построенного из начального состояния {¿о,Хо} .
Целью данной работы является исследование свойств множеств достижимости управляемой системы (1) при заданных ограничениях и построение алгоритма внутреннего эллипсоидального оценивания ее множеств достижимости.
Теорема1. [3] Справедливо следующее равенство X(T; U,E(r,R)) = E(r*,R*) + X(T;U, {0}), r* = $(T,ío)r, R* = Ф(Т,*о)Д(Ф(Г, ío))',
$(í) — фундаментальная матрица уравнения z = A(t)z ( Ф(0) = I), Ф(^ т) = Ф(^Ф-1(т).
Согласно Теореме 1 эллипсоидальная оценка множества достижимости X(T; U, Xo) может быть найдена с помощью формул эллипсоидального оценивания суммы двух эллипсоидов E(r*,R*) и эллипсоида, оценивающего множество достижимости X(T;U, {0}) . Поэтому далее полагаем, что Xo = {0} .
Опишем кратко полученный алгоритм внутреннего эллипсоидального оценивания. Алгоритм. Рассмотрим разбиение {[tj-1, tj ]} : t0<t1< ... <tj-1<tj < ••• < tk=T заданного отрезка [t0, T] , x(tj) € E(yj, Dj) (2), j = 0,..., k .
1. Находим эллипсоид E(a, Q) С E(0, Qo) П Ко в соответствии с Теоремой 2. Теорема2. [4] Справедлива следующая эллипсоидальная оценка
n 1
E(a, Q) С E(0, Qo) П Ко, a = -— ■ Qov(v'QOv)-2,
n + 1
n + 1
Q-1 = (n + 1) ■ (vV) ■ (v'Qov)-1 + ■ Q-1,
где v — единичный вектор внутренней нормали к Ko (||v|| = 1) . Эллипсоид E(a, Q) имеет наибольший объем среди всех эллипсоидов, содержащихся в E(0, Qo) П Ko .
2. Рассматриваем систему (1) с начальным условием Xo = {0} и новым ограничением на управление u(-) € U = E*, где множество E* порождено эллипсоидом E(a, Q)
T
E* = {u(-) € \J y(t)'du(t) < 1 Vy(-) € СП y(t) € E(a, Q) Vt € [to, t1]}.
to
Л е м м а 1. Справедливо следующее включение X(t; U, {0}) С X(t; U, {0}).
3. Строим внутреннюю эллипсоидальную оценку множества достижимости X(t; U, {0}). Здесь, в отличие от задач оценивания состояний динамических систем с управлениями
классического типа, возникает задача внутреннего эллипсоидального оценивания выпуклой оболочки объединения семейства эллипсоидов [3-5].
Т е о р е м а 3. [3] Для системы (1) для любого е > 0 существуют 5 > 0 и конечное множество T¡ = {т1, т2, ..., тт} С [to, t1] такие, что
co( U E(aTi, QTi))CX(t; Ü, {0}) С co( U E(aTi, QTi)) + B(0, е),
Ti€Ts TiGTg
aTi = Ti)a + E(r*, R*), QTi = Ф(T, t*)Q^(T, t¿))'.
Здесь B(0, е) — е -шар в пространстве Rn с центром в начале координат.
Известно, что операции пересечения и объединения для пересекающихся выпуклых множеств являются двойственными операциями в том смысле, что одна из них превращается в другую при переходе от исходных множеств к их полярам [6]. Поэтому перейдем от E(aTi, QTi) к полярам E*(aTi, QTi) и построим эллипсоидальную оценку для
Р| E*(aTi, QTi) С Ё +(a, Q). Здесь используются известные результаты эллипсоидального исчисления полученные академиками А.Б. Куржанским [7] и Ф.Л. Черноусько [8] и их учениками.
Теорема4. Справедлива следующая эллипсоидальная оценка
E-(a-, Q-) = (Ё +(a, ())* С X(t; ¿¡Г, {0}) С X(t; U, {0}).
4. Находим E-(a-,Q-) С E(yi,Di)P|ЕЁ-(a-,Q-) , используя стандартную процедуру внутреннего эллипсоидального оценивания пересечения эллипсоидов [3], [4].
5. Повторяем шаги 2-4 для временных интервалов [tj-i, tj ] , j = 2,...,k, где U порождено эллипсоидом Е-(a-, Q-) , полученным на шаге 4.
В результате применения алгоритма получаем внутреннюю эллипсоидальную оценку трубки траекторий X(t; Ы, Х0), а, следовательно, и X(t; U, Х0).
Разработанный алгоритм оценивания проиллюстрирован модельными примерами.
ЛИТЕРАТУРА
1. Куржанский А.Б. Управление и наблюдение в условиях неопределенности. М.: Наука, 1977.
2. Красовский Н.Н. Теория управления движением. М.: Наука, 1968.
3. Вздорнова О.Г., Филиппова Т.Ф. Внешние эллипсоидальные оценки множеств достижимости дифференциальных импульсных систем // Известия РАН: Теория и системы управления. 2006. № 1. С. 38-47.
4. Matviychuk O.G., Internal Ellipsoidal Estimates of Reachable Set of Impulsive Control Systems under Ellipsoidal State Bounds and with Cone Constraint on the Control // Lecture Notes in Computer Science. 2014. № 8353, P. 125-133.
5. Matviychuk O.G., Estimation Problem for Impulsive Control Systems under Ellipsoidal State Bounds and with Cone Constraint on the Control, // AIP Conf. Proc. 2012. № 1497. P. 3-12.
6. Vazhentsev A.Y., External Ellipsoidal Estimation of the Union of Two Concentric Ellipsoids an Its Application // Computation Mathematics and Modelling. 2004. V. 15. № 2. P. 110-122.
7. Kurzhanski A.B., Valyi I. Ellipsoidal calculus for estimation and control. Laxenburg: IIASA Birkhauser, 1997.
8. Черноусько Ф.Л. Оценивание фазового состояния динамических систем: метод эллипсоидов. М.: Наука, 1988.
БЛАГОДАРНОСТИ: Работа выполнена при финансовой поддержке РФФИ (грант №15-01-02368a) и Программы фундаментальных исследований УрО РАН (проект № 15-16-1-8).
Поступила в редакцию 5 мая 2015 г.
Matviychuk O.G. THE ELLIPSOIDAL ESTIMATION OF TRAJECTORY TUBES OF IMPULSIVE CONTROL SYSTEMS WITH THE CONE CONSTRAINT ON THE CONTROL
The problem of estimating reachable sets of dynamical control systems with impulsive control and uncertainty in initial data is considered. The impulsive control in the dynamical system belongs to the intersection of the special cone with a generalized ellipsoid both taken in the space of functions of bounded variation. Properties of reachable sets of the control system and algorithm of construction of the internal ellipsoidal estimate of trajectory tubes of impulsive systems are studied.
Key words: reachable set; impulse control; ellipsoidal calculus.
Матвийчук Оксана Георгиевна, Институт математики и механики Уральского отделения Российской академии наук, г. Екатеринбург, Российская Федерация, кандидат физико-математических наук, научный сотрудник отдела оптимального управления, e-mail: [email protected]
Matviychuk Oksana Georgievna, Institute for Mathematics and Mechanics of the Ural branch of the Russian Academy of Sciences, Ekaterinburg, the Russian Federation, Candidate of Physics and Mathematics, Researcher of the Optimal Control Department, e-mail: [email protected]