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

В. Ф. Зайцев, А. С. Ложкин

ФУНДАМЕНТАЛЬНЫЕ СИММЕТРИИ ОБЫКНОВЕННЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ

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

V. Zaitsev, A. Lojkin

FUNDAMENTAL SYMMETRIES OF ORDINARY DIFFERENTIAL EQUATIONS

A new approach to the search of symmetries of ordinary differential equations is described. A number of basic theorems are proved, some properties are identified, and possible applications are discussed, with examples provided.

Хорошо известно (см., например, работу [1]), что классический алгоритм поиска допускаемых операторов предполагает решение определяющего уравнения на многообразии, заданном исследуемым уравнением. Например, для уравнения второго порядка

H = у "- F у, у ') = 0 (1)

определяющее уравнение выглядит следующим образом:

dy dy' dy''

= 0, (2)

[H ]

где символ [ H ] означает «на многообразии H и всех его дифференциальных следствиях». Если Ф = Ф(,у,у'), то, выражая у" и у'" из уравнения (1) и его дифференциального следствия Dx ^] = 0, мы получаем вместо уравнения (2)

выражение, зависящее только от первой производной у'. Если задана конкретная зависимость координаты канонического оператора Ф от третьего аргумента (например, для точечного оператора Ф = ^( x, у )-£( x, у) у'), полученное выражение можно «расщепить» по независимой переменной у' и прийти к определяющей системе. Как правило, мы можем найти все решения этой системы, так как число уравнений обычно больше числа неизвестных (т. е. система является переопределенной). Очевидным недостатком этого подхода является необходимость априорного задания типа симметрии. Естественно, для каждого следующего типа задача решается «с нуля».

Тем не менее, существует и другой подход: можно рассматривать инвариантность уравнения не на многообразии его решений, а на всем продолженном

пространстве. По-видимому, этот подход известен с момента введения группового анализа в математическую практику, хотя привести конкретные ссылки весьма затруднительно: все сколько-нибудь значимые результаты применения группового анализа были достигнуты с помощью классического подхода. И это неудивительно — построение эффективного алгоритма применения альтернативного подхода сильно затруднено отсутствием «независимых переменных», по которым можно расщеплять определяющее уравнение.

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

Очевидно, условие инвариантности уравнения (1) на продолженном пространстве можно записать в виде

Опр еделение. Канонический оператор X = Фду, где Ф — решение

уравнения (3), называется фундаментальной симметрией уравнения (1).

Совершенно аналогично определяются фундаментальные симметрии обыкновенного дифференциального уравнения п-го порядка как решения уравнения в полных производных п-го порядка

Решением однородного уравнения (3) является не обычная функция

ства в пространство функций от x). Размерность алгебры фундаментальных симметрий определяется исключительно порядком исследуемого уравнения.

Теорема 1. Размерность алгебры фундаментальных симметрий обыкновенного дифференциального уравнения равно порядку определяющего уравнения, который совпадает с порядком исследуемого уравнения.

Доказательство в виде последовательных утверждений приводится ниже.

Замечание. Для каждого типа локальных симметрий размерность допускаемой алгебры обычных симметрий может быть либо больше, либо равна размерности алгебры фундаментальных симметрий данного типа. Например, для ОДУ 2-го порядка (1) размерность алгебры точечных операторов, согласно теореме Ли, может быть 0, 1, 2, 3 и 8, тогда как размерность алгебры фундаментальных точечных симметрий может быть лишь 0, 1 или 2.

Теорема 2. Фундаментальная точечная симметрия уравнения (1) имеет производящую функцию, не зависящую от производной. Иными словами, опе-

Ф---------+ Ф'----------+ Ф''----------= 0.

ду ду' ду''

(3)

(3а)

