Авторефераты по всем темам  >>  Авторефераты по разным специальностям

На правах рукописи

Муртазина Сария Аширафовна

МЕТОДЫ ИССЛЕДОВАНИЯ ЛОКАЛЬНЫХ БИФУРКАЦИЙ КОРАЗМЕРНОСТИ ДВА В НЕАВТОНОМНЫХ И ДИСКРЕТНЫХ ДИНАМИЧЕСКИХ СИСТЕМАХ

01.01.02 дифференциальные уравнения, динамические системы и оптимальное управление

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата физико-математических наук

Уфа - 2012

Работа выполнена на кафедре дифференциальных уравнений Башкирского государственного университета

Научный консультант: доктор физико-математических наук, профессор Юмагулов Марат Гаязович

Официальные оппоненты: доктор физико-математических наук, профессор Ильясов Явдат Шавкатович, доктор физико-математических наук, профессор Красносельский Александр Маркович

Ведущая организация: Институт математики и механики Уральского отделения РАН г. Екатеринбург

Защита состоится 6 апреля 2012 г. в 15.00 часов на заседании диссертационного совета Д 002.057.01 в Учреждении российской академии наук Институт математики с вычислительным центром Уфимского научного центра РАН по адресу: 450008, г. Уфа, ул. Чернышевского, 112.

С диссертацией можно ознакомиться в библиотеке Учреждения российской академии наук Институт математики с вычислительным центром Уфимского научного центра РАН.

Автореферат разослан 1 марта 2012 года.

Ученый секретарь диссертационного совета кандидат физ.-мат. наук С. В. Попенов

Общая характеристика работы

Актуальность темы. В общей теории динамических систем важное место занимают системы, описываемые дифференциальными уравнениями c периодическими коэффициентами. Такие уравнения возникают при изучении колебательных процессов в механике, физике, биологии, химии, экономике и др. Теория дифференциальных уравнений с периодическими коэффициентами являются хорошо развитой составной частью общей теории дифференциальных уравнений, благодаря работам А.М. Ляпунова, А. Пуанкаре, Н.Н.

Боголюбова, В. Г. Веретенникова, И.Г. Малкина, Ю.А. Митропольского, В.А.

Плисса, М. Розо, Ж. Флоке, Л. Чезари, И. З. Штокало, В.А. Якубовича и многих других математиков.

Одной из наиболее актуальных задач в теории дифференциальных уравнений с периодическими коэффициентами представляется исследование поведения системы в окрестностях стационарных и периодических решений. Здесь особо актуальными представляются исследования поведения системы в предположении, что стационарное или периодическое решение является негиперболическим. В этом случае в системе могут происходить различные бифуркационные явления, возникать новые периодические или квазипериодические решения, становится возможным хаотическое поведение системы. Исследованию этого случая посвяшены работы В.И. Арнольда, М.А. Красносельского, А.М. Красносельского, А.П. Кузнецова, С.П. Кузнецова, В.С. Козякина, Н.А.

Магницкого, Ж.К. Хейла, Л.П. Шильникова и других математиков. Здесь получен ряд важных результатов, связанных с признаками бифуркаций и субфуркаций, построением периодических решений, анализом их устойчивости, исследованием вопросов синхронизации и др. При этом большая часть полученных результатов относится к задачам о бифуркациях коразмерности один. Существенно меньше изучены задачи о локальных бифуркациях коразмерности два и выше. Здесь особо актуальными являются получение достаточных признаков различных сценариев бифуркаций, разработка методов приближенного построения возникающих колебаний, исследование их устойчивости.

Цель работы. Разработка качественных и приближенных методов анализа бифуркационных явлений коразмерности два в динамических системах, описываемых дифференциальными уравнениями c периодическими коэффициентами, получение достаточных признаков основных сценариев локальных бифуркаций, исследование устойчивости возникающих колебаний.

