Эквивалентность линейных дифференциальных уравнений с частными производными и преобразования Эйлера-Дарбу Текст научной статьи по специальности «Математика»

Вычислительные технологии

Том 12, № 4, 2007

ЭКВИВАЛЕНТНОСТЬ ЛИНЕЙНЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ С ЧАСТНЫМИ ПРОИЗВОДНЫМИ И ПРЕОБРАЗОВАНИЯ ЭЙЛЕРА-ДАРБУ*

О. В. КАпцов Институт вычислительного моделирования СО РАН,

Красноярск, Россия e-mail: kaptsov@icm.krasn.ru

An equivalence relation on the class of linear partial differential equations is introduced. The equivalence is defined by means of linear differential operators. We prove that the Euler—Darboux transformations generate the equivalence relation on a certain class of equations. A superposition formula was found for the Euler—Darboux transformations.

The Darboux and Schrodinger equations were considered as applications of the proposed method.

Введение

Задача сведения заданного уравнения к уже изученному представляет несомненный интерес для приложений. Локальные преобразования изучаются в групповом анализе дифференциальных уравнений [1-3]. Два уравнения считаются эквивалентными, если существует точечное преобразование, переводящее одно уравнение в другое. Лучше всего изучены вопросы эквивалентности для линейных обыкновенных дифференциальных уравнений [3]. В работе предлагается дополнить класс преобразований линейных уравнений с частными производными.

Пусть имеется два линейных уравнения с частными производными одного и того же порядка из некоторого класса E:

Ь\п = 0, (0.1)

L2v = 0. (0.2)

Предположим, что существует линейный дифференциальный оператор L такой, что преобразование

v = Lu (0.3)

* Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (грант № 07-01-00489).

© Институт вычислительных технологий Сибирского отделения Российской академии наук, 2007.

переводит решения уравнения (0.1) в решения (0.2). Тогда преобразование (0.3) называют дифференциальной подстановкой. Если найдется еще дифференциальная подстановка и = Му, которая переводит решения уравнения (0.2) в решения (0.1), то уравнения (0.1) и (0.2) будем называть эквивалентными, а операторы С и М — операторами эквивалентности, обратными друг другу.

Таким образом, получается категория EQ, объектами которой являются линейные дифференциальные уравнения, морфизмами — дифференциальные подстановки, а на классе уравнений Е возникает отношение эквивалентности. Из указанных выше условий эквивалентности следуют два соотношения:

Ь2(Си)1Ьуа=о = 0, Ь1(Му)1ь2,и=о = 0.

Классический пример подобного отношения эквивалентности связан с преобразованиями Лапласа [1,4] гиперболических уравнений

иху + а(х, у)их + Ь(х, у)иу + с(х, у)и = 0.

Следует отметить, что в случае совпадения уравнений (0.1) и (0.2) оператор С является симметрией уравнения (0.1) [5]. Если уравнения (0.1) и (0.2) различны, а операторы С, М совпадают, то С будем называть сплетающим оператором.

В данной работе рассматривается специальный класс преобразований линейных уравнений с частными производными. Эти преобразования называются преобразованиями Эйлера—Дарбу (или Э-Д-преобразованиями), поскольку Эйлер первый ввел их в своих работах [6] по построению общих решений уравнений второго порядка, а Дарбу [7] применял к интегрированию обыкновенных дифференциальных уравнений. Доказано, что Э-Д-преобразование имеет обратное, на некотором классе Е2,м линейных уравнений. Указаны формулы суперпозиции Э-Д-преобразований и явный вид преобразованных уравнений. Приведены примеры приложений преобразований Э-Д к уравнениям второго порядка. В частности, описано множество уравнений Эйлера—Пуассона— Дарбу, эквивалентных волновому уравнению с постоянными коэффициентами.

1. Преобразования Эйлера—Дарбу

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

Рассмотрим линейное уравнение с частными производными:

Ьи = Аи + Ви = 0. (1.1)

Предположим, что оператор А — дифференциальный оператор по одной переменной х:

к

А = Х! аг(х)дХ, (1.2)

г=0

а В — дифференциальный оператор по переменным у\,..., уп вида

м

В = £ Ь„(у)а?, (1.3)

