WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!

Pages:     || 2 | 3 |
Сибирский математический журнал Ноябрь декабрь, 2005. Том 46, № 6 УДК 517.929 БИФУРКАЦИОННЫЙ МЕТОД ИССЛЕДОВАНИЯ УСТОЙЧИВОСТИ РЕШЕНИЯ ДИФФЕРЕНЦИАЛЬНОГО УРАВНЕНИЯ С ЗАПАЗДЫВАНИЕМ Ю. Ф. Долгий, С. Н. Нидченко Аннотация: Исследуется устойчивость антисимметрических периодических решений дифференциальных уравнений с запаздыванием. Вводится однопараметрическое семейство периодических решений специальной системы обыкновенных дифференциальных уравнений с переменным периодом. Условия устойчивости антисимметрического периодического решения дифференциального уравнения с запаздыванием формулируются в терминах функции этого периода.

Ключевые слова: дифференциальное уравнение с запаздыванием, периодическое решение, устойчивость.

1. Постановка задачи Рассмотрим скалярное дифференциальное уравнение с постоянным запаздыванием dx(t) = -f(x(t - )), (1.1) dt где f нечетная непрерывно дифференцируемая функция с положительной производной на интервале (-, ), > 0. В работе изучаются вопросы существования и устойчивости периодических решений, удовлетворяющих условию антисимметричности x(t + 2) = -x(t), t (-, +). Для дифференциальных уравнений с запаздыванием изолированные периодические решения исследовались в работах [1, 2], периодические решения с периодом, кратным запаздыванию, в работах [3, 4], а их устойчивость в работах [4, 5].

Используя результаты работы [5], задачу нахождения антисимметрического периодического решения дифференциального уравнения с запаздыванием (1.1) сводим к нахождению решения специальной краевой задачи для системы обыкновенных дифференциальных уравнений:

x = f(x2), x = -f(x1), (1.2) 1 2 x1() = x2(0), x2() = -x1(0). (1.3) Работа выполнена при поддержке программы фундаментальных исследований Президиума РАН Математические методы в нелинейной динамике (№ 15).

© 2005 Долгий Ю. Ф., Нидченко С. Н.

Бифуркационный метод исследования устойчивости 1289 Здесь связь между антисимметрическим периодическим решением x дифференциального уравнения с запаздыванием (1.1) и решением {x1, x2} краевой задачи (1.2), (1.3) определяется формулами x1(t), t [0, ], x(t) = (1.4) x2(t - ), t [, 2].

Система (1.2) продолжается с отрезка [0, ] на всю числовую ось и имеет первый интеграл F (x1) + F (x2) = C = const, x1, x2 (-, ), (1.5) x где F (x) = f(z) dz, x (-, ).

Формула (1.5) определяет расположение интегральных кривых системы (1.2) на фазовой плоскости. Специальным начальным условиям x1(0, µ) = 0, x2(0, µ) = µ, соответствуют замкнутые интегральные кривые, порождающие периодические движения {x1(t, µ), x2(t, µ)}, t R, при 0 < µ <. Периоды T решений системы (1.2) зависят от µ.

Система (1.2) при малых значениях µ является системой Ляпунова [6, с. 153].

Для нахождения асимптотик периодических решений этой системы при малых значениях µ можно воспользоваться методом Ляпунова [6, с. 159].

Утверждение 1.1. Пусть f нечетная трижды непрерывно дифференцируемая функция в окрестности точки x = 0, f (0) > 0. При малых значениях параметра µ периодические решения {x1(t, µ), x2(t, µ)}, t R, периода T определяются формулами xi(t, µ) = yi t, µ, t R, i = 1, 2, T (µ) где yi(s, µ), i = 1, 2, s R, 2-периодические функции. При малых значениях параметра µ функции y1(s, µ), y2(s, µ) на отрезке [0, 2] и функция периодов T задаются асимптотическими разложениями 2 f (0) T (µ) = - µ2 + o(µ2), f (0) 4(f (0)) f (0) y1(s, µ) = sin (s)µ + (- sin3 (s) cos2 (s) - sin5 (s) + sin (s))µ3 + o(µ3), 24f (0) f (0) y2(s, µ) = cos (s)µ + (- sin2 (s) cos3 (s) - cos5 (s) + cos (s))µ3 + o(µ3).

24f (0) Утверждение 1.2. Пусть f нечетная непрерывно дифференцируемая функция с положительной производной на интервале (-, ). Тогда для существования периодического решения дифференциального уравнения с запаздыванием (1.1), удовлетворяющего условию антисимметричности, необходимо и достаточно, чтобы число 4 принадлежало интервалу ( inf T (µ), sup T (µ)).