Методы исследования. В работе использованы методы качественной теории дифференциальных уравнений, теории нелинейных колебаний, нелинейного анализа, приближенного решения операторных уравнений, теории Флоке, малого параметра, метод функционализации параметра, метод Ньютона-Канторовича, метод Розо исследования устойчивости.

Научная новизна определяется впервые проведенными исследованиями, в результате которых разработан математический аппарат для анализа локальных бифуркаций коразмерности два в неавтономных и дискретных динамических системах. При этом получены следующие новые научные результаты:

1. Разработан новый операторный метод исследования бифуркационного поведения основных сценариев бифуркационного поведения двупараметрических неавтономных динамических систем и их дискретных аналогов в окрестностях стационарных решений;

2. Получены новые достаточные признаки локальных бифуркаций коразмерности два неавтономных динамических систем и их дискретных аналогов;

3. Разработаны и обоснованы асимптотические формулы, позволяющие определить основные гармоники вынужденных и субгармонических колебаний неавтономных динамических систем в задачах о бифуркации коразмерности два;

4. Проведен анализ устойчивости вынужденных и субгармонических колебаний, возникающих в неавтономных динамических системах при бифуркиях коразмерности два;

5. Предложены асимптотические формулы в задаче о локализации языков Арнольда неавтономных динамических систем и их дискретных аналогов в основных резонансах.

Теоретическая и практическая значимость. Работа носит теоретический характер. В ней предложены и обоснованы аналитический и приближенный методы исследования задачи об основных сценариях локальных бифуркаций коразмерности два в дифференциальных уравнениях с периодическими коэффициентами. Предлагаемые методы могут быть использованы для анализа бифуркационных явлений в системах, описываемых неавтономными дифференциальными уравнениями с периодическими коэффициентами.

Полученные результаты доведены до расчетных и асимптотических формул.

Апробация работы. Основные результаты диссертации докладывались на региональной научно-технической конференции "Новые программные средства для предприятий Урала"(г. Магнитогорск, декабрь 2004 г.); всероссийской школе-конференции для студентов, аспирантов и молодых ученых "Фундаментальная математика и ее приложения в естествознании"(г. Уфа, октября-3 ноября 2007 г.); научно-практической конференции "Прикладная математика и информационные технологии в науке и образовании"(г. Сибай, 23-24 мая 2008 г.); международной научной конференции "Дифференциальные уравнения и смежные проблемы"(г. Стерлитамак, 24-28 июня 2008 г.);

международной научной конференции "Нелинейные уравнения и комплексный анализ"(г. Уфа, 13-17 декабря 2010 г.); научном семинаре по дифференциальным уравнениям математической физики Института математики с В - УН - РАН (руководители: д.ф.-м.н., профессор Л.А. Калякин и д.ф.м.н., профессор В.Ю. Новокшенов.); научных семинарах кафедры прикладной математики и информационных технологий Сибайского института (филиала) Башкирского государственного университета (руководитель - д.ф.м.н., профессор Юмагулов М.Г.); научном семинаре кафедры дифференциальных уравнений Башкирского государственного университета (руководители: д.ф.-м.н., профессор Султанаев Я.Т и д.ф.-м.н., профессор Юмагулов М.Г.).

Публикации. Основные результаты диссертации опубликованы в работах [1]Ц[14], при этом статьи [1]Ц[4] опубликованы в изданиях, рекомендованных ВАК.

ичный вклад соискателя. Постановки основных задач принадлежат научному руководителю. Основные результаты диссертации получены автором самостоятельно. При выполнении работ [2], [3], [6], [7], [9] и [11], опубликованных в соавторстве, соискатель принимал участие в разработке и обосновании предлагаемых методов исследования. Из результатов этих работ в диссертацию автором включены только результаты, полученные им лично.

Объем и структура диссертации. Работа состоит из введения, трех глав, заключения и библиографии. Общий объем работы составляет 128 страниц. Библиография содержит 92 наименования.

Краткое содержание работы Во введении обоснована актуальность выбранной темы, сформулированы цель и основные задачи исследования, приводится обзор литературных источников, кратко излагается содержание работы.

