УДК 517.925
О ПРИМЕНЕНИИ СИМВОЛЬНЫХ ВЫЧИСЛЕНИЙ ДЛЯ ПОСТРОЕНИЯ ОБОБЩЕННО ПЕРИОДИЧЕСКИХ РЕШЕНИЙ СИСТЕМ ОБЫКНОВЕННЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ С ПОЛИЛИНЕЙНОЙ ПРАВОЙ ЧАСТЬЮ
© М. А. Кириченко, Н. А. Рубанов, Э. П. Агабекян
Ключевые слова: метод последовательных приближений Пикара; построение нелокальных решений класса автономных дифференциальных уравнений; система и аттрактор Лоренца.
На основе метода последовательных приближений Пикара строится решение автономных систем обыкновенных дифференциальных уравнений, правая часть которых представляет собой многомерный многочлен.
1. Введение. Рассмотрим автономную нормальную систему обыкновенных дифференциальных уравнений
X = / (х); (1)
здесь х = (х1,...,хп) - векторная функция действительного переменного £, а / = (/1,...,/п) - действительная векторная функция, каждый элемент которой /г является многочленом переменных х1,... ,хп .
Системы вида (1) давно представляют большой интерес для приложений, поскольку многие модели процессов различной физической, технической и экономической природы описываются подобными системами (см., например, [1] или [2, гл. 3]). Особый интерес системы вида (1) приобрели в последнее время, поскольку все известные системы, предполагаемо содержащие странные аттракторы, оказались именно такими системами [2, с. 109-185].
Стандартные методы численного анализа, используемые для получения решения системы (1), не учитывают конкретный вид ее правой части и, следовательно, не могут быть
оптимальными. Целью настоящей работы является разработка метода построения решений
системы (1), учитывающего то, что каждый элемент /г функции / является многочленом переменных х1,... ,хп .
2. Локальные решения системы. Для построения решения х(Ь) системы (1), удовлетворяющего начальному условию
х(0) = хо, (2)
прежде всего, заменим (1) интегральным уравнением
£
х(Ь) = х0 + J /(х(т)) йт. (3)
о
Обозначим через а и д - некоторые положительное числа, а через Г - компактную часть (п + 1) -мерного евклидова векторного пространства Мп+1 , задаваемую неравенствами
\х — х0\ < а, |£| < д.
Пусть при этом г < д - положительное число, которое будет определено ниже. Наряду с Г введем в рассмотрение более "узкое" компактное множество Гг С Мп+1 , задаваемое неравенствами
\х — х0\ < а, \£\ < г. (4)
Далее, обозначим через Пг - множество непрерывных функций, графики, которых содержатся в Гг . Рассмотрим оператор
£
Лф = хо + У /(ф(т)) йт. (5)
0
Поскольку множество
£ = {х Е Мп : \х — х0\ < а} (6)
компактно, найдется такое положительное число М, что для всех х Е £ выполнено неравенство
\/(х)\< М. (7)
Из условия
а
г < М (8)
следует, что оператор (5) отображает множество Пг в себя (см., например, [3, с. 10]). Поэтому предположим, что число г в (4) выбрано так, что выполнено неравенство (8).
Заметим теперь, что система (1) имеет единственное решение х(Ь) с начальным условием (2), определенное на некотором отрезке [—Т,Т] . Поэтому уравнение (3) также имеет единственное решение х(Ь) , определенное на отрезке [—Т, Т] . Для отыскания этого решения будем использовать метод последовательных приближений Пикара
£
хм+1(Ь) = х0 + ^ /(хм(т)) йт, (9)
0
где начальное приближение примем равным
х1(Ь) = х0. (10)
Несложно заметить, что существует некоторая итерация Лрф = Л-у Л V оператора Л,
р
являющаяся сжатием (см., например, [3, с. 11]). Следовательно, метод последовательных
приближений (9),(10) на отрезке [—г, г] равномерно сходится к решению х(Ь) . Остается
построить последнее.
3. Нелокальные ограниченные решения Пусть теперь х(Ь) - некоторое решение системы (1) с начальным условием (2), определенное для всех значений Ь > 0 и содержащееся при этих значениях Ь в множестве £ , задаваемом равенством (6), в котором а -некоторое надлежащим образом подобранное положительное число. Тогда найдется такое
достаточно большое положительное число М , что при х Е £ выполнено неравенство (7).
а
Для определенных выше чисел а и М зададим число Т = ~^^. При всех значениях Ь > 0 положим
хк(Ь) = х(Ь + (К — 1)Т), К = 1, 2,... (11)
Легко видеть, что каждая из функций семейства (11) является решением системы (1), определенным для всех значений Ь > 0 и содержащимся при этих значениях Ь в множестве £ .
С другой стороны, в силу ограниченности решения х(Ь) число М изначально может быть подобрано так, чтобы при выполнении для всех значений К = 1, 2,... условий
\х — хк(Ь)\ < а, Ь < Т
было также выполнено неравенство (7).
4. Приближенное построение решений. Вновь рассмотрим решение х(Ь) системы (1) с начальным условием (2), определенное для всех значений Ь > 0. Обозначим через
х(0),х(Т),..., х(КТ),... (12)
- положения системы (1) в моменты времени
0,Т,...,КТ,... (13)
Далее, введем в рассмотрение оператор д£ сдвига вдоль решений системы дифференциальных уравнений (1). Имеем х(Ь) = д£х0. В частности, х(Т) = дтх0, и, следовательно,
х(КТ) = дт .. .дт х0, К = 0,1,... (14)
к
Таким образом, если оператор д£ построен, соотношение (14) однозначно задает положения (12) системы (1) в моменты времени (13).
5. Численный эксперимент. В качестве простейшего приложения приведенной схемы вычислений рассмотрим задачу приближенного построения решений системы Лоренца
{х = а (у — х),
у = вх — у — хг, (15)
г = ху — Ьг,
с начальными условиями х(0) = х0, у(0) = у0, г(0) = г0, где а, в и Ь - некоторые положительные числа, играющие роль параметров данной системы. Для простоты ограничимся рассмотрением следующего случая классических значений параметров [1]: а = 10, в = 28, Ь = 8/3.
Положим
с = (х0,у0,г0), дМс = (хМ,у%,г%), N = 1,2,...
Для заданной точности е признаком окончания итеративного процесса построения оператора дМ в силу критерия сходимости Коши служило выполнение неравенства
а вычисления выполнялись с точностью до сорокового знака. Траектории системы (15) строились для широкого спектра начальных условий. Здесь, в частности, рассмотрим случай траектории К , описываемой решением с начальным условием
Хо = -15.720831, Уо = -16.587193, = 36.091132, (16)
взятым в непосредственной близости от аттрактора.
На рис. 1 приведена проекция этой траектории на плоскость хОу .
Таблица 1: Решение системы (15),(16), описывающее типическую траекторию К .
№ Ь х(Ь) у(Ь) г(Ь)
1 0 -15.720831 -16.587193 36.091132
2 17.334 -15.659134 -16.566101 35.954224
3 45.017 -15.652555 -16.566140 35.938103
4 67.104 -15.689783 -16.635305 35.964610
5 86.686 -15.812414 -16.750982 36.164029
х
Рис. 1: Проекция на плоскость хОу дуги типической траектории К , построенной по точкам (14) на отрезке времени [0, 90] .
6. Выводы. Не обнаружено каких-либо убедительных признаков существования циклов в системе Лоренца (см. табл. 1); чтобы полностью убедиться в этом достаточно, например, построить векторное поле системы (15) в приведенных здесь точках. Поэтому поскольку существенного снижения скорости движения точки не зафиксировано, скорее всего, можно считать, что К - компактное минимальное множество, состоящее из незамкнутых рекуррентных траекторий (см., например, [4, с. 402]). Эти траектории, возможно, не являются даже почти периодическими, так как некоторое локальное разбегание типических траекторий в ходе вычислений наблюдалось (см., например, [4, с. 415]).
ЛИТЕРАТУРА
1 . Lorenz E.N. Deterministic Nonperiodic Flow // J. Athmos. Sci. 1963. V. 20. P. 130-141.
2 . Магницкий Н.А., Сидоров С.В. Новые методы хаотической динамики. М.: URSS, 2004.
3 . Шварц Л. Анализ. Т. 2. М.: Мир, 1972.
4 . Немыцкий В.В., Степанов В.В. Качественная теория дифференциальных уравнений. М.: URSS,
2004.
БЛАГОДАРНОСТИ: Работа выполнена при поддержке РФФИ (проекты № 10-07-00136, № 11-07-00098).
Поступила в редакцию 15 мая 2011 г.
Kirichenko M.A., Rubanov N.A., Agabekyan E.P. About application the method of successive approximations of picard for construction non-local solutions of autonomous systems of ordinary differential equations. Based on the method of successive approximations of Picard based solution of autonomous systems of ordinary differential equations, the right side which is a multivariate polynomial.
Key words: the method of successive approximations of Picard; the construction non-local solutions of a class of autonomous differential equations; the system and the Lorenz attractor.
Кириченко Михаил Александрович, Тамбовский государственный технический университет, г. Тамбов, Российская Федерация, аспирант кафедры «Прикладная математика и информатика», e-mail: kirimedia@gmail.com
Рубанов Никита Александрович, Тамбовский государственный технический университет, г. Тамбов, Российская Федерация, аспирант кафедры «Прикладная математика и информатика», e-mail: nikitarubanov@gmail.com
Агабекян Эмиль Паргевович, Тамбоский государственный технический университет, г. Тамбов, Российская Федерация, аспирант кафедры «Экономический анализ», e-mail: emill2007@yandex.ru