0<µ< 0<µ< Этот интервал дополняется точкой t = inf T (µ) (t = sup T (µ)), если 0<µ< 0<µ< inf T (µ) = T (µ1) ( sup T (µ) = T (µ2)) для некоторого µ1 (0, ) (µ2 (0, )).

0<µ< 0<µ< При выполнении условий утверждения 1.2 антисимметрическое периодическое решение x(t, µ), t R, дифференциального уравнения с запаздыванием 1290 Ю. Ф. Долгий, С. Н. Нидченко определяется по формулам (1.4) через периодическое решение системы обыкновенных дифференциальных уравнений (1.2) {x1(t, µ), x2(t, µ)}, t R, для которого T (µ) = 4, 0 < µ <. Таких решений у дифференциального уравнения с запаздыванием (1.1) может быть несколько. Их число определяется количеством положительных решений уравнения T (µ) = 4 при 0 < µ <.

Дифференциальное уравнение с запаздыванием имеет единственное периодическое решение, если функция периода T монотонна.

При исследовании устойчивости антисимметрического периодического решения x(t, µ), t R, дифференциального уравнения с запаздыванием (1.1) рассмотрим уравнение линейного приближения dy(t) = -f (x(t - T (µ)/4, µ))y(t - T (µ)/4) (1.6) dt для уравнения возмущенного движения. В уравнении (1.6) функция a(t, µ) = f (x(t - T (µ)/4, µ)), t R, принимает только положительные значения и периодически зависит от t с периодом T (µ)/2.

В настоящей работе предлагается использовать бифуркационный подход в задаче исследования устойчивости периодического решения. Он связан с переходом от исследования устойчивости дифференциального уравнения с запаздыванием (1.6) к изучению устойчивости однопараметрического семейства дифференциальных уравнений с запаздыванием dy(t) = -f (x(t - T (µ)/4, µ))y(t - T (µ)/4), µ [0, ). (1.7) dt 2. Устойчивость линейной периодической системы дифференциальных уравнений с запаздыванием Проведя в дифференциальном уравнении с запаздыванием (1.7) замену переменных:

T (µ) T (µ) T (µ) t = s, (s) = y s, x(s, µ) = x t, µ, 2 2 находим d(s) T (µ) = - f (x(s - /2, µ))(s - /2), µ [0, ). (2.1) ds Коэффициент полученного дифференциального уравнения при малых µ близок к постоянной, т. е. дифференциальное уравнение с запаздыванием является квазигармоническим.

Утверждение 2.1. Пусть выполняются условия утверждения 1.1. Тогда при малых положительных значениях µ квазигармоническое дифференциаль ное уравнение с запаздыванием (2.1) устойчиво, если T (0) > 0, и неустойчиво, если T (0) < 0.

Доказательство. При µ = 0 дифференциальное уравнение с запаздываd(s) нием (2.1) имеет вид = -(s-/2). Запишем его характеристическое уравds нение = -e-(/2). Используя метод Д-разбиения, можно показать, что оно имеет только два чисто мнимых корня = ±i, а остальные корни имеют отрицательные действительные части. Следовательно, устойчивость дифференциального уравнения с запаздыванием (2.1) при малых положительных µ определяется поведением характеристических показателей этого уравнения, которые Бифуркационный метод исследования устойчивости при µ = 0 совпадают с числами ±i. Паре чисто мнимых корней характеристического уравнения отвечает один полупростой характеристический показатель 0 = i дифференциального уравнения с запаздыванием (2.1). Ему отвечают два линейно независимых решения Флоке: y1(s) = eis, y2(s) = e-is, s R.

Квазигармоническое уравнение с запаздыванием (2.1) имеет -периодическое решение y(s, µ) = (s, µ), s R, которое является решением Флоке этого уравнения с характеристическим показателем = i. В силу двухкратности характеристического показателя 0 квазигармоническое дифференциальное уравнение с запаздыванием (2.1) имеет еще один характеристический показатель (µ) ((0) = 0 = i). Если действительная часть этого характеристического показателя при малых положительных значениях параметра µ больше нуля, то квазигармоническое уравнение неустойчиво. Если действительная часть характеристического показателя (µ) при малых положительных значениях параметра µ меньше нуля, то критический характеристический показатель = i простой, остальные характеристические показатели имеют отрицательные действительные части, а квазигармоническое дифференциальное уравнение с запаздыванием устойчиво. Воспользовавшись результатами утверждения 1.1 и методикой вычисления характеристических показателей для квазигармонических дифференциальных уравнений с запаздыванием [7], находим f (0)T (0) (µ) = i - µ2 + o(µ2).