Основным объектом исследования работы является система дифференциальных уравнений, зависящих от скалярного или векторного параметра с T -периодической по t правой частью:

x = f(x, t, ), x RN, Rk, (1) где функция f(x, t, ) непрерывна по t и непрерывно дифференцируема по x и . Пусть система (1) при всех значениях параметра имеет нулевую точку равновесия x = 0, т.е. f(0, t, ) 0. Уравнение (1) может быть представлено в виде x = A(t, )x + a(x, t, ), (2) где A(t, ) = fx(0, t, ) - матрица Якоби вектор-функции f(x, t, ), вычисленная в точке x = 0, а нелинейность a(x, t, ) равномерно по t и удовлетворяет соотношению a(x, t, ) = O( x ) при x 0; здесь и ниже через обозначена евклидова норма векторов в пространстве RN.

Обозначим через V () матрицу монодромии линейной системы x = A(t, )x. (3) Пусть при = 0 система (3) имеет один или несколько мультипликаторов, равных по модулю 1; тогда точка равновесия x = 0 системы (1) при = 0 является негиперболической. В этом случае значение 0 будем называть точкой бифуркации системы (1).

В частности, при близких к 0 значениях у системы (1) в окрестности точки равновесия x = 0 могут возникать T -периодические решения (вынужденные колебания), kT -периодические решения при k 2 (субгармонические колебания), квазипериодические решения и др. Описанию возможных сценариев бифуркационного поведения дифференциальных уравнений c периодическими коэффициентами, получению достаточных признаков того или иного сценария бифуркации, разработке схем приближенного исследования бифуркаций, исследованию устойчивости возникающих колебаний и посвящена диссертационная работа.

В первой главе приводятся необходимые сведения из теории дифференциальных уравнений с периодическими правыми частями, из теории локальных бифуркаций, описываются основные сценарии таких бифуркаций. Основное содержание главы составляет обоснование метода качественного и приближенного исследования задачи о бифуркациях малых решений операторных уравнений. Глава носит вспомогательный характер. Приведем некоторые необходимые сведения из первой главы, относящиеся к указанному методу.

Рассматривается операторное уравнение x = B()x + b(x, ), x RN, Rm, (4) где матрица B() непрерывно диффиренцируемо зависит от параметра , а нелинейность b(x, ) содержит слагаемые второй и более высоких степеней по x: b(x, ) = O( x ) при x 0.

Пусть e RN - некоторый ненулевой вектор. Значение 0 параметра называют правильной точкой бифуркации уравнения (4) по направлению вектора e, если существуют непрерывные функции = () и x = x() такие, что: (0) = 0, x(0) = 0; x()-e = o() при 0; при этом для каждого 0 вектор x() является решением уравнения (4) при = ().

Правильные точки бифуркации уравнения (4) имеет смысл искать лишь среди тех 0, при которых матрица B(0) имеет собственное значение 1. Рассматриваемые в работе задачи о вынужденных колебаниях таковы, что приводят либо к уравнению вида (4), зависящего от скалярного параметра , и с матрицей B(), имеющей при некотором = 0 простое собственное значение 1, либо к уравнению вида (4), зависящего от двумерного параметра , с матрицей B(), имеющей при некотором = 0 полупростое собственное значение 1 кратности 2. При этом в работе основное внимание уделяется разработке операторной схемы приближенного исследования задачи о точках бифуркации вынужденных колебаний системы (1) для второго случая, т.е. для бифуркации коразмерности 2.

Пусть параметр является двумерным, т.е. = (, ), где и - скалярные параметры. Положим 0 = (0, 0), B0 = B(0, 0). Пусть матрица Bимеет полупростое собственное значение 1 кратности 2. Пусть e, g, e, g - линейно независимые собственные векторы операторов B0 и B0, отвечающие полупростому собственному значению 1. Эти векторы можно выбрать исходя из соотношений:

(e, e) = (g, g) = 1, (e, g) = (g, e) = 0. (5) Положим (B(0)e, e) (B(0)e, e) Q =. (6) (B(0)e, g) (B(0)e, g) Здесь B и B - операторы, полученные дифференцированием оператора B(, ) по и соответственно.

Теорема 1. Пусть det Q = 0. Тогда 0 является правильной точкой бифур кации уравнения (4) по направлению вектора e.

Во второй главе приводятся новые достаточные признаки бифуркаций вынужденных и субгармонических колебаний системы (1), а также получены асимптотические формулы для возникающих колебаний.

Пусть система (3) при = 0 имеет один или несколько мультипликаторов равных по модулю 1, т.е. матрица V (0) имеет собственные значения вида e2i, где 0 . В диссертации основное внимание уделяется рассмотрению случаев, когда матрица V (0) имеет:

S1) простое собственное значение 1;

S2) пару простых собственных значений e2i, где 0 и рациональp но: = - несократимая дробь.

q Во обоих случаях предполагается, что остальные собственные значения матрицы V (0) не равны 1 по модулю.

В зависимости от этих случаев возможны различные сценарии локальных бифуркаций в окрестности состояния равновесия системы (1).

В случае S1) коразмерность бифуркации равна одному; в этом случае естественным будет предположение, что параметр является скалярным. Здесь основным сценарием бифуркации является возникновение у системы (1) в окрестности точки равновесия x = 0 при переходе параметра через 0 (т.е.

ибо при 0, либо при 0) ненулевых T -периодических колебаний малой амплитуды. Данная бифуркация соответствует следующему понятию.

Значение 0 параметра назовем точкой бифуркации вынужденных колебаний системы (1), если каждому > 0 соответствует такое = (), при котором система (1) имеет ненулевое T -периодическое решение x(t, ), причем () 0 и max x(t, ) 0 при 0.

t В случае S2) коразмерность бифуркации равна двум. Поэтому здесь естественным будет предположение, что параметр является двумерным, т.е.

= (, ), где и - скалярные параметры. Здесь в окрестности точки равновесия x = 0 при переходе параметра через 0 могут возникать или исчезать периодические решения различных периодов kT, где k q. Данная бифуркация соответствует следующему понятию.

Значение 0 параметра назовем бифуркацией субгармонических колебаний периода qT cистемы (1), если каждому > 0 соответствует такое = (), при котором система (1) имеет ненулевое qT -периодическое решение x(t, ), причем () 0 и max x(t, ) 0 при 0.

t Рассмотрим сначала случай S1). Обозначим через e и e - собственные векторы, соответствующие простому собственному значению 1 матриц V (0) и V (0) соответственно; здесь V () - транспонированная к V () матрица. Векторы e и e можно выбрать в соответствии с равенствами e = 1, (e, e) = 1. Обозначим через x0(t, ) решение задачи Коши x = A(t, )x, (7) x(0) = e, и положим 0 = [x0(T, )] =, e.

Теорема 2. Пусть выполнено условие S1 и 0 = 0. Тогда 0 является точ кой бифуркации вынужденных колебаний системы (1).

q Рассмотрим теперь случай S2). Положим B() = V (); тогда матрица B(0) (а вместе с ней и транспонированная матрица B(0)) имеет полупростое собственное значение 1 кратности 2. Обозначим через e, g и e, g соответствующие линейно независимые собственные векторы матриц B(0) и B(0). Векторы e, g, e и g будем считать выбранными в соответствии с соотношениями (5). Пусть x0(t, ) - решение задачи Коши (7).

Положим ([x 0(qT, )]=, e) ([x 0(qT, )]=, e) 0 0 =.

([x 0(qT, )]=, g) ([x 0(qT, )]=, g) 0 Теорема 3. Пусть выполнено условие S2 и 0 = 0. Тогда 0 = (0, 0) явля ется точкой бифуркации субгармонических колебаний периода qT системы (1).