Н>о

дг д|а|

где у = (у1,..., уп), дХ = -гг—, да = ^, а = (аь ..., ап) — целочисленный мульти-

дхг у дуГ ,...,даг

индекс. Функции аг(х), Ьа(у) считаются гладкими в соответствующих областях. Класс уравнений вида (1.1) обозначим через Ек,м.

Если Л(х),д(у) — решения уравнений

АЛ = сЛ, (1.4)

Вд + сд = 0, с € Д,

то функция и1 = Лд удовлетворяет уравнению (1.1). Функция и1, порождает преобразование уравнения (1.1).

Теорема 1. Класс уравнений Ек,м обладает следующими свойствами:

1. Если 7 — гладкая функция вида 7 = р(х)д(у) = 0, то преобразование

и —> у = и/7

переводит решения уравнения (1.1) в решения уравнения

Ьу = уЬ(7)/7 + А1у + В1у = 0,

где

к м

аКЫ'-'у ■

ai = X!a1(x)dx bi = X! ba(y)da

i=1 |a|>1

При y = u1 = 0 уравнение Lv = 0 имеет вид

L1v = A1v + B1v = 0. (1.5)

2. Преобразование v —> w = vx переводит решения уравнения (1.5) в решения уравне-

ния

к м

L2w = ^(5x(a1)5X-1w + a^w) + ^ ^d^w = 0. (1.6)

i=1 |а|>1

Доказательство. Для проверки первого свойства заметим сначала, что верны равенства

Lu = L(yv) = vLy + Av + Bv = 0, (1.7)

где

км

Av = ^2 a(x, y)dXv, Bv = ^ 6a(y, y)d^v.

i=1 |a|>1

Коэффициенты a могут зависеть только от x, y и производных от y по x, а коэффициенты ba могут зависеть только от y, y и ее производных по y1,..., yn. Функция Y и ее производные могут входить в коэффициенты a, Ьа лишь линейным образом. Разделив (1.7) на Y, получаем уравнение

Lv = vL(Y) + A1v + B1v = 0,

Y

где операторы А1 , В1 имеют вид

к М

А1 = X!а(х,р(х))дХ, В1 = X! Ь«(у,^(у))дУ'.

1=1 |а|>1

При 7 = и1 приходим к уравнению (1.5). Для доказательства второго свойства достаточно продифференцировать (1.5) по х и ввести новую функцию ,ш = дху. В результате получается уравнение (1.6). Заметим, что все уравнения Ьи = 0, Ь1у = 0, Ь2и> = 0 принадлежат одному классу Ек,м. ^

Следствие. Пусть Ъ — нетривиальное решение уравнения (1.4), г = 0 — гладкая функция от х. Тогда преобразование

Ъ!

г = (их - — и)/г (1.8)

переводит решение уравнения (1.1) в решения уравнения того же класса Ек,м. Действительно, согласно теореме 1 преобразование

г = р(х)д(у)д^ (0)

сохраняет класс уравнений (1.1). Здесь р, д — произвольные гладкие функции, и1 — решение уравнения (1.1), полученное разделением переменных и1 = Ъ(х)д(у). Если положить д = д, р = Ъ/г, то из (1.9) получим (1.8).

Преобразование (1.8) будем называть преобразованием Эйлера—Дарбу, а соответствующий линейный оператор

Ъ 1

дх Ъ

гЪ

— оператором Эйлера—Дарбу первого порядка. В современной литературе преобразование (1.8) и его многомерные аналоги принято называть преобразованием Дарбу [8-11].

Если известно решение Ъ уравнения (1.4) при некотором с, то это позволяет построить преобразование Эйлера—Дарбу. Оказывается, что если известны решения Ъ1,..., Ък уравнения (1.4) при различных с1,..., ск, то легко построить оператор порядка к и соответствующее преобразование, действующее на Ек,м. Для построения таких операторов применим конструкцию, близкую к той, что используется в теории факторизации линейных дифференциальных операторов [8, 12].

Обозначим через оператор Эйлера—Дарбу вида

Ъдх1.

Ъ

Пусть Ъ1,..., Ъм — гладкие, линейно независимые функции от х. Построим последовательность функций и операторов

