Симметрии и законы сохранения для двухкомпонентного дискретного потенциированного уравнения Кортевега-де Фриза Текст научной статьи по специальности «Математика»

ISSN 2074-1871 Уфимский математический журнал. Том 8. № 3 (2016). С. 113-125.

УДК 517.9

СИММЕТРИИ И ЗАКОНЫ СОХРАНЕНИЯ ДЛЯ ДВУХКОМПОНЕНТНОГО ДИСКРЕТНОГО ПОТЕНЦИИРОВАННОГО УРАВНЕНИЯ КОРТЕВЕГА-ДЕ ФРИЗА

М.Н. ПОПЦОВА, И.Т. ХАБИБУЛЛИН

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

Ключевые слова: интегрируемые динамические системы, уравнения на квадратном графе, симметрии, законы сохранения, пара Лакса.

Mathematics Subject Classification: 35Q53, 37K10

1. Введение

Изучению асимптотического поведения системы линейных дифференциальных уравнений вблизи особой точки посвящено значительное количество исследований (см., например, монографию [1]). Асимптотическое представление собственной функции операторов Лакса позволяет эффективно описывать интегралы движения, высшие симметрии и частные решения соответствующей нелинейной динамической системы [2, 3]. Метод формальной асимптотической диагонализации позволил авторам работы [4] установить глубокую связь между интегрируемыми системами и аффинными алгебрами Ли.

Алгоритм решения задачи об асимптотической диагонализации дискретного оператора в окрестности особой точки и ее приложения в теории интегрируемых нелинейных дискретных уравнений подробно обсуждаются в работах [5, 6, 7]. Интересные результаты по неавтономным дискретным динамическим системам получены с использованием метода формальной диагонализации в работах [8, 9]. Альтернативные подходы к проблеме построения асимптотического разложения собственной функции дискретного оператора Лакса предложены в работах [10, 11, 12].

В настоящей работе мы рассматриваем найденное в [13] двухкомпонентное дискретное потенциированное уравнение Кортевега-де Фриза (cdpKdV)

(и - ui,i)(vi,Q - г>од) = р2 - q2, (v - Vi,i)(ui,o - Uo,i) = p2 - q2.

M.N. Poptsova, I.T. Habibullin, Symmetries and conservation laws for a two-component discrete potentiated Korteweg-de Vries equation. © ПопцовА М.Н., ХАБивуллин И.Т. 2016.

Исследование выполнено за счет гранта Российского научного фонда (проект №15-11-20007). Поступила 22 января 2016 г.

Пара Лакса для этого уравнения была предъявлена в [13], отметим также, что в [14] для него построены явные частные решения. Напомним, что однокомпонентное дискретное потенциированное уравнение Кортевега-де Фриза

(«1,1 - u)(ui,0 - uo,i) = 4с2 (1.1)

изучалось ранее во многих работах, начиная с [15, 16]. Оно известно также под названием уравнения H1 из списка Адлера-Бобенко-Суриса (см. [17]). Бесконечная серия законов сохранения для него была получена в работе [18], высшие симметрии построены в статьях [19, 20, 21, 22].

В представленной работе при помощи метода формальной асимтотической диагонали-зации пары Лакса описана бесконечная серия законов сохранения и построены высшие симметрии дискретного двухкомпонентного потенциированного уравнения Кортевега-де Фриза1.

Поясним кратко структуру работы. Во втором разделе описывается класс дискретных линейных уравнений с особыми точками, в окрестности которых строится асимптотическое решение. На уровне примера иллюстритруется способ приведения линейной системы к удобному специальному виду. В третьем разделе обсуждается алгоритм приведения оператора Лакса к квазидиагональному виду. В четвертом излагается известный метод построения производящей функции законов сохранения, использующий условие коммутирования диагонализованных операторов Лакса. В пятом разделе - метод формальной диагонализации применяется к конкретной динамической системе cdpKdV (5.1), для которой предъявляется бесконечная серия законов сохранения. И, наконец, в шестом разделе доказано, что динамическая система cdpKdV (5.1) имеет бесконечную иерархию высших симметрий. Первые две симметрии построены явно, для остальных указан эффективный способ вычисления.

2. Особенности типа полюсл дискретной линейной системы

Рассмотрим линейное дискретное уравнение вида

у(п + 1, А) = f (n,u(n),X)y(n,X), (2.1)

где потенциал f = f (n,u,X) G CNxN зависит от целого n G (-œ, +œ), функционального параметра u = u(n) и комплексного параметра A. Потенциал является мероморфной функцией от А в области Е С C и предполагается, что det f не равен тождественному нулю.

Определим, что мы называем особой точкой. Назовем точку А = Ао особой точкой уравнения (2.1), если хотя бы одна из функций f (n,u,X), f-1(n,u,X) имеет в этой точке полюс. Мы предполагаем, что А0 не зависит от п.

Отметим, что некоторые особые точки можно устранить при помощи преобразования зависимой переменной у(п) = г(п, Х)у(п), которое приводит уравнение (2.1) к тому же виду

