CC BY
0
Math-Net.Ru
А. А. Рогалев, Оценивание множеств решений линейных систем обыкновенных дифференциальных уравнений с возмущениями на основе оператора Коши, Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2023, номер 2, 357-374
001: 10.14498/уБ^и1978
Использование Общероссийского математического портала Math-Net.Ru подразумевает, что вы прочитали и согласны с пользовательским соглашением
http://www.mathnet.ru/rus/agreement
Параметры загрузки:
IP: 109.252.33.182
29 сентября 2024 г., 12:10:17
Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2023. Т. 27, № 2. С. 357-374_
ISSN: 2310-7081 (online), 1991-8615 (print) d https://doi.org/10.14498/vsgtu1978
EDN: MXEZAB
УДК 517.925
Оценивание множеств решений линейных систем обыкновенных дифференциальных уравнений с возмущениями на основе оператора Коши
А. А. Рогалев
Институт космических и информационных технологий
Сибирского федерального университета
Россия, 660074, Красноярск, ул. Академика Киренского, 26 к/1.
Аннотация
Излагается метод численного анализа множеств решений линейных систем обыкновенных дифференциальных уравнений, содержащих возмущения в правой части. Метод определяет экстремальные значения решений, которые составляют множества решений по осям координат или в заданном направлении. Получены оценки на основе использования оператора Коши, записанного символьными формулами вариации произвольных постоянных. Дополнительно реализован контроль отклонения решений при расчете пучка траекторий. Приведены примеры оценивания множеств достижимости систем при воздействии управляющих воздействий и возмущений.
Ключевые слова: возмущающие воздействия; пучок траекторий; символьные алгоритмы.
Получение: 30 декабря 2022 г. / Исправление: 14 апреля 2023 г. / Принятие: 25 мая 2023 г. / Публикация онлайн: 23 июня 2023 г.
Математическое моделирование, численные методы и комплексы программ Научная статья
© Коллектив авторов, 2023 © СамГТУ, 2023 (составление, дизайн, макет)
3 ©® Контент публикуется на условиях лицензии Creative Commons Attribution 4.0 International (https://creativecommons.org/licenses/by/4.0/deed.ru) Образец для цитирования
Рогалев А. А. Оценивание множеств решений линейных систем обыкновенных дифференциальных уравнений с возмущениями на основе оператора Коши // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2023. Т. 27, № 2. С. 357-374. EDN: MXEZAB. DOI: 10.14498/vsgtu1978. Сведения об авторе
Александр Алексеевич Рогалев А https://orcid.org/0000-0003-2176-9639 старший преподаватель; каф. информационных систем; e-mail: gogoba88@mail.ru
Введение. Рассматривается задача оценки множества решений систем обыкновенных дифференциальных уравнений (ОДУ) с возмущениями. Такие задачи возникают при исследовании устойчивости с постоянно действующими возмущениями на конечном интервале времени [1—3], множеств достижимости (МД) управляемых систем [4—11], оценке выживаемости систем [12]. Оценка множеств решений для линейных систем ОДУ с возмущениями актуальна в наше время во многих современных практических сферах, к которым относятся различные системы управления реальными объектами, оценка характеристик объектов в баллистике при массовых вычислениях, управление беспилотными летательными аппаратами и автомобилями [13-14]. Постановка задачи оценки множества решений системы ОДУ с возмущающими воздействиями обусловлена тем, что для рассматриваемой модели отсутствует достаточное описание всех параметров, включающих действующие возмущения (управления) [1, 4, 9-10], а известны только границы изменения возмущений (управлений). В подобных случаях необходимо оценивать все множество решений, соответствующих возмущениям.
Универсального метода получения оценок не существует ввиду многообразия и сложности задач [9]. Существующие методы основаны на различных подходах (например, метод оптимальных двусторонних оценок МД на основе операций над эллипсоидами [4-9], подход работы [10]), в основе которых лежит многократное решение систем ОДУ с измененными значениями параметров, что позволяет учитывать все возможные значения коэффициентов этих систем под воздействием возмущений. В последние годы развиваются исследования в области разработки и применения методов адаптивной интерполяции в задачах с интервальными неопределенностями [15-17], в которых для каждого момента времени строятся кусочно-полиномиальные функции, интерполирующие зависимость решения задачи от точечных значений интервальных параметров с заданной точностью. Авторы метода адаптивной интерполяции применяют й-ё-деревья для хранения информации о найденных значениях численного решения. С помощью этого метода решены несколько практически важных задач [15-17], в которых требуется определить верхние и нижние границы компонент множества решений.
Следует выделить также работы, описывающие подход, основанный на символьной формуле оператора сдвига вдоль траектории и вычислении множеств значений этой формулы по всем множествам управления и управляющим воздействиям [18] и его применение в ряде задач [19-23].
Иногда возникают трудности, связанные с применением многих методов оценки множеств решений с возмущающими и управляющими параметрами [13-14], которые заключаются в расширении этих оценок независимо от поведения самих решений, а также в возможном увеличении ошибок численных расчетов. Вычисление границы МД может иметь большую вычислительную сложность для систем высокой размерности как в линейном, так и нелинейном случаях. Некоторые методы не позволяют установить параметры и управление, ведущие на границу [13]. В настоящее время разработка надежных вычислительных алгоритмов для оценки МД остается актуальной [14].
В данной статье развиваются методы решения задачи оценки множеств решений системы ОДУ с возмущающими воздействиями, основанные на сим-
вольных вычислениях [18]. Предложен и описан новый подход, основанный на символьной формуле Коши и учитывающий особенности интегральной кривой поставленной задачи, то есть свойства решений. В построенных на основе данного подхода численно-символьных алгоритмах вычисление максимальных и минимальных значений границ множеств решений выполняется независимо от предыдущих шагов на каждом шаге по времени. Ошибка вычислений границ множеств решений от шага к шагу не накапливается. Это является преимуществом рассматриваемого метода.
К достоинствам метода относится также использование для оценки множеств достижимости на входе возмущающих или управляющих воздействий как произвольных функций широкого класса.
1. Постановка задачи. Рассматривается следующая начальная задача для системы ОДУ:
^ = А(1)у(1) + и(1), г е кп (1)
у(1о)= уо е Уо, (2)
где уо е Го С Мга, и : [¿о, Т] ^ [-1,1], у : [¿о, Т] ^ Мга, А(1) : [¿о, Т] ^ х Мга, Т > ¿о. Пусть для любого вектора начальных данных существует единственное решение этой задачи у(Ъ) на интервале [Ьо,Т], у(Ъ) принадлежит соответствующему функциональному пространству (является непрерывной, непрерывно-дифференцируемой или абсолютно непрерывной функцией, в общем случае можно сказать — суммируемой функцией). Требуется найти включение (построить границы) множества решений У(Ъ) С Мга такое, что
у(1) е У(I) (3)
на интервале [Ьо,Т] для всех решений у(Ъ) = у(Ь,Ьо,уо) задачи (1) с начальными данными (2) и возмущающим воздействием и(Ь).
Для построения включения (3) множества решений системы ОДУ с начальными данными воспользуемся оператором Коши У(Ьо,Ь) : Мга ^ Мга, сопоставляющим со значением каждого решения у(Ь) системы ОДУ в точке Ь = Ьо значение этого же решения в точке ¿:
Г (1о,1)у(1о)= У(Ь). (4)
Выполняются следующие свойства оператора Коши:
у(ь,ъ) = о, У= У-1(г,ь), у(ъ,т)У(т,г) = У(ь,г).
Если система ОДУ (1) является линейной системой с суммируемой правой частью, то формула (4) записывается в виде следующих форм:
К
ч-Ь
у(Ъ) = у(1, Ь,Уо) = ф(г)(ф-1(ь )уо + £ Ф-1(т )и(т , (5)
у(г) = у(г, ь,уо) = ехр(1А)(вхр-1(1оА)уо + £ вхр-1(тА)и(т. (6)
Р о г а л е в А. А.
Формула (5) справедлива для систем линейных ОДУ с переменными коэффициентами, формула (6) выполняется для систем с постоянными матрицами коэффициентов. Здесь и(т) = (и1(т),..., ип(т)) — вектор внешних воздействий (возмущений); у0 = (у1,0,... ,уп,о) — вектор начальных данных; Ф(£) — матрица фундаментальных решений линейной однородной системы, полученной из системы (1); Ф-1(£) —обратная матрица к матрице фундаментальных решений; ехр(£А) — матрица фундаментальных решений линейной однородной системы с постоянными коэффициентами; ехр-1(£А) —обратная матрица к матрице фундаментальных решений линейной однородной системы с постоянными коэффициентами.
На основе формул (5) или (6) нельзя сразу же получить числовые значения. Необходимо выбрать алгоритм и организовать вычисления; величины
ехр(М), ехр-1 (ЬА), Ф(£), у Ф(£)Ф-1 (¿0) — символьные величины, для определения значений которых нужно реализовать соответствующие алгоритмы.
Представим формулы (5) и (6) как соотношения, связывающие решение и неопределенные воздействия. На их основе необходимо выбрать алгоритм и организовать вычисления, позволяющие получить верхние и нижние оценки множества всех решений. При этом величины, подобные ехр(£А), ехр-1(£А),
Ф(^), У Ф(*)Ф-1(*о), воспринимаются как символьная формула — выражение.
Для представления первой части формул (5), (6) в виде структуры, пригодной для организации вычислений на ЭВМ, необходимо записать матрицу фундаментальных решений как символьную формулу, включающую алгебраические величины (переменные и числа), соединенные знаками арифметических операций, знаками извлечения корня и возведения в степень и скобками, указывающими на последовательность применения арифметических операций. Далее необходимо выполнить такую же запись символьной формулы обратной матрицы фундаментальных решений, умножить символьные формулы матриц между собой и умножить произведение на символьный вектор начальных данных. Матрица фундаментальных решений является матричной функцией, в этом отличие от решения обычной системы ОДУ, решение которого есть векторная функция. Для реализации указанного процесса необходимо разработать алгоритм, не приводящий к значительному росту вычислительных затрат, пригодный к использованию для широкого класса систем вида (1), и реализовать его в виде программного комплекса.
2. Описание численно-символьного алгоритма. Первая часть алгоритма описывает построение символьных формул и организацию вычислений по этим формулам включенных в процесс вычисления оценок областей значений. Следует отметить, что оценивать области значений в этой части формулы Коши не требуется, так как полагаем, что начальные данные и все коэффициенты заданы точно. Далее в описании алгоритма мы не будем ссылаться на номера формул (5) или (6), так как эти действия на всех шагах выполняются как для формул (5), так и для формул (6).
Ш!аг 1. Определяется символьная формула решений фундаментальной матрицы решений матричной системы ОДУ
^Ф Ч / N ж
-Ж =
Для нахождения матричного решения системы (7) используется метод нахождения полиномиального решения линейной функциональной системы либо метод нахождения рационального решения линейной функциональной системы [25], либо метод построения матричной экспоненты, если А — постоянная матрица. В итоге получаем символьные формулы решений. Эти формулы входят как компоненты в первый член формулы Коши Ф(£)Ф-1(£о)^о.
Ш!аг 2. На этом шаге символьная формула Ф(£) преобразуется к символьным формулам Ф(*о), Ф-1(^о)уо, выполняется подстановка символьных переменных и их преобразование. Первый член символьной формулы Коши Ф(£)Ф-1(£о)уо конструируется на основе выражений, полученных после этой подстановки. Размерность произведения матричных и векторных функций равна размерности решения. Под интеграл во втором члене суммы входит также вектор возмущающих воздействий.
Вторая часть шага 2 алгоритма описывает преобразование подынтегрального выражения Ф-1(т)и(т), нахождение первообразной преобразованного выражения, затем подстановку первообразной в символьную формулу Коши и оценивание границ области значений этой формулы.
Это достаточно сложная для реализации часть всего алгоритма. Для внешнего воздействия u(t) известны только границы его значений. Поэтому интегрирование неизвестной функции, зависящей от времени, невыполнимо. Можно преобразовать подынтегральное символьное выражение к виду, для которого определяется первообразная. Затем с использованием этой первообразной, оценивается область значений второй части формулы (1), а значит оценивается область значений формулы (1).
ШШаг 3. На этом шаге оценивается область значений формулы
Ф(*) / Ф-1(т)и(т)dr.
JtQ
В ней изменяется один из сомножителей подынтегрального выражения и(т), \и(т)| ^ 1 для любого т, формулы фундаментальной матрицы решений и обратной матрицы не изменяются.
Чтобы вычислить границы множества значений
N / Ф-1(тМт)^т о
и(т
при и(т) € ["Umin,-umax], в алгоритме определяются экстремальные значения
подынтегрального выражения Ф(г) Ф-1(т)и(т)йт в фиксированных точ-
о
ках времени tk, к = 1,... ,п. Выбор значений возмущающих (управляющих) воздействий как решений задач
max Ф-1(т )и(т) и min Ф-1(т )и(т)
в фиксированные моменты времени эквивалентен тому, что возмущающее воздействие выбирается из условия приближения подынтегрального выражения сверху, затем — из условия приближения снизу. При этом используются методы, в которых применяются только значения целевой функции и не вычисляются производные [26].
Для графического представления границы множества значений и ее оценки строится интерполяционный полином по точкам максимума подынтегрального выражения
Ф(ф-1(т)и(т) = <£(tConst)$>-1(r)u(T),
вычисленным в точках tk; значение t = tConst фиксировано, значения переменной т изменяются в точках tk:
max $-1(ii)u(ii), ..., max Ф-1 (tn)u(tn).
^min 1) ^^max ^min ^n ) ^^max
Для интерполирования максимальных значений применяется кубический интерполяционный сплайн [27], обозначенный здесь как MaxICS(i), to ^ t ^ tn. Кубический интерполяционный сплайн MinICS(i), to ^ t ^ tn, строится аналогично по точкам
min Ф-1^1)и^1), ..., min Ф-1 (tn)u(tn).
^min 11) ^^max ^min ^n ) ^^max
Так как все входящие в подынтегральное выражение функции векторные, данная процедура повторяется для каждой компоненты построенных векторных функций.
Для многих задач возмущающее воздействие [5] может не быть непрерывной, а также дифференцируемой функцией. Возможные варианты выбора возмущающих (управляющих) воздействий: кусочно-гладкая функция, производная которой интегрируема по Лебегу, либо абсолютно непрерывная функция. Поэтому необходимо оценить точность построения интерполяционной или сглаживающей сплайн-функции MaxICS(i), to ^ t ^ tn, вычисленной по точкам максимальных значений подынтегральных выражений, а также сплайн-функции MinICS(i), to ^ t ^ tn, вычисленной по точкам минимальных значений подынтегральных выражений. Компоненты подынтегрального выражения — непрерывные функции, зависящие от t.
Определение. Для любой функции f (t), заданной на отрезке Т, верхняя грань колебаний функции [27] по всем подотрезкам из Т длиной меньше ö является модулем непрерывности этой функции:
u(f,S)=sup{\f (h) - f (t2)\, t1,h e T, \t1 -12\ <ö}.
Сформулируем утверждение, описывающее сходимость интерполяционных полиномов для непрерывной функции на отрезке.
Утверждение. Пусть на отрезке [¿о, Т] заданы непрерывная функция у^) и такая последовательность сеток Ак, что ||Дк|| ^ 0 при к ^ ж. Если кубический сплайн дефекта 2 интерполирует функцию у^) в узлах сетки Дк и его производная в этих узлах равна нулю, то последовательность ф}
сходится к у(Ъ) равномерно на [¿0,Т]. Выполняется неравенство
\т - зДк(4)Ц < 2Ш[ ,у). (8
Доказательство утверждения и оценки (8) проводится методами из [27]. Дополнительным ограничением интерполяционного сплайна является требование равенства производной сплайна нулю в узлах сетки. Поскольку по сути задачи нас более всего интересуют оценки значений подынтегрального выражения в узлах сетки, равенство производной интерполяционного сплайна нулю в добавленных узлах сетки не оказывает существенного воздействия и не усложняет реализацию алгоритма. В противном случае возможно построение сглаживающего сплайна [27].
В выбранных узлах ^, к = 1,... ,п, функции Мах1С8(£) и Мт1С8(£) принимают экстремальные значения, совпадающие со значениями какой-либо траектории задачи (1), (2) в этих узлах. В промежуточных точках интерполяции к значениям интерполирующих функций нужно добавить или отнять оценку ошибки интерполяции (вычисляемая величина) согласно (8). Ниже на шаге 4 описан способ контроля за оценкой отклонения решений в пучке траекторий. Поскольку используется символьная формула точного решения, ошибок приближенного интегрирования системы при реализации этого алгоритма нет. Далее в алгоритме строятся формулы первообразных
1Мах1С8^) = / Мах1С8(т )йт, 1Мт1С8(*) = / Мт1С8(т )йт.
Jto
Ш!аг 4. Для контроля точности алгоритма оценим максимальные отклонения любых двух решений, принадлежащих пучку траекторий
у(1) = ф(1)(ф-1(ь)у0 + £ Ф-1(т)и(т.
Изменяющаяся величина в формуле решения — возмущающее воздействие и(Ь), остальные параметры решений совпадают. Цель шага 4 заключается в демонстрации эквивалентности максимальных отклонений любых двух решений и максимальных отклонений решений, вычисляемых в данном методе по формуле Коши. При этом максимальные отклонения в первом из рассматриваемых примеров не будут превосходить вычисляемые в методе отклонения произвольных значений двух решений.
Преобразуем разность двух решений этого пучка, выходящих из одной начальной точки, но имеющих разные управляющие воздействия:
-1 -1
У1&) - У2$) = Ф(1)Ф-1(Ь)уо + Ф(*) / Ф-1(т)щ(т)б,т-
- Ф(í)Ф-1(íо)yо - Ф(*) / Ф-1(т)и2(т)йт =
■По
= Ф(£)(/ Ф-1 (г)щ(т)йт -I Ф-1(г)и2(т)(1т] =
о о
= Ф&) [ Ф-1(т)(т(т) - и2(т))йт. Ло
Используя свойства интеграла Лебега в этом выражении, оценим покомпонентно модуль разности двух решений:
Ы!) - У2$) < 2ЩЩ ^ 1Ф-1(т)11и1(т) - и2(т)1йт.
о
Далее,
|Ф-1(т)(т(т) - и2(т)| < |Ф-1(т)||т(т) - и2(т)| < 21Ф-1(т)|.
В этом неравенстве учтено, что и1(Ь) € [-1,1], и2(Ь) € [-1,1]. Отсюда для любого £ справедлива оценка |щ(Ь) - и2(¿)| ^ 2.
Это означает, что отклонение любых двух решений пучка траекторий покомпонентно оценивается так:
|У1(1) - У2(Щ < 21Ф(1)1 1Ф-1(т)№. (9)
о
Символьные формулы фундаментальной матрицы решений и обратной к ней матрицы были определены ранее в алгоритме, следовательно, для нахождения отклонений решений необходимо вычислить правую часть неравенства (9). Она применяется для контроля оценок областей решений, полученных в описанном методе.
Метод, представленный в этой статье, оценивает формулу множества точных решений как формулу объединения всех выражений, полученных по формуле Коши:
1
и у(1,1о,уо)= и Фад(ф-1(£о)Уо +[ Ф-1(т)и(т)йЛ. (10) «€[-1,1] «€[-1,1] ^ 1о '
Каждый элемент объединения (10) является точным решением задачи (1), (2). Основная сложность реализации алгоритма заключается в оценке области значений на основе символьных формул этого выражения. Эта сложность преодолена — этап описан на шаге 3.
При этом значение оценок выражения вычисляется сразу в выбранный момент времени. Нет необходимости накапливать (суммировать) множества значений выражений, что возникает при реализации любого численного метода оценки множеств решений. Это свойство является важной характеристикой алгоритма, составленного на основе формулы Коши. Явное достоинство методов, основанных на символьных формулах решений (10), заключается также в отсутствии накопления ошибок со всех временных шагов. В (10) отсутствуют какие-либо промежуточные численные значения, появляющиеся в любом численном методе вида
Ук + ЬЪ(ук,ук-1,...,Ук-т), 0 ^ к ^ И, 0 ^ т ^ к.
Вычисления максимальных и минимальных значений областей решений на каждом шаге времени в методе, описанном в статье, выполняются независимо от предыдущих шагов. Это и позволяет проверить точность символьных методов оценки множеств решений.
ШШаг 5. Печать таблиц значений и построение графиков решений.
Для анализа вычислительной сложности символьного метода на основе оператора Коши отметим следующее. Шаг 3 алгоритма требует пересчета векторной численно-символьной величины Ф(£)Ф(т)и(т) для каждого момента времени £ и определения ее экстремальных значений для всех т, пробегающих узлы сетки от 0 до ¿. Следовательно, асимптотическая зависимость числа операций от количества узлов сетки будет оцениваться как 0(М(Ж - 1)/2), где N — число узлов сетки. Поскольку система является линейной и в методе не используется рекуррентная зависимость решения в каждый момент времени от неопределенных параметров (известна формула решения — формула Коши), в каждом узле сетки выполняется число операций с асимптотикой О(п), где п — размерность системы.
Задача определения экстремумов символьной формулы здесь сводится к линейному программированию на гиперкубе. Сложность этой задачи зависит от числа независимых управляющих воздействий, прилагаемых к системе. Для распространенных практических задач в численных экспериментах подобная оптимизация не превышала сотых долей секунды при использовании СРи массового сегмента.
Данный алгоритм требует хранить в памяти на каждом шаге лишь фундаментальную матрицу решений в символьном виде, текущую формулу сплайна и векторы верхних/нижних оценок в численном виде.
3. Примеры вычисления оценок множеств решений. Задачи оценки множеств решений систем ОДУ появляются в случае, когда на систему влияют различные неопределенные внешние возмущения: неуправляемые изменения параметров, погрешности в измерении начальных условий и коэффициентов системы.
Гарантированный (нестохастический) подход оперирует с множествами, в которых лежат неопределенные переменные, при этом предполагается, что неизвестные возмущения локализованы в известных множествах, а в целом — произвольны.
Рассмотрим управляемую систему, описываемую системой ОДУ четвертого порядка
Лу1 dy2
-тт = У 2, ~тг~ = -У1 + У2 + Щ
аъ аъ (11)
йуз йуА у '
-Ж =у4 ИТ =У1 - У2 -
В системе (11) у1, уз — компоненты вектора состояния, описывающие положение управляемого объекта; у2, у4 — компоненты вектора, описывающие величину скорости управляемого объекта; на возмущения и(Ь) и w(t) наложены ограничения
и(г)^(г) е и = (Н < 1, Н < 1}.
Начальные данные для системы (11) следующие: уг(¿о) = 1, I = 1,2,3,4.
В уравнения движения справа входит вектор равнодействующей всех сил в системе.
Выбор возмущающих или управляющих воздействий системы (11) заранее неизвестен. Поэтому ставится задача определить включение Уи(Ь) для пучка траекторий точных решений системы (11):
Уи^) 5 | и € и, w € и}.
Все члены, входящие в уравнение движения, имеют одинаковые размерности, задача (особенно выраженная линейным оператором) безразмерна.
На основе символьного метода проведены оценки множеств достижимости задачи (11).
На рис. 1, 2 изображены верхние и нижние границы включений области достижимости для первой и третьей компоненты вектора решений, то есть в проекциях на оси у1, у%. На этих рисунках включены только два объекта: верхняя и нижняя границы множества решений в проекции на координатные плоскости.
Оценка множеств решений выполнена также для управляемой системы четвертого порядка [6]:
= 2У2 + У4, = -2У1 + уз + иь
аъ аъ (12)
= У2 + 10щ, — = у1 - 2уз + П2
с возмущающими силами 1и1(Ь)1 ^ 1, 1и2(Ь)1 ^ 1. В [6] построены графики роста объема эллипсоида траекторий системы, входящих в пучок точных решений задачи (12). Полученные численно-символьным методом результаты оценки множеств решений всех компонент вектора состояния задачи (12) совпадают с качественными оценками в [6]. Это подтверждает график на рис. 3.
Рассмотрим задачу оценки множества достижимости при плоском движении самолета, при котором траектории его центра масс расположена в некоторой фиксированной вертикальной плоскости, служащей плоскостью материальной симметрии самолета. Силами, действующим на самолет, являются сила тяги винта, направленная по оси винта и составляющая с хордой крыла постоянный угол, сила тяжести и аэродинамические силы. Совокупность последних может быть приведена к главному вектору. Движение летательного аппарата в вертикальной плоскости на некотором интервале времени можно описать следующей линейной системой обыкновенных дифференциальных уравнений с постоянными коэффициентами [28]:
¿у йх
— = ацу + айв, — = а2зУ + а2^в,
¥ м (13)
ау ав у '
— = а351и\ + азо, -¡- = а4 5и + а4о,
аЪ т
где у, х — координаты местоположения летательного аппарата, V — скорость, в — угол наклона траектории, а^ —динамические коэффициенты. Система (13) получена в [28] из нелинейных дифференциальных уравнений движения летательного аппарата заменой правых частей функциями линейными
Рис. 1. Верхняя и нижняя границы множества достижимости первой компоненты вектора состояний на интервалах [0,4] и [0, 8] для управляемой системы с возмущениями (11)
[Figure 1. The upper and lower bounds of the reachable set of the first component of the state vector are defined on the intervals [0,4] and [0, 8] for the controlled system with
perturbations (11)]
Уз
Рис. 2. Верхняя и нижняя границы множества достижимости третьей компоненты вектора состояний на интервалах [0, 4] и [0, 8] для управляемой системы с возмущениями (11)
[Figure 2. The upper and lower bounds of the reachable set of the third component of the state vector are defined on the intervals [0,4] and [0, 8] for the controlled system with
perturbations (11)]
VI
Рис. 3. Верхняя и нижняя границы множества достижимости первой компоненты вектора состояний на интервале [0,16] для управляемой системы с возмущениями (12)
[Figure 3. The upper and lower bounds of the reachable set of the first component of the state vector are defined on the intervals [0,16] for the controlled system with perturbations (12)]
У
[Figure 4. The upper and lower bounds of the set of attainable values of the first component of the state vector in the vertical takeoff problem within the interval [0,10]]
относительно фазовых координат и управлений. Управляющее воздействие и — значения угла атаки, удовлетворяющие ограничению
^шах ^ и ^ 0,max, &шах — 0.25.
Динамические коэффициенты в работе [28] заданы как величины, имеющие размерности. Они принимают значения а13 = 0, а14 = 1000 м/с, а23 = 1, а24 = -1 м/с, а35 = -50 м/с2, азо = -20 м/с2, = 0.25 1/с, = 0. Требуется построить оценку области достижимости летательного аппарата в вертикальной плоскости Оху для заданного момента времени.
При построении границ множества достижимости использовался символьный метод, описанный в этой статье. Для символьного приближения интеграла в формуле Коши применялось два алгоритма. Один алгоритм расчета основан на аппроксимации управляющего воздействия кусочно-постоянной функцией, второй — на приближении интеграла символьной квадратурной формулой. Динамические коэффициенты, например скорости, принимают большие значения, что означает возможность быстрого роста границ множества значений координат летательного аппарата. Численные значения границ множества решений этой задачи подтверждаются при сравнении со значениями формул точных решений (рис. 4).
Заключение. В статье описан метод оценки множеств решений систем ОДУ с возмущающими воздействиями. Для оценки границ множеств решений строятся формулы общего решения, формулы фундаментальной матрицы решений системы ОДУ и формулы первообразной от функции, в которую входят неточно заданные компоненты (возмущающие воздействия).
Разработан численно-символьный алгоритм построения формул, описывающих точные решения систем ОДУ, организующий упрощения этих формул и вычисления их экстремальных значений. В алгоритм также входят неравенства (9), связывающие максимальные отклонения решений и возмущающие воздействия на шаге 4 алгоритма. Результаты этого этапа нужны для контроля за точностью численно-символьного алгоритма. Построенный численно-символьный алгоритм применим на всем интервале времени, на котором существуют решения всех систем ОДУ, входящих в поставленную задачу оценки множеств решений.
Принципиальное отличие всех известных ранее методов и построенных численно-символьных алгоритмов заключается в том, что вычисление максимальных и минимальных значений границ множеств решений выполняется независимо от предыдущих шагов на каждом шаге по времени. Для известных методов для каждого шага времени к найденным на предыдущем шаге максимальным и минимальным границам множеств решений добавляются значения некоторых выражений, различных для разных методов, что приводит к накоплению ошибок. Ошибка вычислений границ множеств решений символьно-численным методом от шага к шагу не накапливается, такая возможность роста ошибок отсутствует.
Конкурирующие интересы. Конкурирующих интересов не имею. Авторский вклад и ответственность. Я несу полную ответственность за предоставление окончательной версии рукописи в печать. Окончательная версия рукописи мною одобрена.
Финансирование. Работа выполнена в рамках проекта, поддержанного за счет
средств гранта № 20-41-240002, предоставленного РФФИ, Правительством Красноярского края, Краевым фондом науки.
Библиографический список
1. Абгарян К. А. Об устойчивости движения на конечном промежутке времени // Докл. АН СССР, 1968. Т. 183, №3. С. 527-530.
2. Абгарян К. А. Введение в теорию устойчивости движения на конечном интервале времени. М.: Наука, 1991. 160 с.
3. Поляк Б.Т., Хлебников М.В., Щербаков П.С. Линейные матричные неравенства в системах управления с неопределенностью// Автомат. и телемех., 2021. №1. С. 3-54. EDN: VNZYOS. DOI: https://doi.org/10.31857/S0005231021010013.
4. Куржанский А. Б. Управление и наблюдение в условиях неопределенности. М.: Наука, 1977. 392 с. EDN: YOIHKV.
5. Куржанский А. Б., Филиппова Т. Ф. Об описании пучка выживающих траекторий управляемой системы// Диффер. уравн., 1987. Т. 23, №8. С. 1303-1315. EDN: YMJYEV.
6. Черноусько Ф. Л. Оценивание фазового состояния динамических систем. Метод эллипсоидов. М.: Наука, 1988. 320 с. EDN: TTDUCV.
7. Черноусько Ф. Л. Эллипсоидальные аппроксимации множеств достижимости управляемых линейных систем с неопределенной матрицей // ПММ, 1996. Т. 60, №6. С. 940-950.
8. Куржанский А. Б., Месяц А. И. Управление эллипсоидальными траекториями. Теория и вычисления// Ж. вычисл. матем. и матем. физ., 2014. Т. 54, №3. С. 404-414. EDN: RWANON. DOI: https://doi.org/10.7868/S0044466914030120.
9. Chernous'ko F. L. Optimal ellipsoidal estimates of control and uncertain systems (survey) // Applied and Computational Mathematics, 2009. vol.8, no. 2. pp. 135-151. EDN: MWVJOL.
10. Ушаков В. Н., Ершов А. А. Множества достижимости и интегральные воронки зависящих от параметра дифференциальных включений// Докл. РАН. Матем., ин-форм., проц. упр., 2021. Т. 499. С. 49-53. EDN: ZMAGQA. DOI: https://doi.org/10.31857/ S2686954321040159.
11. Ушаков В. Н., Ершов А. А., Ушаков А. В. Управляемые системы, зависящие от параметра: множества достижимости и интегральные воронки// ПММ, 2022. Т. 86, №1. С. 186-205. EDN: RAXDBK. DOI: https://doi.org/10.31857/S0032823522010088.
12. Aubin J.-P., Frankowska H. The value does not exist! A motivation for extremal analysis // Probab. Uncertain. Quant. Risk, 2022. vol.7, no. 3. pp. 195-214. DOI: https://doi.org/ 10.3934/puqr.2022013.
13. Althoff M., Frehse G., Girard A. Set propagation techniques for reachability analysis// Annu. Rev. Control Robot. Auton. Syst., 2021. vol.4. pp. 369-395. DOI: https://doi.org/ 10.1146/annurev-control-071420-081941.
14. Villegas Pico H. N., Alipantis D. C. Reachability analysis of linear dynamic systems with constant, arbitrary, and Lipschitz continuous inputs // Automatica, 2018. vol. 95. pp. 293305. DOI: https://doi.org/10.1016/j.automatica.2018.05.026.
15. Морозов А.Ю., Ревизников Д. Л. Алгоритм адаптивной интерполяции на разреженных сетках для численного интегрирования систем обыкновенных дифференциальных уравнений с интервальными неопределенностями// Диффер. уравн., 2021. Т. 57, №7. С. 976-988. EDN: YOGCJE. DOI: https://doi.org/10.31857/S0374064121070104.
16. Гидаспов В. Ю., Морозов А. Ю., Ревизников Д. Л. Алгоритм адаптивной интерполяции с использованием TT-разложения для моделирования динамических систем с интервальными параметрами// Ж. вычисл. матем. и матем. физ., 2021. Т. 61, №9. С. 1416-1430. EDN: XIHOGZ. DOI: https://doi.org/10.31857/S0044466921090106.
17. Морозов А. Ю., Ревизников Д. Л. Интервальный подход к решению задач параметрической идентификации динамических систем // Диффер. уравн., 2022. Т. 58, №7. С. 962976. EDN: CEMWGM. DOI: https://doi.org/10.31857/S0374064122070081.
18. Новиков В. А., Рогалев А. Н. Построение сходящихся верхних и нижних оценок решений систем обыкновенных дифференциальных уравнений с интервальными начальными данными // Ж. вычисл. матем. и матем. физ., 1993. Т. 33, №2. С. 219-231.
19. Рогалев А. Н., Рогалев А. А. Численные оценки предельных отклонений летательных аппаратов в атмосфере // Вестн. СибГАУ, 2016. Т. 16, №1. С. 104-112. EDN: TRIUWN.
20. Рогалев А. А. Алгоритмы символьных вычислений на основе корневых деревьев для оценки возможностей управления // Сибирский журнал науки и технологий, 2017. Т. 18, №4. С. 810-819. EDN: YNZVUQ.
21. Rogalev A. N., Rogalev A. A., Feodorova N. A. Numerical computations of the safe boundaries of complex technical systems and practical stability // J. Phys.: Conf. Ser., 2019. vol.1399, 033112. EDN: VXANUU. DOI: https://doi.org/10.1088/1742-6596/1399/3/ 033112.
22. Rogalev A. N., Rogalev A. A., Feodorova N. A. Malfunction analysis and safety of mathematical models of technical systems// J. Phys.: Conf. Ser., 2020. vol.1515, 022064. EDN: GIFBHY. DOI: https://doi.org/10.1088/1742-6596/1515/2Z022064.
23. Rogalev A. N. Regularization of inclusions of differential equations solutions based on the kinematics of a vector field in stability problems// J. Phys.: Conf. Ser., 2021. vol.2099, 012045. EDN: WQZHXJ. DOI: https://doi.org/10.1088/1742-6596/2099/1/012045.
24. Смирнов А. В. О билинейной сложности и практических алгоритмах умножения матриц// Ж. вычисл. матем. и матем. физ., 2013. Т. 53, №12. С. 1970-1984. EDN: RLWJJV. DOI: https://doi.org/10.7868/S0044466913120168.
25. Абрамов С. А., Рябенко А. А., Хмельнов Д. Е. Регулярные решения линейных обыкновенных дифференциальных уравнений и усеченные ряды // Ж. вычисл. матем. и матем. физ., 2020. Т. 60, №1. С. 4-17. EDN: AHLUDF. DOI: https://doi.org/ 10.31857/S0044466920010020.
26. Галеев Э. М. Оптимизация: теория, примеры, задачи. М.: УРСС, 2002. 304 с.
27. Ahlberg J. H., Nilson E. N., Walsh J. L. The Theory of Splines and their Applications / Mathematics in Science and Engineering. New York: Academic Press, 1967. xi+284 pp. DOI: https://doi.org/10.1016/s0076-5392(08)x6115-6.
28. Толпегин И. Г. Дифференциально-игровые методы наведения ракет на скоростные маневрирующие цели // Изв. РАРАН, 2003. №1. С. 80-86.
Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki
[J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2023, vol. 27, no. 2, pp. 357-374 ISSN: 2310-7081 (online), 1991-8615 (print) d https://doi.org/10.14498/vsgtu1978
MSC: 65L05, 65L07, 65G40, 34A34, 65L70
The estimation of solutions sets of linear systems of ordinary differential equations with perturbations based on the Cauchy operator
A. A. Rogalev
Institute of Space and Information Technologies, Siberian Federal University, 26 k/1, Academician Kirensky str., Krasnoyarsk, 660074, Russian Federation.
Abstract
The paper outlines a method for numerical analysis of sets of solutions for linear systems of ordinary differential equations that contain perturbations in the right-hand side. The method determines extreme values of the solutions, which comprise the sets of solutions along the coordinate axes or in a specified direction. The estimations are based on using the Cauchy operator, written with symbolic formulas for variations of arbitrary constants. Additionally, control is implemented over the deviation of solutions when calculating a bundle of trajectories. The paper also is devoted to examples of estimating reachability sets of systems under the influence of control and disturbance effects.
Keywords: perturbations, trajectories tube, solutions sets, symbolic algorithms, symbolic formulae.
Received: 30th December, 2022 / Revised: 14th April, 2023 / Accepted: 25th May, 2023 / First online: 23rd June, 2023
Competing interests. I have no conflicting interests.
Authorship contribution and responsibility.. I take full responsibility for submitting the final version of the manuscript for publication. The final version of the manuscript has been approved by me.
Funding. The work was carried out within the framework of a project supported by grant no. 20-41-240002 provided by the RFBR, the Government of the Krasnoyarsk Territory, and the Regional Science Foundation.
Mathematical Modeling, Numerical Methods and Software Complexes Research Article
© Authors, 2023
© Samara State Technical University, 2023 (Compilation, Design, and Layout) 3 ©® The content is published under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/) Please cite this article in press as:
Rogalev A. A. The estimation of solutions sets of linear systems of ordinary differential equations with perturbations based on the Cauchy operator, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2023, vol. 27, no. 2, pp. 357-374. EDN: MXEZAB. DOI: 10.14498/vsgtu1978 (In Russian). Author's Details:
Aleksander A. Rogalev https://orcid.org/0000- 0003- 2176- 9639 Senior Lecturer; Dept. of Information Systems; e-mail: gogoba88@mail.ru
The estimation of solutions sets of linear systems of ODE with perturbations
References
1. Abgaryan K. A. On the stability of motion within a finite interval of time, Sov. Math., Dokl., 1968, vol.9, no. 3, pp. 1110-1112.
2. Abgaryan K. A. Vvedenie v teoriiu ustoichivosti dvizheniia na konechnom intervale vremeni [Introduction to the Theory of Stability of Motion over a Finite Time Interval]. Moscow, Nauka, 1991, 160 pp. (In Russian)
3. Polyak B. T., Khlebnikov M. V., Shcherbakov P. S. Linear matrix inequalities in control systems with uncertainty, Autom. Remote Control, 2021, vol. 82, no. 1, pp. 1-40. EDN: HMIXKI. DOI: https://doi.org/10.1134/S000511792101001X.
4. Kurzhansky A. B. Upravlenie i nabliudenie v usloviiakh neopredelennosti [Management and Supervision in Conditions of Uncertainty]. Moscow, Nauka, 1977, 392 pp. (In Russian). EDN: YOIHKV
5. Kurzhansky A. B., Filippova T. F. On a description of the bundle of surviving trajectories of a control system, Differ. Uravn., 1987, vol.23, no. 8, pp. 1303-1315 (In Russian). EDN: YMJYEV.
6. Chernous'ko F. L. Otsenivanie fazovogo sostoianiia dinamicheskikh sistem. Metod ellip-soidov [Estimation of the Phase State of Dynamic Systems. Ellipsoid Method]. Moscow, Nauka, 1988, 320 pp. (In Russian). EDN: TTDUCV
7. Chernous'ko F. L. Ellipsoidal approximation of the attainability sets of a linear system with an uncertain matrix, J. Appl. Math. Mech., 1996, vol. 60, no. 6, pp. 921-931. DOI:https:// doi.org/10.1016/S0021-8928(96)00114-1.
8. Kurzhanski A. B., Mesyats A. I. Control of ellipsoidal trajectories: Theory and numerical results, Comput. Math. Math. Phys., 2014, vol.54, no. 3, pp. 418-428. EDN: SKQPZP. DOI: https://doi.org/10.1134/S0965542514030117.
9. Chernous'ko F. L. Optimal ellipsoidal estimates of control and uncertain systems (survey), Applied and Computational Mathematics, 2009, vol.8, no. 2, pp. 135-151. EDN: MWVJOL.
10. Ushakov V. N., Ershov A. A. Reachable sets and integral funnels of differential inclusions depending on a parameter, Dokl. Math., 2021, vol.104, no. 1, pp. 200-204. EDN: RZPGMU. DOI: https://doi.org/10.1134/S1064562421040153.
11. Ushakov V. N., Ershov A. A., Ushakov A. V. Control systems depending on a parameter: Reachable sets and integral funnels, Mech. Solids, 2022, vol. 57, no. 7, pp. 1672-1688. DOI: https://doi.org/10.3103/S0025654422070172.
12. Aubin J.-P., Frankowska H. The value does not exist! A motivation for extremal analysis, Probab. Uncertain. Quant. Risk, 2022, vol.7, no. 3, pp. 195-214. DOI: https://doi.org/ 10.3934/puqr.2022013.
13. Althoff M., Frehse G., Girard A. Set propagation techniques for reachability analysis, Annu. Rev. Control Robot. Auton. Syst., 2021, vol.4, pp. 369-395. DOI: https://doi.org/ 10.1146/annurev-control-071420-081941.
14. Villegas Pico H. N., Alipantis D. C. Reachability analysis of linear dynamic systems with constant, arbitrary, and Lipschitz continuous inputs, Automatica, 2018, vol. 95, pp. 293-305. DOI: https://doi.org/10.1016/j.automatica.2018.05.026.
15. Morozov A. Yu., Reviznikov D. L. Adaptive interpolation algorithm on sparse meshes for numerical integration of systems of ordinary differential equations with interval uncertainties, Differ. Equ., 2021, vol.57, no. 7, pp. 947-958. EDN: OECMNR. DOI: https://doi.org/ 10.1134/S0012266121070107.
16. Gidaspov V. Yu., Morozov A. Yu., Reviznikov D. L. Adaptive interpolation algorithm using TT-decomposition for modeling dynamical systems with interval parameters, Comput. Math. Math. Phys., 2021, vol.61, no. 9, pp. 1387-1400. EDN: WLZHVA. DOI: https://doi. org/10.1134/S0965542521090098.
17. Morozov A. Yu., Reviznikov D. L. Interval approach to solving parametric identification problems for dynamical systems, Differ. Equ., 2022, vol. 58, no. 7, pp. 952-965. EDN: PTDLNU. DOI: https://doi.org/10.1134/s0012266122070084.
18. Novikov V. A., Rogalev A. N. Construction of convergent upper and lower bounds of solutions of systems of ordinary differential equations with interval initial data, Comput. Math. Math. Phys., 1993, vol.33, no. 2, pp. 193-203. EDN: XKLCXR.
19. Rogalev A. N., Rogalev A. A. Numerical estimates of the maximum deviations of aircraft in the atmosphere, Vestn. SibGAU, 2016, vol. 16, no. 1, pp. 104-112 (In Russian). EDN: TRIUWN.
20. Rogalev A. A. Algorithms for symbolic calculations based on root trees for assessing management capabilities, Siberian Journal of Science and Technology, 201Т, vol. 1S, no. 4, pp. S10-S19 (In Russian). EDN: YNZVUQ.
21. Rogalev A. N., Rogalev A. A., Feodorova N. A. Numerical computations of the safe boundaries of complex technical systems and practical stability, J. Phys.: Conf. Ser., 2019, vol.1399, 033112. EDN: VXANUU. DOI: https://doi.org/10.1088/1742-6596/1399/3/ 033112.
22. Rogalev A. N., Rogalev A. A., Feodorova N. A. Malfunction analysis and safety of mathematical models of technical systems, J. Phys.: Conf. Ser., 2020, vol. 1515, 022064. EDN: GIFBHY. DOI: https://doi.org/10.1088/1742-6596/1515/2Z022064.
23. Rogalev A. N. Regularization of inclusions of differential equations solutions based on the kinematics of a vector field in stability problems, J. Phys.: Conf. Ser., 2021, vol.2099, 012045. EDN: WQZHXJ. DOI: https://doi.org/10.1088/1742-6596/2099/1/012045.
24. Smirnov A. V. The bilinear complexity and practical algorithms for matrix multiplication, Comput. Math. Math. Phys., 2013, vol.53, no. 12, pp. 1TS1-1T95. EDN: SLOFHJ. DOI: https://doi.org/10.1134/S0965542513120129.
25. Abramov S. A., Ryabenko A. A., Khmelnov D. E. Regular solutions of linear ordinary differential equations and truncated series, Comput. Math. Math. Phys., 2020, vol. 60, no. 1, pp. 1-14. EDN: KIPAWT. DOI: https://doi.org/10.1134/S0965542520010029.
26. Galeev E. M. Optimizatsiia: teoriia, primery, zadachi [Optimization: Theory, Examples, Tasks]. Moscow, URSS, 2002, 304 pp. (In Russian)
27. Ahlberg J. H., Nilson E. N., Walsh J. L. The Theory of Splines and their Applications, Mathematics in Science and Engineering. New York, Academic Press, 196Т, xi+2S4 pp. DOI: https://doi.org/10.1016/s0076-5392(08)x6115-6.
28. Tolpegin I. G. Differential-game methods for targeting missiles at high-speed maneuvering targets, Izv. RARAN, 2003, no. 1, pp. S0-S6 (In Russian).