Р1 = Ъ1, Р2 = СР1Ъ2, ..., рм = £РМ _Г ...^Р1 Ъм, (1.10)

М1 = £рГ , М2 = £р2 МЪ ..., ММ = £рм ММ-1.

Как следует из построения операторов Мк, функции Ъ1,..., Ък удовлетворяют дифференциальному уравнению порядка к

Мк Ъ = 0. (1.11)

Значит, они образуют базис решений уравнения (1.11). Следовательно, действие оператора Мк на произвольную функцию представляется в виде [8, 12]

ка пк пк— 1 (Ъ'1?---?Ъ'к, и) р.4

^МкИ = дх и + ак—ідх и + ... + аои = Т1,/,-----;—^— • (1-12)

^ (Ді, •••, Дк)

Если взять в качестве Д1,...,Дк решения уравнения (1.4), соответствующие различным константам с1, ...,Ск, то получим утверждение.

Лемма 1. Преобразование

W (Ді,...,Дк,и) п 10,

г = Мк и = ----------—— (1.13)

W (Ді,...,Дк)

переводит решение уравнения (1.1) в решение уравнения того же класса ЕК,м.

Операторы Мк отнесем к высшим операторам Эйлера—Дарбу, а преобразование (1.13) — к высшим преобразованиям Эйлера—Дарбу.

Остается открытым вопрос о том, существуют ли операторы порядка к, действующие на Ек,м и не представимые в суперпозиции операторов Эйлера—Дарбу первого порядка.

Заметим, что лемма 1 — аналог известной теоремы Крума [8, 13] из теории линейных обыкновенных уравнений второго порядка. Подобное утверждение имеется в теории сетей [14].

2. Преобразование уравнений класса Е2,м

В данном параграфе будет указан вид уравнений класса Е2,м, полученных действиями преобразования Эйлера—Дарбу, а также найдено обратное преобразование.

Рассмотрим уравнение

^ижж + Сих + Ни = Ви, (2.1)

где С, Н — гладкие функции от х, В — линейный оператор вида (1.3). Уравнения

вида (2.1) относятся к классу Е2,м.

Теорема 2. Преобразование Эйлера—Дарбу (1.8) переводит решения уравнения (2.1) в решения уравнения

+ Сігх + Ніг = Вг, (2.2)

где функции Сі,Ні задаются формулами

г

Сі = С + ^ + 2^—, (2.3)

г

Ні = н + (^г + Сг) + ^,(1п д) + 2^ (1п ду/. (2.4)

г

Здесь штрих означает производную по х, а функция і удовлетворяет уравнению

Пг// + Сі + (Н + с) і = 0, с Є Н. (2.5)

Доказательство. Для краткости обозначим через Аи левую часть (2.1), через Аіг — левую часть (2.2), а преобразование Эйлера—Дарбу запишем в виде

„ их + зи . .

г = Си = —--------, (2.6)

где в = — Л'/Л. В этих обозначениях уравнения (2.1), (2.2) выглядят так:

Аи = Ви, А^ = Вг. (2.7)

Второе уравнение в (2.7), с учетом (2.6), переписывается в форме

А^и = В £и.

Поскольку операторы В и £ коммутируют, то последнее уравнение можно преобразовать к виду

А1£и — £Аи = 0, (2.8)

в силу первого уравнения (2.7). Левая часть (2.8) — многочлен первой степени относи-

тельно мхх,мх,м. Коэффициенты при этих величинах должны быть равны нулю.

Собирая подобные члены при ихх,их, получаем формулы (2.3), (2.4). Используя полученные выражения (2.3), (2.4) и приравнивая к нулю коэффициент при и, приходим к уравнению

*У' + (^' — 2^в + С) в' — ^'в2 + Св' — И' = 0. (2.9)

Прямыми вычислениями убеждаемся в том, что (2.9) — это следствие уравнения

(2.5) и выражения в = — Л'/Л. □

Замечание. Доказанный результат является небольшим обобщением теоремы Эйлера [6]. Можно заметить, следуя [6], что (2.9) допускает первый интеграл

— ^в2 + Св — И = с. (2.10)

Тогда преобразование в = — Л'/Л приводит (2.10) к виду (2.5).