4 + Утверждение доказано.

Устойчивость дифференциального уравнения с запаздыванием (2.1) при больших значениях параметра µ зависит от расположения собственных чисел оператора монодромии, действующего в пространстве C[-, 0]. Оператор монодромии задается формулой (U)() = y( +, ), [-, 0], где y( + ·, ) отрезок решения дифференциального уравнения с запаздыванием (2.1) с начальным моментом t = 0 и начальной функцией.

Задачу нахождения ненулевых собственных чисел C оператора монодромии U можно заменить задачей нахождения ненулевых собственных чисел z C специальной краевой задачи [8, с. 49] 1 = za(, µ)y2, 2 = -za(/2 +, µ)y1, (2.2) y1(-/2) = -zy2(0), y2(-/2) = zy1(0), (2.3) T (µ) где = -z2, a(, µ) = f (x( - /2, µ)), [-/2, /2], µ [0, ), z, C.

Приведем краевую задачу (2.2), (2.3) к специальному виду J = zH(, µ)y, (2.4) y(-/2) = zJy(0). (2.5) Здесь y = (y1, y2), z C, 0 -1 a(/2 +, µ) J =, H(, µ) =, 1 0 0 a(, µ) µ [0, ), [-/2, 0]. Собственные числа z краевой задачи (2.4), (2.5) определяются из характеристического уравнения det(zJY (0, z, µ) - I2) = 0, (2.6) 1292 Ю. Ф. Долгий, С. Н. Нидченко где Y фундаментальная нормированная матрица системы (2.4), Y (-/2, z, µ) = I2, z C, µ [0, ), I2 единичная матрица порядка 2. Используя формулу Лиувилля, находим det(Y (, z, µ)) = exp T r(H(s, µ)) ds = 1, z C, µ [0, ), [-/2, 0].

-/В результате характеристическое уравнение (2.6) преобразуем к виду D(z, µ) = z2 - 2V (z, µ)z + 1 = 0, (2.7) где V (z, µ) = (y12(0, z, µ)-y21(0, z, µ)), Y (, z, µ) = yij(, z, µ) 2, [-/2, 0], z C, µ [0, ). Рассмотрим некоторые свойства фундаментальной нормированной матрицы Y системы (2.4).

Лемма 2.1. Для фундаментальной матрицы Y справедливо равенство Y (-/2 -, µ, z) = JY (, µ, -z)J Y (0, µ, z), (2.8) где [-/2, 0], z C, µ [0, ).

Доказательство. Заменяя аргумент на 1, перепишем систему дифференциальных уравнений (2.2):

1 = za(1, µ)y2, 2 = -za(/2 + 1, µ)y1, 1 [-/2, 0].

Проведя в этой системе дифференциальных уравнений замену переменных:

1 = -/2 -, v1() = y2(-/2 - ), v2() = -y1(-/2 - ), [-/2, 0], получим v1 = -za(, µ)v2, v2 = za(/2 +, µ)v1. (2.9) Произвольное решение {v1, v2} системы дифференциальных уравнений (2.9) представимо в виде (v1(, µ, z), v2(, µ, z)) = J y(-/2 -, µ, z), (2.10) где y(, µ, z) = (y1(, µ, z), y2(, µ, z)), [-/2, 0], z C, µ [0, ), некоторое решение системы дифференциальных уравнений (2.2). Тогда (v1(, µ, z), v2(, µ, z)) = J Y (-/2 -, µ, z)D, где [-/2, 0], z C, µ [0, ), D C2. В силу произвольности решения {v1, v2} векторная постоянная D может принимать любые значения из C2.

Решение {v1, v2} также допускает представление (v1(, µ, z), v2(, µ, z)) = y(, µ, -z) = Y (, µ, -z)C (2.11) при некотором C C2. Справедливо равенство Y (, µ, -z)C = J Y (-/2 -, µ, z)D, [-/2, 0], z C, µ [0, ).

Полагая в этом равенстве = -/2, находим C = J Y (0, µ, z)D, z C, µ [0, ). Следовательно, Y (, µ, -z)J Y (, µ, z)D = J Y (-/2-, µ, z)D, [-/2, 0], z C, µ [0, ).

В силу произвольности D C2 из последнего равенства следует (2.8).

Бифуркационный метод исследования устойчивости Лемма 2.2. Для фундаментальной матрицы Y справедливо равенство Y (-/2 -, µ, z) = SY (, µ, z)Y (0, µ, z), (2.12) 0 где [-/2, 0], z C, µ [0, ), S =.

