На правах рукописи
УДК 517.917
БЫКОВА ТАТЬЯНА СЕРГЕЕВНА
ЛЯПУНОВСКАЯ ПРИВОДИМОСТЬ
ЛИНЕЙНОЙ СИСТЕМЫ С ПОСЛЕДЕЙСТВИЕМ
01.01.02 дифференциальные уравнения
АВТОРЕФЕРАТ
диссертации на соискание ученой степени
кандидата физико-математических наук
Ижевск – 2005
Работа выполнена в ГОУ ВПО Ижевский государственный технический университет.
Научный руководитель: доктор физико-математических наук, профессор Тонков Евгений Леонидович
Официальные оппоненты: доктор физико-математических наук, профессор Долгий Юрий Филиппович доктор физико-математических наук, доцент Попова Светлана Николаевна
Ведущая организация: Институт математики Национальной Академии наук Беларуси 2005 г. в 1400 на заседании
Защита состоится диссертационного совета К 212.275.04 в ГОУ ВПО Удмуртский государственный университет по адресу: 426034, г. Ижевск, ул. Университетская, 1 (корп. 4), ауд. 222. E-mail: [email protected]
С диссертацией можно ознакомиться в библиотеке Удмуртского государственного университета.
Автореферат разослан 2005 г.
Ученый секретарь диссертационного совета к.ф.-м.н., доцент Петров Н.Н.
Актуальность темы. Линейная система с последействием x(t) = dA(t, s)x(t + s), t R = (, ), (1) r может иметь решения x(t), обращающиеся в нуль (с возрастанием t) по истечении конечного промежутка времени, либо не обращающиеся в нуль, но стремящиеся к нулю быстрее любой экспоненциальной функции (см. монографию Дж. Хейла1) ). Это означает, что показаln |x(t)|.
тель Ляпунова (x) = lim этого решения равен.
t t Игнорируя такие решения, мы можем задаться следующим вопросом: будет ли система (1), рассматриваемая только на множестве нетривиальных решений x(t) с конечными показателями Ляпунова (x), асимптотически подобна некоторой системе обыкновенных дифференциальных уравнений? Правда, может оказаться, что пространство решений с конечными показателями Ляпунова (дополненное, конечно, тривиальным решением, показатель Ляпунова которого заведомо равен ) бесконечномерно, а количество различных показателей таких решений по меньшей мере счетно.
Одной из актуальных задач асимптотической теории систем вида (1) является задача выяснения условий, при которых сужение системы (1) на любое конечномерное подпространство решений с конечными показателями Ляпунова, асимптотически подобно некоторой системе обыкновенных дифференциальных уравнений. Это важное свойство подобия обыкновенной системе полезно при изучении асимптотических инвариантов систем с последействием. Неявно оно отмечалось для систем с периодической по t матрицей A(t, s)2)3)4), но в общей ситуации не исследовалось.
Отметим, кроме того, что вопросам асимптотического поведения решений периодических систем с последействием и более общих периодических систем нейтрального типа посвящены работы 1) Хейл Дж. Теория функционально-дифференциальных уравнений. М.:
Мир, 1984. 421 с.
2) Stokes A. Floquet theory for functional-dierential equations // Proc. Nat. Ac.
of Sci., 48:8, 1962. PP. 1330–1334.
3) Шиманов С. Н. Некоторые вопросы теории колебаний систем с запаздыванием // В сб. Пятая летняя матем. школа, Киев, 1968. С. 473–549.
4) Мышкис А. Д. Линейные дифференциальные уравнения с запаздывающим аргументом. М.: Наука, 1972. 352 с.
Ю. Ф. Долгого и С. Н. Шиманова5) Ю. Ф. Долгого6)7), и Ю. Ф. Долгого и В. С. Тарасяна8).
Цель данной работы. Основной целью работы является выяснение условий, при которых сужение системы (1) на всякое конечномерное подпространство существенных (то есть имеющих конечные показатели А. М. Ляпунова) решений асимптотически подобно системе обыкновенных дифференциальных уравнений с ограниченной матрицей коэффициентов. Кроме того, в работе изучаются условия, при которых пространство существенных решений системы (1) конечномерно и исследуется вопрос о равномерной экспоненциальной устойчивости системы (1) и грубости свойства равномерной экспоненциальной устойчивости.
Научная новизна. Построены основы теории асимптотического подобия систем с последействием и систем обыкновенных дифференциальных уравнений.
Теоретическая и практическая ценность. Работа носит теоретический характер. Результаты работы могут быть в дальнейшем применены при исследовании задач управления асимптотическими инвариантами управляемых систем с последействием.
Апробация работы. Основные результаты работы докладывались на городском семинаре по дифференциальным уравнениям и теории управления (Ижевск, 2002 – 2005 годы), на конференции Функционально-дифференциальные уравнения и их приложения, посвященной 80-летию профессора Н. В. Азбелева (Ижевск, 2002), на Международной конференции, посвященной 100-летию А. Н. Колмогорова (Тамбов, ОПУ 2003), на Воронежской весенней математической школе Современные методы теории краевых задач ПонтряДолгий Ю. Ф., Шиманов С. Н. Устойчивость периодической системы дифференциальных уравнений нейтрального типа // Устойчивость и нелинейные колебания. Свердловск. 1982. С.32–39.
6) Долгий Ю. Ф. Асимптотика собственных чисел оператора монодромии для периодических уравнений с запаздыванием // Изв. ВУЗ-ов. Математика. 1994.
№ 11. С.64–72.
7) Долгий Ю. Ф. Устойчивость периодических дифференциальных уравнений с отклоняющимся аргументом. Дис. на соискание степени доктора физ.-мат. наук.
Екатеринбург 1994. 296 с.
8) Долгий Ю. Ф., Тарасян В. С. Конечномерные операторы монодромии для периодических систем дифференциальных уравнений с последействием // Изв.
УрГУ. Математика. 2000. № 18. С. 18–27.
гинские чтения XV, посвященной 65-летию академика В. А. Садовничего (Воронеж, 2004).
Публикации. Основные результаты диссертации опубликованы в шести работах, в том числе в журналах Дифференциальные уравнения и Труды Института математики и механики УрО РАН.
Структура и объем работы. Диссертация состоит из введения, трех глав, восьми параграфов (нумерация параграфов сквозная) и списка литературы. Объем диссертации 108 страниц. Библиографический список содержит 34 наименования.
Во введении описывается общая постановка задачи и излагается краткое содержание работы.
Основная часть диссертации посвящена изучению вопроса об асимптотическом подобии сужений системы (1) на конечномерные подпространства существенных решений, системам обыкновенных дифференциальных уравнений.
На протяжении этой работы предполагается, что интеграл Стилтьеса в (1) рассматривается по переменной s, xt (s) = x(t + s), функция (t, s) A(t, s) ограничена в полосе R[r, 0], равномерно непрерывна по t, имеет ограниченную вариацию по s, A(t, r) 0 и для дробно эти условия, которые далее будем называть естественными, описаны в первом параграфе. В дальнейшем систему (1) будем отождествлять с задающей ее функцией A, а пространство всех систем A, удовлетворяющих естественным условиям, обозначать A.
В качестве пространства начальных функций рассматривается пространство S всех непрерывных функций u : [r, 0] Rn либо = max |u(t)| (в последнем случае S становится банаховым пространством и совпадает с пространством C([r, 0], Rn )).
Для каждого момента времени t0 и любой непрерывной функции u S решением t x(t, t0, u) задачи Коши будем называть непрерывную на [t0 r, ) функцию, совпадающую с функцией u(t t0 ) при t [t0 r, t0 ] и обращающую систему (2) в тождество при t t0 (здесь xt (s) = x(t + s), s [r, 0]).
Всякое решение t x(t, t0, u) задачи (2), (3) порождает движение t xt (·, t0, u) = xt (t0, u) в пространстве S, t t0 (для краткости при t0 = 0 вместо xt (·, 0, u) пишем xt (u)).
В этой работе мы придерживаемся концепции, предложенной Н. Н. Красовским9), состоящей в том, что системе x(t) = v(xt ), обладающей свойством непрерывной зависимости и глобальной правосторонней единственности решения задачи Коши, отвечает некоторая динамическая система с фазовым пространством C([r, 0], Rn ) и потоком на нем t xt, порожденным решениями этой системы. Эта концепция оказалась естественной и очень плодотворной при изучении асимптотического поведения решений системы (1).
Во втором параграфе описана методика преобразования задачи к системе xk+1 (t) = C(t, kr)xk (kr)+ где k = 0, 1,..., функция x0 (t) совпадает при t [r, 0] с функцией u(t), а q(t, xk (·)) = 9) Красовский Н. Н. Некоторые задачи теории устойчивости движения (линейные системы). М.: Гос. издат. физ.-матем. лит. 1959. 211 с.
образом, при каждом k = 0, 1,... и всех t [kr, (k + 1)r] имеет место равенство x(t, u(·)) = x(k+1)r (s), где s [r, 0], t = s + (k + 1)r.
Показано, что эта методика не позволяет улавливать негрубые свойства системы (4), а также свойства устойчивости, равномерные относительно сдвигов системы. Однако, равенство (6) можно применять при получении оценок экспоненциального роста решений системы (4)10). Действительно, обозначим тогда найдется константа такая, что при всех t [kr, (k + 1)r] выполнено неравенство |q(t, xk (·))| xk. Пусть далее, Из (6) следует неравенство из которого, в свою очередь, получаем (при каждом k) неравенства и некоторых констант N и, то k N exp(r). Поэтому непосредственно из (7) для всех k = 1, 2,... получаем неравенства Таким образом, имеет место следующее утверждение:
для любого > 0 найдется такая константа N, что для всякой начальной функции u S и всех t 0 решение x(t, u) задачи (4), (5) допускает экспоненциальную оценку 10) Азбелев Н. В., Максимов В. П., Рахматуллина Л. Ф. Введение в теорию функционально–дифференциальных уравнений. М.: Наука, 1991. 280 с.
Поэтому, если < 0, то тривиальное решение системы (4) экспоненциально устойчиво.
В третьем параграфе исследуется вопрос о равномерной экспоненциальной устойчивости системы (1). Этот вопрос для систем, не обладающих свойством периодичности, не исследован.
О п р е д е л е н и е. Будем говорить, что система (1) удовлетворяющая естественным условиям, C-равномерно экспоненциально устойчива, если найдутся такие константы > 0 и M > 0, что для всякого движения t xt (·), порожденного системой (1), для любого t0 0 и всех t t0 выполнено неравенство Т е о р е м а 1. Система (1), удовлетворяющая естественным условиям, C -равномерно экспоненциально устойчива в том и только в том случае, если показатель Боля системы (1) удовлетворяет неравенству B0 (A) < 0.
Аналогичным образом вводится понятие L2 -равномерной экспоненциальной устойчивости.
Далее, пусть A0 подпространство всех систем из A, обладающих свойством продолжаемости влево (см. утверждение в) формулируемой ниже теоремы 4).
Т е о р е м а 2. Свойство C-равномерной экспоненциальной устойчивости на пространстве A0 с метрикой является грубым свойством.
Доказательство этой теоремы использует формулируемую ниже лемму, представляющую самостоятельный интерес.
Л е м м а 1. Показатель Боля B0 (A) системы A A устойчив вверх, то есть каждому > 0 отвечает такое > 0, что для любой системы A + B A0, где B удовлетворяет естественным место неравенство Если на пространстве A всех систем вида (1) определить метрику равенством то утверждения аналогичные выше сформулированным, останутся справедливыми без условия продолжаемости влево, то есть на всем пространстве A. Приведем эти утверждения.
Л е м м а 2. Показатель Боля B0 (A) системы A A устойчив вверх, то есть каждому > 0 отвечает такое > 0, что для любого возмущения B(t, s), удовлетворяющего естественным условиям и неравенству sup var |B(t, s)|, имеет место неравенt 0 s[r,0] ство (8).
Т е о р е м а 3. В пространстве A с метрикой, определенной равенством (9), свойство C-равномерной экспоненциальной устойчивости, является грубым свойством.
В четвертом параграфе введено понятие L2 -показателя Ляпунова решения системы (1) и сформулировано основное утверждение диссертации.
Пусть t0 = 0. Для xt (u) = xt (·, 0, u) определим L2 -показатель Ляпунова Тогда множество S = {u S : (u) = } образует линейное подпространство в S. Пусть S+ прямое дополнение подпространства S до пространства S, то есть S = S+ S. Тогда для всех ненулевых функций u S+ выполнено неравенство (u) >.
Пусть u1 (·),..., up (·) фиксированный набор p линейно независимых функций из S+. Линейное подпространство в S+, порожденное этим набором, обозначим Sp и всякой начальной функции u Sp поставим в соответствие движение t xt (u), отвечающее решению задачи (2), (3) при t0 = 0. Таким образом, построено движение t xt (Sp ) = Sp пространства Sp. Мы будем говорить, что это движение порождено сужением системы A на подпространство Sp. Такое сужение обозначим (A, Sp ). Наряду с системой (A, Sp ) рассмотрим систему обыкновенных дифференциальных уравнений с непрерывной на полуоси R+ матричной функцией t B(t). Будем далее отождествлять систему (10) с задающей ее матрицей B и называть системой B. По аналогии с подпространством Sp, введем t в рассмотрение линейное пространство Rp размерности p с базисом y 1 (t),..., y p (t), образующем столбцы матрицы Коши Y (t, ) системы B при = 0.
Пусть L(Sp, Rp ) пространство линейных операторов, действуюt t щих из Sp в Rp с нормой · L2 Rp.
Функцию t L(t) L(Sp, Rp ) будем называть обобщенным ляt t пуновским преобразованием систем (A, Sp ) и B, если: 1) функция t L(t) непрерывна на R+ = [0, ); 2) при t 0 оператор L(t) является гомеоморфизмом пространств Sp и Rp и 3) выполнено неравенt t ство sup L(t) L2 Rp + L1 (t) Rp L2 <. Будем говорить также, что система (A, Sp ) приводима обобщенным ляпуновским преобразованием L к системе B, или что системы (A, Sp ) и B асимптотически подобны.
В диссертации показано, что система (A, Sp ) имеет не более p различных L2 -показателей Ляпунова 1 (A),..., p (A) и что асимптотически подобные системы (A, Sp ) и B сохраняют показатели Ляпунова:
для всякого показателя i (A) системы (A, Sp ) найдется такое решеln |yi (t)| норма в R. Верно и обратное утверждение.
Основным утверждением диссертации является следующая теорема.
Т е о р е м а 4. Пусть (A, Sp ) сужение системы A на подпространство Sp и Sp S+. Тогда:
а) найдутся система B с непрерывной (p p)-матрицей B(t) и обобщенное ляпуновское преобразование L, приводящее систему (A, Sp ) к системе B;
б) в множестве {B} всех систем, асимптотически подобных системе (A, Sp ), найдется система с непрерывной на R+, верхней треугольной матрицей B(t);
в) если, в дополнение к сказанному, всякое решение системы (A, Sp ) продолжаемо влево, то есть найдется константа > такая, что для каждого u Sp, любого [r, 0] и всех t R+ выполнено неравенство xt+ (·, u) 2 xt (·, u) 2, то в множестве {B} всех систем, асимптотически подобных системе (A, Sp ), най- дется система B с ограниченной на полуоси R+ матрицей B(t) (и следовательно, с огрниченной на R+ верхней треугольной матрицей B(t));
г) если A(t + T, s) = A(t, s) для всех (t, s) R[r, 0], то найдутся система B с вещественнозначной непрерывной T-периодической матрицей B(t) и T-периодическое по t обобщенное ляпуновское преобразование L, приводящее систему (A, Sp ) к системе B.
д) в множестве обобщенных ляпуновских преобразований, приводящих систему (A, Sp ) к системе B с непрерывной на R+, верхней треугольной матрицей B(t), найдется ортогональное (L(t)L(t)=Ip ) обобщенное ляпуновское преобразование.
В диссертации описан алгоритм построения обобщенного ляпуновского преобразования, приводящего сужение (A, Sp ) системы A к системе B с треугольной матрицей B. Из этого алгоритма следует, что увеличение размерности треугольной системы обыкновенных дифференциальных уравнений на единицу за счет пополнения пространства началных условий не влечет за собой больших вычислительных затрат, так как при этом новая система содержит старую в качестве подсистемы.
Отметим, что без дополнительных условий не удается доказать ограниченность матрицы B(t) системы B, асимптотически подобной системе (A, Sp ) (весьма важное условие при исследовании асимптотического поведения решений системы B). Однако, матрица B(t) будет ограниченной, если выполнено, так называемое, условие продолжаемости влево (см. утверждение в) теоремы 4). Условие продолжаемости влево трудно проверяемо, но например, для периодических систем это условие и не требуется (см. утверждение г)).
Вероятно, этот факт имеет место и для более широкого класса систем, а именно для систем с рекуррентной по t матрицей A. В диссертации не доказан факт продолжаемости влево решений рекуррентных систем, но показано, что при некотором дополнительном условии, более слабом, чем условие продолжаемости влево, утверждение об ограниченности матрицы B(t) остается верным и для случая рекуррентных систем.
В пятом параграфе приводится доказательство теоремы 4.
В шестом параграфе выделен некоторый класс систем вида (1), имеющих бесконечномерное пространство решений, но подпространство решений S+ которых конечномерно и подсчитана размерность подпространства S+.
Рассмотрим систему Здесь x(t) Rn, y(t) Rm, функции A : R M(n), B : R M(n, m), D : R M(m) непрерывны и ограничены на прямой R, а функция C : R[r, 0] M(n, m) удовлетворяет естественным условиям.
Т е о р е м а 5. Для любой системы вида (11) размерность пространства S+ равна n + m.
Далее получен ряд следствий теоремы 5, касающихся асимптотического поведения решений системы Здесь x(t), y(t) R, скалярные функции (t), (t) и (t) непрерывны и ограничены на R, а функция g : R [r, 0] R удовлетворяет естественным условиям.
В седьмом параграфе введено определение рекуррентрой системы с последействием и исследуютя некоторые свойства рекурентных систем.
Функцию (t, s) A(t, s) (или, что эквивалентно, систему A A), удовлетворяющую естественным условиям, будем называть рекуррентной (по переменной t), если для любых > 0 и T > 0 множество относительно плотно на прямой R (то есть найдется такая константа В восьмом параграфе получен аналог теоремы В. М. Миллионщикова11) для линейных систем уравнений с последействием.
При каждом фиксированном s [r, 0] сдвиг функции t A(t, s) на константу обозначим A (t, s) = A(t +, s). Пусть далее, R(A) замыкание множества {A (t, s) : R} сдвигов функции A в локально открытой топологии. Это означает, что A R(A) в том и только в том случае, если для некоторой последовательности {i } i= и любых > 0 и T > 0 найдется такой номер i0, что для всех i i выполнено неравенство Зафиксируем подпространство Sp S+ и для каждой A R(A) полный набор L2 -показатель Ляпунова системы (A, Sp ) обозначим 1 (A),..., p (A). Будем считать, что 1 (A)... p (A).
Формулируемую ниже теорему можно рассматривать, как частичное распространение теоремы В. М. Миллионщикова на линейные системы уравнений с последействием.
Т е о р е м а 6. Пусть Sp S+, система A A рекуррентна и для всех A R(A) и некоторой константы > выполнено 11) Миллионщиков В. М. О связи между устойчивостью характеристических показателей и почти приводимостью систем с почти периодическими коэффициентами. // Дифференц. уравнения. 1967. Т. 3, № 12. С. 2127–2134.
неравенство 1 (A). Тогда найдутся система B с непрерывной и ограниченной на R верхней треугольной матрицей B(t) и обобщенное ляпуновское преобразование L, приводящее сужение (A, Sp ) системы A на подпространство Sp к системе B.
1. Быкова Т. С., Тонков Е. Л. О ляпуновской приводимости системы с последействием // Известия Института математики и информатики УдГУ. Ижевск, 2002. Вып. 2(25). С. 27–30.
2. Быкова Т. С., Тонков Е. Л. Ляпуновская приводимость линейной системы с последействием // Дифференциальные уравнения.
2003. Т. 39, № 6. С. 731–737.
3. Быкова Т. С. Ляпуновская приводимость системы с последействием // Вестник Тамбовского университета. Тамбов, 2003. Т. 8, вып. 3. С. 355–356.
4. Быкова Т. С., Тонков Е. Л. Распространение теоремы Перрона– Миллионщикова о триангуляции на линейные системы с последействием // Вестник Удмуртского университета. Серия Математика. Ижевск, 2004. № 1. С. 51–66.
5. Быкова Т. С. О ляпуновской приводимости систем с последействием // Современные методы теории краевых задач. Материалы Воронежской весенней математической школы Понтрягинские чтения XV. Воронеж, 2004. С. 41–42.
6. Быкова Т. С., Тонков Е. Л. Приводимость линейной системы с последействием // Труды Института математики и механики УрО РАН. Екатеринбург, 2005. Т. 11, № 1. С. 53–64.
Отпечатано с оригинал-макета заказчика Подписано в печать..2005. Формат 60x84/16.
Типография ГОУВПО Удмуртский государственный университет.
426034, Ижевск, ул. Университетская, 1, корп. 4.