Вернемся к высшим преобразованиям Эйлера—Дарбу.

Теорема 3. Пусть ^1,... , — решения уравнений (2.5), соответствующие раз-

личным константам с1,... ,с&. Тогда преобразование (1.13) переводит решения уравнения (2.1) в решения уравнения

+ Ск + Икг Вг. (2.11)

При этом коэффициенты С&, И задаются формулами

С = С + к*1', И = И + кС' + к(к 2— 1) ^'' + ^' (1п Ж )' + 2^ (1п W)'', (2.12)

где Ж — вронскиан функций Л1,... , Л,.

Доказательство. Выражение для С& получается по индукции, последовательным применением формулы (2.3) при г = 1. Используя (2.4) и конструкцию (1.10) функций р1,...,р&, легко видеть, что индукционное построение коэффициентов И, приводит к выражениям

И, = И + кС' + к(к2— 1) ^ + ^'(1пР1.. .р*)' + 2^(1пР1.. .№)''. (2.13)

Необходимо найти произведение р1 ...рк. Так как, согласно (1.10) и (1.12), имеют место соотношения

Ж(ЛЬ . . . ,1,1+1)

Р*+1 = МЛ+1

Ж(ЛЬ . . . , Л*)

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

Л Ж(Л1,Л2) Ж(кЬ...,кк) ЖЛ ^

Р1.. .рй = к---- -----... ——-------------г = Ж(кь... , кк).

к Ж(кх,..., кк_х)

Остается подставить последнее выражение в (2.13) и получить вторую формулу в (2.12). □

Замечание. Поскольку функция к* (г = 1,...,к) удовлетворяет уравнению (2.5) с параметром с = с*, оператор М, зависит от с1,...,с,. Если некоторые параметры с;,... , с^+т стремятся к одной величине Ь, то в соответствии с методом Даламбера [15] решениями уравнения (1.11) будут функции к;, дък;, ...,5тк^. Чтобы лемма 1 и теорема 3 оставались справедливыми, необходимо в соответствующих формулах заменить к;+1,... , к;+т на дък,,... , дтк.

Перейдем теперь к построению преобразования Эйлера—Дарбу, переводящего решения (2.2) в решения (2.1). Такое преобразование назовем обратным к (1.8).

Теорема 4. Пусть задано преобразование Эйлера—Дарбу, переводящее решения уравнения (2.1) в решение уравнения (2.2). Тогда обратное преобразование Эйлера—

Дарбу

и = — — к1 Л (2.14)

г1 к1

переводит решения уравнения (2.2) в решения (2.1), если

= 1 Ь =1 ( )

Г1 г^, 1 кг ехр(/ С/^ж) 5

Доказательство. Согласно (2.3), первые условия того, что преобразования Эйлера— Дарбу (1.8) и (2.14) являются взаимно обратными, имеют вид

С1 = С + ^' + 2^ (1п г)', (2.16)

С = С1 + ^' + 2^ (1п Г1)'. (2.17)

Исключая С и С1 из (2.16) и (2.17), получаем

^' + ^ (1п г)' + ^ (1п г1)' = 0.

Это соотношение позволяет найти

с

г1 = —=;, с € Д. г^

В дальнейшем будем считать, что с =1. Используя (2.4), записываем вторые два условия того, что (1.8) и (2.14) взаимно обратны:

И1 = И + (^г + Сг) + ^' (1п к)' + 2^ (1п к)'', (2.18)

г

и = И1 + ^ + С1г1) + ^' (1п к1)' + 2^ (1п к1)''. (2.19)

г1

Покажем теперь, как получается выражение для к1. Функция к должна удовлетворять уравнению (2.5). Предположим, что к1 удовлетворяет аналогичному уравнению

^к? + С1к1 + (И + с)к1 = 0. (2.20)

Впоследствии проверим справедливость этого предположения. Подставим в (2.18) и (2.19) функции С1,г1. В результате получим два линейных уравнения относительно И, И1 . Исключив из этих уравнений И, И1 , приходим к выражению