Во второй главе работы также приведены асимптотические формулы для существующих в условиях теорем 2 и 3 бифурцирующих решений. Полученные в работе асимптотические формулы используют квадратичную и кубическую нелинейности, входящие в правые части основных дифференциальных уравнений. Здесь ограничимся приведением асимптотических формул для уравнения (2), когда нелинейность a(x, t, ) представима в виде a(x, t, ) = a2(x, t, ) + a3(x, t, ), (8) где a2(x, t, ) содержит квадратичные по x слагаемые, а нелинейность a3(x, t, ) удовлетворяет соотношению a3(x, t, ) = O( x ) при x 0 равномерно по t и .

Теорема 4. Существующие в условиях теоремы 2 бифурцирующие решения x(t, ) системы (1) и соответствующие значения параметра () представимы в виде:

x(t, ) = e(t) + 2e1(t) + e11(t, ), () = 0 + 1 + 11(). (9) Здесь e(t) = x0(t, 0); для функции e1(t) и числа 1 в диссертации также получены расчетные формулы, а функции 11(), e11(t, ) удовлетворяют соотношениям 11() = O(2), max e11(t, ) = O(3) при 0.

t Теорема 5. Существующие в условиях теоремы 3 бифурцирующие решения x(t, ) системы (1) и соответствующие значения параметра () = ((), ()) представимы в виде:

x(t, ) = e(t) + 2e1(t) + e11(t, ), (10) () = 0 + 1 + 11(), () = 0 + 1 + 11().

Здесь e(t) = x0(t, 0); для функции e1(t), чисел 1, 1 в диссертации также получены расчетные формулы, а функции 11(), 11(), e11(t, ) удовлетворяют соотношениям 11() = O(2), 11() = O(2), max e11(t, ) = O(3) t при 0.

Полученные во второй главе асимптотические формулы вида (10) позволяют получить приближенные представления языков Арнольда - множеств соответствующих рационально синхронизированным соотношениям параметров и . Функции = () и = () на плоскости параметров (, ) описывают некоторую кривую, начинающуюся в точке (0, 0). Эту кривую называют кривой синхронизации. Она получена в предположении, что в качестве направления бифуркации выбран один из собственных веторов, отвечающим полупростому собственному значению 1 кратности 2 матрицы B(0), а именно вектор e. Для различных собственных векторов получаются разные кривые синхронизации. Совокупность этих кривых и образует язык Арнольда.

Основной задачей исследования в третьей главе является анализ устойчивости бифурцирующих решений, возникающих в условиях теорем 2 и 3.

Ограничимся рассмотрением случая S1). Здесь основным объектом исследования является система (2), в котором матрица A(t, ) представима в виде A(t, ) = A0 + ( - 0)A1(t) + A2(t, ).

Предполагается, что постоянная матрица A0 имеет простое собственное значение 0, матрицы A1(t) и A2(t, ) являются T -периодическими, причем max A2(t, ) = O(| - 0|2), при | - 0| 0.

t Предлагаемая схема анализа устойчивости использует методы М. Розо.

исследования устойчивости решений дифференциальных уравнений с периодическими коэффициентами.

Определим число T = ( [1A1() + a 2x(e(), , 0)]de, e).

Теорема 6. Пусть нелинейность a(x, t, ) имеет вид (8). Тогда при всех малых значениях > 0 бифурцирующие решения x(t, ) системы (2), существующие в условиях теоремы 2, асимптотически устойчивы, если < 0, и неустойчивы, если > 0.

Здесь a 2x(x, t, ) - матрица Якоби квадратичной нелинейности a2(x, t, ).

Автор выражает благодарность своему научному руководителю М.Г. Юмагулову за неоценимую помощь в работе.

ПУБЛИКАЦИИ АВТОРА ПО ТЕМЕ ДИССЕРТАЦИИ В рецензируемых изданиях из списка ВАК 1. Муртазина С.А. Бифуркации вынужденных колебаний в системах автоматического регулирования // Вестник Тамбовского университета: общие проблемы управления и их приложения (ОПУ-2009): материалы международной конференции Колмогоровские чтения, Тамбов, 2009. T.14.