Ф^), а функционал Ф(,у,у', ...) (действующий из продолженного простран-

ратор фундаментальной точечной симметрии уравнения (1) имеет вид X = Ф(,у)ду .

Доказательство. Докажем, что функция Ф зависит от x и у, т. е. Ф = Ф( у) . Предположим противное: пусть Ф зависит от производной по у(k)

некоторого конечного порядка k < да . При подстановке в условии инвариантности (3) эта производная в третьем слагаемом перейдет в производную порядка k+2 и не сократится ни с каким другим слагаемым, т.е. условие инвариантности

не будет выполнено ни при какой функции Ф(x,у,у',..., у^)) . Например,

пусть имеется зависимость от первой производной у', тогда после подстановки Ф(, у, у') в выражение (3) получим

-ФFy -

дФ дФ дФ

-^+^- у +^т у '

V

дx ду ду' у дx дxдy ду

д2Ф, 2 дФ

■ + 2----+ 2-----у' у"+—- (у") +-------у= 0.

ду дxдy' дуду' у ду' ду'

Поскольку переменные x, у, у\ у", у'" рассматриваются как независимые, можно «расщепить» это условие по степеням переменной у'". Приравняв нулю коэффициент при первой степени у'", получим

дФ=0,

ду'

т. е. функция Ф зависит только от x и у.

Теорема доказана.

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

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

Теорема 3. Линейное однородное уравнение в полных производных первого порядка имеет одно и только одно фундаментальное решение.

Доказательство. Рассмотрим уравнение

-дH ^}дH л Ф + Ф'-----= 0 .

ду ду'

Очевидно, это уравнение можно проинтегрировать как уравнение с разделяющимися переменными. Имеем

дН

Ф' = ду Ф дН ’

ду'

то есть

дН

■ ду

Ф = е

ду'

Таким образом, для уравнений первого порядка теорема 1 верна.

Далее, если известна одна фундаментальная симметрия уравнения (1), то мы можем найти вторую [2]. Пусть производящая функция известной фундаментальной симметрии — функция Ф1 (x, у, у',...) . Будем искать Ф2 в виде

Ф2 = Ф1^, где и - новая неизвестная функция, которая будет удовлетворять уравнению

дН

Фх— [и ] +

1 ду''

Ф

дН

1 ду'

дН

+2 о [ф 1^) о [и ]=0

Выполнив замену Ж = Бх [и], получим относительно Ж линейное однородное уравнение уже первого порядка, решЬё которое получим

Ж = Ф1 е

дН

т1 ду!

* ж_

ду"

то есть

ф 2 = фб;1

Ф-2е

дН

ду_

дН

ду"_

о

Здесь Бх — символ оператора полной производной, в уравнениях (2), (3) для краткости этот оператор обозначен штрихом ('), Б-1 — символ оператора,

обратного к оператору полной производной. Этот оператор называют еще полным интегралом. Мы будем в ряде случаев использовать знак интеграла, при необходимости указывая, что интеграл — полный. Заметим, что теперь для доказательства теоремы 1 нам остается доказать лишь, что любое однородное ли-

нейное уравнение в полных производных (За) имеет хотя бы одно нетривиальное решение.

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

в нуль на [ Н ]:

, дН дН дН -А / /а\

Ф^~+ф'+фV; = ЕА,(у' ...))[Н], (4)

ду ду ду к=о

где Ак — произвольные функции из бесконечно-продолженного пространства. Совершенно аналогично классическому методу Лагранжа ищем решение в

виде

¥ = С1 (y,у', ...)ф1 (y,у\ ...) + С2 (x,^у \ ...)ф2 (^y,у\ ...),

где Ф1, Ф 2 — фундаментальные решения исходного однородного уравнения, а

С1, С2 — произвольные константы, которые мы рассматриваем как функции от

всех переменных продолженного пространства.

Вычисляем последовательно производные функции ¥, накладывая дополнительные условия на первые производные функций С1, С2. Получим

¥' = С1Ф '1 + С2Ф'2 + С \ Ф1 + С\ Ф 2, С \ Ф1 + С\ Ф 2 = 0,

¥'' = С1Ф '\ + С2Ф ''2+С \ Ф \ + С '2 Ф '2,

ад

СФ', + С '2 Ф'2 = 1 Ак (, у, у ',...)) [Н ].

к =0

Очевидно, что ¥ является решением неоднородного уравнения (4), если функции С1, С2 удовлетворяют линейной системе

С \ Ф1 + С '2 Ф 2 = о,

ш (5)

С ',Ф',+ с, Ф'2 =1 А, (, у, у',...)) [Н ]

