CC BY
39
УДК 517.929
О РОБАСТНОЙ УСТОЙЧИВОСТИ НЕАВТОНОМНЫХ ФУНКЦИОНАЛЬНО-ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
© В.В. Малыгина
Ключевые слова: функционально-дифференциальное уравнение; робастная устойчивость; эффективные признаки.
Для линейного неавтономного дифференциального уравнения с несколькими запаздываниями получен эффективный признак робастной устойчивости. В качестве примера рассмотрено уравнение с двумя запаздываниями, для которого построена область ро-бастной устойчивости в пространстве параметров исходной задачи.
1. Робастная устойчивость для функционально-дифференциальных уравнений. Одно из противоречий, которое приходится преодолевать при адаптации математической модели к приложениям — это невозможность точного задания параметров модели. Как правило, параметры задаются «интервальными оценками»: известно, что параметр принадлежит заданному интервалу, и ответ следует дать в терминах границ этих интервалов. При исследовании устойчивости дифференциальных (в частности, функционально-дифференциальных) уравнений вопросы такого рода характеризуются специальным термином «робастная устойчивость» (от англ. robust — крепкий, грубый).
Отметим, что задача робастной устойчивости содержательна уже для автономных функционально-дифференциальных уравнений [1, 2], т. к. постоянство коэффициентов и запаздываний еще не означает, что они точно заданы. В ситуации, когда рассматриваемый объект — неавтономный, задача робастной устойчивости, сохраняя свою актуальность, становится значительно сложнее.
Первый результат, который можно классифицировать как признак робастной устойчивости для неавтономных функционально-дифференциальных уравнений — знаменитая «3/2-теорема» А.Д. Мышкиса.
Рассмотрим дифференциальное уравнение с одним переменным запаздыванием
X(t) + a(t)x(t - r(t)) = 0, t € R+. (1)
Т е о р е м а 1 (А.Д. Мышкис, 1951, [3]). Пусть 0 < r(t) < r, m < a(t) < M .
• Если m> 0, Mr < 3/2, то уравнение (1) асимптотически устойчиво.
• Если m ^ 0, Mr ^ 3/2, то уравнение (1) устойчиво по Ляпунову.
Постоянная 3/2 в теореме Мышкиса является точной: построенные в той же работе [3] примеры показывают, что ее нельзя увеличить ни на какую, сколь угодно малую величину.
Теорема 1 в дальнейшем многократно переоткрывалась и уточнялась (см.[4-6]); в частности, было замечено, что в условиях первой части теоремы уравнение (1) является экспоненциально устойчивым, а в условиях второй части — равномерно устойчивым.
Обращает на себя внимание то, что в теореме 1 роль границ интервалов, задающих параметры уравнения (1), существенно различна. Значения постоянных M и r определяют величину области устойчивости, и чем точнее они заданы, тем больше область в пространстве параметров, попадание в которую гарантирует устойчивость уравнения (1). Напротив, величина границы m не существенна: в первой части теоремы важна отделенность коэффициента a от нуля, во второй части требуется лишь его неотрицательность. Более того,
если рассмотреть частный случай уравнения (1), положив т = М, то постоянная 3/2 по-прежнему останется точной границей области устойчивости. Повысить ее удается, только если сделать постоянным также и запаздывание: при г(*) = г > 0 область экспоненциальной устойчивости увеличится до интервала вида 0 < Мг < п/2 , а область равномерной устойчивости — до отрезка 0 ^ Мг ^ п/2 .
Причины такого явления, по-видимому, объясняются спецификой уравнений с последействием: наличие запаздывающего аргумента изначально делает ситуацию «несимметричной».
2. Обобщение на случай нескольких запаздываний. Рассмотрим скалярное функционально-дифференциальное уравнение вида
N
Х(*) + ^ ак(*)х(* - гк(*))=0, * € М+, (2)
к=0
где функции ак локально суммируемы, функции Гк измеримы. Решение уравнения (1) принято считать [7, с. 9] принадлежащим классу абсолютно непрерывных на каждом конечном отрезке функций. При отрицательных значениях аргумента полагаем решение доопределенным некоторой начальной функцией.
Пусть параметры уравнения (2) при всех * € М+ подчинены неравенствам: ак ^ ак(*) ^ ^ Ак , 0 ^ Гк (*) ^ Гк , где ак ,Ак,Гк —заданные положительные постоянные. Поставим задачу найти условия устойчивости уравнения (2) в терминах этих констант.
Поставим в соответствие уравнению (2) так называемое ЬввЬ-уравненив:
N
Рк(*)у(* - Гк)=0, * € М+. (3)
к=0
При любых £ ^ 0 полагаем у(£) = 1, функции Рк определяем по следующему правилу: 'Рк (*) = Ак , если * ^ г к + *0 и рк (*) = ак , если * > Гк + *0 , где *0 - первый нуль решения test-уравнения. Случай ¿0 = ^ не исключается.
Теорема2. Пусть I — первый минимум решения ЬввЬ-уравнения.
• Если у(1) > —1, то уравнение (2) экспоненциально устойчиво.
• Если у(1) ^ —1, то уравнение (2) 'равномерно устойчиво.
3. Пример построения области робастной устойчивости. Покажем, как работает теорема 2 применительно к конкретным объектам. Рассмотрим уравнение вида
х(*) + а(*)х(*) + ь(*)х(* — г(*)) = 0, * € м+. (4)
Не нарушая общности можно считать, что 0 ^ г(*) ^ 1. Коэффициенты уравнения (4) подчиним неравенствам 0 < а ^ а(*) ^ А и 0 < Ь ^ Ь(*) ^ В .
Test-урaвнение для (4) легко строится «по шагам».
Пусть * € [0,1] .В малой окрестности нуля решение уравнения (3) имеет вид у(*) = = (1 + В/А)е-М — В/А, откуда легко определяется его первый нуль: *0 = 1/А 1п(1 + + А/В) . Следовательно, на отрезке [0,^] вид функции у не меняется. Далее, выстраивая функции '1 и '2 в соответствии с определением, находим решение на отрезке [*0, 1] , где оно становится отрицательным, но продолжает монотонно убывать, и, наконец, на отрезке [1,1 + , где решение достигает первого минимума * = I. Первый минимум определяется из уравнения
=(1-0 =
1- e-
t A - a A + B
A\ a/A 1 + B
1/A-a
В соответствии с теоремой 2, кривая у(1) = —1 есть граница области устойчивости, которая для уравнения (4) имеет вид:
Bi Л + A
As V + B
e-As —
Bq
As(1 - s)
i+A
+
Bq1/s
A(I^)
i+A
+
A2 s
= 1,
(5)
где s = a/A ,а q = (l - e-AsA+# (1 + A)S)
s/1-s
На рисунке 1 приведены графики кривых (5) при в , меняющемся от 0 до 1 с шагом 1/10. При в = 0 и в = 1 кривые (5) формально не определены, но, как несложно убедиться, имеют пределы при в ^ 0 и в ^ 1 и, следовательно, могут быть в этих точках доопределены по непрерывности.
Области устойчивости уравнения (4) — множество точек плоскости, лежащих между осью ОА и кривыми (5).
Интересно отметить, что нижняя граница коэффициента Ь, как и в теореме 1, не входит в формулу для границ области устойчивости, тогда как обе границы коэффициента а оказываются востребованными. Граница, которая получается в результате предельного перехода в ^ 1, совпадает, как и следовало ожидать, с точной границей устойчивости уравнения (4) в ситуации а = а(£) = А , Ь = Ь(£) = В (см. [8]).
Рис. 1
s
ЛИТЕРАТУРА
1. Харитонов В.Л. Критерий устойчивости одного семейства квазиполиномов запаздывающего типа // Автоматика и телемеханика. 1991. Вып. 2. С. 73-82.
2. Харитонов В.Л. Семейство устойчивых квазиполиномов // Автоматика и телемеханика. 1991. Вып. 7. С. 75-88.
3. Мышкис А.Д. Линейные дифференциальные уравнения с запаздывающим аргументом. Москва; Ленинград: Гостехиздат, 1951.
4. Yorke J.A. Asymptotic stability for one-dimensional delay equations //J. Different. Equat. 1970. V. 7. № 1. P. 189-202.
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: mavera@list.ru
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: mavera@list.ru
УДК 517.977
ЭЛЛИПСОИДАЛЬНОЕ ОЦЕНИВАНИЕ ТРУБОК ТРАЕКТОРИЙ УПРАВЛЯЕМЫХ ИМПУЛЬСНЫХ СИСТЕМ С КОНИЧЕСКИМ ОГРАНИЧЕНИЕМ НА УПРАВЛЕНИЕ
© О.Г. Матвийчук
Ключевые слова: множества достижимости; импульсное управление; эллипсоидальное исчисление.
В работе рассматриваются задача оценивания множеств достижимости управляемых динамических систем, описываемых дифференциальными уравнениями с импульсным управлением и неполной информацией о начальных данных. Предполагается, что импульсные управления в динамической системе находятся в пересечении специального конуса и обобщенного эллипсоида в пространстве функций ограниченной вариации. Целью данной работы является исследование свойств множеств достижимости управляемой системы заданного вида и построение алгоритма внутреннего эллипсоидального оценивания траекторий динамических систем импульсного типа.
