Вычисление оптимального гарантированного результата в системах с запаздыванием в управлении Текст научной статьи по специальности «Математика»

УДК 517.917

© М. И. Гомоюнов, Н. Ю. Лукоянов

ВЫЧИСЛЕНИЕ ОПТИМАЛЬНОГО ГАРАНТИРОВАННОГО РЕЗУЛЬТАТА В СИСТЕМАХ С ЗАПАЗДЫВАНИЕМ В УПРАВЛЕНИИ1

Дана процедура вычисления оптимального гарантированного результата в системе с запаздыванием в управлении при показателе качества, представляющем собой евклидову норму совокупности отклонений движения системы в заданные моменты времени от начала координат.

Ключевые слова: оптимальное управление, дифференциальные игры, запаздывание.

§ 1. Постановка задачи

В рамках теоретико-игрового подхода [1,2] рассматривается задача о вычислении величины оптимального гарантированного результата управления в системе

x(t) = A(t)x(t) + B(t)u(t) + Bh(t)u(t — h) + C(t)v(t), to ^ t* ^ t ^ §, (1)

x € Rn, u € P С Rr, v € Q С Rs, h = const > 0, при начальных условиях

x(i*) = x* € Rn, Utt (■) = {uit (£) = u(i* + £), £ € [—h, 0)} = p*() € U, (2)

и показателе качества