у(п + 1) = f(n,u,X)y(n)

с новым потенциалом f (n, u, X) = г-1(п + 1, X)f (п, и, X)r(n, X).

В качестве простого иллюстративного примера рассмотрим уравнение (2.1) с потенциалом

t = ( t a£) • *< = *<м-

Это уравнение имеет две особые точки А = œ, А = 0. Обе особые точки удаляются преобразованием у(п) = Хпу(п). Действительно, f = X-1f = {gij}.

1Мы благодарим авторов работы [14], обративших наше внимание на эту задачу.

Менее тривиальный пример доставляет хорошо известное линейное уравнение, соответствующее в контексте интегрируемости уравнению (1.1). Его потенциал имеет вид

/(п,и,Х)=( +!} , п /ч V (2.2)

■> к ' ' ' \ -А 2 - и(п)и(п + !) и(п) J к 7

Здесь и(п) - произвольная функция. Уравнение (2.1), (2.2) имеет две особые точки Л = 0, Л = то. Особая точка Л = 0 удаляется преобразованием

( Х~п 0

у(п) = г(п,Х)у(п), г(п,Х)=[ п л-п_ 1

0 А"

Действительно, в данном случае новый потенциал приобретает форму

-и(п + !)Л !

?(п\)=( -и(п +!)Л ! ^ 1 (П,л)={ -! - и(п)и(п + !)А2 и(п)Х )

и имеет только одну особую точку Л = то.

Ключевым шагом алгоритма диагонализации является сведение исходной системы (2.1) к специальному виду

у(п + ! , Л) = Р(п,Х)гУ(п,Х) (2.3)

в окрестности особой точки Л = Ао. Здесь функции Р(п,Х) и Р-1(п,Х) - аналитические в окрестности А0, главные миноры матрицы Р(п,Х) удовлетворяют условиям (3.8) ниже, а матрица Z имеет диагональный вид (3.7).

Если нам удается привести систему к указанному виду, тогда коэффициенты асимптотического ряда эффективно вычисляются. По известным коэффициентам ряда легко строятся законы сохранения и симметрии. Обсудим на примере уравнения

у(п +!,Х) = ¡(п,и(п),Х)у(п,Х), ! = ( >! ^^ (2.4)

один из способов приведения исходной системы к виду (2.3). Легко видеть, что уравнение (2.4) имеет единственную особую точку А = то. Первоначально представим £ в виде произведения f (п,и(п),Х) = а(п,и(п), X)Zft(п,и(п),Х) трех матриц - диагональной Z и треугольных а и [3:

-=(а". «(0,))=(0 л0.)=(! "(1л-1 )•

Особо отметим, что а(п,Х) и [3(п,Х) являются аналитическими и невырожденными в окрестности А = то. Далее после замены гф = уравнение (2.4) примет требуемый вид

ф(п + !,А) = Р (п,и(п),Х)гф(п,Х).

Здесь Р определяется по формуле Р (п, и(п), X) = ¡3 (п +!,и(п + !), Х)а(п, и(п), X), а именно

, , ( ! — и(п + !)А°2 -и(п)и(п + !)А°2 \ Р (П,и(п),Х)=[^ \о. > ^(П) )•

При этом Р(А) аналитична на бесконечности и главные миноры матрицы Р(то) отличны от нуля.

3. Асимптотическая дилгонАлиздция дискретного оператора в окрестности особой точки

Предположим, что /(п,и,Х) имеет полюс при А = А0. Тогда / раскладывается в ряд Лорана в окрестности этой точки:

/ (п, и, Х) = (Х - Хо)-к Ып) + (Х - Ао Г^^ф) + ••• ,к > !• (3.1)

Цель этого раздела состоит в обсуждении достаточных условий «диагонализуемости» уравнения (2.1) в окрестности точки Л = Ао. Согласно работе [6] уравнение (2.1) диа-гонализуется, если существуют формальные ряды

R(n, X) = R(o) + R(i) (X - Xa) + Rm (X - Xq)2 + h(n, X) = h(Q) + h/i)(X - Xq) + h/2)(X - Xq)2 +

(3.2)

(3.3)

с матричными коэффициентами Е СмхМ, где Н^) V] диагональная (блочно-

диагональная) матрица, такие, что формальная замена зависимой переменной у = Р<£> приводит уравнение (2.1) к виду

ух = кг^. (3.4)

Здесь Z = (X — Х0)л, ё, - диагональная матрица с целыми элементами. Предполагается, что det R(0) = 0, det Н(0) = 0. Из формул (3.2)-(3.4) следует, что уравнение (2.1) обладает формальным решением, заданным следующим асимптотическим разложением

y(n, X) = R(n, A)eE—0logh(s'x)Zn

c «амплитудой» А = R(n, X) и «фазой» ф = nlog Z + YTsZn0 l°g h(s, X).

Предположим, что потенциал f (n,X) можно представить в следующем виде

f (n, u(n), X) = a(n, u(n), X)Z¡3(n, u(n), X),

где a(n,X) и @(n,X) аналитические и невырожденные в окрестности А = Ао, Z нальная матрица вида

Z

( (X - ХоУ1 Ei 0

0 (А - Xq)72Е2

0 0

\

(3.5)

(3.6) диаго-

(3.7)

V 0 0 ... (А — Ао)™Ем /

где Ез единичные матрицы размера ej х ej, показатели степени попарно различны 71 <72 < ... < ^(и. Положим Р(п, и,Х) = [3(п + 1,и(п + 1), Х)а(п, и(п), X) и допустим, что главные миноры матрицы Р(п,и, Х0) удовлетворяют следующим условиям:

det Р(n,u, Xq) = 0 for j = ei,ei + e2, ei + e2 + e3,

,N.

(3.8)

Теорема 1. Предположим, что X = А0 особая точка и потенциал f (п,и(п), X) удовлетворяет условиям (3.6), (3.8) в окрестности А0 ив области изменения и(п). Тогда существуют формальные ряды, "диагонализующие" уравнение (2.1) в окрестности X = А0, т.е. формальная замена у = Яр сводит (2.1) к виду (3.4), где к имеет следующую блочно-диагональную структуру

/Ни 0 ... 0 \

0 Н22 ... 0

h

(3.9)

\ 0 0 ... hrr /

Здесь hjj квадратные матрицы размера ej х ej. Коэффициенты R(i) и h^) зависят от, конечного, зависящего от г набора переменных из бесконечного множества {u(k)}k=-(X(X).

Доказательство теоремы 1 приведено в работе [7]. Необходимо указать, что в процессе доказательства теоремы строится формальный ряд Т = ßR, удовлетворяющий уравнению

Dn(T )h = Р (п,и(п),Х)Т, Т = ZTZ-1. (3.10)

Следствие 1. Линейное уравнение (2.1), переписанное в следующей специальной форме

ф(п + 1, X) = Р(п, и(п),Х)гф(п, X) сводится к блочно-диагональному виду

р(п + 1, X) = h(n, X)Zip(n, X)

преобразованием ф(п,Х) = Т(п,Х)^(п,Х), если Р = Оп(@)а, / = aZfl и выполняются условия (3.6), (3.8).

4. Асимптотическая дилгонАлизлция оператора Лакса

и законы сохранения

Рассмотрим динамическую систему

Р(0т0пу,Вту,Впу,у) = 0, (4.1)

где искомый объект - векторнозначная функция V = у(п,т) с координатами (п,т), ] = !,..., N, зависящими от целых п, т. Операторы сдвига Ип и Ит действуют по правилам Опу(п,т) = у(п + !,т) и Оту(п,т) = у(п,т + !). Мы предполагаем, что (4.1) является условием совместности линейных уравнений

у(п +!,т) = Р(п,т, [ь], Х)гу(п,т), (4 2)

у(п,т + !) = К(п,т, [и], Х)у(п,т). ( . )

Выражение [и] указывает на то, что функции Р и К зависят от переменной V и конечного числа ее сдвигов И^ь, О^у. Введем дискретные операторы Ь = D^оlPZ и М = О-^Я. Тогда условие совместности системы (4.2) можно записать в следующем виде:

[Ь, М] = 0, где [Ь, М] = ЬМ - МЬ. (4.3)

Отметим, что первое уравнение (4.2) имеет вид (2.3). Предположим, что Р(п,т, [и], А) удовлетворяет условиям теоремы 1 (Р аналитична в окрестности А = А0 для любых целых п, т и для всех значений и = [и] из некоторой области, а главные миноры (3.8) отличны от нуля в этой области). Мы также предполагаем, что функция К(п, т, [и], А) мероморфна в окрестности точки А = \0, когда и принимает значения в рассматриваемой области.

Из теоремы 1 следует, что дискретный оператор Ь = П"1PZ сводится к квазидиагональной форме Ь0 = Dо1hZ преобразованием Ь ^ Т-1ЬТ = Ь0, где Т(п, X) = ^¿>0 -А0)г. Из (4.3) следует, что [Ь0, М0] = 0, где

М0 := Т-1 МТ. По построению и в силу предположения, сделанного выше, коэффициент Б в формуле М0 = О—.Б есть формальный ряд вида

^ = (А - \0)к Е7=0 ^ - \0)-\

Теорема 2. Коэффициенты Бг ряда Б имеют ту же блочно-диагональную структуру, что и матрица h.

В силу блочно-диагональной структуры Б коммутирует с Z, и мы находим, что

^ = От(К)8. (4.4)

Переходя к блочному представлению Б = (Б^}, h = (^} в равенстве (4.4), получаем Оп(Зц)^г = Отфц)8ц. Теперь ясно, что уравнение

(Ип - !)\ogdet Бц = (Бт - !)\ogdet г = !, 2,..., N0 (4.5)

генерирует бесконечную последовательность законов сохранения для уравнения (4.1). Так как функция Б = Л = Бц не равна нулю тождественно, логарифмы в (4.5) корректно определены.

5. Законы сохранения двухкомпонЕнтного дискретного

потЕнциировднного уравнения Кортевега-де Фриза

Рассмотрим двухкомпонентное дискретное потенциированное уравнение Кортевега-де Фриза (еёрКёУ)

(^п,т '^п+1,т+1)(^п+1,т Vп,т+1) $ & , (5 1)

)(ип+1,т - ип,т+1) = б2 - а2. ' ;

(5.2)

В данном параграфе мы опишем бесконечную последовательность законов сохранения и построим высшие симметрии для системы (5.1) при помощи метода асимптотической диагонализации операторов Лакса. Пара Лакса для (5.1) построена в работе [13] и задается системой уравнений

Уп+1,т — /Уп ,m, Уп,т+1 — 9 Уn,m, где потенциалы и записываются следующим образом

( 0 —ип+1,т 0

ип+1,т 0 1

0 А ип+1,т11п,т 0 1)п

\ А ип,тип+1,т 0

0 а2 —

\2 г2 л_1

а — о — А — ип,тип,т+1

1 0

0 ип,т 0

ип,т+1 0 1 \

0 1 0

А ип,т+1 ип,т 0

0 0 /

Представим потенциалы / и д в виде / = РП и д = СП соответсвенно, где ( —ип+1_т 0

Р

\

Ц"п,+ 1,т 0

— А 1 — ип+1,тУ, 0

V п+1 п,т 0

-А"1

С =

( V

— /

0

— А"1 + а2 — б2 — ип>т+1 ьп 0

Уп+1

0 0

1 0

,п,п 0

0

1 0

/

1 0

0

0 1 0

иг,

\

П

—А 1 + о2 — Ь2 — ип,т+1ип,т

/0100 \ 10 0 0 0001 0010

Пусть Я - кольцо матриц Хпхп, удовлетворяющих условию аХа-1 = X, где а = diag(1, —1,1, —1). Легко видеть, что матрицы Р и С принадлежат группе 0 обратимых элементов кольца Я.

Преобразуем систему (5.2) при помощи замены = уп>т к системе уравнений

(5.3)

с новыми потенциалами

р _ П-(п+т+1)рПп+т+1 _ I Рп+1,т

= = —\-11 — Рп+1,тР',

(

п,т уп,т

Рп

)

С_(п+т+1) п+т+1_

рп,т+1

—А 11 + (а2 — 82)1 — рп,т+1 р

п,т 1-"п,т

Рп

Здесь I обозначает единичный блок diag(1,1), переменная Рп,т определена следующей формулой:

0

Рп,т

^-(п+т)( ип,т 0 \ V 0 ^п,т )

^ п+т

(5.4)

где

Е =

01 10

Теперь систему (5.1) можно рассматривать как условие коммутирования двух дискретных операторов С = Б-1 Р и М = Б^С.

Приведем уравнение п+\,т = к специальному виду (2.3). Для этого представим

потенциал Р в виде произведения Р = aZP трех матриц - блочной нижнетреугольной а, блочно-диагольной Z и блочной верхнетреугольной @:

~ ( 1 0 \

у Яп,т + Л Рп+1,т Рп+1,т )

% = I ^ \-1

0_п,т + А Рп+1,т Рп+1,\ ( 1 0 А Ъ = ( -Рп+1,ш I \

^ 0 а- ч ), р = ^ 0 I )

Тогда замена ф = приводит первое уравнение (5.3) к виду гфп+1,т = Рфп,т, где Р = Пп+1(/3)а, а именно

р(Д) = ( Яп,т - 'Рп+1,т —Pn+1,m \ + д-1 ( Рп+1,т 0 \ \ йп,т -Рп+1,т ) V Рп+1,т 0 /

Из последней записи видно, что Р Е'Я,. При этом миноры

¿еЬ(дп,т - рп+1,т), ¿еЬР(то) (5.5)

матрицы Р(то) отличны от нуля. Напомним, что (Ц.п,т и рп+1,т есть 2 х 2 блоки. В силу теоремы 1 существуют формальные ряды

Т = Т0 + Т.\-1 + ••• , h = ¡10 + ¡г.Х-1 + •••

такие, что оператор Ь0 := Т-1ЬТ, где Ь = 1PZ будет диагональным оператором вида Ь0 = 1hZ. Найдем ряды Т и 1г из уравнения

вп(т% = РГ, Т = гтг-1. (5.6)

Поскольку

/ Г>1,1 Л±1,2

"1 = ( ^1,1 ^1,2 | \ ^ 1Г^2,1 Т2,2 )

то справедливы равенства

ХТ.,2 = Т 1,2, А ±2,1 = Т2,1, Тг,г = Ti,i, % =1, 2. Подставляя сюда формальные ряды

т = т0 + \~1т. + ••• , Т = Т0 + \~.Т. + ••• ,

получим, что

ХТ0 ,1,2 + Т. ,1,2 ,1,2 + • • • = т 0 ,1,2 + \-1т 1 ,1,2 + 2Т 2 ,1,2 + • • •

А .^0,2,1 + А 2Т.,2,1 + Л 3 Т2,2,1 + • • • = Т0,2,1 + А 1Т 1,2,1 + А 2Т 2,2,1 +

^0,г,г + А 1Тщ + • • • = Т0,',

Т0,г,г + А 1Тщ + • • • = Т0,г,г + А 1Тщ + • • • , г = !, 2.

Из последних равенств следует, что Т0.,2 = 0, Т0,2,. = 0, т.е. матрицы Т и Т являются блочными нижнетреугольной и верхнетреугольной соответсвенно, и справедливы равенства

Тр+1,1,2 = Тр,1,2, Тр,2,1 = Tp+l,2,l, Тр,г,г = Tp,i,i, % = 1, 2. (5.7)

Вернемся к уравнению (5.6) и перепишем его в виде

Вп(% + А"1Т1 + • • • )(Ы + % + ••• ) = (Р0 + А"1Р.)(Т0 + .Т1 + •••) (5.8)

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

Вп(То)1го = РаТа, (5.9)

ВП(Т1 )ка + Оп{Та)К1 - Р0Т1 = Р^а, (5.10)

к-1

Пп(Тк)%а + Оп(Т0)кк - Р0Тк = РгРк-г - ^ Вп(Тк-3)~к3, к> 2. (5.11)

3 = 1

Уравнение (5.9) есть задача Гаусса о разложении матрицы Ра в произведение трех матриц в группе 0 - блочной нижнетреугольной Оп{Та), блочно-диагональной ка и блочной

-1

верхнетреугольной Та . Разрешимость этой задачи гарантируется условием регулярности (5.5). Единсвенность решения этой задачи обеспечиваем, задав диагональные блоки матрицы Та равными единичным матрицам размера 2 х 2. Поскольку задача Гаусса решается в группе 0, то мы получаем, что матрицы Та, ка и Та принадлежат группе 0.

Далее, полагаем, что диагональные блоки матриц Тк и Тк для всех к > 0 являются

нулевыми. Разложим каждую из матриц Тк и Тк в сумму нижней блочно-треугольной и верхней блочно-треугольной матриц с нулевыми диагональными блоками:

Тк = ТкЬ + Тки, Тк = Ткь + Тки. (5.12)

Матрицы Т\и и Тц находятся легко, действительно, используя (5.7), имеем Т\д,2 = Та,1,2, Т\,2,1 = Та,2,1. Т.е. эти элементы уже найдены на предыдущем шаге. Для нахождения неизвестных Т\^ и Т1и воспользуемся уравнением (5.10), записанным в следующем виде

кф-1 + ВП(Т0-1Т1Ь) - каТа Тшк-1 = Нь

~ ~ ~ ~ ~ ~ л -1

где правая часть Н1 = Оп(Т0 )Р1Так- - Оп(Т- Т1и) + каТа Т 1ьк- содержит уже известные матрицы. При этом Н1 Е Я. Для того чтобы найти неизвестные к1, Т^, Тщ, нужно разложить Н1 в сумму трех слагаемых: блочно-диагональной к1к-1, блочной нижнетреугольной Оп^-1^^ и блочной верхнетреугольной -каТ° Тщк-1. Ясно, что так как задача решается в кольце Я, то и искомые слагаемые будут принадлежать Я, а значит

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

ккк- + °П{Т0 ТкЬ) - каТа Ткик- = Нк,

где член Нк содержит слагаемые, найденные на предыдущем шаге. Таким образом, доказано, что все коэффициенты ряда Т принадлежат кольцу Я.

Выпишем первые элементы формальных рядов Т и к в явном виде:

(

0 I

Т — ( р„_х,т т ) +

Рп + 1,т Рп — 1 ,т

0

I / Рп + 1 ,т (Т?п + 2 ,т Т?п,т

+ I 1^ + 1 т Г> I Л +

Рп+1 ,т

(Рп + 1,т-Рп-1 ,т)2(Рп,т-Рп-2,т)

1

h

(

Рп

Рп+2,',

0

Рп+2,п

+

Рп + 1,т

_1_

Рп + 1 ,т Рп -1,т 0

)

+

0

Рп + 2 ,тРп + 3 ,т Рп + 2 ,тРп + 1 ,т +Рп,тРп +1 ,гг 1,т(Рп + 3,тоРп + 1,т)(Рп + 2,тоРп,т)2

Г)2

Рп +

)

А" 1 + ••

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

Оператор М = О^С диагонализуется следующим образом: М0 = О^Б, где Б = (Т-1Т. Выпишем первые слагаемые ряда Б в явном виде:

5

(

Рп,т + 1 Рп + 1,т 0

+

(

0

62о<Т2+Р„,т(Рп,т + 1оРп + 1,т) Рп+1 ,т

)

+

1

Рп + 1 ,т Рп-1,Т1 0

0

ё2 оО +Рп,т + 1Рп,т+Рп + 2 ,т,Рп + 1 ,т —Рп + 1 ,тРп,п Рп + 1 -т('Рп,'т0Р>п + 2,т )

А"1 + ••

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

(Ип - !)\ОЕ-

!

(Оп - !)

Рп+1,т Рп,т+1 Рп+1,т Рп,т+1

(О т !)1оё(Р п,т Рп+2,т),

!

(

]п+1

- Рп- 1,т)($2 - Я2)

(От - !)

(Рп+1 - Рп )(Рп,т - Рп+2,т)

(Оп - !)

Ь2 - О2 + рп

( п,т+1 Рп+1,т ) + Рп+2,тРп+1,т Рп+1,т ( п,т Рп+2,т )(62 - а2 + Рп,т(Рп,т+1 - Рп+1,ш))\

Рп+2,тРп+3 ,т Рп+2,тРп+1,т + Рп+1,тРп

(Вт - !)

Рп+1,тРп+2,т ( п,т Рп+2,т )(Рп+1,т - Рп+3,т)'

Переходя к исходным переменным и и V, получим

!

(Бп -

(щ,1 - и.,0)(Ь0,1 - г>1,0)

(Бт - !) 1о§(м - щ,0)(у - Ь2,0),

(Пп - !)

П0,1 - и.,0

+

0,1 - 1,0

(р2 - Ч2)(и-1,0 - Щ,0) (р2 - Я2)(у-1,0 - ^1,0)

(V - У2,0)(и-1,0 - Щ,0) (и - Ч2,0)(у-1,0 - ^1,0)

(От - !)

(и-1,0 - Щ,0)(у - У2,0) (V-1,0 - У1,0)(и - 42,0)

(Оп - !)

(р2 - д2 + ии + У.,042,0) (р2 - д2 + VII + и.,0У2,0)

щ,0р (р2 - д2 + и и)

и1,0(г(р2 - д2 + ьр)

(От - !)

Уз,0

+

!

+

Щ,0

+

!

У1,0Р&1,0 и2,0&1,0 Щ,0&Р1,0 Щ,0Р1,0

Здесь используются обозначения р = и - и2,0, а = V - у2,0, р = и0,. - и.,0, и = - ь.^.

Теперь перейдем к построению симметрий системы (5.1). Схема построения высших симметрий дискретной динамической системы при помощи алгоритма диагонализации подробно изложена в работе [6].

52-а2

Вернемся к первому из линейных уравнений, образующих пару Лакса (5.3) для системы (5.1), и перепишем его, используя дискретный оператор С = О^Ё, в виде

р = Ср.

Как было показано выше, замена переменных ф = [р приводит это уравнение к специальному виду

ф = Ьф

с оператором Ь = Do1РZ. Далее были найдены формальные ряды

те те

Т = ^Т(к)Л-к, к = ^ к(к)Л-к, (5.13)

к=а к=а

«диагонализующие» указанное уравнение в окрестности Л = то, т.е. такие, что оператор Ьа := Т-1ЬТ является оператором с (блочно-) диагональными коэффициентами Ьа = Doo1~кZ.

Следуя работе [6], построим формальный ряд

Ва = ^ (Ва)(к)Л к,

к=-М

с коэффициентами (Ва)(к), имеющими ту же блочно-диагональную структуру, что и элементы ряда к, и не зависящими от п, и такой, что [Ьа, Ва] = 0. В работе [6] доказано, что в этом случае формальный ряд В' = ^Гк=-М В(к)Л-к, заданный формулой В' = ТВаТ-1,

коммутирует с оператором Ь. Тогда формальный ряд В = ^ск=~М В(к)Л-к, заданный формулой

В = [-1ТВа([-1Т )-1, (5.14)

коммутирует с оператором С.

Теорема 3. Коэффициенты В(к) ряда (5.14) лежат в кольце Я, т.е. для любого к удовлетворяют соотношению аВ(к)а-1 = В(к), где а = diag(1, -1,1, -1).

Выше было показано, что коэффициенты Т(к) ряда Т для любого к лежат в кольце Я,

[ Е К - по построению. Таким образом, доказательство Теоремы 3 следует из формулы (5.14).

Возьмем Ва в виде

Ва = В(-м)Л-М, В(-М) = ^(1,1, -1, -1), (5.15)

тогда прямым вычислением доказывается, что для любого к коэффициенты В(к) ряда В будут удовлетворять условию

к _ т~к

22 = -В11.

Здесь через Вк3 обозначены 2 х 2 блоки матрицы В(к):

эк ~~к

В(к) =| В12 | . (5.16)

В1 1 В1 2 В2 1 В2 2

(к) = \ Бк

21 В22

Далее разложим ряд В в сумму В = А +(В -А), где А = ^ ^ А(к)Лк = ^ В(к)Л-к. Тогда

[С, В] = [С, А] + [С,В - А] = 0.

Потенциал Ё первого из уравнений пары Лакса (5.3) является рациональной функцией

вида Ё = Л 1Ё1 + Ёа, где

Ё =( 0 0 ^ Ёо =( -Р'п+1,™ } ^

У 1 0 ) \ Рп+1,тРп,т Рп ,т /

Выясним какой вид имеет коммутатор [С, А]:

[С, А] = [Б-1Г, А]

С другой стороны,

[С, А] = -[С, в -А]

м

м

[Б —.(Х~1Р1 + р0),^2 А(к)Хк ] = [В-1Р.,А(1)] + акХ

к=1

к=1

[Б —.(Х-1Р1 + р0),^В(к)Х-к ] = - [О-.р0,В(0)]^^2 а'кХ~

к=0

к=0

Из последних равенств следует, что

[с, А] = н,

где

к = -[в~1р0,В(0)] = [в-1Р.,А(1)]. Из второго равенства (5.17) следует, что

вп(в) = Р.А(1) - вп(А(1))Р. = ( а^^ ) .

(5.17)

(5.18)

Таким образом, Оп(В) является (блочной) нижнетреугольной матрицей. Используя этот факт, из первого равенства (5.17) получим, что Ип(К) имеет вид

Оп(В) = -р0В(0) + Вп(В(0))р0 = (

Д11 0

-рп+1 В22 + Рп В.. В22

)

(5.19)

Положим 1Р

<и п

К или <_Р0

- Вп(В). Это соотношение определяет дифференциально-разностное уравнение, действительно, левая часть имеет вид

(

_<Рп + 1,т

Лрп+1,т ~ + ~ <рп,т

< Рп,т + Рп+1,т

0

<Рп,т

<м

)

и правая в силу (5.18) и (5.19)

Вп(В)

(

-Вп (А.2)

IуРп+1,тА12 + Рп,тВп(^А12)^

А

12

Следовательно,

¿рп

А1 . А12.

(5.20)

Ниже мы приведем в явном виде две высших симметрии. Выберем затравочный ряд В0 в виде В0 = diag(1,1, -1, - 1)А. Теперь формальный ряд В будет иметь вид

В = р-1ТВ0Т-1 (3 = В(-1)Х + В(0) + В(1)Х~1 + ••• .

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

Рп + 1 ,т +р п — 1 ,т а _ i рп + 1,т—рп-1,'

^ = I 2рп-1>трп+1>1

Рп + 1 ,т Рп — 1 ,т

Рп + 1,т Рп — 1 ,т Рп + 1,т +р п — 1 ,т Рп + 1,т Рп — 1 ,т

)

к

0

1

2

Используя (5.20), выписываем симметрии

dPn,m 2

dt Pn+l,m Pn-l,m

переходя к исходным переменным u и v, получаем симметрии

du 2 d ; 2

n, m n, m

dt Vn+l,m Vn—l,m dt Un+i,m Un—i,m

системы (5.1).

Далее для того чтобы найти симметрию следующего порядка, выберем затравочный ряд В0 в виде В0 = diag(1,1, — 1, — 1)Л2. Теперь искомый формальный ряд будет иметь вид

в = г1тВоТ-1/з = В{-2)Л2 + В-)Л + •••.

Тогда формальный ряд А будет иметь вид А = А(2)Л2 + А(^Л, где А(2) = В(-2) = А,

Aw = B-D =( an а12 ) , (1) ( 1) \ a2i — an J

где

2(Pn+2 ,mPn+l,m - Pn , m Pn+1 , m + Pn — l,m'Pn,m - Pn—l,mPn—2,m)

an = -

ai2

(Pn+2 , m - 'Pn , m

)(p

n, m - Pn 2, m )(Pn+l , m - Pn—l,m)2

2(Pn+2 ,m Pn—2,m)

2l

СPn+2 , m - - Pn , m )(Pn , m - Pn -2, m )(Pn+l , m - Pn— l,m)2' 2(Pn+l,mPn+2,m - Pn+l,mPn,m + Pn—l,mPn,m - Pn—l,mPn—2,m)

(.Pn+2 , m - Pn , m )(Pn+l , m - - Pn l, m )2(Pn , m - Pn—2,m)

Используя (5.20), получаем

dpn , m 2(Pn+2 , m - Pn—2,m)

dt (Pn+2 , m - - Pn , m )(Pn , m - - Pn 2, m )(Pn+l , m - Pn—l,m)2

или в исходных переменных, записываем симметрии

dun,m 2(un+2,m - un—2,m)

dt (un+2,m un,m)(un,m un- — 2,m)( Vn+l,m Un—l,m)

d n, m 2(Pn+2,m - n- 2 , m )

dt ( Vn+2,m Vn,m)( Vn,m Vn- -2,m)(un+l,m un—l,m)

системы (5.1)

СПИСОК ЛИТЕРАТУРЫ

1. W. Wasow Asymptotic Expansions for Ordinary Differential Equations Dover: Dover Books on Advanced Mathematics, 1987. 374 p.

2. Захаров В.Е., Манаков С.В., Новиков С.П., Питаевский Л.П. Теория солитонов: метод обратной задачи М.: Наука, 1980. 320 с.

3. Тахтаджян Л.А., Фаддеев Л.Д. Гамильтонов подход в теории солитонов М.: Наука, 1986. 528 с.

4. V.G. Drinfeld, V.V. Sokolov Lie algebras and equations of Korteweg-de Vries type // Journal of Soviet Mathematics. 1985. V. 30. No. 2. P. 1975-2036.

5. Хабибуллин И.Т. Дискретная система Захарова-Шабата и интегрируемые уравнения // Зап. научн. сем. ЛОМИ. 1985. Т. 146. С. 137-146. Англоязычная версия: Journal of Soviet Mathematics. 1988. V. 40. No. 1. P. 108-115.

6. Хабибуллин И.Т., Янгубаева М.В. Формальная диагонализация дискретного оператора Лакса и законы сохранения и симметрии динамических систем // ТМФ. 2013. Т. 177. № 3. С. 441467.

7. I.T. Habibullin, M.N. Poptsova Asymptotic diagonalization of the Discrete Lax pair around singularities and conservation laws for dynamical systems J. Phys. A: Math. Theor. 48:11 (2015) 115203 (37pp)

8. R.N. Garifullin, R.I. Yamilov Integrable discrete nonautonomous quad-equations as Baecklund autotransformations for known Volterra and Toda type semidiscrete equations // Journal of Physics: Conference Series 621 (2015) 012005 (18pp).

9. R.N. Garifullin, I.T. Habibullin, R.I. Yamilov Peculiar symmetry structure of some known discrete nonautonomous equations //J. Phys. A: Math. Theor. 48 (2015) 235201 (27pp).

10. A.V. Mikhailov Darboux transformations and symmetries of partial difference equations // Geometric Structures in Integrable Systems, International Workshop, Moscow (Russia), October 30 - November 2, 2012.

11. Da-jun Zhang, Jun-wei Cheng, Ying-ying Sun Deriving conservation laws for ABS lattice equations from Lax pairs //J. Phys. A: Math. Theor. 2013. V. 46. No. 26. id 265202.

12. Jun-wei Cheng, D-j Zhang Conservation laws of some lattice equations // Front. Math. China. 2013. V. 8. No.5. P. 1001-1016.

13. T. Bridgman, W. Hereman, G.R.W. Quispel and P.H. van der Kamp Symbolic computation of Lax Pairs of partial difference equations using consistency around the cube // Foundations of Computational Mathematics. 2013. V. 13. No. 4. P. 517-544.

14. Fu Wei, Zhang Da-Jun, Zhou Ru-Guang A Class of Two-Component Adler-Bobenko-Suris Lattice Equations // Chinese Phisics Letters. 2014. V. 31. No. 9. id 090202.

15. R. Hirota and S. Tsujimoto Conserved quantities of a class of nonlinear difference-difference equations // J. Phys. Soc. Jpn. 1995. V. 64. No. 9. P. 3125-3127.

16. F. Nijhoff and H. Capel The discrete Korteweg-de Vries equation // Acta Applicandae Mathematica. 1995. V. 39. No. 1-3. P. 133-158.

17. V.E. Adler, A.I. Bobenko, Yu.B. Suris Classification of integrable equations on quad-graphs. The consistency approach // Commun. Math. Phys. 2003. V. 233. No. 3. P. 513-543. arXiv: nlin.SI/0202024.

18. A.G. Rasin, J. Schiff Infinitely many conservation laws for the discrete KdV equation //J. Phys. A: Math. Theor. 2009. V. 42. No. 17. 175205 (16 pp.)

19. D. Levi, M. Petrera Continuous symmetries of the lattice potential KdV equation //J. Phys. A: Math. Theor. 2007. V. 40. No. 15. P. 4141-4159.

20. D. Levi, M. Petrera, C. Scimiterna The lattice Schwarzian KdV equation and its symmetries // J. Phys. A: Math. Theor. 2007. V. 40. No. 42. P. 12753-12761.

21. O.G. Rasin and P.E. Hydon Symmetries of Integrable Difference Equations on the Quad-Graph // Stud. Appl. Math. 2007. V. 119. No. 3. P. 253—269.

22. A. Tongas, D. Tsoubelis, P. Xenitidis Affine linear and D4 symmetric lattice equations: symmetry analysis and reductions //J. Phys. A: Math. Theor. 2007. V. 40. No. 44. P. 13353-13384.

Мария Николаевна Попцова, Институт математики с ВЦ УНЦ РАН, ул. Чернышевского, 112, 450008, г. Уфа, Россия E-mail: mnpoptsova@gmail.com

Исмагил Талгатович Хабибуллин, Институт математики с ВЦ УНЦ РАН, ул. Чернышевского, 112, 450008, г. Уфа, Россия E-mail: habibullinismagil@gmail.com