Исследование вопроса о существовании периодических решений Текст научной статьи по специальности «Математика»

7. Хартман Ф. Обыкновенные дифференциальные уравнения / Ф. Хартман. М. : Мир, 1970. 720 с.

8. Шестаков А. А. Об асимптотическом поведении решений нелинейной системы дифференциальных уравнений / А. А. Шестаков, А. У. Пайвин // ДАН. 1948. Т. 58, № 5.

С. 495 498.

9. Hartman P. Asymptotic Integrations of linear Differential Equations / P. Hartman, A. Wintner // Amer. J. Math. 1955, 77. P. 45 87.

10. Perron O. Üeber lineare Differentialgleichungen, dei denen die unab-hangig Variable reel list, I und II / O. Perron // Journ. für die Reine und Ange-wandte Mathematik. 1913. Bd. 142. S. 254 270.

11. Perron O. Über stabilität und asymptotishen Verchalten der Integrate von Differentialgbichungs systemen / O. Perron // Math. Zeitschrift. 1929. Bd. 29. S. 129 160.

12. Perron O. Über eine Matrixtransformation / O. Perron // Math. Zeitschrift. 1930. Bd. 32. S. 465 473.

13. Petrowsky I. Üeber das Verhalten der Integralkurven eines Systems gewhnlicher Differentialgleichungen in der Nähe eines singulären Punktes / I. Petrowsky // Мат. сб. 1934.

Т. 41. № 1. С. 107 155.

14. Petrowsky I. Nachtag zu meiner Arbeit «Über das Verhalten der Inte-gralkurven eines systems gewöhnlicher Differential gleichungen in der Nähe eines singulären Punktes» / I. Petrowsky // Мат. сб. 1935. Т. 42, № 3. С. 403.

Поступила 13.02.2012.

УДК 517.912

ИССЛЕДОВАНИЕ ВОПРОСА О СУЩЕСТВОВАНИИ ПЕРИОДИЧЕСКИХ РЕШЕНИЙ*

А. В. Зубов, О. А. Зубова, О. А. Иванова, К. А. Пешехонов

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

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

Постановка задачи. Сначала рассмот-

рим систему обыкновенных дифференциальных уравнений:

= Л (х1, ■■■, хк, гЪ ■■■, 2п-к),

2) = д}- ■■■, Ч, ■■■, ■■■, 2п-к) > (1) 5 = 1, ■■■, к, ) = 1, ■■■, п - к

Условия существования и единственности решений системы (1) будем считать выполненными. Пусть правые части системы (1) являются периодическими функциями

* Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (РФФИ) (проект № 10-08-000624).

© Зубов А. В., Зубова О. А., Иванова О. А., Пешехонов К. А., 2012

30 ВЕСТНИК Мордовского университета | 2012 | № 2

¿1, ..., с периодом 2р. Пусть решения

системы (1) являются уходящими по ..., гп_ь т. е.:

1) существует неограниченная последовательность значений параметра {Ьп} : Ьп ^ да,

выполнено х5 (Ьп ) ^ х5 е (-да, +да);

2) для любой неограниченной последовательности {Ьп} : Ьп ^ да выполнено

п^да

г 1 (£и ) ^да-

Рассмотрим решение системы (1) с начальными данными:

Х0, ] = 1, ..

х0, г\

к, X

Хо + 2яжу, [0,2 л],

Зх л ■ 1 Зх л —- + Ё —

Ы к=1

Л

Ё Ск]х] + Ё + дк (ь)

Ё + Ё Ву-г,- + ^ (t), г - 1,

которая, являясь начальным условием задачи Коши для системы (2), определяет 2л-пе-риодическое решение (3).

Введем в рассмотрение векторы

X = (Х1,

¿и

(4)

ф = (/1 ■■■, /т gl, ■■■, дт)

и рассмотрим линейную систему

X = РХ,

где Р — блочная матрица

р-(С 3-

Пусть У(£) — фундаментальная матрица системы (4). Разобьем ее на блоки размеров, соответствующих размерам матриц А, В, С, В:

где т, целые. Если в системе (1) сделать замену переменных = X] + 2лгП] то правые

части не изменятся, система перейдет «сама в себя»:

Х5 = / (X, х), X] = д (X, х),

5 = 1, ..., к, ] = 1, ..., п - к.

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

У -

У11 У12 У21 У22

(2)

Теорема 1. Если матрица Уц(2л) — Е является неособенной, то у системы (1) существует единственное 2л-периодическое решение (4).

Доказательство. Введем обозначение

У (Ь, т) = У (Ь) У-1 (т). Приступим к изучению поведения решений системы дифференциальных уравнений характеристик системы (1) по отношению к координатам

Х1, ..., хп. Решение систем дифференциальных уравнений характеристик системы (1) имеет вид:

/=1 /-1

Будем считать, что матрицы А размерности п х п, В размерности п х т, С размерности т х п, В размерности т х т, векторные

функции / = (/1, ¡п), д = (gl, дт) являются 2л-периодическими функциями независимой переменной £. Наша задача будет заключаться в установлении условий, при которых у системы (1) существует 2л-перио-дическое по £ решение

х(^) (3)

и в построении этого решения. Наша задача будет решена, если будет найдена кривая Хо = Хо(2), или, в параметрической форме,

Х0 = Х(£), 2 = Z(t),

X = УХ о + I У (£, т) ф (т) йх.

(5)

Отсюда

(6)

X = УиХо + У^о + / (Уц (¿, т) / (т) +

¿о

+ У12 (¿, т) д (т)) йт,

1 = У21Хо + У22^) + I У21 ((¿, т) / (т) +

¿о

+ У22 (¿, т) д (т)) йт.

Положим £ = 2л, Х(2л) = Хо. Тогда из последних соотношений получим

Хо = (Е - Уп (2л))

(7)

У12 (2л) 2о + 2\ (Уп (^ т) / + У^д) dт *о

Подставим (7) во вторую группу формул (6) и выразим 2о как функцию

Ь, 2 : 2о = X (Ь, 2). Теперь, исключив ¿о из выражения (6), получим функцию

Хо = X (Ь, 2), (8)

которая, будучи заданной в качестве начального условия задачи Коши при £ = о, определяет 2л-периодическое решение системы (2). Чтобы убедиться в этом, достаточно повторить весь ход рассуждений в обратном порядке. На этом и закончим доказательство этой теоремы.

Рассмотрим систему дифференциальных уравнений

X = Г (£, X, т), (9)

где вещественная и непрерывная по совокупности своих аргументов векторная функция Г

задана при Ь е (-да, +да), X е Еп, т е [0, т] ■ Независимая переменная £ фактически входит в правую часть системы (9), которая является 2л-периодической функцией "X е Еп.

Пусть, кроме того, выполнены условия теоремы существования и единственности решений системы (9) в любой конечной области

G с Еп при Ь е (-да, +да) ■

Рассмотрим вопрос о существовании периодических решений у системы (9) при т = о, т. е. рассмотрим систему

X = ^ (£, X), (1о)

где функции Г удовлетворяют всем условиям, наложенным на правую часть системы (9). Для решения поставленной задачи попытаемся построить в Еп динамическую систему, определяемую системой дифференциальных уравнений, решения которой сходятся к точке Хо, обладавшей тем свойством,

что решение X (ь, Х°,Ь0) системы (1о) будет

2л-периодическим [2].

В случае успешного построения указанной дифференциальной системы мы получим не только критерии существования периодических решений, но и методы построения этих решений, основанные на численном ин-

тегрировании указанной дифференциальной системы.

Введем следующую функцию переменного 2 е Еп:

V (2) = 2 - X (2л,2,0). (11)

Функция У(2) определена в любой точке 2 е Еп, это следует из условий, наложенных на систему (1о). Потребуем, чтобы функция У(2) удовлетворяла на решениях искомой дифференциальной системы дифференциальному уравнению

dУ

dx

= -У.

(12)

Здесь и далее т — параметр искомой дифференциальной системы. Из этого следует:

ВХ (2л, 2,0) ^ й2 йх

(13)

В2

= -2 + X (2л, 2,0). Умножим это соотношение слева на матрицу (Е - 0X102) ; получим искомую дифференциальную систему:

dх

(Е - Р (2)) 1 (-2 + X (2л, 2,0)),

