Введение 3
Постановка задачи 5
Обзор литературвх 5
Глава 1. Матрица Ляпунова. Приближённое ввхчисление и оценка
точности приближения 7
1.1. Система дифференциалвно-разностнвхх уравнений ... 7
1.2. Матрица Ляпунова как решение интегралвного уравнения 8
1.3. Устранение диагоналвнвхх разрвхвов ядра уравнения . . 12
1.4. Метод квадратур для построения приближённого реше¬
ния интегралвного уравнения 13
1.5. Оценка точности приближения методом Анселоне. ... 14
1.6. Реализация в системе Matlab 19
Глава 2. Применение матриц Ляпунова для исследования устой¬
чивости систем с запаздвхванием 21
2.1. Критерий устойчивости 21
2.2. Реализация в системе Matlab 23
Выводы 26
Заключение 27
Список литературы 28
В данной работе предложен метод приближённого вычисления матрицы Ляпунова для дифференциально-разностных систем с запаздывающим аргументом. С помощью систем уравнений данного вида описываются многие химические, физические и биологические процессы, происходящие в мире. Благодаря исследованию уравнений с запаздыванием, мы имеем возможности получать более точную информацию о процессах и предсказывать поведение объектов и систем в будущем. Необходимость использования запаздываний продиктована тем, что скорость изменения в системах зависит не только от состояния в настоящий момент времени, но и от происходившего в уже прошедший промежуток времени. При решении некоторых задач запаздыванием пренебрегают, и в ход идёт теория обыкновенных дифференциальных уравнений. Однако нужно понимать, при неучитывании запаздываний можно получить результат, никоим образом не соответствующий реальности. Поэтому, при детальном изучении процессов, целесообразно использовать функции с запаздывающим аргументом.
Понятие матрицы Ляпунова играет важную роль в теории дифференциально-разностных уравнений, так как поиск решений уравнений данного типа - трудоёмкая задача, хотя бы потому, что собственных чисел в системе бесконечно много. Все их не найти, а значит, проверить, что они лежат в левой комплексной полуплоскости непросто. Поэтому для изучения свойств систем используется метод функционалов Ляпунова-Красовского полного типа и матрица Ляпунова, являющаяся их основным элементом. В настоящее время матрица Ляпунова применима к таким задачам, как исследование систем на устойчивость, вычисление H2нормы придаточной матрицы. А сами функционалы полного типа Ляпунова-Красовского используются для исследования робастной устойчивости, построения экспоненциальных оценок и стабилизирующих управлений, вычисления значений функционалов качества в задачах оптимального управления и для других задач.
Тем не менее, в настоящее время не существует общего метода построения матриц Ляпунова. В случае одного запаздывания и нескольких кратных запаздываний задача построения сводится к решению линейных систем обыкновенных дифференциальных уравнений с линейными граничными условиями. Для некратных запаздываний применимы эвристические численные методы без оценки точности приближения. Также существуют методы с оценкой точности приближения для экспоненциально устойчивых систем.
Постановка задачи
Основная цель данной работы состоит в построении метода приближённого вычисления матрицы Ляпунова для дифференциально-разностных уравнений с некратными запаздываниями, а также в получении оценки точности приближения. Для этого задача построения матрицы Ляпунова сводится к решению уравнения Фредгольма II рода. Для решения поставленной задачи построения оценки приближения квадратурных формул применяется метод Анселоне.
Вторая задача касается применения матриц Ляпунова для исследования устойчивости линейных систем с запаздыванием. Для её решения применяется недавно полученная теорема, дающая необходимое и достаточное условие устойчивости систем с запаздыванием, выраженное в терминах матрицы Ляпунова. Основная задача состоит в написании программы с акцентом на уменьшении времени вычисления.
Обзор литературы
В последнее время появляется всё больше работ, посвящённых теории и приложениям дифференциальных уравнений с отклоняющимся аргументом, среди них работа [1] занимает особое место. В этой книге можно найти хороший обзор по теории систем с запаздываниями. Более новые результаты, касающиеся исследования устойчивости и связанных с нею характеристик для систем с запаздыванием можно найти в [10],[2],[11]. В книге Харитонова В.Л. [10] вводится понятие матрицы Ляпунова для систем с запаздываниями и описываются различные задачи, решаемые с помощью этой матрицы. В частности, матрица Ляпунова позволяет построить функционал Ляпунова-Красовского полного типа, а с помощью функционала исследовать робастную устойчивость, строить экспоненциальные оценки решений, вычислять значения интегральных критериев качества. В статье [2] приводится необходимое и достаточное условие устойчивости линейных систем с запаздываниями, выраженное в терминах матрицы Ляпунова.
На сегодняшний день существует несколько методов построения матрицы Ляпунова. В случае кратных запаздываний задача построения матрицы Ляпунова может быть сведена к решению линейной системы обыкновенных дифференциальных уравнений с линейными граничными условиями [5], [10]. Этот метод называется полуаналитическим.
Для общего случая существует ряд эвристических методов [6], [8], [9], [4] то есть численных методов без оценки точности приближения. Такие методы основаны на поиске решения в виде полиномиальной или кусочно-полиномиальной функции, коэффициенты которой находятся из условия минимизации невязки, которая характеризует точность выполнения свойств из определения матрицы Ляпунова.
Кроме того, недавно был представлен метод [3], позволяющий построить приближение матрицы Ляпунова с оценкой точности приближения. Но стоит отметить, что этот метод работает только при условии, что система экспоненциально устойчива. Следовательно, он не применим в случае, когда с помощью матрицы Ляпунова мы пытаемся исследовать устойчивость.
В выпускной квалификационной работе [16] было показано как задача построения матрицы Ляпунова может быть сведена к интегральному уравнению Фредгольма II рода. Полученное уравнение решалось численно методом квадратур. Стоит отметить, что попыток оценить точность приближения не предпринималось.
В настоящей работе построен метод приближённого вычисления матрицы Ляпунова для дифференциально-разностных систем с запаздывающим аргументом, что позволяет исследовать системах дифференциальных уравнений с запаздываниями, и была построена оценка точности приближения данного эвристического метода. При этом, не накладывалось условие экспоненциальной устойчивости линейных дифференциально-разностных систем с запаздывающим аргументом.
Исследована недавно полученная теорема о построении области устойчивости линейных систем с запаздывающим аргументом. На практике был реализован алгоритм, основанный на данной теореме, и получены конкретные значения. В результате чего, можно понять, что можно улучшить в дальнейшем в данной задаче. Например, уменьшить значения размерности матрицы, построения которой требует данный алгоритм проверки на устойчивость системы с запаздываниями. Что позволит затрачивать значительно меньше времени на выполнение программы.
1. Беллман Р., Кук К. Дифференциально-разностные уравнения. М.: Мир, 1967.
2. Egorov A.V., A Finite Necessary and Sufficient Stability Condition for Linear Retarded Type Systems // Труды конференции "IEEE 55th Conference on Decision and Control (CDC)”, Лас Вегас, США. 2016. с. 3155-3160.
3. Egorov A.V., Kharitonov V.L. Approximation of delay Lyapunov matrices // International Journal of Control, 2018. Vol. c. 2588-2596.
4. Егоров А.В., Вычисление матриц Ляпунова для систем с запаздыванием // XII Всероссийское совещание по проблемам управления, Москва. 16-19 июня 2014. с. 1292-1303.
5. Garcia-Lozano EL, Kharitonov V.L. Lyapunov matrices for time delay system with commensurate delays // 2nd Symposium on System, Structure and Control, Oaxaca, 2004. c. 102-106.
6. Garcia-Lozano EL, Kharitonov V.L. Numerical computation of time delay Lyapunov matrices // Department of Automatic Control, Mexico. 2006.
7. Graham A. Kronecker Products and Matrix Calculus with Applications. Ellis Horwood, Chichester, UK, 1981. 21-30 p.
8. Huesca E. , Mondie S. , Santos O. (2009). Polynomial approximations of the Lyapunov matrix of a class of time delay systems // In E. Petre (Ed.), Proceedings of the 8th IFAC workshop on time delay systems, Sinaia, Romania, 1-3 September 2009. c. 261-266. IFAC.
9. Jarlebring E., Vanbiervliet J., Michiels W. (2011). Characterizing and computing the H2 norm of time-delay systems by solving the delay Lyapunov equation // IEEE Transactions on Automatic Control, 56 (4).
10. Kharitonov V.L. Time-delay systems: Lyapunov functionals and matrices. Birkhauser, Basel, 2013.
11. Kharitonov V.L., Zhabko A.P. Lyapunov-Krasovskii approach to the robust stability analysis of time-delay systems // Automatica 39, 2003. c. 15-20.
12. Mondie S., Egorov A.V. Some necessary conditions for the exponential stability of one delay systems// Труды конференции "8th International Conference on Electrical Engineering", 2011. c. 106.
13. Манжиров А.В., Полянин А.Д. Справочник по интегральным уравнениям. М.: Факториал Пресс, 2000.
14. Олемской И.В., Численные методы. Паств II. Учебное пособие. СПб, 2012.
15. Хатсон В., Пим Дж. Приложения функционального анализа и теории операторов. М.: Мир, 1983. 215-224 с.
16. Чернышева Л. Выпускная квалификационная работа. Матрица Ляпунова как решение уравнения Фредголвма. СПб, 2016. 10-16 с.