Об успокоении системы управления нейтрального типа на временном графезвезде с запаздыванием, пропорциональным времени
- Авторы: Леднов А.П.1,2,3
-
Учреждения:
- Саратовский национальный исследовательский государственный университет им. Н.Г. Чернышевского
- Московский центр фундаментальной и прикладной математики
- Московский государственный университет им. М.В. Ломоносова
- Выпуск: Том 71, № 4 (2025)
- Страницы: 642-654
- Раздел: Статьи
- URL: https://journals.rudn.ru/CMFD/article/view/48172
- DOI: https://doi.org/10.22363/2413-3639-2025-71-4-642-654
- EDN: https://elibrary.ru/MGDOHS
- ID: 48172
Цитировать
Полный текст
Аннотация
На временно´м графе-звезде рассматривается задача об оптимальном успокоении системы управления для обобщенного уравнения пантографа, представляющего собой уравнение нейтрального типа с запаздыванием, пропорциональным времени. Запаздывание в системе распространяется через внутреннюю вершину графа. Исследуется вариационная задача минимизации функционала энергии с учетом вероятностей сценариев, соответствующих различным ребрам. Установлено, что оптимальная траектория удовлетворяет условиям типа Кирхгофа во внутренней вершине. Доказана эквивалентность вариационной задачи некоторой краевой задаче для функционально-дифференциальных уравнений второго порядка на графе и установлена однозначная разрешимость обеих задач.
Полный текст
1. Введение Дифференциальные операторы на графах, часто называемые квантовыми графами, активно изучаются с прошлого века в связи с моделированием различных процессов, протекающих в сложных системах, представимых в виде пространственных сетей [7, 10, 17, 20, 21]. Для таких моделей помимо условий непрерывности в вершинах характерны также условия Кирхгофа. Для задания на графах функционально-дифференциальных операторов с запаздыванием, С.А. Бутериным в работе [11] была предложена концепция глобального запаздывания. Последнее означает, что запаздывание распространяется через внутренние вершины графа. Другими словами, решение уравнения на входящем ребре служит начальной функцией для уравнений на © А.П. Леднов, 2025 This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License https://creativecommons.org/licenses/by-nc/4.0/legalcode 642 исходящих ребрах. Глобальное запаздывание стало альтернативой локально нелокальному случаю, рассмотренному в [24], когда уравнение на каждом ребре имеет свой собственный параметр запаздывания и может быть решено отдельно от уравнений на остальных ребрах. Использование концепции глобального запаздывания позволило перенести на графы класс задач об успокоении управляемых систем с последействием. Впервые задача этого типа была поставлена и исследована на интервале Н.Н. Красовским [6] для системы управления с постоянным запаздыванием, описываемой уравнением запаздывающего типа. Эта задача получила дальнейшее развитие в работах А.Л. Скубачевского [9, 23] и позже в работах других авторов (см. [1] и литературу там), где рассматриваемая система управления имеет нейтральный тип, т. е. содержит запаздывание и в главных членах. Это существенно усложняет задачу и, в частности, приводит к понятию обобщенного решения соответствующей краевой задачи для оптимальной траектории. С.А. Бутерин в работах [3, 12] распространил на графы задачу об успокоении систем управления с постоянным запаздыванием. В [12] рассмотрен случай уравнения первого порядка запаздывающего типа, а в [3] - общий случай нестационарной управляемой системы, уравнения которой относятся к нейтральному типу и имеют произвольный порядок. Рассмотрение указанной задачи на графах, в свою очередь, привело к концепции временн´ого графа, ребра которого, в отличие от пространственной сети, отождествляются с промежутками времени, а каждая внутренняя вершина понимается как точка разветвления процесса, дающая несколько различных сценариев дальнейшего его протекания. В [3-5, 12] показано, что на временны´х графах также могут возникать условия Кирхгофа. А именно, им будет удовлетворять траектория течения процесса, являющаяся оптимальной с учетом сразу всех сценариев. Кроме того, в [5, 12] была предложена стохастическая интерпретация системы управления на временн´ом дереве. В частности, к системе на дереве приведет замена коэффициентов в уравнении на интервале дискретными случайными процессами с дискретным временем. В работах [18, 19] на графы была перенесена задача об успокоении системы управления, описываемой так называемым уравнением пантографа [16]. В данном случае запаздывание не постоянно, а является пропорциональным времени сжатием. Рассматривалась система управления, задаваемая классическим уравнением пантографа , (1.1) где b,c ∈ R, q > 1. Уравнение вида (1.1) широко применяется в прикладных задачах. Так, например, это уравнение используется при моделировании динамики контактного провода электроснабжения подвижного состава [22]. Для q ∈ (0,1) оно возникает в астрофизике при описании поглощения света межзвездной материей [2], а также в биологии при моделировании процесса роста клеток [13]. В данной работе мы переходим к рассмотрению системы управления для обобщенного уравнения пантографа , (1.2) где a,b,c ∈ R, q > 1. Уравнение вида (1.2) изучалось в [14, 15] и ряде других работ. Были получены различные представления решения и показано, что разрешимость задачи (1.2) зависит от коэффициента a и от класса гладкости решений. В частности, при, существует единственное решение в C∞[0,+∞); при этом в зависимости от значения a могут существовать и другие C1-решения, не принадлежащие C∞[0,+∞). На интервале задача об успокоении системы управления для обобщенного уравнения пантографа была рассмотрена Л.Е. Россовским в работе [8], где исследовалась следующая система управления нейтрального типа: , (1.3) (1.4) где a,b,c ∈ R, q > 1, а u(t) - управляющее воздействие, которое является вещественнозначной функцией; состояние системы в начальный момент времени задается условием (1.4). Задача управления формулируется следующим образом: требуется найти u(t), приводящее систему (1.3), (1.4) в равновесие y(t) = 0 при t T для некоторого T > 0. Для этого достаточно найти u(t) ∈ L2(0,T), приводящее систему в состояние (1.5) а затем сбросить управление, положив u(t) ≡ 0 при t > T. При этом из всех возможных управлений ищется управление, обладающее минимальной энергией . В результате получается вариационная задача о минимизации квадратичного функционала (1.6) на множестве функций y (t) ∈ W21[0,T], удовлетворяющих краевым условиям (1.4), (1.5). Исследование вариационной задачи (1.4)-(1.6) включает сведение ее к эквивалентной краевой задаче для функционально-дифференциального уравнения второго порядка с растяжением и сжатием аргумента. В частности, установлено, что если функция y (t) ∈ W21[0,T], удовлетворяющая условиям (1.4), (1.5), минимизирует функционал (1.6), то она является решением краевой задачи для уравнения при краевых условиях (1.4) и (1.5). При этом, поскольку задача (1.4), (1.5), (1.7) может не иметь решения в, ее решение является обобщенным в смысле выполнения условия . Обратное утверждение также верно: если y (t) ∈ W21[0,T] является обобщенным решением задачи (1.4), (1.5), (1.7), то y доставляет минимум функционалу (1.6). Следующая теорема при предположении, обеспечивающем коэрцитивность функционала J(y), устанавливает существование и единственность обобщенного решения краевой задачи (1.4), (1.5), (1.7) и, стало быть, однозначную разрешимость вариационной задачи (1.4)-(1.6). Теорема 1.1 (см. [8]). Пусть . Тогда задача (1.4), (1.5), (1.7), имеет единственное обобщенное решение y ∈ W21[0,T]. В следующем разделе дается постановка вариационной задачи на графе-звезде. Далее, в третьем разделе, следуя общей стратегии для интервала, устанавливается эквивалентность вариационной задачи некоторой краевой задаче для функционально-дифференциальных уравнений второго порядка на графе. В заключительном разделе доказывается однозначная разрешимость обеих задач. 2. Постановка вариационной задачи на графе-звезде Рассмотрим изображенный на рис. 1 граф типа звезды Γm, состоящий из m ребер. Как обычно, под функцией y на графе Γm будем понимать кортеж y = [y1,...,ym], в котором компонента yj определена на ребре ej, т. е. yj = yj (t), t ∈ [0,Tj]. Пусть до момента времени t = T1 > 0, ассоциированного с единственной внутренней вершиной v1 графа Γm, наша система управления с запаздыванием, пропорциональным времени, на Γm описывается уравнением , (2.1) заданным на ребре e1 графа Γm, с начальным условием y1(0) = y0. (2.2) При t = T1, т. е. в вершине v1, система разветвляется на m - 1 независимых параллельных процессов, описываемых уравнениями , (2.3) но имеющих общую историю, определяемую уравнением (2.1) с начальным условием (2.2) и условиями - (2.4) а также условиями непрерывности в вершине v1, которые в данном случае согласуются с (2.4) при t → -0: (2.5) Как и в случае интервала, мы предполагаем, что q > 1, y0 ∈ R и все aj, bj, cj ∈ R. Рис. 1. Граф Γm Fig. 1. Graph Γm В (2.3) j-ое уравнение задано на ребре ej графа Γm, представляющем собой, вообще говоря, бесконечный луч, выходящий из вершины v1. Условия (2.4) означают, что запаздывание распространяется через вершину v1. Пример 2.1. Пусть m = 2, a := a1 = a2, b := b1 = b2, c := c1 = c2, и , , Тогда система управления (2.1)-(2.5) принимает вид (1.3), (1.4). Предположим для определенности, что Tj > (q - 1)T1 при всех Для успокоения системы (2.1)-(2.5) сразу при всех сценариях нужно привести ее в состояние (2.6) выбрав подходящие управления Тогда, положив uj (t) ≡ 0 при будем иметь yj(t) = 0 при тех же t и j. Другими словами, система будет приведена в равновесие на каждом ребре, выходящем из вершины v1. Поскольку такие uj (t) не единственны, будем искать их, минимизируя усилия. Кроме того, аналогично тому, как это было сделано в [12] для случая постоянного запаздывания, мы можем регулировать степень участия каждого в соответствующем функционале энергии, выбирая определенный положительный вес αj. Таким образом, приходим к вариационной задаче (2.7) при условиях (2.2), (2.4)-(2.6), где фиксированы. Для выбора весов можно применить вероятностный подход в соответствии с интерпретацией системы управления на временн´ом графе, предложенной в [12]. А именно, для этого нужно положить α1 = 1, а в качестве взять вероятности сценариев, задаваемых соответствующими уравнениями в (2.3). Тогда α2 + ... + αm = 1. Последнее тождество также обеспечивает соответствие случаю интервала, если уравнения в (2.3) не зависят от j, т. е. являются искусственными копиями единственного возможного сценария (см. [12, пример 2]). Заметим, что условия (2.4) никаких ограничений на функцию y = [y1,...,ym] не накладывают. Поэтому условимся, что взятие J (y), равно как и при от какой бы то ни было функции y на Γm автоматически подразумевает применение условий (2.4). Для краткости также введем обозначение. 3. Сведение к краевой задаче Рассмотрим вещественное гильбертово пространство со скалярным произведением где - скалярное произведение в - скалярное произведение в . Обозначим через W замкнутое подпространство W2 (Γm), состоящее из всех наборов [y1,...,ym], удовлетворяющих условиями (2.5), (2.6) и y1 (0) = 0. m Также введем пространство . Лемма 3.1. Если y ∈ W21 (Γm) является решением вариационной задачи (2.2), (2.4)-(2.7), то . (3.1) Обратно, если для некоторого y ∈ W21 (Γm) выполняется (2.2), (2.5), (2.6) и (3.1), то y является решением задачи (2.2), (2.4)-(2.7). Доказательство. Пусть y ∈ W21 (Γm) - решение задачи (2.2), (2.4)-(2.7). Тогда для произвольной фиксированной функции w ∈ W сумма y + sw принадлежит W21 (Γm) при любом s ∈ R и удовлетворяет условиям (2.2), (2.5), (2.6). Нетрудно видеть, что J (y + sw) = J (y) + 2sB (y,w) + s2J (w). Так как для всех s ∈ R, то выполняется (3.1). Обратно, пусть y ∈ W21 (Γm) удовлетворяет условиям (2.2), (2.5) и (2.6). Тогда (3.1) влечет . Таким образом, y доставляет минимум функционалу (2.7) при условиях (2.2), (2.5) и (2.6). Преобразуем (3.1), сделав замену переменных в членах, содержащих и . В результате выражение для B (y,w) примет вид . Применяя (2.4) к w = [w1,...,wm] ∈ W, можно представить Тогда перепишем (3.1) в эквивалентном виде где (3.3) при ν ∈ {0,1} и Обозначим через B краевую задачу для функционально-дифференциальных уравнений второго порядка (3.4) при условиях (2.2), (2.4)-(2.6) и условии типа Кирхгофа , (3.5) где, определены в (3.3), а выражения имеют вид (3.6) тогда как (3.7) Определение 3.1. назовем обобщенным решением краевой задачи B, если а функции удовлетворяют уравнениям (3.4) и условиям (2.2), (2.4)-(2.6), (3.5). Имеет место следующее утверждение. Лемма 3.2. Если y ∈ W21 (Γm) удовлетворяет условиям (2.2), (2.5), (2.6) и (3.1), то y является обобщенным решением краевой задачи B. Обратно, любое обобщенное решение краевой задачи B подчиняется условию (3.1). Доказательство. Пусть y ∈ W21 (Γm) и удовлетворяет условиям (2.2), (2.5), (2.6) и (3.1). Учитывая, что (3.1) эквивалентно (3.2) и применяя лемму 2 из [12] к (3.2) вместе с (3.7), получаем, что и выполняется (3.5). Кроме того, интегрируя по частям в (3.2) и используя (2.5), (2.6), будем иметь . (3.8) В силу (3.5), а также произвольности wj, из (3.8) получаем (3.4). Обратно, пусть y - обобщенное решение задачи B. Тогда левое равенство в (3.8) дает (3.1). Объединив леммы 3.1 и 3.2, получаем следующий результат. Теорема 3.1. Функция y ∈ W21 (Γm) является решением вариационной задачи (2.2), (2.4)- (2.7) тогда и только тогда, когда y является обобщенным решением краевой задачи B. 4. Однозначная разрешимость В данном разделе устанавливается однозначная разрешимость краевой задачи B, а согласно теореме 3.1 -и вариационной задачи (2.2), (2.4)-(2.7). Введем обозначения где автоматически предполагается (2.4). Лемма 4.1. Существуют C0 и C1 такие, что , . (4.3) Доказательство. Пусть w ∈ W. Используя (4.1), (4.2) и неравенство , (4.4) для n = 2, получаем , . Учитывая (2.4) применительно к w, для ν = 0,1 вычисляем В частности, это дает вторую оценку в (4.3). Наконец, аналогично лемме 5 в [12] нетрудно показать, что порождает норму в W, эквивалентную норме . Таким образом, используя оценку , получаем первое неравенство в (4.3). Лемма 4.2. Пусть. Тогда существует C2 > 0 такое, что . (4.5) Доказательство. Предположим от противного, что найдутся при s ∈ N, такие что и . (4.6) Используя первое выражение в (4.1) и (4.6), для всех s ∈ N получаем . (4.7) Применим неравенство Коши-Буняковского к среднему интегралу: . Подставляя эту оценку в (4.7), имеем . (4.8) С другой стороны, согласно (4.1), для всех таких что, почти всюду на (0,Tj) справедливы оценки . (4.9) Используя (4.1), (4.4), (4.6) и (4.9), при всех получаем оценки если aj = 0, , если, где. Учитывая, что из полученных неравенств следует при s → ∞ для Более того, если , то для таких j дополнительно имеем при s → ∞. Последнее, вместе с тождествами и оценками (4.8), позволяет заключить, что при s → ∞. В итоге имеем при s → ∞, где . Поскольку порождает норму в W, эквивалентную норме , то и при s → ∞. Последнее противоречит . Лемма 4.3. Пусть. Тогда существует C3 > 0 такое, что . Доказательство. Снова предположим от противного, что найдутся w(s) ∈ W, s ∈ N, такие, что , но теперь . (4.10) Из неравенства , совместно с оценками (4.3) и (4.5) получаем . (4.11) В силу компактности вложения W21 (Γm) в L2 (Γm) найдется подпоследовательность {w(nk)}k∈N, фундаментальная в L2 (Γm). Тогда неравенство (4.11) дает . Кроме того, используя (4.4) при n = 2 и (4.10), имеем . Таким образом, последовательность {w(sk)}k∈N является фундаментальной в W и сходится к некоторой функции w(0) ∈ W. В силу леммы 4.1 сходимость w(sk) к w(0) в W влечет сходимость . Следовательно, в силу (4.10) будем иметь , т. е. w(0) ∈ W удовлетворяет уравнениям Если для некоторого выполняется aj = cj = 0, то уравнение для w(0),j сводится к обычному дифференциальному уравнению, решение которого, с учетом принадлежности w(0) классу W, дает w(0),j(t) = 0 при 0 < t < Tj. Рассмотрим при которых . Поскольку при соответствующие уравнения на интервале принимают вид или Отсюда с учетом следует, что при . Продолжая дальше влево аналогичным образом, получаем w(0),j(t) = 0 Учитывая, что w(0),j(t) = w(0),1(t+T1) при (q-1 -1)T1 < t < 0, из чего следует w(0),1(t) = 0 для q-1T1 < t < T1, аналогичными рассуждениями можем получить w(0),1(t) = 0 для 0 < t < T1. В итоге имеем w(0)(t) = 0, что противоречит . Следующая теорема является основным результатом данного раздела. Теорема 4.1. Пусть . Тогда краевая задача B имеет единственное обобщенное решение y ∈ W21 (Γm). Кроме того, существует C такое, что . (4.12) Доказательство. Рассмотрим функцию Φ = [Φ1,...,Φm] ∈ W21 (Γm) такую, что В силу леммы 3.2, для того чтобы функция y ∈ W21 (Γm), удовлетворяющая условиям (2.2), (2.5), (2.6), была решением краевой задачи B, необходимо и достаточно, чтобы выполнялось условие (3.1). Другими словами, y является решением краевой задачи B, если и только если x := y - Φ ∈ W и B (Φ,w) + B (x,w) = 0 ∀w ∈ W. (4.13) Так как B (w,w) = J (w), то в силу леммы 4.3 билинейная форма (·,·) := B (·,·) является W скалярным произведением в W. Кроме того, справедлива оценка , (4.14) где. Таким образом, по теореме Рисса об общем виде линейного непрерывного функционала в гильбертовом пространстве существует единственная функция x ∈ W такая, что выполняется (4.13). Следовательно, краевая задача B имеет единственное решение y = Φ + x. Кроме того, согласно (4.13) и (4.14) имеем , что позволяет получить оценку (4.12). Замечание 4.1. Утверждение теоремы 4.1 может быть дополнено следующим результатом из [18]. А именно: в случае, когда система (2.1)-(2.5) имеет запаздывающий тип (т. е. для всех), краевая задача B имеет единственное решение; при этом y ∈ W21(Γm) ∩ W22(Γm) и также справедлива априорная оценка (4.12).Об авторах
А. П. Леднов
Саратовский национальный исследовательский государственный университет им. Н.Г. Чернышевского; Московский центр фундаментальной и прикладной математики; Московский государственный университет им. М.В. Ломоносова
Автор, ответственный за переписку.
Email: lednovalexsandr@gmail.com
ORCID iD: 0009-0002-7088-8693
Саратов, Россия; Москва, Россия
Список литературы
- Адхамова А.Ш., Скубачевский А.Л. Задача об успокоении системы управления с последействием с различным числом входов и выходов// Соврем. мат. Фундам. направл.- 2024.- 70, № 2.-С. 189- 200.- doi: 10.22363/2413-3639-2024-70-2-189-200.
- Амбарцумян В.А. К теории флуктуаций яркости в Млечном пути// Докл. АН СССР. - 1944.- 44.- С. 244-247.
- Бутерин С.А. Об успокоении системы управления произвольного порядка с глобальным последействием на дереве// Мат. заметки.- 2024.- 115, № 6.- С. 825-848.-doi: 10.4213/mzm14223.
- Бутерин С.А. Интегро-дифференциальная система управления на временном графе// Материалы Межд. конф. «Динамические системы: устойчивость, управление, дифференциальные игры». - Екатеринбург: ИММ УрО РАН, 2024.-С. 71-75.
- Бутерин С.А. Об управляемой системе на бесконечном временном дереве// Мат. заметки.- 2025.- 117, № 3.-С. 462-467.-doi: 10.4213/mzm14540.
- Красовский Н.Н. Теория управления движением. -М.: Наука, 1968.
- Покорный В.В., Пенкин О.М., Прядиев В.Л., Боровских А.В., Лазарев К.П., Шабров С.А. Дифференциальные уравнения на геометрических графах.-М.: Физматлит, 2005.
- Россовский Л.Е. Эллиптические функционально-дифференциальные уравнения со сжатием и растяжением аргументов неизвестной функции// Соврем. мат. Фундам. направл.-2014.- 54.- С. 3-138.
- Скубачевский А.Л. К задаче об успокоении системы управления с последействием// Докл. РАН. - 1994.-335, № 2.-С. 157-160.
- Berkolaiko G., Kuchment P. Introduction to Quantum Graphs.- Providence: Am. Math. Soc., 2013.
- Buterin S. Functional-differential operators on geometrical graphs with global delay and inverse spectral problems// Result. Math.- 2023.- 78.-79.-doi: 10.1007/s00025-023-01850-5.
- Buterin S. On damping a control system with global aftereffect on quantum graphs: Stochastic interpretation// Math. Methods Appl. Sci. -2025.- 48.- С. 4310-4331.- doi: 10.1002/mma.10549.
- Hall A., Wake G. A functional differential equation arising in modelling of cell growth// J. Aust. Math. Soc. Ser. B. Appl. Math.- 1989.- 30.- С. 424-435.-doi: 10.1017/S0334270000006366.
- Iserles A. On the generalized pantograph functional-differential equation// Eur. J. Appl. Math.- 1993.- С. 1-38.-doi: 10.1017/S0956792500000966.
- Iserles A., Liu Y. On neutral functional-differential equations with proportional delays// J. Math. Anal. Appl. - 1997.- 207, № 1. -С. 73-95.-doi: 10.1006/jmaa.1997.5262.
- Kato T., McLeod J. Functional-differential equation y˙ = ay(λt) + by(t)// Bull. Am. Math. Soc.- 1971.- 77, № 6.- С. 891-937.-doi: 10.1090/S0002-9904-1971-12805-7.
- Langese J., Leugering G., Schmidt J. Modelling, analysis and control of dynamic elastic multi-link structures.-Boston: Birkha¨user, 1994.
- Lednov A.P. On damping a control system on a star graph with global time-proportional delay// ArXiv.- 2025.-2503.02522.
- Lednov A.P. On damping a delay control system with global contraction on a temporal tree// ArXiv.- 2025.-2509.02608.
- Montrol E. Quantum theory on a network// J. Math. Phys.- 1970.-11, № 2.-С. 635-648.-doi: 10.1063/1.1665178.
- Nicaise S. Some results on spectral theory over networks, applied to nerve impulse transmission// Math. Model. -1987.- 9, № 6.-С. 437-449.
- Ockendon J., Tayler A. The dynamics of a current collection system for an electric locomotive// Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. -1971.- 322.-С. 447-468.-doi: 10.1098/rspa.1971.0078.
- Skubachevskii A. Elliptic functional differential equations and applications.- Basel: Birkh¨auser, 1997.
- Wang F., Yang C.-F. Traces for Sturm-Liouville operators with constant delays on a star graph// Results Math. -2021.-76.-220.-doi: 10.1007/s00025-021-01529-9.
Дополнительные файлы




