Известия Института математики и информатики УдГУ_2016. Вып. 1 (47)
УДК 517.958, 530.145.6 © Л. И. Родина, А. X. Хаммади
ХАРАКТЕРИСТИКИ ИНВАРИАНТНОСТИ МНОЖЕСТВА ДОСТИЖИМОСТИ УПРАВЛЯЕМОЙ СИСТЕМЫ1
Изучаются характеристики, связанные с инвариантностью или слабой инвариантностью заданного множества М = {(Ь,х) е [0, х К" : х е М(£)} относительно управляемой системы х = /(Ь,х,и) на конечном промежутке времени. Одной из таких характеристик является относительная частота Ьщ^ т+щ(Е,М) поглощения множества достижимости Б(Ь,Х) данной системы множеством М та отрезке [т,т + $], равная отношению меры Лебега тех Ь из [т,т + $], при которых Б(Ь,Х) С М(Ь), к длине данного отрезка. Другая характеристика, (Б,М) = М Т+,](Б,М) отображает свойство равномерности пребывания множества достижимо-
сти Б(Ь,Х) в мпожестве М та отрезке заданной длины Доказаны теоремы об оценке и вычислении этих характеристик для различных многозначных функций М(Ь) и Б(Ь,Х). В частности, получены равенства для нахождения &едт(Б, М) для функции М(Ь), периодической с периодом Т и функции Б(Ь, X), которая при всех Ь ^ 0 удовлетворяет включению Б(Ь + Т,Х) С Б(Ь,Х). Рассмотрены примеры вычисления и оценок данных характеристик.
Ключевые слова: управляемые системы, дифференциальные включения, множество достижимости.
Введение
Задача исследования инвариантности множеств относительно различных управляемых систем и дифференциальных включений является одной из важнейших задач математической теории управления и теории дифференциальных игр. Данной тематике посвящены работы H.H. Красовского и А. И. Субботина, А. Б. Куржанского и Т. Ф. Филипповой, Ж. П. Обена, Е. Л. Тонкова и Е. А. Панасенко, В.Н. Ушакова, Ф. Хартмана и многих других авторов.
Первый результат в этой области опубликован М. Нагумо в 1942 году, которым было изучено свойство слабой инвариантности заданного множества относительно дифференциального уравнения. Эти исследования продолжил Ф. Хартман [1], сформулировав необходимые и достаточные условия слабой инвариантности для системы дифференциальных уравнений. Большое внимание изучению вопросов инвариантности и «выживаемости» (как называют в иностранной литературе слабую инвариантность) уделял Ж. П. Обен [2], который получил условия выживаемости множества K С Rn относительно управляемой системы
x = f (x, u), U £ U(x)
в предположении, что для каждого x £ Rn множество f (x,U(x)) выпукло и компактно.
K
любой начальной точки xo £ K существует хотя бы одно решение x(t) этой системы, начинающееся в xo и выживаемое в том смысле, что x(t) £ K для всех t ^ 0.
Е. Л. Тонков и Е. А. Панасенко [3, 4] сформулировали условия, при которых заданное множество обладает свойствами положительной инвариантности, инвариантности, устойчивой или асимптотически устойчивой инвариантности относительно нестационарного дифференциального включения. Наличие инвариантного множества позволяет им рассматривать сужение включения на это множество и исследовать в нем различные экстремальные движения. А. Б. Куржанский исследовал структуру слабо инвариантных множеств гибридных систем, движение которых осуществляется путем мгновенного переключения с одной из «стандартных систем» на другую.
Отметим, что свойство слабой инвариантности находится в тесной связи с задачами о сближении управляемой системы с компактным целевым множеством, исследуемыми H.H. Кра-совским, А. И. Субботиным, В.Н. Ушаковым и многими другими авторами. При решении этих
1 Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (проект № 16-0Ю0346-а) и Министерства образования и науки Российской Федерации в рамках базовой части госзадания (проект № 2003).
задач возникает вопрос о том, в какой степени заданное множество не является инвариантным относительно дифференциального включения, порожденного управляемой системой (см. работы В. Н. Ушакова и его учеников [5,6]). Один из возможных подходов к решению этого вопроса состоит в применении инфинитезимального представления свойства инвариантности и вычисления дефекта инвариантности, который оценивает степень несогласованности множества и динамики системы с точки зрения понятия инвариантности.
В работах [7-11] также исследуются множества, не являющиеся инвариантными в «классическом» смысле; для таких множеств вводится естественное расширение понятия инвариантности, которое названо статистической инвариантностью. Пусть D(t, X) — множество достижимости управляемой системы
x = f (t,x,u), (t,x,u) G [0, +x>) x Rra x Rm, (0.1)
в момент времени t из начального множества X. Множество
M = {(t,x) G [0, +го) x Rra : x G M (t)}
называется статистически инвариантным относительно системы (0.1), если относительная частота пребывания множества достижимости D(t, X) в множестве M равна единице.
В данной работе изучаются характеристики, связанные с инвариантностью или слабой инвариантностью заданного множества M = {(t,x) G [0, +œ>) x Rra : x G M (t)} относительно управляемой системы (0.1) на конечном промежутке времени. Одной из таких характеристик является относительная частота freq^+^D, M) поглощения множества достижимости D(t, X) данной системы множеством M на отрезке [т, т+$], равная отношению меры Лебега тех t из [т,т + $], при которых D(t,X ) С M (t), к длине данного отрезка. Другая характеристика,
freq^(D,M) = inf freq[r r+^](D, M)
T ^ U
отображает свойство равномерности пребывания множества достижимости D(t, X) в множестве M на отрезке заданной длины Доказаны теоремы об оценке и вычислении этих характеристик для различных многозначных функций M (t) и D(t,X ). В частности, получены равенства для нахождения freqr (D,M ) для функции M (^периодической с периодом T и функции D(t,X), которая при всех t ^ 0 удовлетворяет включению D(t + T,X) С D(t,X). Рассмотрены примеры вычисления и оценок данных характеристик.
§ 1. Основные определения и свойства характеристик
Рассматривается управляемая система
x = f (t,x,u), (t, x, u) G [0, +œ) x Rra x Rm, (1.1)
где функция / (Ь, х, и) непрерывна по совокупности переменных, управление и содержится в компактном множестве и(Ь, х) С Мт, и функция и(Ь, х) полунепрерывна сверху в метрике Хаусдорфа при всех (Ь, х) е [0, х Rn.
Определение 1.1. Допустимым процессом, управляемой системы (1.1) назовем функцию Ь — (и(Ь),х(Ь)) е Rn х Rm, которая удовлетворяет следующим условиям:
1) управление и(Ь) определено для всех Ь ^ 0, ограничено и измеримо по Лебегу;
2) решение х(Ь) системы дифференциальных уравнений х = /(Ь,х,и(Ь)) в смысле Каратео-дори определено для всех Ь ^ 0;
3) имеет место включение и(Ь) е и(Ь,х(Ь)).
Рассмотрим отвечающее системе (1.1) дифференциальное включение
х е ^(Ь,х), (Ь,х) е [0, +то) х Rn, (1.2)
где для каждой фиксированной точки (Ь, х) е [0, х Rn множество F(Ь, х) состоит из всех предельных значений функции /(Ьг,х1,и(Ь1 ,х^)) при (Ьг,хг) — (Ь,х).
Пусть D(t,X) — множество достижимости системы (1.1) в момент времени t из начального множества X, то есть множество, состоящее из всех значений в момент t решений p(t, x) включения (1.2), когда начальное условие ^>(0,x) = x пробегает все множество X. Предполагаем, что для каждого X множество достижимости D(t, X) существует для всех t ^ 0. Это означает, что для каждой точки x G X существует решение ф(t,x) включения (1.2), удовлетворяющее начальному условию ^>(0,x) = x и продолжаемое на полуось R+ = [0, +œ>).
Пусть множество M = I (t, x) G [0, +то) x Rra : x G M (t)} задано функ цией t ^ M (t), непрерывной в метрике Хаусдорфа, и для каждого t G [0, +œ>) множество M (t) непусто и компактно. Для определения характеристик множества достижимости введем в рассмотрение множество a(r,ê,X) = {t G [т,т + tf] : D(t,X ) С M (t)}.
Определение 1.2 (см. [7,9]). Относительной част,от,ой поглощения множества достижимости D(t,X) системы (1.1) множеством M называется следующий предел
г . , mes a(0,tf,X) mes{t G [0,tf]: D(t,X) С M(t)} freqp.lí = lim -2—'- = lim -ï-——'—-1.3
Q^œ tf Q^œ tf
mes
„ ^ . -i—mes а(0, i?, X) „ ,. , mes a(0, tf, X)
îreq (D,M)= lim-—J-, freq*(£>, M) = lim-—'- (1.4)
!?-)■ oo tf ^oo tf
называются, соответственно верхней и нижней относительной част,от,ой поглощения множе-D(t, X) M.
Определение 1.3 (см. [10]). Относительной част,от,ой поглощения множества достижимости D(t,X) системы (1.1) множеством M на отрезке [т,т + tf] называется характеристика
,п . mes а(т^,и) mes {t G [т,т + tf] : D(t,X ) С M (t)} freq [T,T+0](D,M) =--- =---.
Важно рассматривать относительную частоту freq[r r+Q](D, M) для любого момента времени т ^ 0, поэтому естественно для заданного tf > 0 определить характеристику
. mes {t G [т,т + tf] : D(t,X ) С M (t)} freq*(A M) = inf freq[T>T+0](D, M) = inf-1-Lj-^ V' WJ.
Эта характеристика отличается от пределов (1.3), (1.4) тем, что она отражает свойство рав-
D(t, X) M
данной длины. Величина freqQ(D,M) равна минимальной доли времени, в течение которого D(t, X) M tf.
Теорема 1.1. Имеют место следующие свойства:
1) для любо го tf > 0 выполнено неравенство
freqQ(D, M) < freq,(D, M) ; (1.5)
2) если функ ции t ^ D(t, X ) и t ^ M (t) периодические с пер иодом, T > 0, m,о предел, freq(D, M) существует и
г mes {t G [0,T] : D(t,X ) С M (t)} freqT(A M) = freq(_D, M) =-1-1 J ^ ! ~-K-LL ;
3) если функ ция t ^ M (t) периодическая с пер иодом, T > 0 и для, вс ex t ^ 0 имеет место включение D(t + T,X) С D(t,X), то
freqT(A M) = ■«s{t6[0.r]:g(f,X)CAf(i)} (Le)
Доказательство. Для доказательства первого утверждения отметим, что неравенство (1.5) следует из равенства
к-1
Е шв8{Ь е [г0, (г + 1)0] : Б(Ь, X) С М(Ь)} йге^ДМ) = Ит ^-—-
к^-оо
и неравенства
те8{^€[г1?,(г + 1)1?]: Р&Х) С М{Ь)} ^ -^-^ йе^ДМ),
которое верно для любых г = 0,1,... и 0 > 0.
Докажем второе утверждение. Поскольку функции Ь — 0(Ь, X) и Ь — М(Ь) периодичны с общим периодом Т > 0 и 0 = Т, то выполнены следующие равенства
&е<1т(АМ) = ш * [т,т + Ц гДЛ*)£*(,)}
т > U T
mes {t G [0,T] : D(t,X ) С M (t)}
= freq(D,M ).
T
Последнее равенство получено в [9, с. 53-54].
Докажем третье утверждение. Пусть функция t ^ M (t) периодическая с периодом T и D(t + T, X) С D(t, X) для всех t ^ 0. Тогда если для некоторого t* ^ 0 имеет место включение D(t*,X) С M (t* ), то
D(t* + T, X) С D(t*,X) С M (t* ) = M (t* + T). Поэтому для любого т ^ 0 выполнено неравенство
mes {t G [0, т] : D(t, X ) С M (t)} < mes {t G [Г,т + T] : D(t, X ) С M (t)}. (1.7)
Пусть т ^ T, тогда
mes {t G [0, т] : D(t, X ) С M (t)} =
= mes {t G [0, T] : D(t, X ) С M (t)} + mes {t G [T, т] : D(t, X ) С M (t)},
mes {t G [T, т + T]: D(t, X ) С M (t)} =
= mes {t G [T, т] : D(t, X ) С M (t)} + mes {t G [т,т + T] : D(t, X ) С M (t)}.
Следовательно, из (1.7) получаем
mes {t G [0,T] : D(t,X ) С M (t)} < mes {t G [т,т + T] : D(t, X ) С M (t)}. (1.8)
Далее, если т < T, то отрезки [0,T] и [т,т + T] имеют непустое пересечение — отрезок [r,T}. В этом случае неравенство (1.8) получаем из неравенств
mes {t G [0,T] : D(t,X ) С M (t)} =
= mes {t G [0, т] : D(t, X ) С M (t)} + mes {t G [т, T] : D(t, X ) С M (t)} < < mes {t G [T, т + T]: D(t, X ) С M (t)} + mes {t G [т, T] : D(t, X ) С M (t)} =
= mes {t G [т,т + T]: D(t, X ) С M (t)}.
Таким образом, (1.8) верно для всех т ^ 0 и поэтому
inf mes {t G [т,т + T]: D(t, X ) С M (t)} =mes{t G [0, T]: D(t, X ) С M (t)}.
т ^ U
Последнее равенство равносильно (1.6). □
Следствие 1.1. Если M (t) = M для всех t ^ 0 и D(t + §,X ) С D(t, X ) для некоторого > 0 и всех t ^ 0, то
&е<1„(0,M) = £[»,<>] :Щ,Х) ЯМ)
Доказательство. Поскольку M (t + ê) = M (t) для всех t ^ 0, то данное утверждение
□
§ 2. Об оценке и вычислении относительных частот для некоторого класса многозначных функций
Рассмотрим множество M = {(t,x) Е [0, +œ>) х Rn : x Е M(t)}, заданное непрерывной многозначной функцией M(t), и многозначные функции D(t), D(t), также непрерывные в метрике Хаусдорфа. В частности, множество D(t) = D(t, X) может является множеством достижимости управляемой системы (1.1), а множество D(t) = D(t,X) — множеством достижимости управляемой системы
X = f(t,x,u), (t,x,u) Е [0, +ж) х Rn х Rm,
в момент времени t из некоторого начального множества X.
Предполагаем, что для каждого t ^ 0 множества M (t), D(t) и D(t) непустые, компактные и функции M (t), D(t) периодические с пер иодом T > 0 (под T-периодической функцией будем понимать функцию d(t), удовлетворяющую при всех t ^ 0 равенству d(t+T) = d(t)). Обозначим через dist(A, B) расстояние по Хаусдорфу между множествами A и B в пространстве Rn.
В следующей теореме получена оценка характеристики freqT(D, M) и приведены условия, при которых можно найти ее значение.
Теорема 2.1. Пусть функции M (t), D(t) периодические с периодом T > 0 и выполнено равенство lim dist(D(t), D(t)) = 0. Тогда имеют место следующие свойства,:
1) < freqT(ДМ) =
2) если D(t) Ç D(t) для всex t ^ 0, то freqT(D, M) = freqT(D, M). Доказательство. Для каждого т ^ 0 определим функцию
R(t ) = mes {t Е [т,т + T] : D(t) Ç M (t)} - mes {t Е [0, T] : D(t) Ç M (t)}. Для доказательства первого утверждения теоремы нужно показать, что inf R(t) ^ 0. По свой-
т
ствам меры Лебега (так как из а-аддитивности меры следует ее непрерывность [12, с. 274]) для любого е > 0 найдется такое значение £о = £о(^), что имеет место неравенство
mes {t Е [0,T] : D(t) Ç M£0 (t) \ M (t)} < е.
Далее, поскольку h(t) = dist(D(t), D(t)) 0 при t +œ>, для ео существует момент времени t0 такой, что h(t) ^ ео для всех t ^ to. Следовательно, для всех т ^ to выполнены неравенства
mes {t Е [т,т + T] : D(t) Ç M (t)} < mes {t Е [т,т + T] : D(t) Ç Mh(t) (t)} <
< mes {t Е [t, t + T]: D(t) Ç M£0 (t)}. (2.1)
Из (2.1) в силу периодичности функций M (t), D(t) получаем, что для всех т ^ to
R(t ) < mes{t Е [т,т + T]: D(t) Ç M£0 (t)} -mes{t Е [0,T]: D(t) Ç M (t)} =
= mes {t Е [0, T] : D(t) Ç M£0 (t)} - mes {t Е [0, T] : D(t) Ç M (t)} =
= mes {t Е [0, T] : D(t) Ç M£0 (t) \ M (t)} < е. Таким образом, имеет место неравенство inf R(t) ^ 0, из которого в силу определения функ-
т ^0
ции R(t) и относительных частот следует неравенство
freqT(D, M) ^ freqT(D, M).
Докажем второе утверждение теоремы. Если D(t) Ç D(t) для всех t ^ 0, то при всех т ^ 0 справедливо неравенство
mes {t Е [т,т + T] : D(t) Ç M (t)} ^ mes {t Е [т,т + T] : D(t) Ç M (t)} =
= mes {t Е [0, T] : D(t) Ç M (t)},
поэтому R(t) ^ 0 для всех т ^ 0 и inf R(t) ^ 0. Учитывая доказанное выше неравенство
т
inf R(t) ^ 0, получаем, что inf R(t) = 0, следовательно freqT(D,M) = freqT(D,M). □
т^0 " т^0
Рассмотрим множество M G R2 вида
M = {(t,x) G [0, +ж) х (—ж,с]}, с G R.
Ниже приведено следствие из теорем 1.1 и 2.1.
Следствие 2.1. Пусть функции p : R+ ^ R и p : R ^ R непрерывны, p(t) периодическая с периодом T > 0 и ^ lim (p(t) — p(t)) = 0. Тогда имеют место следующие свойства:
1W , , °°, л mes{¿ G [0, Т] :<p(t) ^ с}
1) freqT(<¿>, (-ж, с]) ^ freqT(<¿>, (-ж, с]) =---;
2) если p(t) ^ p(t) для всex t ^ 0, то freqT (pp, (—ж, с]) = freqT(р, (—ж, с]);
3) есл и p(t) > p(t) и p(t + T) < p(t) для вс ex t ^ 0, то
mes{¿ G [0, T] freqT(^, (-oo, c]) =---.
Доказательство. Отметим, что функции p(t) и p(t), удовлетворяющие условию следствия, ограничены, поэтому существует постоянная E > 0 такая, что |p(t)| ^ Ей \p(t)\ ^ E для всех t ^ 0. Для каждого t ^ 0 рассмотрим множества D(t) = [—Е, p(t)] и D(t) = [—Е, p(t)]. Тогда выполнены равенства
dist(D(t), D (t)) = |p(t) — p(t)|, freqT(p, (—ж,с]) =freqT (p, —Е,с]) = freqT (D, —Е,с}),
и неравенство p(t) ^ p(t) равносильно неравенству D(t) С D(t). Поэтому данное утверждение
□
Замечание 1. Если функция g(t) = p(t) — p(t) — убывающая на [0, +ж), то для любого T > 0 неравенство g(t + T) < g(t) (равносильное неравенству p(t + T) < p(t)) выполнено для всех t ^ 0. Обратное утверждение неверно. Например, функция g(t) = sin t + e-t не является убывающей, но g(t + 2п) < g(t) для всех t G [0, +ж) (для всех t G R).
Пример 2.1. Рассмотрим управляемую систему
x = —x + (cos t + 1)u, (t, x) G [0, +ж) x R,
где u G U = [—1,1], и найдем множество достижимости D(t,X) данной системы из начального множества X = [x\,x0] : D(t,X) = [xi(t, x\),x2(t, x°)] , где
Если X = X =
xl(t,x°1) = e~t(Kx01 + ^j + -1,
x2(t,x°2) = е~*(х°2 - + ^ sin(í + 1) + 1, t> 0.
3 3
2' 2 J
, то многозначная функция
t ^ D(t, X) = sin(í + - 1, ^ sin(í + £ ) + 1
периодическая с периодом 2n.
Пусть Ш = х) € [0, +оо) х —^ тогда freq27r(-D, М) = ^ в силу теоремы 1.1. Если
X С X, то В(Ь + 2п,Х) С Б(Ь,Х), поэтому (см. следствие 1.1)
шев{Ь е [0, 2п] : Б(Ь,Х) С М}
(Б,М) =
2п
Если X С Х, то Б(Ь,Х) С Б(Ь,Х) для всех Ь ^ 0. Далее, Иш Б(Ь,Х),Б(Ь,Х)) = 0, поэтому в силу теоремы 2.1 выполнено равенство
1гщ2ж{0,М) = 1гщ2ж{Ъ,М) =
§ 3. Характеристики инвариантности решений, возникающие в физических, химических, биологических и экономических процессах
Пусть г(Ь) — решение задачи Коши
г = г), г(0) = г0.
Во многих прикладных задачах величина г(Ь) не может принимать отрицательные значения, например, в физических процессах неотрицательными являются энергии частиц, в химических — концентрации реагирующих веществ, в биологических — размер популяции, в экономических — величины производств и цены на продукцию (соответствующие примеры приведены, в частности, в работах [13-16]). Поэтому для исследования этих задач введем следующие характеристики, определенные для любого с > 0 :
те8{* € [г, г + 0] : < с} й^Дг, (-оо,с|) = ---,
0 0 (3 1) те8{*€ [0,0] ^ с} (—оо, с ) = ит ------
(если последний предел существует).
Отметим, что если по содержанию задачи величина г(Ь) может быть отрицательной, то можно определить характеристики (3.1) для любого с е К. Также представляет интерес исследование характеристики
, шев{Ь е [т,т + 0] : г(Ь) е [сис]} йге^ (г, [а, с2]) = ---,
равной минимальной доли времени, в течение которого решение г(Ь) находится в множестве [с\, С2] на каждом временном отрезке длины 0.
Для управляемых процессов, которые имеют периодический характер, будем рассматривать характеристики (г, (—го,с]) и freq^(г, [с1,с2]) при 0 = Т, где Т > 0 — период данного процесса.
Пример 3.1. Получим равенства для вычисления характеристики freqy [г, (—го,с]), где г(Ь) — численность популяции, динамика которой задана задачей Коши
г = (е(Ь) — а(Ь)г)г, г(0) = го. (3.2)
Здесь го ^ 0, е(Ь) и а(Ь) — непрерывные периодические функции с периодом Т > 0. Предполагаем, что существуют постоянные а1, а^ и е\, £2 такие, что
0 < а1 < а(Ь) < а2, 0 < £1 < £(Ь) < £2. (3.3)
Несложно показать, что для любого начального размера популяции го > 0 с течением времени кривая численности популяции г(Ь) стремится к периодической кривой, которая задается следующей функцией:
z0 exp^ e(r)drj exp^ e(r)dr^J — 1 *(*) =-71--> W zo = -pr Jo
zoJ a(s)exp^J e(r)dr^Jds + 1 J a(s)exp^J e(r)dr^jds
/o vo
Действительно, решение задачи Коши (3.2) имеет вид:
rt
(4 o
Тогда z(t) — z(t) = h(t) ■ g(t), где
z0 ехщ / е(т)dr\ z{t) = J
z0 j a(s)exp^ J e(r)dr^Jds + 1
z z exp( e(r)dr) Kt) = —t-^-» 9® = —t-^--
z0 J a(s)exp^у е(т)dr^Jds + 1 z0 j a(s)exp^ J е(т)dr^Jds + 1
Из неравенств (3.3) несложно получить, что функция g(t) ограничена, а функция h(t) стремится к нулю при t ^ (если zo > 0). Поэтому, если zo > 0, то
lim (z(t) — z(t)) = 0. (3.4)
Из (3.3) также следует, что z0 > 0.
Пусть z0 £ (0, z0), тогда ^(t) = z(t) — z(t) ^ 0 для всех t ^ 0. Из (3.4) в силу следствия 2.1 получаем равенства
freqT(z, (-оо,с|) = freq(z, (-oo, c\) =---.
Пусть zo > zo. Несложно проверить, что неравенство h(t + T) ^ h(t) выполнено для всех t ^ 0. Учитывая периодичность и положительность функции g(t), получаем
p(t + T) = h(t + T)g(t + T) = h(t + T)g(t) < h(t)g(t) = <p(t)
для всех t ^ 0. Тогда в силу следствия 2.1 выполнено равенство
mes{te [0,Г] :z(t) ^ с} freqT(z, (-oo, cjj =---.
Список литературы
1. Хартман Ф. Обыкновенные дифференциальные уравнения. М.: Мир, 1970. 720 с.
2. Aubin J.-P. Viability theory. Boston-Basel-Berlin: Birkhauser, 1991. 543 p.
3. Панасенко E.A., Тонков Е.Л. Функции Ляпунова и положительно инвариантные множества дифференциальных включений // Дифференциальные уравнения. 2007. Т. 43. № 6. С. 8594360.
4. Панасенко Е.А., Тонков Е.Л. Инвариантные и устойчиво инвариантные множества дифференциальных включений // Труды Математического института им. В.А. Стеклова. 2008. Т. 262. С. 202-221.
5. Ушаков В.Н., Малев Я. А. К вопросу о дефекте стабильности множеств в игровой задаче о сближении // Труды Института математики и механики УрО РАН. 2010. Т. 16. № 1. С. 199-222.
6. Ушаков В.И., Зимовец A.A. Дефект инвариантности множеств относительно дифференциального включения // Вестник Удмуртского университета. Математика. Механика. Компьютерные науки. 2011. Вып. 2. С. 98-111.
7. Родина Л.И., Тонкое Е.Л. Статистические характеристики множества достижимости управляемой системы, неблуждаемость и минимальный центр притяжения // Нелинейная динамика. 2009. Т. 5. № 2. С. 265-288.
8. Родина Л.II. Оценка статистических характеристик множества достижимости управляемых систем // Известия вузов. Математика. 2013. № 11. С. 20-32.
9. Родина Л.И. Инвариантные и статистически слабо инвариантные множества управляемых систем // Известия Института математики и информатики УдГУ. 2012. Вып. 2 (40). С. 3-164.
10. Родина Л.И., Хаммади А.Х. Характеристики множества достижимости, связанные с инвариантностью управляемой системы на конечном промежутке времени / / Вестник Удмуртского университета. Математика. Механика. Компьютерные науки. 2013. Вып. 1. С. 35-48.
11. Родина Л.И., Хаммади А.Х. Статистические характеристики множества достижимости управляемых систем со случайными коэффициентами // Известия вузов. Математика. 2014. № 11. С. 50-63.
12. Колмогоров А.И., Фомин С.В. Элементы теории функций и функционального анализа. М.: Наука, 1976. 544 с.
13. Давыдов А.А., Пастрес Р., Петренко И.А. Оптимальное распределение выброса загрязнения в одномерный поток // Труды Института математики и механики УрО РАН. 2010. Т. 16. № 5. С. 30-35.
14. Кузенков О.А., Рябова Е.А. Математическое моделирование процессов отбора. Нижний Новгород: Издательство Нижегородского госуниверситета, 2007. 324 с.
15. Ризниченко Г.Ю. Лекции по математическим моделям в биологии. Часть 1. Ижевск: НИЦ «Регулярная и хаотическая динамика», 2002. 232 с.
16. Недорезов Л.В. Курс лекций по математической экологии. Новосибирск: Сибирский хронограф, 1997. 161 с.
Поступила в редакцию 31.03.2016
Родина Людмила Ивановна, д.ф.-м.н., зав. кафедрой математического анализа, Удмуртский государственный университет, 426034, Россия, г. Ижевск, ул. Университетская, 1. E-mail: LRodina67@mail.ru
Хаммади Алаа Хуссейн, аспирант, кафедра математического анализа, Удмуртский государственный университет, 426034, Россия, г. Ижевск, ул. Университетская, 1. E-mail: alaairaqmath@yahoo.com
L.I. Rodina, A. H. Hammady
Characteristics of invariancy for the attainability set of a control system
Keywords: control systems, differential inclusions, attainability set. MSC: 34H05, 34H99, 93C10
We study characteristics associated with invariancy or weak invariancy of a given set M = {(t,x) € [0, x Rn : x € M(t)} with respect to a control system X = f (t,x,u) on a finite time interval. One of such characteristics is relative frequency freq^ T+^](D,M) of containing the attainability set D(t,X) of this system in the set M on a segment [т, т + $]. This characteristic is equal to the quotient of the Lebegues measure of those t from [т, т + at which D(t, X) С M(t) to the length of the given segment. Other characteristic, freq^(D, M) = inf freq[T (D, M)
displays uniformity of containing the attainability set D(t, X) in the set M on a segment of the fixed length We prove theorems about estimation and calculation of these characteristics for various multivalued functions M(t) and D(t,X). In particular, we receive equalities for freqT(D,M) if the function M(t) is periodic with a period T and the function D(t, X) satisfies the inclusion D(t + T, X) С D(t, X) for all t ^ 0. We consider examples of calculation and estimations of these characteristics.
REFERENCES
1. Hartman Ph. Ordinary differential equations, New York-London-Sydney: John Wiley and Sons, 1964. Translated under the title Obyknovennye differentsial'nye uravneniya, Moscow: Mir, 1970, 720 p.
2. Aubin J.-P. Viability theory, Boston-Basel-Berlin: Birkhauser, 1991, 543 p.
3. Panasenko E.A., Tonkov E.L. Lyapunov functions and the positive invariant sets of differential inclusions, Differ. Uravn., 2007, vol. 43, no. 6, pp. 859-860 (in Russian).
4. Panasenko E.A., Tonkov E.L. Invariant and stably invariant sets for differential inclusions, Proceedings of the Steklov Institute of Mathematics, 2008, vol. 262, no. 1, pp. 194-212.
5. Ushakov V.N., Malev A.G. On the question of the stability defect of sets in an approach game problem, Proceedings of the Steklov Institute of Mathematics, 2011, vol. 272, no. 1, pp. 229-254.
6. Ushakov V.N., Zimovets A.A. Invariance defect of sets with respect to differential inclusion, Vestn. Udmurt. Univ. Mat. Mekh. Komp. Nauki, 2011, no. 2, pp. 98-111 (in Russian).
7. Rodina L.I., Tonkov E.L. Statistical characteristics of attainable set of controllable system, non-wandering, and minimal attraction center, Nelin. Dinam., 2009, vol. 5, no. 2, pp. 265-288 (in Russian).
8. Rodina L.I. Estimation of statistical characteristics of attainability sets of controllable systems, Russian Mathematics, 2013, vol. 57, issue 11, pp. 17-27.
9. Rodina L.I. Invariant and statistically weakly invariant sets of control systems, Izv. Inst. Mat. Inform. Udmurt. Gos. Univ., 2012, no. 2 (40), pp. 3-164 (in Russian).
10. Rodina L.I., Hammady A.H. The characteristics of attainability set connected with invariancy of control systems on the finite time interval, Vestn. Udmurt. Univ. Mat. Mekh. Komp'yut. Nauki, 2013, no. 1, pp. 35-48 (in Russian).
11. Rodina L.I., Hammady A.H. Statistical characteristics of attainability set of controllable systems with random coefficients, Russian Mathematics, 2014, vol. 58, no. 11, pp. 43-53.
12. Kolmogorov A.N., Fomin S.V. Elementy teorii funktsii i funktsional'nogo analiza (Elements of functions theory and functional analysis), Moscow: Nauka, 1976, 544 p.
13. Davydov A.A., Pastres R., Petrenko I.A. Optimization of the spatial distribution of pollution emission in 1D flow, Tr. Inst. Mat. Mekh. Ural. Otd. Ross. Akad. Nauk, 2010, vol. 16, no. 5, pp. 30-35 (in Russian).
14. Kuzenkov O.A., Ryabova E.A. Matematicheskoe modelirovanie protsessov otbora (Mathematical modeling of processes of selection), Nizhnii Novgorod: Nizhnii Novgorod State University, 2007, 324 p.
15. Riznichenko G.Yu. Lektsii po matematicheskim modelyam v biologii. Chast' 1 (Lectures on mathematical models in biology. Part 1), Izhevsk: Regular and Chaotic Dynamics, 2002, 232 p.
16. Nedorezov L.V. Kurs lektsii po matematicheskoi ekologii (Course of lectures on mathematical ecology), Novosibirsk: Sibirskii khronograf, 1997, 161 p.
Received 31.03.2016
Rodina Lyudmila Ivanovna, Doctor of Physics and Mathematics, Head of the Department of Mathematical
Analysis, Udmurt State University, ul. Universitetskaya, 1, Izhevsk, 426034, Russia.
E-mail: LRodina67@mail.ru
Hammady Alaa Hussein, Post-Graduate Student, Department of Mathematical Analysis, Udmurt State University, ul. Universitetskaya, 1, Izhevsk, 426034, Russia.
E-mail: alaairaqmath@yahoo.com