Вып. 4. С. 773 - 775.

2. Вышинский А. А., Ибрагимова Л. С., Муртазина С. А., Юмагулов М. Г.

Операторный метод приближенного исследования правильной бифуркации в многопараметрических динамических системах // Уфимский математический журнал, Уфа, 2010. Т 2, № 4. С. 3 - 26.

3. Юмагулов М.Г., Вышинский А.А., Нуров И.Д., Муртазина С.А. Операторный метод исследования локальных бифуркаций многопараметрических динамических систем // Вестник Санкт-Петербургского госуниверситета: прикладная математика, информатика, процессы управления, С.-Петербург, 2009. Серия 10. Вып. 2. С. 146 - 155.

4. Муртазина С.А. Исследование устойчивости вынужденных колебаний в многопараметрических системах // Вестник Тамбовского университета:

общие проблемы управления и их приложения (ОПУ-2011): материалы международной конференции Колмогоровские чтения - V, Тамбов, 2011.

T.16. Вып. 4. С. 1135 - 1137.

В других изданиях 5. Муртазина. С.А. Метод малого параметра в задачах приближенного построения малых автоколебаний. //Новые программные средства для предприятий Урала: сборник трудов региональной научно-технической конференции, Магнитогорск, 2004. Вып. 3. С. 199 - 201.

6. Вышинский А.А. Муртазина С.А. Приближенное исследование бифуркации малых решений операторных уравнений // Новые программные средства для предприятий Урала: сборник научных трудов, Магнитогорск, 2006. Вып. 5. С. 100 - 102.

7. Вышинский А. А., Муртазина С. А. Асимптотические формулы в задаче о бифуркациях нелинейных колебаний // Фундаментальная математика и ее приложения в естествознании: материалы Всероссийской школыконференции для студентов, аспирантов и молодых ученых, Уфа, 2007.

С. 20 - 21.

8. Муртазина С.А. Алгоритмы приближенного построения периодических колебаний в двупараметрических системах автоматического управления //Создание и внедрение корпоративных информационных систем (КИС) на промышленных предприятиях РФ: сборник трудов Международной научo-технической конференции, Магнитогорск, 2007. Вып. 2. С. 300 - 301.

9. Юмагулов. М.Г., Муртазина С.А. Коразмерность бифуркации векторных полей //Прикладная математика и информационные технологии в науке и образовании: материалы научно-практической конференции, Сибай, 2008. C. 102 - 109.

10. Муртазина С.А. Асимптотические формулы в задаче о бифуркации вынужденных колебаний в многопараметрических системах управления // Прикладная математика и информационные технологии в науке и образовании: материалы научно-практической конференции, Сибай, 2008. C.

59 - 65.

11. Юмагулов. М.Г., Муртазина С.А. Бифуркация вынужденных колебаний в многопараметрических системах управления // Дифференциальные уравнения и смежные проблемы: труды международной научной конференции, Стерлитамак, 2008. Т.3. C. 51 - 55.

12. Муртазина С.А. Признаки бифуркации вынужденных колебаний в двупараметрических системах // Уральский регион Республики Башкортостан: человек, природа, общество: материалы региональной научнопрактической конференции, Сибай, 2009. C. 365 - 369.

13. Муртазина С.А. Бифуркация субгармонических колебаний в многопараметрических динамических системах // Дифференциальные уравнения и их приложения: труды Всероссийской научной конференции с международным участием, Стерлитамак, 2011. C. 113 - 115.

14. Муртазина С.А. Расчет устойчивости вынужденных колебаний многопараметрических динамических систем // Математическое и программное обеспечение систем в промышленной и социальной сферах: международный сборник научных трудов, Магнитогорск, 2011. Ч. 1. С. 198 - 203.

Авторефераты по всем темам  >>  Авторефераты по разным специальностям