к=0

Так как Ф1, Ф 2 линейно независимы, то вронскиан системы (5) не равен 0. Из первого уравнения системы находим

Ф

С 1 = —С 1 2

С1 - С 2 Ф1,

из второго

2 А (уу,у[Я]

С’г =

к=0

Ф,

Интегрирование дает функцию С2:

С = О-'

2 Ак (уу,у',...)) [Я]

к=0

Ф -

( Ф '2 Ф \ ^

Ф- Ф'

Соответственно для С1 получим:

С' = -О-

2 Ак (у у, у',...)) [Я ]

к=0

Ф,

Конкретные классы симметрий можно выделить и найти, наложив дополнительные условия на найденное общее решение уравнения (4) до перехода на многообразие (1). Заметим, что предложенный выше алгоритм нахождения всех симметрий уравнения применим и для операторов, в которых Ф' задана в виде

бесконечного ряда, т. е. с нелокальной Ф .

Теперь завершим доказательство теоремы 1, для чего покажем, что любое уравнение (За) имеет хотя бы одно нетривиальное решение в виде формального степенного ряда. Из-за громоздкости выкладок мы приведем лишь несколько первых шагов алгоритма на примере уравнения (3). Для уравнения (1) определяющее уравнение (3) можно записать в виде

Ф''- ¥уФ'- ¥у Ф = 0.

Пусть функция ¥ может быть разложена в ряд Тейлора по третьему аргументу:

¥ (у у, у ') = 2 /т (х у)у')" .

т=0

Тогда решение — функция Ф — можно разложить в обобщенный степенной ряд

Ф = 2Фк ( у, у',..., у|к-,))у|к).

к =0

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

ад

ЕкУ‘+(2(Ф, -F^,)-(, -FyD^, -Fy%)] = 0.

к=0

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

ФоУ "+ (2(Фо - Fy 'Фо ) У (Dx Фо - FyDxФо - Fy Фо ) .

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

Фо,, - /¡Фо, - /оФо = о.

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

Фо ( У ) = С1 (У ) Фо1 (Х) + С2 (У ) Фо2 (Х) .

Подставив полученное решение в следующие коэффициенты ряда и выделив коэффициент при у' в первой степени, находим C1 (у) и С2 (у) . Дальнейшие

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

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

Огромные возможности применения данного алгоритма порождает вопрос

о возможности редукции простой симметрии в фундаментальную. Справедлива следующая теорема.

Теорема 4. Любая симметрия обыкновенного дифференциального уравнения 2-го порядка (1), заданная производящей функцией Ф(х, у), является фундаментальной симметрией равносильного ему уравнения /Н=0, при этом множитель / равен СФ-1, где С = const.

Доказательство. Пусть функция Ф(х, у) определяет канонический оператор, являющийся симметрией уравнения (1), тогда она удовлетворяет определяющему уравнению:

. дН ^,дН . „дН,

Ф----+ Ф'----+ Ф''-- „ = 0.

ду ду1 дуy_F

Домножим И на f (х, у, у') . Тогда определяющее уравнение запишется в следующем виде:

f (Ф ( + 2фхуУ Чфуу (у f +ф уУ ") + (ф, +фуУ ') Ну ^ +

+(ф,+ф,.у')Н/,+ФHу.f+фи/, _ 0.

Расщепляем это уравнение по степеням у". Коэффициент при первой степени дает уравнение

f ф у +(ф., +ф уУ ') /у,+ф/у _ 0,

решая которое, находим:

у _ф z (х ) + |^фх dy ].

В определяющем уравнении переходим к новым переменным:

1 С ф

z _— У - I —2 dy, у _ t; фу J ф2 у

после преобразований получаем:

f Фt +ф/ _ о,

откуда следует требуемое равенство: f _ Сф-1, С = const.

Теорема доказана.

Таким образом, можно сделать вывод, что любая точечная симметрия, удовлетворяющая условию ^_ 0, может быть переведена в фундаментальную симметрию эквивалентного уравнения домножением исходного уравнения на множительf Назовем такой множитель фундаментализирующим.

Рассмотрим пример.

(у ')2

Пусть дано уравнение у" _-—— + ху' с допускаемым оператором

у

X _ Уд у.

Таким образом, мы имеем И _ уу-1 (у')2 - ху'; это уравнение обладает точечной симметрией с порождающей функцией ф _ у. По теореме 2 получаем, что фундаментализирующий множитель f _ Су-1, С _ const. Проверим, что это действительно так.

После домножения на фундаментализирующий множитель определяющее уравнение приобретает следующий вид:

/у *+ /у'(-2у-1 у'- х) + '(У''- У-1 (у')2 - ху') +

+/у"' (у1)2 + /уУ (У'- У-1 (у f - ху') _ 0.

Расщепляя это уравнение по степеням у", получаем

f + yfy + У ' fy'_ 0, (6)

откуда

dy' dy

У “ 7 ’

следовательно, у' _ yz (х). Переходим к новым переменным: z _ —, y=t.

В новых переменных уравнение (6) приобретает следующий вид:

f + f _ 0.

Отсюда получаем f _ Ct-1 или f _ Су-.

Приведем еще один пример.

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

Возьмем уравнение

И _У"+ У^1 (у')2 -3ху _ 0,

имеющее точечную симметрию ф _ уду.

По теореме 4 можем превратить эту симметрию в фундаментальную, умножив уравнение И на множитель f _ Су-1, где С = const. Действительно, тогда

X /и _ у (у-2 у ( 2 у-1 (у')2) + у' (2у-2 у') + у'' у-| _ 0,

т. е. симметрия является фундаментальной.

Найдем вторую симметрию:

дН “if ^

ф2 _ С2ф0|ф-2£ дУ dx _ С2у|y~2exp(-J2у-1у'dx^dx _ С2у|y~4dx.

Весь класс симметрий находим с помощью метода Лагранжа, используя формулы для нахождения С1 и С2 :

с, = -1

Г V

У'+ У

У

У

у

с=1А

У + У

Л

Г1 4 У' I У~4дх

У ■’ у

дх,

где

к=0

Ху ')2

У "+Ш— 3хУ

У

Обратимся вновь к уравнению (1). Как мы уже говорили, зная одну фундаментальную симметрию этого уравнения, без труда можно найти вторую, которая будет иметь следующий вид:

Г -Г-д*

ф2 = ф, Гф-2е 5дУ' дх .

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

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

¥ = п(х, у )ехр (ГС(*, У, У ')дх),

то она обладает первым дифференциальным инвариантом, а следовательно, уравнение может быть факторизовано. Справедлива следующая теорема.

Теорема 5. Для того чтобы вторая фундаментальная симметрия обыкновенного дифференциального уравнения второго порядка могла быть приведена к экспоненциальной форме, необходимо, чтобы это уравнение имело следующий вид:

н = у / (х у )(у ' )2 - £ (у У )у - Н (у у ).

Доказательство. Для того чтобы вторая фундаментальная симметрия могла быть приведена к экспоненциальной форме, необходимо, чтобы подынтегральное выражение в формуле для ф2 было равно полной производной по х для некоторой функции ¥ :

-ГдНА

ф-2е -1 дУ' = Вх [¥]. (7)

Пусть ¥ представима в виде

Г--дх

¥ = ¥0е д' .

Тогда подстановка этого выражения в равенство (7) дает уравнение для определения ¥0:

ф-2 = Б [¥ ]-¥ —

1 х^ ^ 0 ду’

Очевидно, что ¥0 может зависеть только от х и у. Исходя из предположения, что фундаментальные симметрии зависят только от тех же переменных, получаем

дн

ду'

откуда

¥ 0 =¥ 0

д2 Н

д(у ')2

или

¥п

¥ 0

д2 Н

д(у')!

Интегрируя это выражение, получаем, что

д2 Н

Чу)

может зависеть только

от х и у. Следовательно, Н должно иметь следующий вид:

н = У"- / (х У )(У' )2 - £ (х У )У'- Н (х У).

Теорема доказана.

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

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

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

2. Зайцев В. Ф., Ложкин А. С. О разрешимости обыкновенных дифференциальных уравнений второго порядка // Труды математического центра имени Н. И. Лобачевского. Казань, 2004. Т. 24. С. 48-52.

У