Введение 3
Постановка задачи 6
Глава 1. Линейные дифференциальные системы с постоянным запаздыванием 10
1.1. Фундаментальная матрица и её свойства 11
1.2. Функционалы полного типа 12
1.3. Матрицы Ляпунова 15
1.4. Алгоритм проверки положительной определенности квадратичного функционала 16
1.5. Построение матрицы Ляпунова и анализ положительной определенности функционала на практике 22
Глава 2. Линейные дифференциальные системы с линейно возрастающим запаздыванием 25
2.1. Анализ устойчивости 27
2.2. Фундаментальная матрица и ее свойства 28
2.3. Функционал с заданной производной 29
2.4. Матрица Ляпунова 32
Заключение 34
Список литературы 35
В настоящее время область приложений дифференциальных уравнений с отклоняющимся аргументом заметно расширяется. Системы таких уравнений позволяют моделировать динамические процессы в различных областях физики, техники, биологии, медицины, экономики, экологии.
Дифференциальные уравнения с отклоняющимся аргументом — дифференциальные уравнения, в которые неизвестная функция входит при различных аргументах. Такие уравнения можно разделить на три типа: уравнения с запаздывающим аргументом, уравнения с опережающим аргументом и уравнения нейтрального типа. С основами теории дифференциальных уравнений с отклоняющимся аргументом можно познакомиться в работах [1, 2].
В настоящей работе рассматриваются линейные дифференциальные уравнения с запаздывающим аргументом, описывающие процессы, в которых скорость изменения в любой момент времени зависит от предыдущих состояний, а не только от текущего состояния.
Для дифференциальных уравнений существует несколько методов анализа устойчивости. Два основных подхода были разработаны Ляпуновым А. М., они известны как первый и второй (прямой) методы Ляпунова.
Под вторым методом Ляпунова понимают совокупность приемов и средств исследования устойчивости решений систем дифференциальных уравнений при помощи специальных функций Ляпунова. Для систем вида x = Ax, где x 2 Rn, A — постоянная матрица, известен критерий: система экспоненциально устойчива тогда и только тогда, когда существует положительно-определенная квадратичная форма (функция Ляпунова), производная которой вдоль решений данной системы представляет собой отрицательно-определенную квадратичную форму.
Существуют обобщения второго метода Ляпунова на системы с запаздыванием. Данная работа будет основана на одном из таких методов — методе Красовского Н. Н.. В силу того, что состояние линейной системы с запаздыванием — это сегмент траектории, а не значение в текущий момент времени, то функции Ляпунова Красовский заменяет функционалами, зависящими от сегмента траектории. Именно эти функционалы, используемые для анализа экспоненциальной устойчивости, в дальнейшем получили название функционалов Ляпунова-Красовского.
Для систем с запаздыванием доказаны теоремы, аналогичные теоремам Ляпунова. Эти теоремы основаны на анализе расположения нулей характеристического квазиполинома на комплексной плоскости. Отрицательность действительных частей всех его нулей — необходимое и достаточное условие экспоненциальной устойчивости системы [1].
Существует аналог критерия Ляпунова для линейных стационарных систем с запаздыванием. Он известен как теорема Красовского: пусть существует положительно-определенный функционал, допускающий квадратичные оценки, производная которого вдоль решений системы представляет собой отрицательно-определенный функционал, тогда система является экспоненциально устойчивой.
Ставится вопрос о нахождении функционала, позволяющего исследовать экспоненциальную устойчивость системы с запаздыванием. Существует ли такой функционал? Если существует, как его построить? Является ли он единственным, и как определить его положительную определенность?
Известны способы построения положительно-определенных функционалов с заданной отрицательно-определенной производной, которые, согласно теореме Красовского, пригодны для анализа устойчивости [5]. Один из таких функционалов — функционал полного типа, который допускает квадратичную оценку снизу в случае экспоненциальной устойчивости системы.
Важным элементом, определяющим такие функционалы, является матрица Ляпунова. Для нее установлено необходимое и достаточное условие единственности [6]. В данной работе остановимся на вычислении матрицы Ляпунова для линейных дифференциальных уравнений с постоянным запаздыванием. По определению эта матрица является решением специальной системы уравнений, в которую входят некоторое условие симметрии, дифференциально-разностное уравнение и граничное условие. Причем матрица Ляпунова существует и единственна.
В данной работе рассмотрены методы проверки положительной определенности функционалов, которые пригодны для анализа экспоненциальной устойчивости для линейных дифференциальных систем с постоянным запаздыванием. Существует несколько таких методов, мы же сосредоточились на подходе, приведенном в [3].
Заметим, что не во всех случаях можно считать запаздывание постоянным и ограниченным, поэтому появляется необходимость в рассмотрении дифференциальных уравнений с переменным запаздыванием. На рисунке 1 представлены модели, которые описываются уравнениями с линейно возрастающим запаздыванием. Такие уравнения требуют особого внимания, так как в отличие от уравнений с ограниченным переменным запаздыванием, предыстория их решения, влияющая на последующую динамику, неограниченно возрастает при увеличении времени.
Рис. 1: Динамические модели с линейно возрастающим запаздыванием
Первая модель (слева) — модель смесительного бака, который наполняется разными жидкостями с помощью двух труб. Жидкости втекают с массовыми скоростями F1(t) и F2(t). В обоих потоках содержится некоторое растворимое вещество с постоянными концентрациями C-1 и C2. А из нижней части бака вытекает полученная смесь с массовой скоростью F(t) и концентрацией C(t).
Вторая модель (справа) — модель обменивающихся между собой информацией и удаляющихся друг от друга объектов. В данном случае присутствует транспортное запаздывание — время, необходимое для передачи информации между объектами, причем это время линейно возрастает при увеличении расстояния между объектами.
Существуют модели, при описании которых появляется не только линейно возрастающее запаздывание, но и постоянное. Примером такой модели является модель динамики транспортного потока на КАД («Ring road traffic»). Эта модель подробно рассмотрена в работе [4], но без учета постоянного запаздывания. В данной работе рассмотрена дифференциальная система, описывающая модель динамики транспортного потока на КАД с учетом постоянного запаздывания.
Вопрос распространения метода Ляпунова на линейные системы с постоянным запаздыванием рассмотрен в [5, 6], а распространение метода на системы с линейно возрастающим запаздыванием — в [9]. В настоящей работе представлено распространение метода функционалов на линейные системы с постоянным и линейно возрастающим запаздываниями.
Данная работа посвящена анализу положительной определенности линейных стационарных систем с постоянным запаздыванием. В ней использован подход, объединяющий метод Разумихина и метод функционалов Ляпунова-Красовского. Получены условия экспоненциальной устойчивости, основанные на оценке положительной определенности функционалов Ляпунова-Красовского на множестве функций, которые удовлетворяют аналогу условия Разумихина. Реализован алгоритм вычисления матрицы Ляпунова и проверки положительной определенности функционалов.
Также в данной работе рассмотрены линейные дифференциальные уравнения с постоянным и линейно возрастающим запаздываниями. Для них построен квадратичный функционал Ляпунова - Красовского, который может использоваться для анализа устойчивости и неустойчивости систем, если найдено решение системы уравнений, определяющих матрицу Ляпунова.
[1] Эльсгольц Л. Э., Норкин С. Б. Введение в теорию дифференциальных уравнений с отклоняющимся аргументом. М.: Наука, 1971. 296 с.
[2] Беллман Р., Кук К. Дифференциально-разностные уравнения / Пер. с англ. Под ред. Л. Э. Эльсгольца. М.: Мир, 1967. 548 с.
[3] Жабко А. П., Медведева И. В. Алгебраический подход к анализу устойчивости дифференциально-разностных систем // Вестник Санкт-Петербургского университета. Серия 10: Прикладная математика. Информатика. Процессы управления. 2011. № 1. С. 9-20.
[4] Zhabko, A., Chizhova, O., Zaranik, U. Stability Analysis of the Linear Time Delay Systems with Linearly Increasing Delay // Cybernetics and Physics. VOL. 5, NO. 2, 2016 , 67-72.
[5] Харитонов В. Л. Функционалы Ляпунова с заданной производной. I. Функционалы полного типа // Вестник Санкт-Петербургского университета. Серия 10: Прикладная математика. Информатика. Процессы управления. 2005. № 1. С. 110-117.
[6] Харитонов В. Л. Функционалы Ляпунова с заданной производной. II. Матрицы Ляпунова // Вестник Санкт-Петербургского университета. Серия 10: Прикладная математика. Информатика. Процессы управления. 2005. № 2. С. 200-209.
[7] Зубов В. И. К теории линейных стационарных систем с запаздывающим аргументом // Известия ВУЗов. Математика. 1958. № 6. С. 86-95.
[8] Зубов В. И. Устойчивость движения. М.: Высшая школа, 1973. 272 с.
[9] Меденников И. П. Прямой метод анализа устойчивости систем с линейно возрастающим запаздыванием // Вестник Санкт-Петербургского университета. Серия 10: Прикладная математика. Информатика. Процессы управления. 2014. № 3. С. 125-140.
[10] Жабко А. П. ,Чижова О.Н. Гибридный метод анализа устойчивости линейных дифференциально-разностных систем с линейно возрастающим запаздыванием // Вестник ТГУ, 2015. т.20. вып.4. С. 843-850.
[11] Разумихин Б. С. Об устойчивости систем с запаздыванием // Прикладная математика и механика. 1956. Т. 20. Вып. 4. С. 500-512.
[12] Красовский Н. Н. О применении второго метода А.М.Ляпунова для уравнений с запаздыванием времени. // Прикладная математика и механика. 1956. Т. 20. С. 315-327.