7 = ^||£)[Ma;(tM)||2 + ... + \\DWx(tW)\\2. (3)

Здесь x — фазовый вектор, X(i) = dx(t)/dt, u — управляющее воздействие, v — воздействие

помехи; to и § — зафиксированы; P и Q — компакты; матрицы-функции A(i), B(t), Bh(t) и C(t)

кусочно-непрерывны; h — запаздывание; U — множество измеримых (по Борелю) функций из [—h, 0) в P; DM — постоянные (p[i хп)-матрицы (1 ^ p[i ^ n), t[i] € [t0, §] — оценочные моменты, г = 1, iV; ¿1г+11 = ¿W + h, i = 1, N — 1, ¿W < to + h, = $; /¿(t) = min {г = 1, iV : ^ t},

h* = h(t*); || ■ || — евклидова норма. Цель оптимизации - минимизация показателя (3).

Под стратегией управления понимаем произвольную функцию U = U(t, x,p(-),e) € P, где t € [to, §), x € Rn, p(-) € U и e > 0 — параметр точности [1]. Стратегия действует в дискретной по времени схеме на базе разбиения As = {tj : т\ = t*, 0 < Tj+1 — Tj ^ 5, j= 1, k, \ = $}, формируя управляющее воздействие согласно правилу

u(t) = U (Tj ,x(Tj ),uTj (0,e) t € [Tj,Tj+1), uTj (•) = {uTj (£) = u(Tj + £ € [ — h, 0)}. (4)

Пусть S = S(t*,x*,p*(),U(■), Дг,e) — множество троек {x(-), u(-), v(-)} таких, что v(-) — измеримая функция из [t*, §) в Q, u(-) — удовлетворяющая (2) измеримая функция из [t* — h, §) в P, которая на [t*, §) формируется в согласии с (4), x(-) — удовлетворяющая (2) абсолютно непрерывная функция из [t*,§] в Rn, которая вместе с u(-) и v(-) удовлетворяет уравнению (1) при почти всех t € [t*,§]. Оптимальным гарантированным результатом управления будет

Г°(i*, ж*,£>*(•)) = inf lim lim sup {7 I {x(-),u(-),v(-)} € S').

“v 7 u(0e^oг^о д 1 J

В случае, когда нет запаздывания в управлении, эффективные процедуры вычисления оптимального гарантированного результата в подобных задачах даны, например, в [1,3].

§2. Процедура вычисления оптимального гарантированного результата

Пусть X(t, т) — матрица Коши для уравнения x(t) = A(t)x(t). Пусть разбиение Дг содержит все точки разрыва матриц-функций A(t), B(t), Bh(t + h), С it) и моменты времени i = h*,N. Для j = 1,к обозначим h{Tj — 0) = h{rj); а также h{Tj + 0) = h{Tj) + 1, если Tj = и

h(Tj + 0) = h(Tj), если Tj < t[h(Tj)]. Положим h(Tk+1 + 0) = h(Tk+1 — 0) = h(Tk+1).

хРабота выполнена в рамках программы Президиума РАН «Динамические системы и теория управления», при финансовой поддержке УрО РАН (проект 12-П-1-1002), а также при поддержке РФФИ (грант № 12—01— 00290—а) и Программы государственной поддержки ведущих научных школ (НШ-5927.2012.1).

Попятно по шагам разбиения построим последовательность множеств С/ (¿*, т/ ± 0) С

Мр[ ( J 11 х К" и функций ± 0,1, т), (I, т) € т? ± 0), } = 1, к + 1.

При j = к + 1 полагаем

С&+1(£*,Тк+1 ± 0) = {(1,т) : ||1|| ^ 1, т = 0}, <^+1^*,Тк+1 ± 0,1, т) = 0.

Пусть 1 ^ j ^ к и уже построены С-/+1(£*, т/+1 ± 0) и ^/+1 (^*, т?+1 ± 0,1, т). Тогда определяем

С/(£*,т.,' +0) = С/+1(£*,т/+1 - 0) (£*,т.,' +0,1,т) = {^'(^ 0}*а,{и>Т. +0),

где символ |^(-)}с означает выпуклую сверху оболочку функции ^(1, т) на множестве С и

•0j(t*, l, m) = +i(t*, Tj+i—0, l, m)+ i maxmi^/l,D[h(rj +0)]X(t[h(rj +0)], t) (B(t)u+C(t)v)\ +

JT. veQ «еР у \ /

+ ^m,X($, t) ^(B(t) + X(t,t + h)Bh(t + h))u + C(t)vj^dT, (l, m) € Gj(t*,Tj + 0). Далее, если h(Tj — 0) = h(Tj +0), то полагаем

Gj (t*, Tj — 0) = Gj(t*, Tj + 0), (t*, Tj — 0, l, m) = ^j(t*, Tj + 0, l, m).

Jh(T7 -°)]

Если же h(Tj — 0) + 1 = h(Tj + 0), то, рассматривая для (l,m) € Rp J x Rn множество Yj(l,m) = {(v*, (l*,m*)) € [0,1] x Gj(t*,Tj + 0) : ||l||2 ^ 1 — v^, m = v*(m* + XT (#(r? +0)],^)D[h(TJ +0)1T l*)} , определяем

Gj(t*,Tj — 0) = {(l,m): Yj-(l,m) = 0}, ^j(t*,Tj — 0,l,m)= max v*^j(t*,Tj + 0,l*,m*).

(v* ,1* ,m*) eYj (1,m)

Рассмотрим величину

e(t*,x* ,p* (-);Дг) = max f^i (t* ,Ti — 0, l,m) + (l,D[h*] (X (t[h*],t*)x* +

(¡,m)eGi(t*,ri-0) V N 4

/* t[h*] — h г t* .

+ / X (t[h*] ,t + h)Bh (t + h)u(T)dr)) + (m,X (#,t* )x* + / X (#,t + h)Bh (t + Л,)и(т)^т)).

t* — h t* — h

Теорема 1. Для любого числа ( > 0 найдется число ¿(Z) > 0, такое, что, каковы

бы ни были 0 < £ ^ ¿(Z) и разбиение Д<$, будет справедливо неравенство

|r«(t*,x*,p*(-)) — e(t*,x*,p*(-); ДЙ) | ^ Z.

Список литературы

1. Krasovskii A.N., Krasovskii N.N. Control under Lack of Information. Berlin etc.: Birkhauser, 1995. 322 p.

2. Осипов Ю.С., Пименов В.Г. К теории дифференциальных игр в системах с последействием // Прикладная

математика и механика. 1978. Т. 42. Вып. 6. С. 963-977.

3. Лукоянов Н.Ю. К вопросу вычисления цены дифференциальной игры для позиционного функционала // Прикладная математика и механика. 1998. Т. 62. Вып. 2. С. 188-198.

Поступила в редакцию 14.02.2012

M. I. Gomoyunov, N. Yu. Lukoyanov

Calculating the optimal guaranteed result in systems with control delays

The procedure of calculating the optimal guaranteed result is obtained for systems with control delays and the quality index, which is the Euclidian norm of the set of deviations of the system motion at the given instants from the origin. Keywords: optimal control, differential games, delay.

Mathematical Subject Classifications: 49N70, 93B52, 49J35

Гомоюнов Михаил Игоревич, программист, Институт математики и механики УрО РАН, 620990, Россия, г. Екатеринбург, ул. С. Ковалевской, 16. E-mail: [email protected]

Лукоянов Николай Юрьевич, д.ф.-м.н., зав. сектором, Институт математики и механики УрО РАН, 620990, Россия, г. Екатеринбург, ул. С. Ковалевской, 16. E-mail: [email protected]

Gomoyunov Mikhail Igorevich, Programmer, Institute of Mathematics and Mechanics, Ural Branch of the Russian Academy of Sciences, ul. S. Kovalevskoi, 16, Yekaterinburg, 620990, Russia

Lukoyanov Nikolai Yur’evich, Doctor of Physics and Mathematics, Head of Department, Institute of Mathematics and Mechanics, Ural Branch of the Russian Academy of Sciences, ul. S. Kovalevskoi, 16, Yekaterinburg, 620990, Russia