1 Доказательство. Проведя в системе дифференциальных уравнений (2.2) замену переменных: 1 = -/2 -, v1() = y2(-/2 - ), v2() = y1(-/2 - ), [-/2, 0], получим новую систему дифференциальных уравнений v1 = za(, µ)v2, v2 = -za(/2 +, µ)v1.

Произвольное решение {v1(, µ, z), v2(, µ, z)}, [-/2, 0], z C, µ [0, ) новой системы дифференциальных уравнений представимо в виде (v1(, µ, z), v2(, µ, z)) = Y (, µ, z)C, где [-/2, 0], z C, µ [0, ), C C2. Но это же решение можно представить в виде (v1(, µ, z), v2(, µ, z)) = SY (-/2 -, µ, z)D, где [-/2, 0], z C, µ [0, ), D C2. Тогда Y (, µ, z)C = SY (-/2 -, µ, z)D для всех [-/2, 0], z C, µ [0, ). Из этого равенства при = -/вытекает, что C = SY (0, µ, z)D. Тогда Y (, µ, z)SY (0, µ, z)D = SY (-/2, µ, z)D, [-/2, 0], z C, µ [0, ).

В силу произвольности D C2 из последнего равенства следует (2.12).

Лемма 2.3. Пусть f нечетная непрерывно дифференцируемая функция с положительной производной на интервале (-, ). Тогда функция V удовлетворяет условиям (а) V (1, µ) = 1, µ [0, );

(б) V (-z, µ) = -V (z, µ), z C, µ [0, ).

Доказательство. Можно показать, что уравнение (2.1) имеет антисимdx(s,µ) метрическое 2-периодическое решение (s, µ) =, s R, µ [0, ). Ему ds отвечает собственное число оператора монодромии = -1, так как для антисимметрического периодического решения дифференциального уравнения с запаздыванием (1.1) имеем тождество (/2 +, µ) -(, µ), [-/2, 0].

Этому числу соответствуют два собственных числа z = ±1 краевой задачи (2.4), (2.5). Из уравнения (2.7) следует, что функция V удовлетворяет условиям V (±1, µ) = ±1, µ [0, ). Первое утверждение леммы доказано.

Покажем, что y12(0, µ, z) = -y21(0, µ, z), µ [0, ), z C. По лемме 2.имеем Y (0, µ, z)SY (0, µ, z) = S, µ [0, ), z C, или, в развернутой форме, a11(µ, z) a12(µ, z) 0 =, a21(µ, z) a22(µ, z) 1 где a11(µ, z) = y11(0, µ, z)(y21(0, µ, z) + y12(0, µ, z)), a12(µ, z) = y12(0, µ, z) + y11(0, µ, z)y22(0, µ, z), a21(µ, z) = y11(0, µ, z)y22(0, µ, z) + y21(0, µ, z), 1294 Ю. Ф. Долгий, С. Н. Нидченко a22(µ, z) = (y12(0, µ, z) + y21(0, µ, z))y22(0, µ, z), µ [0, ), z C.

Пусть y12(0, µ, z) + y21(0, µ, z)) = 0 при некоторых µ [0, ), z C. Тогда из матричного равенства вытекает, что y11(0, µ, z) = y22(0, µ, z) = 0, y12(0, µ, z) = y21(0, µ, z) = ±1.

В этом случае det Y (0, µ, z) = -1; противоречие. Следовательно, y11(0, µ, z) y12(0, µ, z) Y (0, µ, z) =, µ [0, ), z C.

-y12(0, µ, z) y22(0, µ, z) -Теперь согласно лемме 2.1 Y (0, -z) = J Y (0, z)J, или, в развернутом виде, y11(0, µ, -z) y12(0, µ, -z) y11(0, µ, z) -y12(0, µ, z) =, -y12(0, µ, -z) y22(0, µ, -z) y12(0, µ, z) y22(0, µ, z) где µ [0, ), z C. Из матричного равенства находим y12(0, µ, -z) = -y12(0, µ, z), µ [0, ), z C.

Следовательно, V (-z, µ) = y12(0, µ, -z) = -y12(0, µ, z) = -V (z, µ), µ [0, ), z C.

Лемма доказана.

При изменении параметра µ корни характеристического уравнения (2.7) перемещаются по комплексной плоскости симметрично относительно мнимой оси. Поэтому при изучении их движений можно ограничиться рассмотрением правой полуплоскости. Из леммы 2.3 следует, что характеристическое уравнение (2.7) имеет корни z = ±1, для каждого действительного корня z1 = существует действительный корень -z1 и для каждого комплексного корня zсуществуют комплексные корни -z2, z2, z2.

Pages:     || 2 | 3 |



© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.