Г'гкк' [(гкк1)' Г + (гкк1)С] + 5 = 0, (2.21)

где 5 не содержит Г'. Явный вид 5 не приводится из-за его громоздкости. Требуя, чтобы выражение в квадратных скобках обратилось в ноль, приходим к соотношению

С

(1п гкк1)' = — —.

Г

Полагая теперь

к, 1

гк exp(J С/Г^ж)’

прямыми вычислениями можно проверить, что данная функция удовлетворяет (2.20) и (2.21). □

Следствие. Операторы Эйлера—Дарбу порождают отношение эквивалентности на классе уравнений Е2,м.

3. Преобразование подкласса Е2,м

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

ихх + С(ж)их = Ви, (3.1)

где В — оператор (1.3), С — некоторая гладкая функция от ж. В этом параграфе будут указаны преобразования, действующие на /2, найдены обратные преобразования и получен аналог теоремы 3. Интерес к уравнениям типа (3.1) связан с их приложениями в механике жидкости и теории упругости.

Как следует из (2.4), преобразование Эйлера—Дарбу

г = 1 (и* — ки) (3.2)

переводит решения уравнения (3.1) в решения уравнения

+ С^ж)^ = Вг, (3.3)

если функция г удовлетворяет уравнению

г'' + Сг' + (С' + 2(1п к)'') г = 0. (3.4)

Прямой проверкой убеждаемся в том, что одним из решений (3.4) является функция

к

г = — к

Знак минус выбран из соображений удобства. Второе решение (3.4) более громоздкое и здесь не используется. Таким образом, справедливо следующее утверждение.

Лемма 2. Преобразование Эйлера—Дарбу

к

г = — ^их + и (3.5)

переводит решения уравнения (3.1) в решения (3.3), если коэффициент С1 задается

формулой

С1 = С + 2(1п(к'/к))', (3.6)

а функция к удовлетворяет уравнению

к'' + Ск' + ск = 0, с € Д. (3.7)

Следствием теоремы 4 является лемма 3.

Лемма 3. Если к — нетривиальное решение уравнения (3.7), то преобразование

(3.5) имеет обратное:

к1 +

и = —ту + г, к1

где

к. = — 1

Ъ7 exp(/ С^ж)

Таким образом, преобразования Эйлера—Дарбу порождают преобразование эквивалентности на /2.

Перейдем к построению высших операторов Эйлера—Дарбу. Введем линейный оператор

= — Т7дх + 1)

Ъ

соответствующий преобразованию (3.5). Предположим, что известны решения Ъ^ ... , Ъ^ уравнения (3.7) при некоторых значениях С1. Как в разд. 2, определим рекуррентным способом функции и операторы

р1 Ъ1) р2 ^рі Ъ2) • • • , ^рм—1 •••^Р1 Ъп)

1 ^рі) 2 ^рг 1) • • • ) N ^рм N—1 •

Каждый оператор М& имеет вид

Мк = 5^ + ^—1$!: 1 + ... + 1,

где а — функции от ж.

Если к1,... , к^ — линейно независимы, то при любом к < N они составляют фундаментальную систему решений уравнения

Мк и = 0. (3.8)

Это уравнение эквивалентно следующему:

Ш(и, к1,... , кк) = 0,

где Ш — вронскиан. Поскольку коэффициент при и в (3.8) равен единице, разлагая вронскиан по первому столбцу и разделив его на Ш (к1,... , кк), получаем равенство

Ш(и,кь ... ,кк)

Мк и = 4 -----—^.

Ш (к/1,...,кк)

Значит, высшее преобразование Эйлера—Дарбу в данном случае имеет вид

г Ш(и,к1,...,кк) (3 9)

Ш(к1,...,кк) . ( . )

Лемма 4. Преобразование (3.9) переводит решения уравнения (3.1) в решения уравнения

+ Ск

где задается формулой

W (Ъ1)...,Ък)

С + 2дх і1" W(Ъь.^Ък)) • (3-10)

Доказательство. Как следует из леммы 2, функции и Сг+1 связаны соотно-

шением

Сі+1 = Сі + 2дх(1п Ъ^/Ъг).

Значит, функция представляется в виде

с+2а*(1п ЪЫ) • (з-11)

Нам остается вычислить произведение

р = Ъ1 ••• Ък Ъ1... Ък'

По построению оператора Мк справедливо представление

Ъ^п , Л / Ъ

