CC BY
38
2004
Доклады БГУИР
январь- март
№ 1
УДК 517.925
АНАЛИТИЧЕСКИЕ СВОЙСТВА РЕШЕНИЙ УРАВНЕНИЙ ПЕНЛЕВЕ-ТИПА И НЕКОТОРЫЕ ИХ ПРИЛОЖЕНИЯ
В В. ЦЕГЕЛЬНИК
Белорусский государственный университет информатики и радиоэлектроники П. Бровки, 6, Минск, 220013, Беларусь
Поступила в редакцию 18 ноября 2003
Для каждого из шести уравнений Пенлеве получены новые формулы, определяющие общий вид ассоциированных с ними полиномиальных гамильтонианов. Показано существование ассоциированных с данными уравнениями гамильтонианов, отличных от полиномиальных.
Ключевые слова: уравнения Пенлеве-типа, гамильтониан, солитон, преобразования Беклун-да.
Введение
На протяжении четырех последних десятилетий теория волн в нелинейных средах и аналитическая теория уравнений с частными производными пережили — параллельно и при взаимном влиянии — значительный подъем.
Этот подъем был стимулирован нуждами физической науки, в различных областях которой — физике плазмы, нелинейной оптике, физике ферромагнетиков — в начале шестидесятых годов стали систематически возникать проблемы взаимодействия волн большой амплитуды. Опыт показал, что при всем огромном многообразии существенно нелинейных задач основную роль здесь играет небольшое число универсальных математических моделей. Многие из этих моделей обладают уникальным свойством интегрируемости, позволяющим глубоко и эффективно провести их исследование. Для исследования вопроса интегрируемости систем Гарднером, Грином, Крускалом и Миурой в 1967 г. был предложен метод обратной задачи рассеяния (ОЗР) [1]. Суть математического открытия, сделанного указанными авторами, состояла в том, что они обнаружили связь уравнений Кортевега-де Фриза (KdV) с линейным уравнением Шредингера на прямой и предложили метод точного решения нелинейных уравнений с частными производными (НУЧП), получивший название метода ОЗР.
Следует отметить, что открытию метода ОЗР предшествовали глубокие исследования по теории солитонов. Так, например, еще в XVIII в. было установлено, что уравнение KdV и уравнение sin-Gordon имеют локализованные точные решения — солитоны. Уже тогда для уравнения sin-Gordon, изучавшегося в связи с его применением в теории поверхностей отрицательной кривизны, был открыт способ "размножения" точных решений — преобразование Бек-лунда.
Применение метода ОЗР позволило установить точную интегрируемость большого числа как НУЧП, так и обыкновенных дифференциальных уравнений (ОДУ), имеющих важные приложения. В настоящее время имеются десятки содержательных монографий и обзоров (не претендуя на полноту списка отметим наиболее известные и значимые, на наш взгляд, работы [2-8]), связанных с методом ОЗР и его приложениями. Поскольку в определении интегри-
руемых систем в существующей литературе имеются разночтения, то более подробная информация по данному вопросу содержится в [9-11].
Одна из важных особенностей, связанных с открытием метода ОЗР заключается в том, что исследование НУЧП, допускающих применение метода ОЗР, позволило установить тесную связь между последними и ОДУ, решения которых не имеют подвижных критических особых точек , в частности, шестью уравнениями Пенлеве
= 6w2 + г, (Рг)
w'' = 2w3 + ™ + а, (Р2)
,2 ' 1 £
w" = — - — +1 (аw2 +в) + YW3 + —, (Р3)
w г г —
—' = —--^^ + 4гм^1 + 2(г2 -а> + —, (Р4)
2w 2 w
„ 3w -1 ,2 w' ^ -1 )2 р, у 5w(w +1)
w =-w--+ --(аw + —) + — w + —!--, (Р5)
2w( w -1) г г w г w -1
'1 ( 1 1 1 ^ ,2 ( 1 1 1 ^ '
w =—\ — +-+-I w - I — +-+-к +
21 w w -1 w - г I I г г -1 w - г I
+ w(w -1)^ - г)
г2(г -1)2
рг , у(г -1) , 5г(г -1)
а +—т +--2 +--т"
w ^ -1 )2 ^ - г)2
( Р6 )
с произвольными параметрами а, в, у, 5 , решения которых в общем случае не сводятся к известным трансцендентным функциям. Указанное обстоятельство способствовало тому, что уравнения Пенлеве в последние десятилетия являются объектом пристального внимания благо -даря не только обнаружению их многочисленных приложений к различным задачам теоретической и математической физики, но и благодаря их богатой внутренней симметрии. Именно этот факт позволил П. Бутру еще задолго до открытия метода ОЗР описать некоторые асимптотические свойства решений уравнений Пенлеве, а Р. Фуксу, Л. Шлезингеру и Р. Гарнье найти связь трансцендентов Пенлеве с изомонодромными деформациями дифференциальных операторов. Н.П. Еругин, А.И. Яблонский, Н.А. Лукашевич, В.И. Громак и др. нашли классы точных и од-нопараметрических семейств решений и описали трансформационные свойства уравнений Пенлеве также без учета связи между НУЧП, интегрируемыми с помощью метода ОЗР.
История появлений уравнений Пенлеве, возникающих при их решении задачи и полученные результаты достаточно подробно изложены в работах [12-15].
Говоря о приложениях уравнений Пенлеве, необходимо указать на использование этих функций для описания определенных переходных и автомодельных режимов, а также подчеркнуть тот факт, что трансценденты Пенлеве играют в нелинейной теоретической физике ту же роль, что и классические специальные функции (типа Эйри, Бесселя, Вебера-Эрмита) для линейных дифференциальных уравнений с частными производными [16-19].
Несмотря на то что уравнения (Р2)-(Р6) получили название неприводимых (трансцендентных), длительное время существовали разные точки зрения на проблему сведения этих уравнений к линейным алгебраическим дифференциальным уравнениям. В настоящее время доказано, что уравнения Пенлеве не сводятся к линейным уравнениям [20]. Следовательно, ре-
Это свойство решений называют Р-свойством, а уравнения с Р-свойством решений Пенлеве-типа или Р-типа.
— уравнениями
шения каждого из уравнений (Р2)-(Р6) в общем случае представляют новые трансцендентные функции. Однако при определенных значениях параметров их решения могут быть выражены через элементарные функции (например, (Р3) при а = у = 0 или в = 8 = 0; (Р5) при у = 8 = 0).
Для уравнений (Р2)-( Р6) получены рекуррентные соотношения, или преобразования Беклунда, как их обычно называют в этом случае. При помощи преобразований Беклунда для уравнений Пенлеве можно построить иерархии решений уравнений (Р2)-( Р6), которые являются либо рациональными, либо алгебраическими, либо функциями, связанными (через повторное дифференцирование и умножение) с классическими трансцендентными функциями: Эйри, Бесселя, Вебера-Эрмита, Уиттекера и гипергеометрическими функциями.
Гамильтонианы, ассоциированные с уравнениями Пенлеве и их свойства
Еще одно важное свойство (наряду с другими, отмечавшимися ранее) шести трасцен-дентных уравнений Пенлеве состоит в том, что каждое из данных уравнений допускает представление в виде системы Гамильтона
йм = дИ^ йи = дИ^ (1)
йг ди ' йг дм
где И= И(г,м,и) (] = 1, 6) — полиномиальный относительно м, и гамильтониан, соответствующий уравнению (Р^).
Представление уравнений Пенлеве в виде системы (1) впервые рассматривалось в [21] и тесно связано с проблемой деформации уравнений, сохраняющих группу монодромии [22].
К. Окамото использовал [23] представление уравнений Пенлеве в виде системы (1) для построения уравнений второго порядка Р-типа с иррациональной правой частью. С этой целью им были построены уравнения, определяющие функции вида г) = ф( г) И}(г),и(г)),
где И^ — гамильтониан, отвечающий уравнению (Р^), а функция ф(г) определена следующим образом: ф(г) = 1 при ] = 1,2, 4; ф(г) = г для ] = 3, 5 ; ф(г) = г -1), если ] = 6. Полученные уравнения, которые будем называть ^-уравнениями, так же, как и уравнения (Р[) — (Р6), имеют многочисленные приложения в нелинейной физике, в частности, в квантовой теории поля [24-26], а свойства их решений во многом схожи со свойствами решений уравнений Пенлеве.
В данном разделе для каждого из уравнений Пенлеве будут получены новые формулы, определяющие общий вид ассоциированных с ними полиномиальных гамильтонионов, включающие простейшие (ранее известные). Будут указаны также другие, ассоциированные с уравнениями (Р[)-(Р6) гамильтонианы, и рассмотрены их свойства. 1. Полиномиальный гамильтониан вида
и2 3
И, =--+ 2м + zw.
1 2
ассоциированный с уравнением (Р[), впервые получен Мальмквистом в работе [21]. Легко видеть, что функция
и2
И, =--2м3 - (2)
12
также определяет гамильтониан, ассоциированный с уравнением (Р). Укажем формулу, определяющую общий вид полиномиального гамильтониана, ассоциированного с уравнением (Р1 ). Справедлива [27]
Теорема 1. Пусть Рт (г, w), дп (г, w) — многочлены относительно w степени т, п е 2о соответственно с аналитическими по г коэффициентами такие, что
дРт (г, *) _дQn (г, *)
дг дн>
Тогда выражение
(3)
Я, = [" + Р.2г.'')Г + д^) - 2\\>1 - „ (4)
определяет гамильтониан, ассоциированный с уравнением (Р1) , с точностью до произвольной (зависящей от г) функции и преобразования " ^ р", Я ^ рЯ, где р — отличный от нуля параметр.
Рассмотрим эквивалентную уравнению (Р1) систему w' = " + S(w), "' = -[" + '(w) + 6w2 + г , (5)
отвечающую гамильтониану ТТ [" + S(w)]2 0 3
Я3 = --^—— - 2w - zw , (6)
2
где S(w) — аналитическая функция w . Вводя функцию Их(г) = Я3 (г, w(г),"(г)) , находим, что
И'(г) = И'(г) = -" - S(w) .
Тогда из (5) с учетом того, что Я 3 = И1 (г) находим уравнение И2 + 4И1'3 + 2(гИ[ - И ) = 0, (7)
w'2 3
определяющее функцию И(г) . Так как в силу (5) И1(г) = --2w - ^ , то уравнение (7)
является уравнением Р-типа. Отметим, что (7) с точностью до обозначений совпадает с уравнением [23] для простейшего гамильтониана Я1, ассоциированного с (Р1) .
Из формулы (6) следует существование рациональных (по переменной w) гамильтонианов, ассоциированных с уравнением (Р1) .
Примером отличного от полиномиального (по переменной w) гамильтониана является гамильтониан
2
Я =(-- 2wi - ^ + lnw , (8)
4 2
ассоциированный с уравнением (Р1) . Характерной особенностью гамильтониана (8) является
то, что определяющее функцию И1(г) = Я4 (г,w(г),"(г)) дифференциальное уравнение не является (в отличие от уравнения (7)) уравнением Р-типа.
Существование гамильтонианов (6), (8) объясняется тем, что справедлива [27]
Теорема 2. Пусть К( ),Ь( ) — аналитические функции переменных г, н такие, что Ц г, н) = I-ан. Тогда выражение
1 дг
[и + К( 2 )]2 , . „ 3
Н5 = --' 77 + Ь(г, н) - 2н - гн
5 2
определяет гамильтониан, ассоциированный с уравнением (Рх) . Преобразование и ^ ри, Н5 ^ рН5 с параметром р Ф 0 определяет гамильтониан Н5 также ассоциированный с уравнением (Р1) .
2. Нетрудно убедиться в том, что уравнение (/2) можно представить в виде эквивалентной ему системы Гамильтона с гамильтонианом
Н„ = [и + Р (2г,Н)]' -«V + *>[„ + /,(,,„)] + (9)
+ Оп(г,н) - (а + т)н, в2 = 1,
где Рт (г, н), Qn (г, н) — многочлены относительно н, степени т,п е Z0 соответственно с аналитическими по 2 коэффициентами, удовлетворяющие условию (3).
3. Ассоциированный с уравнением (Р3) полиномиальный гамильтониан в случае, когда у = 0, а5 Ф 0 (без ограничения общности можно считать а = -5 = 1) имеет вид
Нш = - гн2[и + Р (г, м>)]2 + 1-еР [и + рт(2,н)] + в[и + Рт( г,н)] + вн + Qn(z,н) . (10) 2 г г
Два ненулевых параметра у и 5 в уравнении (Р3) также всегда можно фиксировать, полагая а = -5 = 1. Ассоциированный с (Р3) в данном случае полиномиальный гамильтониан имеет вид
~ш= 2чн [и + Рт(-.,н>] 2 +вн2 [и + рМ] - (1 )[и + Р.(г,н)1н +
г
г
+ а [и + Рт(г,н)] + Qm(z, н) - в( ав +2) н, а2 = 1. (11)
4. Уравнение (Р4) можно представить в виде эквивалентной ему системы Гамильтона с гамильтонианом
Н1у = 2вн[и + Рт(г,н)]2 -вн2[и + Рт(г,н)] - 2вгн[и + Рт(г,н)] -- а [и + Рт(г,н)] + Ьн + Qm(z,w), а2 + 2р = 0, 2Ь = 2в-1 + . (12)
5. Для уравнения (Р5) в случае 5 Ф 0 (без ограничения общности можно считать 25 = -1) ассоциированный с ним полиномиальный гамильтониан имеет вид
Ну = ^[н(н -1 )2у12 - а(н -1 )2 V - Ь(н -1 )ну1 + вгну1 + $(н -1) + Qm(г,н), V = и + Рт(г,н), а2 + 2р = 0, Ь = -^-1, ^ = 4[(а + Ь)2 - 2а ] . (13)
в
6. Формула ассоциированного с уравнением (Р^) полиномиального гамильтониана имеет вид
~ 1 2
Ип =-[и(и -1)(и - 1)У2 - В0(и -1)(и - 1)У2 - Г{Н>(и - 1)У2 + (т2 +1 )и(и - 1)у2 -
г(г -1)
- к2 - (е.- 1)2 (и - 2) - (2 - 1)90 Г2 - ^ ] + дт(г> и), (14)
У2 = и + Рт(г,и), Во2 + 2в = О, Г!2 = 2у, г22 = 1 - 25, (В.-1)2 = 25, к = 1 -В0 - Г + Г2 .
Отметим, что в формулах (10)-(14) многочлены Рт(1,и), связаны услови-
ем (3).
Некоторые классы уравнений нелинейной физики, связанные с уравнениями Пенлеве-типа
1. В работах [9, 28, 29] дано описание и проведена классификация систем вида и = ихх + 1(и,У,их,Ух), - V = ^ + g(u,v,ux,vx), (15)
интегрируемых методом обратной задачи рассеяния и обладающих бесконечной серией законов сохранения. Из таких систем отметим нелинейные системы интегрируемых уравнений типа Шредингера (список содержит более 40 систем), среди которых можно выделить несколько основных (опорных) систем, связанных с помощью подходящих преобразований с остальными системами из списка (см. [28], [29]).
Поскольку системы (15) интегрируемы с помощью метода обратной задачи рассеяния, то они (согласно гипотез [6]) могут иметь автомодельные решения, выражающиеся через решения уравнений Р-типа.
Пусть в (15) / = 2иих + vx, g = -2(ш)х . Тогда имеем систему
и = ихх + 2иих + ^ , V = + 2(uv)x, (16)
моделирующую процесс распространения слабо нелинейных волн на поверхности воды в приближении длинных волн. С помощью преобразования
и(х,г) = (-г)'172 р(г), 1 = х(-г)-17 2 (17)
от системы (16) перейдем к системе обыкновенных дифференциальных уравнений
2 (р + 1р' ) = р' + д' + 2рр', д +1 щ' = -дщ + 2(рд)'. (18)
Систему (18) можно записать в виде
21р = р' + д + р2 - 2 а,д +11д' =-д" + 2(рд)', (19)
где а - произвольная постоянная. Уравнение, определяющее функцию р, имеет вид
1 3 а
р" - 6 р2р' - ^ 12р' + 31рр' + ар ' + 2 р2 - ^ 1р - — = 0. (20)
Выполняя в (20) преобразование р ^ р\, 1 ^ 2 т, придем к уравнению
1 т 3 2 1 2 г 3 , а , 2З
8а-4А А -2т А + 2тР1А + ^А + А -2^-а =
Полагая в (21) Р1 = у + 2т относительно функции у, получаем уравнение
у " = 6у2у ' + 12туу' + 4(т2 - а)у' + 4у2 + 4ту, (22)
первым интегралом которого является уравнение (Р4) с произвольными параметрами а, в и независимой переменной т . Таким образом, с учетом (17) имеет место [30]
Теорема 3. Пусть у = у(т) — решение уравнения (Р4) при произвольных значениях параметров а, в . Тогда функции
и(х,0 = у( т ) + 2т
v(x,г) = (4г)~
^ + у2 + 2 ту + 2( 1 -а) ат
, 2т = 1, 1 = х(-г)
-1/2
являются решениями системы (16).
Наряду с системой (16) рассмотрим систему связанных уравнений Шредингера
Рг = Рх, + р2д> - д = дхх+д2р. (23)
Известно, что если функции р,д — решения этой системы, то функции
Рх
и = , v = рд р
определяют [28] решение системы (16). Преобразованием
р(х,г) = ехр[хг+3 г3 ] р( т), (24)
д(х,г) = ¡42ехр[-хг-2г3]д(т), т = х + г2, /2 =-1
система (23) редуцируется к системе обыкновенных дифференциальных уравнений
а2р „ 2 а2д 2 __
= 2р2д + тр, —2 = 2рд2 +тд, (25)
ат ат
которая при условии р = д сводится к уравнению (Р2) в случае а = 0 , а при равенстве д = 0 - к уравнению Эйри с независимой переменной т . Имеет место [30]
Теорема 4. Пусть и = и(т) — решение уравнения (Р2) при произвольном значении параметра а . Тогда функции
р = С ехр[г| и(т)ат], (26)
д = (2Сехр[-е|и(т)От](и2 + в- т),
в = -1, С — произвольная постоянная (С Ф 0), являются решениями системы (25).
Используя соотношения (24), (26) можно получить явные формулы автомодельных решений системы (23), выражающихся через решения уравнения (Р2) .
2. Для исследования сложного динамического поведения нелинейных систем чаще всего применяются методы численного моделирования и физические эксперименты. Однако эти
методы и эксперименты не всегда обеспечивают полную адекватность описания нелинейных систем, если они не дополнены аналитическими методами.
Система
х = -ух + ьу + уг, у = -уу - ьх + хг, у = Т - ху (27)
с постоянными параметрами V,ь,Т (точка означает производную по независимой переменной г) является предложенной [31, 32]) для описания временного изменения магнитного поля Земли моделью двудискового динамо, в которой при определенных положительных значениях параметров возможно наличие хаотических движений. Очевидно, что явные аналитические решения системы (27) в этом случае невозможны и фактически ожидается, что они будут сложными функциями времени с подвижными особыми точками в плоскости комплексного переменного г.
Один из возможных вопросов, который возникает после сказанного выше, сводится к следующему: существуют ли значения параметров V, ь, Т , для которых общее решение системы (27) может быть выражено через известные аналитические функции, подвижными особыми точками которых могут быть только полюсы. Если отвлечься от физического смысла параметров системы (27), то (ниже переменную г считаем комплексной) ответ на поставленный вопрос положителен. Так, например, в работе [33] показано, что решения системы (27) в случае ь = Т = 0 выражаются через решения третьего уравнения Пенлеве в одном частном случае. Аналитические свойства решений системы (27) исследовались с использованием метода Пен-леве-анализа. При этом было отмечено (подобная ситуация имеет место и для других трехмерных динамических систем с квадратичными нелинейностями, которые сводятся к известным моделям, обладающим хаотическим поведением), что наличие у системы (27) свойства Пенлеве гарантирует отсутствие у нее решений с хаотическим поведением. Данное обстоятельство при исследовании конкретной системы позволяет, с одной стороны, сузить область поиска параметров, при которых ее решения обладают хаотическим поведением, а с другой стороны, дает возможность получить условия для параметров (обеспечивающие существование) и строить [34 - 37] специальные классы точных решений, обладающих свойством Пенлеве.
Для полноты изложения отметим, что ветвление решений принято связывать с хаотической динамикой и отсутствием интегралов - законов сохранения [38].
Замечание. Метод Пенлеве-анализа ( в основе которого лежит преложенная С.В. Ковалевской идея исследования интегралов уравнений движения твердого тела около неподвижной точки) первоначально был использован для исследования аналитических свойств решений известной системы Лоренца
X = а(у - х), у = - у + рх - хг, г = -Ьг + ху, (28)
2
имеющей в случае а = 10, р = 28, Ь = — странный аттрактор. Сигур получил 4 набора параметров [33], при которых система (28) либо вырождается в линейную систему, либо ее решения выражаются через эллиптические функции, либо через функции-решения второго или третьего уравнения Пенлеве.
Заключение
Современный интерес к уравнениям Пенлеве-типа объясняется тем, что последние (как уже отмечалось выше) возникают как редукции многих вполне интегрируемых уравнений в частных производных. Обобщенные симметрии таких уравнений дают возможность записать высшие аналоги ОДУ. Таким образом, можно получить иерархии ОДУ из вполне интегрируемых уравнений в частных производных. В связи с этим остается открытой проблема классификации решений высших аналогов уравнений Пенлеве.
ANALYTICAL PROPERTIES OF SOLUTIONS OF PAINLEVE' TYPE EQUATIONS AND THEIR APPLICATIONS
V.V. TSEGELNIK Abstract
For each of the six Painleve' equations were received new formulas which define the general form of with them polinomial namiltonians associated. There have been shown the existence of associated with given equations hamiltonians different from polinomial.
Литература
1. Gardner C.S., Greene I.M., KruskalM.D., Miura R.M. //Phys. Rev. Lett. 1967. Vol. 19. P. 1095-1097.
2. Теория солитонов: метод обратной задачи / В.Е. Захаров, С.В. Манаков, С.П. Новиков, Л.П. Питаев-ский. М., 1980.
3. Солитоны / Под ред. Р.Буллаф, Ф. Кодри. М., 1983.
4. Калоджеро Ф., Дегасперис А. Спектральные преобразования и солитоны. М., 1985.
5. Тахтаджян Л.А., Фаддеев Л.Д. Гамильтонов подход в теории солитонов. М., 1986.
6. Абловиц М., СигурX Солитоны и метод обратной задачи. М., 1987.
7. Солитоны и нелинейные волновые уравнения / Под ред. Р. Додд, Дж. Эйлбек, Дж. Гиббон, Х. Моррис. М., 1988.
8. Ньюэлл А. Солитоны в математике и физике. М., 1989.
9. Интегрируемость и кинетические уравнения для солитонов: Сб. науч. тр. / Ред. В.Г. Барьях-тар, В.Е. Захаров, В.Н. Черноусенко. Киев, 1990.
10. Маханьков В.Г., Рыбаков В.Ю., Санюк В.И. // Успехи физических наук. 1994. Т.164, №2. С. 121-148.
11. KruskalM.D., Clarkson P.A. // Stud. Appl.Math. 1992. Vol.86. P. 87-165.
12. Громак В.И., Лукашевич Н.А. Аналитические свойства решений уравнений Пенлеве. Мн., 1990.
13. The Paileve' property. One century later. CRM series in Mathematical Physics Ed. R.Conte. New-York, 1999.
14. Gromak V.I., Laine I., Shimomura Sh. Painleve' differential equations in the complex plane. De Gruyter Studies in Mathematics. Berlin - New-York, 2002.
15. Цегельник В.В. Уравнения Пенлеве-типа: аналитические свойства решений и их приложения: Авто-реф. дис. ... д-ра физ.-мат. наук. Мн.: Ин-т математики НАН Беларуси. 2002.
16. AblowitzM.I., Segur H. // Stud. Appl. Math. 1977. Vol. 57. Р. 43-44.
17. Новиков В.Ю. // Дифференц. уравнения. 1985. Т. 21, №11. С. 1915-1926.
18. Манаков С.В. // Письма в журнал эксперим. и теорет. физики. 1982. Т. 35, вып. 5.С. 193-195.
19. Creamer D.B., Thacker H.B., Wilkinson D. // Phys. Rev.D. 1981. Vol. 23, №12. P. 3081-3084.
20. Громак В.И. // Тр. Ин-та матем. НАН Беларуси. 2000. Т.6. С. 65-74.
21. MalmquistI. // Ark. Mat. Astr. Fys. 1922-23. Vol.17. P. 1-89.
22. Fuchs R. // Comptes Rendus. 1905. Vol. 141. P. 555-558.
23. Okamoto K. // Proc. Japan Acad. 1980. Vol. 56A. P. 264-268; Р. 367-371.
24. Сато М., Дзимбо М., Мива Т. Голономные квантовые поля. М., 1983.
25. Jimbo M, Miwa T, Ueno K. // Physica D. 1981. Vol. 2. P. 306-352.
26. Sato M, Miwa T, Jimbo M. // Proc. Japan Acad. 1977. Vol.53A. P. 6-10.
27. Цегельник В.В. // Докл. НАН Беларуси. 2001. Т. 45, №5. С. 45-47.
28. Михайлов А.В., Шабат А.Б., Ямилов Р.И. // Успехи матем. наук. 1987. Т.42, вып.4. С. 3-53.
29. MikhailovA.V., ShabatA.B., YamilovR.I. // Commun. Math. Phys. 1988. Vol. 115. P. 1-89.
30. Цегельник В.В. // Диф. уравнения. 1994. Т.30, №6. С. 992-997.
31. Rikitake T. // Proc. of the Cambridge Phylos. Society. 1958. Vol. 54, Part 1. P. 89-105.
32. Нелинейные волны. Сб. статей / Отв. ред. А.В. Гапонов-Грехов. М., 1979.
33. Bountis T.C., Ramani A., Grammaticos B, Dorizzi B. // Physica. 1984. Vol. 128A, №1-2. P. 268-288.
34. ГромакВ.И., ЦегельникВ.В. // Диф. уравнения. 1991. Т. 27, №3. С. 396-402.
35. Цегельник В.В. // Изв. вузов. Прикладная нелинейная динамика. 1998. Т.6, №5. С. 84-88.
36. Цегельник В.В. // Диф. уравнения. 1999. Т.35, №7. С. 1003-1004.
37. Tsegelnik V.V. //Regular and chaotic dynamics. 1999. Vol. 4, №4. P. 77-80.
38. Козлов В.В. // Прикладная математика и механика. 1998. Т. 62, №1. С. 3-11.
