Современный этап развития теории дифференциальных уравнений характеризуется как стремлением к анализу и переосмыслению огромного количества накопленного научно-практического материала, так и созданием новых, более общих и содержательных точек зрения, разработкой новых качественных и приближенных методов исследования дифференциальных уравнений, направленных на решение сложных задач теории и практики. Указанным направлениям исследования посвящена обширная литература (см., например, [4], [10], [11], [39], [44], [46], [53] и имеющуюся там библиографию).
Современные методы качественного анализа дифференциальных уравнений берут свое начало в работах А. Пуанкаре [45], А. М. Ляпунова [29], А. А. Андронова [1] и др. Эти методы условно можно разбить на две группы. Первая группа методов ориентирована на изучение нелокальной теории дифференциальных уравнений, когда решения изучают либо во всем фазовом пространстве, либо в существенной ее части. Вторая группа ориентирована на изучение локальной теории дифференциальных уравнений, когда решения изучают, например, в окрестностях особых точек, периодических решений, инвариантных многообразий и т. п.
Методы нелокальной теории развивались по ряду направлений. Упомянем лишь некоторые из большого числа методов, показавших свою эффективность при решении многих задач теории и практики. Были разработаны методы инвариантных многообразий, методы построения вполне непрерывных операторов, неподвижные точки которых определяют периодические решения систем дифференциальных уравнений, методы теории абсолютной устойчивости, с наибольшей полнотой разработанные в нелокальных проблемах теории управления, и многие другие (см., например, [20], [34], [36], [40], [44], [55]).
При изучении локальной теории дифференциальных уравнений возникают две принципиально различные ситуации.
Первая ситуация связана с тем, что линеаризованное уравнение в окрестности особой точки является невырожденным. Поведение решений таких уравнений хорошо изучено, здесь разработан ряд эффективных методов, таких как метод малого параметра, метод усреднения, метод аналитического продолжения (Н.Н.Боголюбов А. Н. Крылов, А.А.Митро-польский [6], [27], И. Г. Малкин [31], М. Розо [48], И. З. Штокало [56] и другие).
Вторая ситуация возникает в задачах исследования поведения решений дифференциальных уравнений, когда линеаризованное уравнение является вырожденным. При этом, как правило, дифференциальные уравнения содержат различные параметры. Возникающие здесь задачи приводят к необходимости исследования эволюции поведения системы в окрестностях особых точек в зависимости от значений параметров. Типичными здесь являются такие эффекты как ветвление решений и различные бифуркации (положений равновесия, периодических или почти периодических колебаний и т. п.). Существенный вклад в развитие теории бифуркаций и теории ветвления решений нелинейных уравнений внесли В. И. Арнольд [2], [3], Р. И. Богданов [5], М. М. Вайнберг, В. А. Треногин [8], Н. К. Гаврилов [7], Дж. Гукенхеймер, Ф. Холмс [9], Ж. Йосс, Д. Джосеф [61], Ю. С. Ильяшенко [13], Ю. А. Кузнецов [75], Дж. Марсден, М. Мак-Кракен [32], М. Т. Терехин [50], Л. П. Шильников [57], А. Н. Шошиайшвили [62] и.
ДР.
Важный раздел локальной теории дифференциальных уравнений составляют задачи о рождении малых циклов из положений равновесия автономных и неавтономных систем при изменении параметров системы. Этой проблематике, восходя ей к классическим работам А. Пуанкаре [45], А. А. Андронова [1] и Е. Хопфа [70] посвящена обширная литература (см., например, [2], [12], [24], [25], [32], [54], [58], [60], [65]-[68], [71], [72], [78]). Методы исследования бифуркаций Андронова-Хопфа в системах с гладкими нелинейностями основаны на использовании аналитической теории, теорем о центральном многообразии, нормальных форм. Негладкая ситуация впервые изучена М. А. Красносельским и В. С. Козякиным [17], которые использовали специальный метод функционализации параметров.
В задачах исследования поведения решений дифференциальных уравнений в окрестностях особых точек широкое распространение получил подход, основанный на применении методов функционального анализа, алгебры и геометрии. Этот подход, который можно называть операторным, показал свою эффективность в работах В. И. Арнольда [2], П.П.Заб-рейко, М. А. Красносельского [21], Э. М. Мухамадиева [38], Е.Н.Розенвас-сера [49], М. М. Вайнберга, В. А. Треногина [8] и многих других математиков. На основе разработанных методов удалось решить ряд важных для теории и практики задач, в частности, классифицировать основные типы локальных бифуркаций в нелинейных динамических системах, получить эффективные признаки различных ветвлений и бифуркаций, провести анализ устойчивости решений, предложить методы построения решений и др.
Многие вопросы современной теории дифференциальных уравнений и многочисленные приложения требуют дальнейшего развития операторных методов. Здесь особо актуальны следующие основные направления исследований. Первое связано с разработкой методов, приводят, не только к признакам ветвления или бифуркации решений, но и к возможности приближенного построения решений, получения асимптотических (по параметрам) формул, проведения анализа устойчивости решений и т. д.
Второе направление связано с разработкой методов, учитывающих специфику данного дифференциального уравнения для различных классов динамических систем, в частности, для задачи о возникновении вынужденных и свободных колебаний, для дифференциальных уравнений теории управления и др.
Третье направление относится к приложениям, в частности, к разработке алгоритмов и программ численного исследования поведения решений дифференциальных уравнений. Сложное поведение решений дифференциальных уравнений в окрестностях особых точек часто выдвигает на первый план именно компьютерное моделирование системы. В этой связи особый интерес вызывает разработка операторных методов, которые могут быть доведены до алгоритмов и программ численного исследования системы.
Основной целью диссертации является:
— Разработка новых общих операторных методов исследования поведения решений нелинейных дифференциальных уравнений в окрестностях особых точек, приводящих к качественным и количественным характеристикам решений.
— Конструирование семейств операторных уравнений, определяющих основные сценарии поведения решений широкого класса дифференциальных уравнений в окрестностях вырожденных особых точек.
— Разработка итерационных процедур приближенного исследования операторных уравнений. Получение новых признаков бифуркаций периодических и почти периодических решений нелинейных дифференциальных уравнений, основанных на анализе эффективно вычисляемых спектральных характеристик особых точек.
— Исследование сценариев бифуркационного поведения систем с медленно меняющимися и слабоосциллирующими параметрами в окрестностях особых точек.
— Разработка пакета программ компьютерного моделирования поведения решений дифференциальных уравнений, основанных на итерационных процедурах численного построения решений эквивалентных операторных уравнений.
Диссертация в значительной степени носит теоретический характер. В ней на основе функционально-операторного подхода изучаются вопросы локальной теории нелинейных дифференциальных уравнений, зависящих от параметров, разработаны методы исследования периодических колебаний в окрестностях особых точек, их приближенное построение, анализ устойчивости, исследованы основные сценарии бифуркационного поведения динамических систем, предложена и обоснована новая итерационная процедура численного исследования бифуркации автоколебаний.
Полученные результаты доведены до расчетных формул, составлены и отлажены соответствующие программы. Предложенная итерационная процедура позволяет эффективно строить бифурцирующие решения, а также позволяет в новых условиях обнаруживать возникновение периодических колебаний при изменении параметров. Полученные результаты важны в задачах локальной теории дифференциальных уравнений, в задачах приближенного построения периодических колебаний нелинейных динамических систем теории управления. Предложенные методы итерационного построения решений могут быть использованы при составлении алгоритмов и программ численного исследования колебаний.
Результаты диссертации докладывались на Международной конференции «Дифференциальные уравнения и смежные вопросы», посвященной памяти И. Г. Петровского (Москва, ИПУ, 2007 г.) — на Международной конференции «Дифференциальные уравнения, теория функций и приложения», посвященной 100-летию И. Н. Векуа (Новосибирск, 2007 г.) — на Международной конференции «Дифференциальные уравнения и их приложения» (Душанбе, 2002 г.) — на Международной конференции «7 th International Pure Mathematics conference» (Пакистан, Исламабад, 2006 г.) — на Международной конференции «International Congress on Ghiyath Al-Din Jamshid Kashani (ICGK, 2000 r.)» (Kashan, I.R. Iran) — на Второй Международной конференции по проблемам управления, (Москва, ИПУ РАН 1999 г.) — на Второй и Третьей Всероссийских научных конференциях «Проектирование инженерных и научных приложений в среде MAT-LAB «(Москва, ИПУ РАН, 2004 г. и Санкт-Петербург, СПбГУ, 2007 г.) — на Второй Всероссийской научной конференции «Математическое моделирование и краевые задачи» (Самара, 2005 г.) — на семинаре ICTP, Триест (Италия, 2005 г., руководитель — профессор Д. Ли) — на семинаре в Институте проблем управления Российской академии наук (Москва, 2002 г., руководитель — профессор Н.А.Бобылев) — на семинарах Сибайского института Башгосуниверситета (2004;2007 гг., руководитель — профессор М.Г.Юмагулов) — на семинарах в Институте математики АН Республики Таджикистан (2000;2007 гг.).
По теме диссертации опубликовано более 20 научных статей. Список основных публикаций приведен в конце диссертации.
Работа состоит из введения и четырех глав. Главы разбиты на параграфы.
Список литературы
содержит 103 наименований.
СОДЕРЖАНИЕ РАБОТЫ.
Первая глава (§§ 1.1−1.4) носит вспомогательный характер. В § 1.1 рассматриваются задачи о периодических решениях для различных классов дифференциальных уравнений и указаны способы перехода к эквивалентным интегральным уравнениям. В § 1.2 приводятся вспомогательные сведения из качественной теории дифференциальных уравнений: вопросы структурной устойчивости, краткие сведения из теории бифуркаций. В § 1.3 приводятся основные положения метода функционализации параметра, предложенного М. А. Красносельским [21] для решения задач со связными континуумами неподвижных точек. В § 1.4 приводятся основные положения метода Ныотона-Капторовича [14] в той форме, в которой он используется в работе.
Вторая глава (§§ 2.1−2.4) содержит основные результаты диссертации, связанные с разработкой качественных и приближенных методов исследования свободных колебаний дифференциальных уравнений в окрестностях вырожденных особых точек.
В § 2.1 рассматривается автономная система х' = А (Л)ж + а (ж, А), ж е (1) в которой N ^ 2 и а (х, Л) — нелинейная вектор-функция равномерно по Л удовлетворяющая соотношению ||а (ж, Л)|| = °(||ж||) при ||ж|| —> 0. Предполагается, что матрица А (Х) при Л = Ао имеет одну простую пару чисто мнимых собственных значений то > 0. В этом случае в системе (1) в окрестности решения х = 0 могут возникать устойчивые или неустойчивые предельные циклы — это явление обычно связывают со следующим понятием.
Число Ао называют точкой бифуркации Андронова-Хопфа системы (1), если найдется последовательность Ап такая, что при каждом, А = Ап система (1) имеет ненулевое периодическое решение х = хп{Ь) некоторого периода Тп, при этом Ап —Ао и гпах хп{1) —> 0 при п —> оо.
В § 2.1 приводятся основные положения операторного метода исследования задачи о бифуркации Андронова-Хопфа системы (1), приводящего к новым признакам бифуркации, итерационной процедуре построения бифурцирующих решений и анализу их устойчивости.
В § 2.2 и § 2.3 изучаются задачи о свободных колебаниях дифференциальных уравнений теории управления.
В начале § 2.2 приводятся вспомогательные сведения из теории автоматического управления, такие как линейное звено, интегральное звено первого порядка, звено с дробно-рационалыюй передаточной функцией, импульсно-частотная характеристика и т. п. Рассматриваются также нелинейные дифференциальные уравнения теории автоматического управления.
Рассматривается нелинейное дифференциальное уравнение.
2) где.
Цр, А) = рп + а1(АУ1−1 + а2(Х)рп~2 +. + ап (А) М (р, А) = Ь0(Х)рт + б1(А)р—1 +. + Ът{А).
— многочлены с нерерывно зависящими от скалярного параметра Л вещественными коэффициентами, п > т > 0- функция /(ж, Л) предполагается непрерывной по совокупности переменных, причем представимой в виде fix, А) = с (Х)х + ip (x, X), где функции с (А) и (р (х, Х) непрерывны по своим переменным и равномерно по Л выполнено условие ¡-<£>(а-, А)(= °(|ж|), —^ 0. Пусть при некотором Л — Ло выполнены условия: Uli Уравнение.
L (p, А0) — с (А0)М (р, А0) = 0 имеет два чисто мнимых корня р — ic^o", с^о > 0, и не имеет корней вида где к = 0, 2,3, • • ¦- U12 Имеют место соотношения.
L (±uu0ki, А0) ф 0, к = 0,±1,±2, .
В этих условиях значение Ао будет бифуркационным для уравнения (2), т. е. в ней могут возникать периодические колебания с периодом близпп 2тг ким к То = —.
LU о.
В § 2.2 приводятся основные положения операторного метода исследования задачи о бифуркации Андронова-Хопфа уравнения (2), приводящего к новым признакам бифуркации, итерационной процедуре построения бифурцирующих решений и анализу их устойчивости. В § 2.2 проводится также анализ устойчивости бифурцирующих решений уравнения (2).
В § 2.3 приводится метод импульсно-частотных характеристик (ИЧХ) решения задачи о бифуркации для дифференциальных уравнений вида (2).
В § 2.4 приводятся доказательства основных утверждений второй главы.
Третья глава (§§ 3.1−3.3) посвящена рассмотрению двумерных динамических систем. Для таких систем изучаются условия возникновения периодических решений при потере устойчивости стационарных состояний (задача о бифуркации Андронова-Хопфа). Приводится новое доказательство теоремы о признаках бифуркации, основанное на методах теории вращения векторных полей.
В § 3.1 изучаются дифференциальные уравнения х" + а (Х)х' + Ъ{)х = 1р (х, А), (3) в котором коэффициенты а (А) и Ь (А) и функция <�р (х,) зависят от скалярного параметра Л. Предполагаются, что:
— равномерно по параметру Л выполнено соотношение р (х, Х) = о (|а?|), х 0 ;
— при некотором Л = Ао выполнены соотношения: а (А0) = 0, а'(Ло) ф О, Ь (Х0)>0, при этом к = 2,3,.
В этих условиях значение Ао будет бифуркационным для уравнения (3). Явление бифуркации здесь изучается методами теории вращения векторных полей.
В § 3.2 рассматривается уравнения Льенара вида х" + ai (A)z' + а2(Х)х = f (x, ж — А), (4) правая часть которого представима в виде 1 f{x, у, X) = Ь1(Х)ф{х-Х)у + b2{X) J ф{рX)dp, о при некоторой непрерывной функции ф (х, А) такой, что ф{О, А) = 0.
Для уравнения (4) получены необходимые и достаточные условия бифуркации периодических решений, при этом для построения решений разработана итерационная процедура.
В § 3.3 приводятся иллюстративные примеры исследования двумерных динамических систем.
В четвертой главе (§§ 4.1−4.5) основным объектом исследования снова является дифференциальное уравнение (1). Исследуется ситуация, когда параметр, А слабо осциллирует по периодическому закону в окрестности точки бифуркации Ао:
А = fs{t) = А0 + 6.
< 1. В этом случае уравнение (1) принимает вид: x' = A[fs (t)]x + a[xJ6(t)], (5) в котором 5 является параметром.
В § 4.1 предполагается, что число 0 является простым собственным значением матрицы Л (Ао). Показано, что в этом случае бифуркация двукратного равновесия системы (1) преобразуется в бифуркацию вынужденных колебаний системы (5) с сохранением основных характеристик бифуркации.
В § 4.2 система (1) исследуется в предположении, что значение До является точкой бифуркации Андронова-Хопфа. Другими словами, здесь изучается вопрос о том, во что может преобразоваться бифуркация Андро-нова-Хопфа в ситуации, когда се параметры слабо осциллируют в окрестности точки бифуркации. Показано, что явление бифуркации Андронова-Хопфа является в естественном смысле устойчивым по отношению к малым периодическим возмущениям параметров системы: она преобразуется в бифуркацию почти периодических колебаний.
В § 4.3 приводятся примеры численного исследования дифференциальных уравнений, а в § 4.4 приводятся доказательства основных утверждений четвертой главы.
Заключительный § 4.5 носит характер приложения. В нем приводится описание алгоритмов и программ численного исследования бифуркации, разработанных в соответствии с предложенными в диссертации методами исследования.
4.2.1 Основные результаты.
В силу приведенных выше предположений матрица ^4(А) при А, близких к Ао, имеет пару простых собственных значений вида /¿-(А) = о-(А)±го-(А), где а (А) и а>(А) — непрерывно дифференцируемые функции такие, что а (Ао) = 0 и а-(До) = сооТак как собственные значения матрицы.
Ао = Л (Ао) и транспонированной к ней матрицы Ао простые, то найдутся такие пары линейно независимых векторов е, д? Яп и е*, д*? что.
А0е = —щд, А0д = ш0е, Л? е* = си0д*, А^* = -и>0е*;
4.11).
М = |Ы| = 1, (е, е*) = {д, д*) = 1, (е, д*) = (д, е*) = 0. Пусть, наряду с 111, выполнено также условие.
112. 7 ^ 0, где 7 = (Л'е, е*) + (А'д, д*) и А' = А'{А0).
Тогда [25] значение Ао параметра, А является точкой бифуркации Андро-нова-Хопфа для системы (4.7). Возникающие при, А = А (е) периодические решения системы (4.7) имеют период Т = Т (е), близкий к.
2тг числу То = —, а именно: |Т{е) — То| < е. Отметим также [26], что число.
Шо.
7 связано с числом а'(Хо) равенством 7 = 2<�У (Ао). Положим = (Р (т)(1т.
Теорема 4.7. Пусть 70 < 0 ('усро > 0) — тогда х = 0 является асимптотически устойчивым (неустойчивым) решением уравнения (4-Ю) при малых 5 > 0.
Доказательство этой теоремы приводится в п. 4.4.
Отметим [25], что в задаче о бифуркации Андронова-Хопфа для системы (4.7) решение х = 0 становится неустойчивым при 7(А — Ао) > 0- т. е. если, например, 7 > 0, то решение х = 0 неустойчиво при, А > Ао.
Условия теоремы 4.7 гарантируют устойчивость системы (4.8) в окрестности точки бифуркации Ло, если даже ее параметр Л периодически попадает в область неустойчивости, лишь бы выполнялось соотношение 7</?о < 0. Например, если 7 > 0, то при (ро < 0 решение х = 0 является устойчивым, хотя параметр Л при этом может периодически удовлетворять неравенству Л > Ло, т. е. переходить границы устойчивости исходной автономной системы (4.7). Отмеченный факт в естественном смысле соответствует эффекту затягивания потери устойчивости стационарных состояний динамических систем [41].
Пусть выполнены условия теоремы 4.7. Тогда 7 ф 0 и, следовательно, для системы (4.7) значение Л = Ло является точкой бифуркации Андро-нова-Хопфа, т. е. при переходе через значение Ло у системы (4.7) возникают малые периодические ненулевые решения. Возникает естественный вопрос о том, сохранится ли подобное поведение для возмущенной системы (4.10), в которой параметром является 5- значение 5 = 0 для нее, очевидно, является критическим.
Число 5 = 0 называют [19] точкой бифуркации почти периодических колебаний системы (4.10), если каждому е > 0 соответствует такое 5 =.
5(e) 6 [0,е), при котором система (4.10) имеет нестационарное почти периодическое решение x (t, e), причем sup ||ж (£, е)|| < е. t.
Теорема 4.8. Пусть усро ф 0. Тогда значение 5 = 0 является точкой бифуркации почти периодических колебаний системы (4−10).
Теорема 4.8 допускает усиление в следующем направлении.
Теорема 4.9. В условиях теоремы 4−8 каждому е > 0 соответствует такое 5 = S (s) G [0,s) — при котором система (4-Ю) имеет нестационарное почти периодическое решение x (t:e), представимое в виде x (t, e) = х0 (i, e) + h (t,?), где xq (t, е) — Т-периодическая функция некоторого периода Т = Т (е).
27 Г такого, что Т (г) —> То = —, a h (t, е) — это почти периодическая.
CdQ функция, удовлетворяются соотношению sup ||Л (?, е)|| = o (max \xo{t, ?" ||) t 1 при е —> 0. Функция хо (t,£) — это бифурцируюш}ее решение задачи о бифуркации Андронова-Хопфа для системы х' = [А (Л0) + 6(s).
Доказательства этих теорем приводятся в п. 4.4. Таким образом, теорема 4.8 показывает, что бифуркация Андронова-Хопфа исходной автономной системы (4.7) преобразуется в бифуркацию почти периодических колебаний системы (4.10). При этом согласно теореме 4.9 основной вклад в формирование возникающих почти периодических колебаний x (t, e) системы (4.10) вносят бифурцирующие периодические решения хо (t, e) автономной системы (4.7) с периодом Т (е), близким к То, а слабая осцилляция параметра, А приводит лишь к малому дрожанию этих решений.
4.3 Примеры численного исследования.
Приведем в качестве иллюстрации некоторые результаты компьютерного моделирования локальных бифуркаций со слабой осцилляцией параметров.
4.3.1 Пример 1: уравнение Ван-дер-Поля.
Рассмотрим уравнение Ван-дер-Поля вида (см. также п. 3.3.1).
При Л < 0 решение х = 0 уравнения (4.13) является устойчивым, а при Л > 0 оно становится неустойчивым и в результате у уравнения (4.13) возникают малые ненулевые периодические решения с периодом близким.
Пусть параметр Л уравнения (4.13) слабо осциллирует по закону Л = где + р) = > (?), а именно рассмотрим уравнение х" + (Зх2 — Х) х' + ж = 0.
4.13) к То = 2тт. х" + [Зх2 — 0,1ф)}х +х = 0.
4.14).
Рисунок 4(а, б, в) а) б) в).
На рисунке 4 приведены результаты компьютерного моделирования уравнения (4.14) для различных функций </?(?). Во всех рассмотренных случаях начальные условия решений выбраны одинаково: ж (0) = 0,1 и ж'(0) = 0. В случае (a).
27г-периодической функции, что соответствует теореме 4.9.
Для сравнения на рисунке 4 (в) приведена ситуация, когда условия теореми 4.8 и 4.9 не выполнены, а именно ip (t) = cos t. В этом случае, как показывает рисунок в окрестности нуля также возникают устойчивые почти периодические колебания, однако их форма далека от от периодической.
4.3.2 Пример 2.
Рассмотрим теперь дифференциальное уравнение х' = х-х2 (4.15).
Стационарными решениями уравнения (4.15) являются функции х — 0 и х = А, т. е. А = О является точкой бифуркации двукратного равновесия. Наряду с (4.15) рассмотрим уравнение вида х' = гп cos2 (t)x — х2, (4.16) полученную из (4.15) путем слабой осцилляции параметра Л. Несложно проверить условия теорем 4.2, 4.3 и 4.5:
А0 = 0, <р0 = тгфО, А (А) = А, ¿-(Ао) = 0, А'(А0) = 1, до = е0 = 1, к0тт0, ?2 = ~2тг ф 0, ке2 = -2Л < 0.
Следовательно при переходе через, А = 0 происходит бифуркация вынужденных устойчивых колебаний уравнения (4.16). Этот факт проиллюстрирован на рисунке 5.
Здесь: верхний график — это решение уравнения (4.15), а нижнийэто решение уравнения (4.16), при этом началные условия одинаковы.
4.4 Доказательства основных утверждений. 4.4.1 Вспомогательные утверждения.
При доказательстве основных утверждений настоящей главы используются некоторые установленные в [59] результаты, связанные с бифуркацией малых решений операторных уравнений. Приведем их в краткой форме.
Рассмотрим операторное уравнение вида х = В (Х)х + а (х, Х), (4.17) где В (Л): Е —>• Е — вполне непрерывный оператора (ж, А) удовлетворяет соотношению || а (ж, А) ||= о (|| х ||), || х ||-> 0.
Число Ао назовают точкой бифуркации малых решений уравнения (4.17) если существуют последовательности Ап —у А0 и хп —у 0, хп ф 0,.
Пусть выполнены следующие условия VI и У2.
VI. Число 1 является простым собственным значением оператора В (Ао).
Обозначим через ео — собственный вектор матрицы Во, отвечающий нулевому собственному значению. Транспонированная матрица Вд также имеет простые нулевые собственные значения. Пусть до — соответствующий собственный вектор матрицы Будем считать, что ео и до выбраны из условий ||ео|| = 1 и (ео, до) = 1.
2. Имеет место соотношение (В'(Хо)ео, до) ф 0.
Теорема 4.10. Пусть выполнены условия VI и V2. Тогда Ао является точкой бифуркации малых решений уравнения (4−17).
Пусть нелинейность а (х, А) имеет вид, а (¡-с, А) = а^ж. А) + оз (ж, А) + Ь (х, А).
Теорема 4.11. Пусть выполнены условия VI и V2. Пусть (^(ео, Ао), до) Ф 0. Тогда бифуркация в уравнении (4−17) является двусторонней. При этом уравнение (4−17) имеет в точности одно семейство бифурциру-ющих решений х (Х).
Теорема 4.12. Пусть выполнены условия VI и V2. Пусть [а2{ео, Ао), до) О и (аз (ео, Ао),.
4.4.2 Доказательство теоремы 4.2.
Задача о бифуркации вынужденных колебаний уравнения (4.3) равносильна задаче о бифуркации малых решений операторного уравнения х = У (Т,<5)х + и (х, 5), (4.18) где У (Т, — матрица монодромии линейного уравнения х' = а оператор и (х, 5) имеет вид т у (х, 5) = J Х (Т, 5) Х~1(г, 6) а[Х (г, 5) х, 6](1т, о где Х (1−7 6) — фундаментальная матрица решений линейного уравнения ММФ.
Матрица монодромии линейной части системы (4.3) имеет вид: А[5ШсН.
3) = е° - (4.19) тогда К (Т, 0) = еТА ()^ Пусть ц (5) — простое непрерывно зависящее от 6 собственное значение оператора У (Т, 5) такое, что ?/(0) = 1.
В соответствии с теоремой 4.10 для доказательства теоремы 4.2 достаточно установить, что 1? <�т (У (0)), и /?'(0) ф 0. Выполнение первого из этих требований вытекает из условия II1. Далее, из [26] получим ц (6) = 1 + <�ро (Х (о)ео, до) 5 + о (5) 175.
Тогда.
0) = ?>о (А'(А0)е0, уо) Ф 0.
Теорема 4.2 доказана.
4.4.3 Доказательства теорем 4.3 и 4.4.
Справедливость утверждения теоремы относительно бифуркации двукратного равновесия системы (4.1) установлена в [59].
Докажем утверждение о бифуркации вынужденных колебаний системы (4.3). Эта задача приводит к уравнению (4.18), в котором У (Т, ё) оператор (4.19), а оператор у (х, ё) представим в виде т у (х, 5) = J У{Т, 5) У-т, ё) а (У (т, 8) х, 5)<1т. о.
В соответствии с теоремой 4.11 для доказательства теоремы 4.3 достаточно установить, что (у2(ео, 0), до) Ф 0. Действительно, т е2 = Ыео, 0),$<,) = I {У{Т, 0) У-т, 0) а2(У (т, 0) е0,0), д0) ¿-т = о т I (е-тА^а2(е0,о), до) ¿-т. о.
Так как 1 — собственное значение оператора е~тЛ (х° то верно равенство е~тА (Хо^ео = ео. Пространство Яп может быть представлено в виде Яп — Щ ф Н°, где Щ одномерное подпространство содержащее вектор ео, а — это собственное подпространство отвечающее остальному спектору оператора А (Ао). Ниже используется тот факт, что из и € Н°, следует (и, д0) = 0. Далее <12(60, А0) = Р0а2(е0,о) + Р°а2(ео, А0), где Р0х = (х, до) ео и Р° = I — Pq. Поэтому е-г^а2(ео, А0), ро) = (е-тА (Ло)(^оа2(е0, А0) + Р%(е0) А0))^0) = (е-тЛ^Р0а2(ео, Хо), до) = (е'^^Ыео,), д0) е0, до) = (02(^0, Ао), .
4.4.4 Доказательства теорем 4.5 и 4.6.
Характер устойчивости решения х (е) системы (4.1) определяется собственными значениями матрицы.
B (e) = A[(e)] + a'2x[x©, 4t)}.
Так как <22(2, А) является квадратичной, то матрица а2х (х, Х) линейно зависит от х. Поэтому а'2х (х (е), Х (е)) = а'2х (?ео+е2е1 + о (е2), A0+eAi + o (s)) = ?a'2x (eQ, А0) + о (е). Тогда.
В (е) = Л (А0) + eAiA’fAo) + ?а'2х{е0, А0) + о (е).
Согласно [26] собственное значение fi (e) матрицы В (е:), близкое к нулю, представляется в виде li (e) = (AiA'(Ao)eo + «2а-(ео, А0) е0, д0) е + о (е). (4.20).
Так как а'2х (х, Х) х = 2а2(х, А), то равенство (4.20) примет вид = (А1Л/(Ао)ео + 2а2(ео, Ао))ро)? + о (£).
Отсюда получим х (е) = (а2(е0, Ло)> Ро)^ + о (е), из которого и следует утверждение первой части теоремы.
Для доказательства второй части заметим, что ?1 = T? q. Согласно [20] бифурцирующие решения системы (4.3) устойчивы, если k5? i < 0. Т.к. к > 0, 6 > 0, и < 0, то получим требуемый результат. Теорема 4.5 доказана.
Доказательство теоремы 4.6 проводится по той же схеме, что и доказательство теоремы 4.5.
4.4.5 Доказательство теоремы 4.7. Вспомогательные утверждения.
Приведем сначала вспомогательные утверждения. Рассмотрим линейную систему где A2(t, 5) — это Т — периодическая матрица, ограниченная по t и S. Согласно[48] существует замена переменных у = Кх, где К = K (t) -некоторая невырожденная ограниченная Т-периодическая матрица, переводящая систему (4.21) к виду х' = [А (А0) + 6.
4.21) dy dt [Л (А0) + 6S + S2R (t, 6)]y,.
4.22) где ?) — Х-периодическая ограниченная матрица, а постоянная матрица 5 удовлетворяет равенству т т.
I Х-1ЗХ (г)(т)с1т = I Х~ф[г)А'(о)Х{г)(1г] (4.23) о о здесь Х{Ь) = ел° и А0 = Л (А0).
Из результатов работы [26] следует.
Лемма 4.1. Пусть матрица Л (Ао) имеет простые собственные значения ±а-ог. Тогда матрица ^4(Ао) + 63 имеет пару простых собственных значений рь{5) = ?(6) +гт{6), где ?(5), т{5) — непрерывно дифференцируемые функции причем ?(0) = 0, т (0) — шо и.
Лемма 4.2. Имеет место равенство.
5е, е*) + (5д,/) = ^ (4.24).
Завершение доказательства теоремы 4.7.
Теорема 4.7 будет доказана, если показать, что нулевое решение линейного уравнения х' = А[А0 + 6<�р{Ь)]х (4.25) является асимптотически устойчивым (неустойчивым) при малых 5 > О, если матрица В (5) = Л (Ао) + 55 при таких 5 будет устойчивой (неустойчивой).
В силу леммы 4.1 устойчивость матрицы В (5) определяется знаком числа ?'(0): если ?'(0) < 0, то матрица В (6) устойчива, если же ?'(0) > 0, то матрица В (6) неустойчива. Так как по лемме 4.2 имеем ?'(0) = то отсюда получим утверждение теоремы.
Уо7.
2 Т '.
4.4.6 Доказательства теорем 4.8 и 4.9.
При доказательстве этих утверждений без специальных ссылок используются некоторые понятия и факты из теории почти периодических колебаний [19].
Для системы (4.12) значение До является точкой бифуркации Андро-нова-Хопфапоэтому каждому е > 0 соответствует такое 5 = 5(е) Е [0, бг), при котором система (4.12) имеет нестационарное периодическое решение е), причем тах ||жо (£, е)11 < еДля функций 5(е) и хо^, е) могут быть получены [26] асимптотические формулы по вспомогательному параметру е, при этом 5(е) = 0{е2) и хо^, е) = 0{е) при е —> 0. Ниже без ограничения общности будем считать, что 6(е) > 0.
Произведем в системе (4.10) (в которой положим 5 = £(е)) указанную выше замену переменных по формуле у — Кх, где К = К (Ь, 8) в результате перейдем к системе у' = [А (А0) + + + Ка[К~1у, Д0 + 5ф)). (4.26).
Из леммы 4.1 следует, что дифференциальный оператор Ь (8)у = у' — [А (До) + 8Б]у при малых 5 > 0 является регулярным. Поэтому задача о почти периодических колебаниях системы (4.26) равносильна нелинейному интегральному уравнению.
3/(0 = I оо где г, — функция Грина регулярного оператора Ь (5) и у, *, = 5) у + Ка[К-1у, Л0 + 5ф)].
Для завершения доказательства на основе принципа сжатых отображений можно показать, что последнее уравнение в шаре Т{хо, р) бана-хового пространства В (ЯМ) почти периодических функций с центром в точке радиуса р{е) — о (е:) имеет единственное почти периодическое решение.
Заключение
.
В диссертационной работе получены следующие основные результаты:
1. Разработаны теоретические положения о новом операторном методе исследования решений нелинейных дифференциальных уравнений в окрестностях вырожденных особых точек.
2. Предложены новые семейства операторных уравнений в задачах о бифуркации периодических и почти периодических колебаний нелинейных дифференциальных уравнений, решения которых изолированы в соответствующих пространствах и могут быть построены итерационными процедурами.
3. Получены новые признаки бифуркации периодических и почти периодических решений нелинейных систем, основанные на анализе эффективно вычисляемых спектральных характеристик особых точек.
4. Исследованы общие свойства функций Грина в задаче о периодических решениях нелинейных дифференциальных уравнений теории управления, позволившие получить новый метод построения периодических решений дифференциальных уравнений.
5. Предложена новая процедура исследования основных сценариев бифуркационного поведения систем с медленно меняющимися и слабоосцилирующими параметрами.
6. Разработаны алгоритмы и программы численного исследования бифуркации периодических и почти периодических колебаний нелинейных дифференциальных уравнений.