Следовательно, коэффициент при старшей производной равен (—1)к/Р. С другой стороны, как показано выше, выражение Мки равно Ш(и, к1,... , кк)/Ш(к1,... , кк). Поэтому коэффициент при дХи имеет вид

(—1)к+1Ш (к1,...,кк)

1 ; ш (к1,...,кк )•

Таким образом, искомое Р равно

Ш (к1 ,...,кк)

Ш(к1,... ,кк)’

и лемма доказана. □

Для линейного обыкновенного уравнения второго порядка

у'' + д(ж)у' + су = 0, с € Д,

формулы, аналогичные (3.9) и (3.10), были получены в работе [16].

4. Применения преобразований Эйлера—Дарбу

Рассмотрим уравнение Эйлера—Пуассона—Дарбу

ИХХ + С(х)их иуу • (4-1)

Оно возникает в исследованиях по газовой динамике и теории упругости [17, 18]. Хорошо известны его решения, найденные Эйлером, соответствующие функции

2п

С — —, п € X. х

Групповой анализ произвольного линейного уравнения второго порядка с двумя независимыми переменными выполнен в [1]. В частности, доказано, что уравнения

с

ихх + их -- , с Е ^,

х

с разными константами с не являются эквивалентными относительно точечных преобразований. Заметим, что относительно преобразований Эйлера—Дарбу эти уравнения при с — 2п, где п — любое целое число, эквивалентны простейшему волновому уравнению [6]:

Ихх 0

Все линейные дифференциальные подстановки первого порядка, соответствующие указанной функции С, получены в [19]. Г. Дарбу [7] несколько расширил список уравнений, приводящихся к (4.1) точечными преобразованиями и допускающих общее явное решение.

Используем лемму 4 для интегрирования уравнений (4.1). Для ее применения необходимо знать решение уравнения (4.1) с некоторой функцией С и решения обыкновенного дифференциального уравнения

Ъ" + С(х)Ъ/ + сЪ — 0, с € Д. (4.2)

Рассмотрим наиболее простой случай С — 0. В зависимости от константы с, решениями уравнения (4.2) являются функции

Й1Х + 61, а2 сов(пх) + 62 в1п(пх), а3 еЬ(пх) + 63 вЬ(пх), (4.3)

где йг,6г € Д, п2 — |с|. На основании леммы 4 можно утверждать, что уравнение (4.1) с коэффициентом

С -2ЭЛ 1п 'Г>'

W(Ъь... ,Ъ*)

имеет решение

— W (И°А .•••.Ъ.fc) (4 4)

где Ъ1,... , Ък — функции вида (4.3) такие, что W(Ъ/ь ... , Ък) — 0, и0 — и(х+у) +V(х-у), а и, V — произвольные гладкие функции.

Последовательным применением преобразования Эйлера—Дарбу можно получать решения, включающие рациональные функции. Например, решением уравнения