где

Р (2 ) =

БХ (£, 2,0)

Б2

£=2л

Матрица Р(2) удовлетворяет следующему дифференциальному соотношению:

DF (Х,£)

ВХ

Р,

(14)

где матрица Якоби DF (X,t)/DX вычисляется при X = X (Ь, 2,0).

Таким образом, если матрица Р(2) не имеет собственных чисел, вещественная часть которых равна единице, то правая часть полученной дифференциальной системы определена. Матрица Р(2) находится как решение задачи Коши для линейной системы (14) с начальным условием Р = Е при £ = о. Значение матрицы Р(2) берется в момент £ = 2л.

Для того чтобы функция У(2) экспоненциально стремилась к нулю в силу построенной дифференциальной системы при х ^ да, достаточно, чтобы решение 2 = 2 (т, 2о), где

32

ВЕСТНИК Мордовского университета | 2о12 | № 2

е Еп построенной дифференциальной системы было неограниченно продолжаемо в положительном направлении [1].

Если требование (12) можно заменить на следующее:

^=-ку, йх

где k — некоторая положительно определенная матрица, и матрицу k выбрать в виде

к = (Е - Р (2))(Е - Р (2)) , то система дифференциальных уравнений примет вид

¿2 ¿х В

(^ X0, о)

ния последовательности (2к). Для простоты возьмем за основу метод Эйлера интегрирования систем дифференциальных уравнений.

1. Выбирая некоторое положительное целое число Ы, оценку величины которого можно получить из (12) и погрешности метода Эйлера, строим последовательность Хд = 20,

Хк+1 = Х

2 к

N

к, 20,0 I, к = 1,

N.

= (Е - Р (2))* (-2 + X (2п, 2,0)). (15)

2. Вычисляем значение Р(2д) следую-

2 п

случае продолжаемости решения ^ = 2(т,2д) этой системы в положительном направлении будем двигаться по этому решению, например, с помощью численного метода типа Рунге — Кутты. Поскольку решение 2 = 2(х, 2д) сходится к значению X0, определяющему периодическое решение системы (10), причем траектория попадает в е-окрестность точки X0 за конечное время (оценку можно получить из (12), то можно доказать, что путем выбора параметров численного метода можно построить последовательность {2к}, сходящуюся к X0 [4].

Все вышеизложенное можно резюмировать в следующей теореме.

Теорема 2. Если существует точка 20 е Еп такая, что решение 2 = 2(х, 2^) системы (15) неограниченно продолжимо в положительном направлении, то это решение при х ^ да, сходится к значению X0, определяющему периодическое решение X (t, X системы (10).

Перейдем теперь к алгоритму построе-

щим образом:

. рр X) п

Р0 =

Р+1 = Р +

N

РХ

tk

2 п

N

к, к = 0,

N -1,

Р (20) = р

N.

21 = 20 + ^(Е

3. Вычисляем 21 по формуле

Р (2о))* (-2о + X(2п,2о,0)).

Таким образом, строится последовательность {2к} . Критерием окончания будет близость к нулю функции V, критерием выбора шага интегрирования будет неравенство

IV (2к+1)\\ < IV (2к)||.

Выводы. В данной статье вопрос о существовании периодического решения системы (10) сведен к вопросу о существовании продолжаемого решения системы (13). На основании многих известных достаточных признаков продолжаемости [3] можно предложить несколько достаточных условий существования периодического решения системы (10).

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Зубов А. В. Стабилизация и управление в динамических системах / А. В. Зубов. СПб. : СПбГУ, 2007. 132 с.

2. Зубов А. В. Динамическая безопасность управляемых систем / А. В. Зубов, Н. В. Зубов. СПб. : Изд-во НИИ химии СПбГУ, 2009. 172 с.

3. Зубов Н. В. Автоматизация проектирования устойчивости и надежности колебательных систем / Н. В. Зубов. СПб. : Мобильность-плюс, 2010. 355 с.

4. Зубов А. В. Теория устойчивости и теория квазипериодических систем / А. В. Зубов, Н. В. Зубов, С. А. Стрекопытов. СПб. : Мобильность-плюс, 2010. 206 с.

Поступила 22.02.2012.