(2 6х2 \

ихх + I 3 + 6 I Пх — , 6 Е

V х х | 6 у

является функция

х3 + 6 и — —— (и" + V") - х (и + V') + и + V. 3х

Если

то функция

—2(3х8 + бкх5 — 10тх3 + 30к2х2 + 5кт)

С — ------т~с.---^---------------г-------------------, к,т € Д,

(х6 + 5кх3 — 5к2 + 5тх)(х3 + к)

^ _ (х6 + 5Кх3 — 5к + 5тх)(х3 + к) ,",ттШ , Т^,Л гг„ Т^Л ,

И — 15(х7 + 2кх4 + к2х) (х(и + V ) — и — V ) +

х3 I к

+(и'' + V'') — х (и' + V') + и + V 3х

удовлетворяет уравнению (4.1). Вероятно, указанные решения можно получить из (4.4) некоторым предельным переходом. Следует отметить, что все указанные уравнения принадлежат одному классу эквивалентности, порожденному волновым уравнением.

Пусть теперь п не обязательно целое число. Легко видеть, что уравнение

2п

ихх х Их иуу, п Е (4.5)

обладает частным решением

и — с1 (х2 — у2)п + с2, с1, с2 € Д.

Для того чтобы применить преобразование Эйлера—Дарбу к (4.5), нужно решить уравнение

2п

Ъ"------ЪЪ + сЪ — 0, с € Д.

х

Решения последнего уравнения выражаются через функции Бесселя. Используя лемму 4, можно построить класс эквивалентных уравнений, частные решения которых выражаются через функции Бесселя.

Многомерное уравнение Дарбу

2т

иц — Ди---^и* — 0,

где Ди — п-мерный оператор Лапласа, при целом т эквивалентно волновому уравнению

и^ — Ди — 0.

Эквивалентность устанавливается последовательным применением преобразования Эйлера— Дарбу по переменной £.

Рассмотрим п-мерное стационарное уравнение Шредингера

Ди + С(х1,... , хп)и — 0. (4.б)

Предположим, что

П

С(х)^ Су (х,)

.7 = 1

и для каждой функции Gj справедливо представление

Gj = 2dXj (ln W(hi-----hjj)) ,

где функции hj(xj) задаются формулами (4.3) и все вронскианы W(hj,... , ) не равны

нулю. Тогда, действуя последовательно преобразованиями Эйлера—Дарбу по каждой переменной Xj, можно показать, что уравнение (4.6) эквивалентно n-мерному уравнению Лапласа.

Применение преобразований Эйлера—Дарбу к линейным уравнениям с частными производными позволяет существенно расширить списки уравнений [20], допускающих частные или общие решения.

Список литературы

[1] Овсянников Л.В. Групповой анализ дифференциальных уравнений. М.: Наука, 1978.

[2] Ибрагимов Н.Х. Группы преобразований в математической физике. М.: Наука, 1983.

[3] Olver P.J. Equivalence, Invariants and Symmetry. Cambridge University Press, 1995.

[4] Байков В.А., Жибер А.В. Уравнения математической физики. М.; Ижевск: Институт компьютерных исследований, 2003.

[5] Миллер У. Симметрия и разделение переменных. М.: Мир, 1981.

[6] Эйлер Л. Интегральное исчисление. М.: ГИФМЛ, 1958. Т. 3.

[7] Darboux G. Legons sur la theorie generale des surfaces et les applications geometriques du calcul infinitesimal. P.: Gauthier-Villars, 1915. Vol. 2.

[8] Matveev V.B., Salle M.A. Darboux trasformations and solitons. N. Y.: Springer-Verlag, 1991.

[9] Rogers C., Sohief W.K. Backlund and Darboux transformations: geometry and modern applications in soliton theory. Cambridge University Press, 2002.

[10] Berest Y., Veselov A. On the structure of singularities of intereble Schrodinger operators // Lett. Math. Phys. 2000. N 52. P. 103-111.

[11] Chalykh O.A., Oblomkov A.A. Harmonic oscillator and Darboux transformations in many dimensions // Phys. Lett. A. 2000. Vol. 267. P. 256-264.

[12] Айнс Э.Л. Обыкновенные дифференциальные уравнения. Харьков: ГНТИ, 1939.

[13] Crum М.М. Associated Sturm-Liouville systems // Quart. J. Math. 1955. Vol. 6.P. 121-127.

[14] Eisenhart L.P. Transformations of surfaces. N.Y.: Chelsea, 1962.

[15] Гурса Э. Курс математического анализа. М.; Л.: ГТТИ, 1933. Т. 2. Ч. 2.

[16] Yurova A.A., Yurov A.V., Rudnev М. Darboux trasformation for classical acoustic spectral problem // Intern. J. Math. and Math. Sciences. 2003. Vol. 49. P. 3123-3142.

[17] Овсянников Л.В. Лекции по основам газовой динамики. М.: Наука, 1981.

[18] Аксенов А.В. Симметрии и соотношения между решениями класса уравнений Эйлера— Пуассона—Дарбу // Докл. РАН. 2001. Т. 381, № 2. С. 176-179.

[19] Аксенов А.В. Линейные дифференциальные соотношения между решениями класса уранений Эйлера—Пуассона—Дарбу // Изв. РАН. Механика твердого тела. 2001. Т. 63, № 1. С. 15-20.

[20] Полянин А.Д. Справочник по линейным уравнениям математической физики. М.: Наука, 2001.

Поступила в редакцию 18 апреля 2007 г.