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

УДК 517.9

ЭЛЕМЕНТАРНАЯ ТЕОРИЯ ИСПОЛЬЗОВАНИЯ ИНВАРИАНТОВ ДЛЯ РЕШЕНИЯ МАТЕМАТИЧЕСКИХ УРАВНЕНИЙ

© 2008 А.Д. Полянин1

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

Ключевые слова: инвариант, алгебраическое уравнение, дифференциальное уравнение, групповой анализ.

1 Полянин Андрей Дмитриевич ([email protected]), Учреждение Российской академии наук, Институт проблем механики им. А.Ю. Ишлинского РАН, 119526, Россия, г. Москва, пр-т Вернадского, д. 101, корп. 1.

1. Симметрии. Общая схема использования инвариантов для решения математических уравнений

1.1. Симметрии. Преобразования, сохраняющие вид уравнений. Инварианты

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

Под симметриями математических уравнений будут пониматься преобразования, сохраняющие вид уравнений. Ниже приведены примеры конкретных уравнений, которые сохраняют свой вид при некоторых простых преобразованиях.

Пример 1: Рассмотрим биквадратное уравнение

х4 + ах2 + 1 = 0. (1.1)

Замена х= — х приводит к точно такому же уравнению

Х + a х2 + 1 = 0

(т.е. уравнение (1.1) сохраняет вид при преобразовании х = — х).

Два других преобразования х = ±1/Х также сохраняют вид уравнения (1.1),

4

поскольку после умножения на х получим

х4 + ах2 + 1 = 0.

Пример 2: Вид дифференциального уравнения

y'L — = 0 (1.2)

не изменится, если сделать любое преобразование х = х + a, y = y (а — любое число);

х = х, y = y + b (b — любое число); (1.3)

х = х, y = cy (c — любое число, не равное нулю),

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

хх х— Хх = 0.

Далее будет показано, что преобразования, сохраняющие вид уравнений, позволяют "размножать" решения уравнений.

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

Пример 3: Преобразование одновременного сдвига по обоим осям

х = х + а, у = у + а, где а — любое число, имеет инвариант

I = У - х = X - X.

Если х — независимая переменная, а у — зависимая переменная, то другими инвариантами этого преобразования являются производные

12 = УХ = ух 13 = У'Хх = ~хуХ X ••• •

Пример 4: Преобразование одинакового изменения масштабов длины по обоим осям

х = ах, у = ах, где а ф 0 — любое число, имеет инвариант

1=1 = 1 х X

Если х — независимая переменная, а у — зависимая переменная, то имеются также более сложные инварианты, которые зависят от производных и сохраняют вид при действии данного преобразования, например,

12 = ух = Ъ 13 = ху"ххх = х х X ••• •

1.2. Общая схема использования инвариантов для решения математических уравнений

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

Принципиальная схема исследования математических уравнений

После выполнения указанных действий уравнение часто упрощается или приводится к разрешимому виду. Важно отметить, что данная схема с успехом может применяться для самых различных типов математических уравнений (см. далее разделы 2-4).

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

2. Алгебраические уравнения и системы уравнений

2.1. Уравнения, содержащие четные степени

Рассмотрим алгебраическое уравнение

а2пх2п + а2п-2 х2"-2 + а2п-4 х2п-4 + ••• + а4 х4 + а2 х2 + а0 = 0, (2.1)

содержащие только четные степени. Биквадратное уравнение является частным случаем уравнения (2.1) при п = 2. Замена

х = - х (2.2)

приводит точно к такому же уравнению для х (говорят, что уравнение (2.1) инвариантно относительно преобразования (2.2)). Отсюда следует, что если уравнение (2.1) имеет решение х = хх, то оно имеет также другое решение х = -хх.

Возводя в квадрат обе части (2.2), получим простейшую алгебраическую функцию, которая сохраняет вид при преобразовании (2.2):

:2 = X2.

(2.3)

Эта функция и является инвариантом преобразования (2.2). Если выбирать инвариант (2.3) за новую переменную, г = х2, то уравнение (2.1) порядка 2п преобразуется к уравнению порядка п:

а2пгп + а2п-2 тп-1 + «2п-4 тп-2 + ••• + а^т1 + а2г + а0 = 0.

Таким образом в данном случае переход от исходной переменной х к инварианту преобразования (2.3) г = х2 позволяет упростить рассматриваемое уравнение (понизить его порядок в два раза).

2.2. Возвратные уравнения

Возвратное алгебраическое уравнение четной степени имеет вид

а0х2п + а1х2п-1 + а2х2п-2 + ••• + а2х2 + а1х + а0 = 0 (а0 ф 0). (2.4)

Левая часть этого уравнения называется возвратным многочленом. Замена

1 , л

х = ^ (2.5)

преобразует уравнение (2.4) в точно такое же уравнение Отсюда следует, что если уравнение (2.4) имеет корень х = Ь, то оно имеет также другое решение х = 1/Ь.

Простейшее возвратное уравнение — квадратное уравнение

а0 х2 + а1 х + а0 = 0.

Делим его на х2:

а§{х Л—) + а\ =0.

Результат удобно представить в форме уравнения первого порядка

а0т + а1 = 0,

где

11

г = х + - = х +- (2.6)

х х

— простейший инвариант преобразования (2.5).

Теорема 1: В общем случае возвратное уравнение порядка 2п допускает понижение порядка с помощью подстановки (2.6). В результате получается алгебраическое уравнение порядка п.

Пример 1: Рассмотрим возвратное уравнение четвертого порядка ах4 + Ьх3 + сх2 + Ьх + а = 0.

Делим его на х2:

а(х2 + + + -) + с = 0. (2.7)

Учитываем, что

(х + -)2 = х2 + 2 + \ => х2 + \ = I2 - 2.

х' г2 г2

X'

В результате уравнение (2.7) путем замены (2.6) (основанной на инварианте преобразования (2.5)) сводится к квадратному уравнению

аг2 + Ъг + с - 2а = 0.

Теорема 2: В общем случае возвратное алгебраическое уравнение нечетного порядка

Р2и+\( х) = 0,

где Р2п+1 (х) = аох2п+1 + а\х2п + а2х2п-1 +-----+ а2х2 + а\х + ао имеет корень х = -1,

а его левая часть допускает представление

Р2п+1( х) = (х + 1)02п( х),

где Q2n(х) — возвратный многочлен степени 2п.

Пример 2: Возвратное уравнение третьего порядка

а0 х3 + а1 х2 + а1 х + а0 = 0

может быть представлено в виде

(х + 1)[а0х2 + (а1 - а0)х + а0] = 0.

Из теоремы 2 следует, что возвратное уравнение порядка 2п + 1 после деления на (х +1) и введения новой переменной (2.6) сводится к уравнению порядка п.

2.3. Обобщенные возвратные уравнения

Рассмотрим алгебраическое уравнение

а0х2п + а1 х2п-1 + ••• + ап-1 хп+1 + апхп+ (2 8)

+Лап-1 хп-1 + Х2ап-2хп-2 + ••• + \п-1а1 х + \па0 = 0 (а0 ф 0). (.)

Первые п+1 членов этого уравнения (которые расположены на первой строке) совпадают с соответствующими членами возвратного уравнения (2.4), а оставшиеся члены (которые расположены на второй строке) отличаются множителями вида ~кт от соответствующих членов возвратного уравнения. Нетрудно проверить, что преобразование

х = (2.9)

сохраняет вид уравнения (2.8), а простейший инвариант преобразования (2.9) записывается так:

1-1 , г = х+ - = х+- (2.10)

х х

Введение новой переменной по формуле (2.10) приводит (2.8) к уравнению п-го порядка.

Пример 3: Рассмотрим уравнение четвертой степени ах4 + Ьх3 + сх2 - Ьх + а = 0, которое является частным случаем уравнения (2.8) при п = 2, X = -1. Замена

1

г = х — х

приводит к квадратному уравнению

ат2 + Ьт + 2а + с = 0.

2.4. Системы алгебраических уравнений,

симметричные относительно перестановки аргументов

Многочлен Р(х, у) от двух переменных называется симметрическим, если он не меняется при перестановке аргументов: Р(х, у) = Р(у, х).

Замечание: В терминах преобразований симметрический многочлен определяется как многочлен сохраняющий вид при преобразовании х = у, у = х.

Простейшие симметричные многочлены

и = х + у, w = ху (2.11)

называются элементарными. Эти многочлены являются простейшими алгебраическими инвариантами при перестановке аргументов.

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

Р( х, у) = 0, 0(х, у) = 0,

где Р и Q — симметричные многочлены, полезно в качестве новых переменных использовать элементарные симметрические многочлены (2.11). Подобные системы обладают следующим свойством: если система имеет решение х = х0, у = у0, то она имеет также решение х = у0, у = х0.

Пример 4: Рассмотрим нелинейную систему алгебраических уравнений

х2 + аху + у2 = Ь,

х4 + схх2у2 + у4 = й. ( . )

Система (2.12) не меняется при перестановке аргументов.

(2.13)

Переходя в (2.12) от х, у к переменным (2.11) и учитывая формулы х2 + у2 = (х + у)2 - 2ху = и2 - 2w,

х4 + у4 = (х2 + у2)2 - 2х2у2 = (и2 - 2w)2 - 2w2 = и4 - 4и^ + 2w2,

получим

и2 - (а - 2^ = Ь, и4 - 4и^ + (с + 2^2 = й.

Исключая из уравнений и, приходим к квадратному уравнению

(а2 + с - 2^2 - 2aЬw + Ь2 - й = 0.

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

22

х2+у2 = а' (214)

х3 + у3 = Ь.

Переходя к переменным (2.11) и учитывая равенство х3 + у3 = (х + у)3 - 3ху(х + у), получим

и2 - 2w = а, и3 - 3uw = Ь.

Исключив w, приходим к кубическому уравнению

и3 - 3аи + 2Ь = 0. (2.15)

Отметим, что прямое исключение из системы (2.14) переменной у приводит к существенно более сложному чем (2.15) уравнению шестого порядка:

(а - х2)3 = (Ь - х3)2 2х6 - 3ах4 - 2Ьх3 + 3а2х2 + Ь2 - а3 = 0.

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

Краткие выводы

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

3. Обыкновенные дифференциальные уравнения

3.1. Преобразования, сохраняющие вид уравнения, и их инварианты

Говорят, что обыкновенное дифференциальное уравнение

F (х, у, у'х,...,у(п}) = 0 (3.1)

инвариантно относительно (обратимого) преобразования

х = ф(х, у), у = у(х, у), (3.2)

если в результате подстановки (3.2) в (3.1) получим точно такое же уравнение

F (х, у, Ууп)) = 0. (3.3)

Функция F в уравнениях (3.1) и (3.3) одинакова.

Преобразования, сохраняющие вид уравнения, позволяют "размножать" его решения. Действительно, пусть известно частное решение

У = *(х) (3.4)

уравнения (3.1). Поскольку уравнение (3.1) после перехода к новым переменным (3.2) сохраняет такой же вид, то преобразованное уравнение (3.3) имеет решение

У = 8( х). (3.5)

Возвращаясь в (3.5) к старым переменным по формулам (3.2) (их предварительно надо разрешить относительно х и у) получим решение уравнения (3.1), которое, в случае общего положения будет отличаться от исходного решения (3.4).

Пример 1: Уравнение третьего порядка

ух'х - ух = 0 (3.6)

имеет частное решение

У = ех.

Преобразование

х = х + а, у = у + Ъ (3.7)

сохраняет вид этого уравнения, поэтому преобразованное уравнение третьего порядка Ухху- Уу = 0 имеет решение у = еу. Подставляя сюда старые переменные, полученные обращением формул (3.7), имеем новое решение уравнения (3.6):

у = Аех + Ъ, А = е~а,

содержащее две произвольные постоянные А и Ъ.

Функция I(х, у, у'х) (отличная от постоянной) называется инвариантом преобразования (3.2), если она сохраняется при этом преобразовании

1(х, у, ух) = I(х, У, У).

Замечание: Если I = 1(х, у, у'х) — инвариант преобразования (3.2), то ¥(!), где ¥ — произвольная функция, тоже является инвариантом данного преобразования.

3.2. Процедура понижения порядка уравнений при п ^ 2 (приведение уравнений к разрешимому виду п = 1)

Конкретизируем описанную в разделе 1 общую схему использования инвариантов для анализа математических уравнений применительно к обыкновенным дифференциальным уравнениям. На первом этапе ищется преобразование

х = ф( х, у; а), у = х, у; а), (3.8)

которое сохраняет вид уравнения (3.1). Преобразование (3.8) обязательно должно зависеть от одного свободного параметра а е [а1, а2] (само исходное уравнение (3.1) от этого параметра не зависит).

На втором этапе рассматриваемой схемы для уравнений второго и более высоких порядков (при п ^ 2) строятся два функционально-независимых инварианта преобразования (3.8):

¡1 = Ь(х, у), ¡2 = ¡2(х, у, у'х). (3.9)

На третьем этапе инварианты (3.9) выбираются в качестве новых переменных для уравнения (3.1), т.е. делается преобразование

и = ¡2, т = ¡1, и = и(т).

Указанная процедура приводит к понижению порядка уравнения на единицу.

Для уравнений первого порядка (при п = 1) на третьем этапе в (3.1) надо сделать замену

т = ¡1 , т = т( х).

В результате уравнение приводится к разрешимому виду (к уравнению с разделяющимися переменными).

3.3. Используемые преобразования. Процедура определения инвариантов

Далее в качестве преобразований (3.8) будем использовать только простейшие преобразования вида

х =х + А, у = у + В преобразование сдвига; х = Ах, у = В у преобразование масштабирования

и их комбинации

х = А1У + В!, у = А2 у + В2. (3.10)

В этом случае между производными имеют место линейные соотношения

' А2~1 " А2~>> (п) А2"(п) /о11\

Ух = -гЪ Ухх = ~аТх-х Ух = ■ (З-11)

А 1 А 1 А 1

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

Справедливо следующее утверждение. Пусть преобразование (3.10) сохраняет вид некоторого уравнения, имеющего частное решение (3.4). Тогда это уравнение имеет также решение

/ х - В1 ч у = В2 + Агк{—^-).

Первый инвариант 11 получаем путем исключения параметра а из равенств (3.10), а второй инвариант 12 — путем исключения параметра а из одного из соотношений (3.10) и первого соотношения (3.11).

3.4. Анализ конкретных обыкновенных дифференциальных уравнений

Пример 2: Уравнение второго порядка, не зависящее явно от у:

ух = Ях, ух). (3.12)

Это уравнение не меняется при произвольном сдвиге по зависимой переменной: у у + а (соответствует у = у + а), где а — свободный параметр. При этом из трех величин х, у, у'х две остаются неизменными:

х ух.

Это инварианты уравнения (3.12), т.е. 11 = х, 12 = ух. Их выбираем за новые переменные:

и = у'х, г = х, и = и(г).

В результате получим уравнение первого порядка: и'х = Р(х, и).

Пример 3: Автономное уравнение второго порядка, не зависящее явно от х:

уТ = р (у, ух). (3.13)

Это уравнение не меняется при произвольном сдвиге по независимой переменной: х х + а, где а — свободный параметр. При этом из трех величин х, у, у'х две остаются неизменными:

у> ух.

Это инварианты уравнения (3.13), т.е. 11 = у, 12 = у'х. Их выбираем за новые переменные:

и = у'х, г = у, и = и(г). В результате получим уравнение первого порядка: ии'у = Р(у, и).

Пример 4: Нелинейное уравнение второго порядка

у;х = ур{х^) (3.14)

не меняется при растяжении зависимой переменной: у ау. При этом две комбинации из трех величин х, у, у'х остаются остаются неизменными:

х, —.

у

Это инварианты уравнения (3.14). Их выбираем за новые переменные:

и — —, 1 = х, и = и(г). у

Дифференцируя и по х, имеем

х у у у

Исключая с помощью этого соотношения у'хх в (3.14), приходим к уравнению первого порядка

и!х = (х, и) - и2.

Замечание: При _Р(х, и) = #(х)+/(х)и исходное уравнение (3.14) является общим линейным однородным уравнением второго порядка. Указанным преобразованием оно сводится к уравнению первого порядка с квадратичной нелинейностью.

Пример 5: Рассмотрим нелинейное уравнение второго порядка

уу'хх - (ух)2 = куъеХх. (3.15)

Ищем инвариантное преобразование в виде

х = х + Ь, у = ау. (3.16)

Подставив (3.16) в (3.15), а затем сократив на а, получим

XXх- (у)2 = аекЪку3екх.

Требуя совпадения с уравнением (3.15), получим соотношение для определения параметра Ь:

аекЪ = 1 => Ъ = -\\па. (3.17)

к

При этом параметр а остается свободным.

Подставим (3.17) в (3.16), а затем из второго соотношения исключим параметр а с помощью первого соотношения. Имеем

у = уе(х-х) уе^ = уекх.

Таким образом получен первый инвариант преобразования (3.16)-(3.17) в виде

¡1 = уекх. (3.18)

Чтобы найти второй инвариант, вычислим производную

ух = аух

Исключая здесь параметр а с помощью второго соотношения (3.16), получим второй инвариант

у' у— = ^ = 12. (3.19)

у у

Для понижения порядка исходного уравнения надо принять за новые переменные инварианты (3.18)-(3.19):

г = еХху, и = у, и = и(г). (3.20)

С одной стороны

= (3.21)

у ку' у

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

и'х = игт!х = (ХеХху + екху'х)и'г = (к + е1ху^)и'г = (Хг + ги)и'г. (3.22)

Приравнивая (3.21)-(3.22), имеем

у"^ у"

— - и1 = (Хг + гиК => — = и + (М + ги)и'.

уу

Подставив в (3.3), после элементарных преобразований получим уравнение первого порядка с разделяющимися переменными:

(X + и)и'г = к.

Замечание: Аналогичными свойствами обладает более общее нелинейное уравнение второго порядка

Преобразование (3.20) приводит его к уравнению первого порядка:

и2 + (Хг + ги)и[ = Р (г, и). Пример 6: Рассмотрим теперь нелинейное уравнение первого порядка

, _ ^/ахл: + |3и + ух

Ух ~

/ - Л«1*+Р1У + У1\ ,3 23)

К а2 х + в2у + 72 >

Левая часть этого уравнения не меняется при преобразованиях вида

х = ах + Ъ, у = а у + с, (3.24)

где а, Ъ, с — произвольные постоянные. Подставим (3.24) в аргумент правой части уравнения (3.23):

а!Х+Р1У +У1 _ аСсдх+взО + а^ + Рхс + у! а2х + $2у + у2 а(а2 х + у) + с^Ь + (32с + у2'

Чтобы уравнение (3.23) было инвариантным относительно преобразования (3.24) в (3.25) надо положить

aib + Pic + Yi = ayi, агЪ + + уг = ауг- (3.26)

Эти соотношения можно рассматривать как систему двух линейных алгебраических уравнений первого порядка для определения коэффициентов Ъ и с. При этом коэффициент а остается произвольным. Таким образом уравнение (3.23) инвариантно относительно преобразования (3.24), (3.26), а первым инвариантом этого преобразования является аргумент правой части уравнения (3.23). Поэтому в уравнении (3.23) надо сделать замену

+ где * = (3.27)

аг х + вгУ + Y2

Разрешим (3.27) относительно у, а затем продифференцируем по х. Заменив затем у'х на f(z) (это следует из (3.23) и (3.27)), после элементарных преобразований приходим к уравнению с разделяющимися переменными

(a2|3i - ai|32)x + P1Y2 ~ P2Y1 / _ я ч °-iz - ai (Piz-PO2 Zx~m+ fcz-fr'

В таблице 1 представлены некоторые обыкновенные дифференциальные уравнения второго порядка, допускающие понижение порядка путем использования простейших инвариантных преобразований. Для уравнений первого порядка (когда функция F(u, v, w) не зависит от третьего аргумента) уравнения, указанные в таблице 1, решаются путем перехода от у к новой зависимой переменной z = Ii(x,у), где Ii —первый инвариант.

Приведенные в таблице 1 результаты легко обобщаются на нелинейные уравнения произвольного порядка.

Краткие выводы

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

Таблица 1

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

Уравнение Инвариантное преобразование Первый инвариант Второй инвариант

F(x,y'x,y'D = о У = У + а 1\ - X Ь= У

F(y,?M = 0 х = х + а к =У Ь=У

F(ax + ßy + у, y'x, yxx) = 0 х = х + aß, у = у-аа 1\ = ах + (Зу + у ь = у

У у" Fix,—, — =0 У У 1 У = ау 1\ - X У

F(y,xy'x,x1y;x) = 0 х - ах к = у Ь = ху'х

F^yJ^ik J = 0 х = х — Ina, к-У = ау 1х = екху 1г = ~ У

F(ekxy, ekxy'x, ekxyxx) = 0 - 1, х = х - — та, к- ' у = ау 1х = екху к = екху'х

F(xeky, xy'x, x2y';x) = 0 X — & Х^ у= У~\Хпа 1х = хеку 12 = Ху'х

F(xky,xk+1yx,xk+2y^x) = 0 х - ах, у - аку 1\ = **у к = **+1у'х

F x"y"\Ухх =0 y ' У 1 .v а" .V. у а ' у к = хпут 2 У

4. Дифференциальные уравнения с частными производными

4.1. Преобразования, сохраняющие вид уравнения, и их инварианты

Будем рассматривать обратимые преобразования

х = ф(х,г, w), г = х(, w), w = х(х$, w). (4.1)

Говорят, что уравнение с частными производными / дw дw д^ д2w д2w \

инвариантно относительно преобразования (4.1), если в результате подстановки (4.1) в (4.2) получим точно такое же уравнение

/-- - д w д w д2 w д2 w д2 w \ „ ...

F[x,t,w,—,—, ) = 0. (4.3)

^ д х dt д? дхд dt2 ' V 7

Функция Р в уравнениях (4.2) и (4.3) одинакова.

Величина 1(х, г, м>) называется инвариантом преобразования (4.1), если она сохраняется при этом преобразовании

I(х, г, м>) = 1( х,г,

Замечание. Если I = !(х, г, w) — инвариант преобразования (4.1), то ¥(!), где ¥ — произвольная функция, тоже является инвариантом данного преобразования.

4.2. Процедура построения точных решений

Конкретизируем описанную в разделе 1 общую схему использования инвариантов для анализа математических уравнений применительно к дифференциальным уравнениям с частными производными. На первом этапе ищется преобразование

х = ф(х,г , w; С), г = ху w; С), w = х( х[, w; С), (4.4)

которое сохраняет вид уравнения (4.2). Преобразование (4.4) обязательно должно зависеть от одного свободного параметра С е [Сь С2] (причем само уравнение (4.2) от этого параметра не зависит).

На втором этапе строятся два функционально-независимых инварианта преобразования (4.4):

I! = Ь(х,у, w), I2 = Ь(х,у, w). (4.5)

На третьем этапе инварианты (4.5) выбираются в качестве новых переменных и точное решение уравнения (4.2) ищется в виде

I2 = Ф(!0, (4.6)

где Ф(г) — функция, удовлетворяющая обыкновенному дифференциальному уравнению, которое получается путем подстановки (4.6) в (4.2) (соответствующее решение называется инвариантным).

Как и ранее, в качестве преобразований (4.4) будем рассматривать только простейшие преобразования

х = х + а, г = г + Ъ, w = w + с (преобразование сдвига); (4.7)

х = ах, г = Ъ, w = cw (преобразование масштабирования) (4.8)

и их комбинации.

4.3. Примеры построения точных решений нелинейных уравнений математической физики

Пример 1: Рассмотрим нелинейное уравнение нестационарной теплопроводности

дw д г дw 1

аГ = (4'9)

Это уравнение не меняется при произвольном сдвиге по зависимым переменным

г = г + Ск, х = X + СХ, w = W, (4.10)

где С — свободный параметр, а к и Х —некоторые заданные числа (которые могут выбираться произвольно). Исключая из первых двух соотношений (4.10) параметр С находим один из инвариантов

х- х г-г -

—-— = —— => кх - м = к х - м = 1\. Хк

Вторым инвариантом здесь является /2 = w = w. В силу (4.6) уравнение (4.9) допускает решение вида

w = Ф(г), г = кх - Хг. (4.11)

Подставив (4.11) в (4.9), приходим к обыкновенному дифференциальному уравнению

к2[/ (Ф)Ф']' + ХФ' = 0.

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

J ХФ + С1

где С1 и С2 — произвольные постоянные.

Замечание 1: Решения вида (4.11) называются решениями типа бегущей волны [6,7].

Замечание 2: Решения типа бегущей волны допускают уравнения общего вида, которые не зависят явно от независимых переменных:

/ дw дw д2w д2w д2w ч

Подставляя (4.11) в (4.12), получим автономное обыкновенное дифференциальное уравнение относительно функции Ф(г):

Р(Ф, кФ', -ХФ', к2Ф'', -кХФ'', Х2Ф'', ...) = 0.

Пример 2: Рассмотрим уравнение нестационарной теплопроводности с нелинейным источником степенного вида

дw д2w

Ищем инвариантное преобразование масштабирования. Подставив (4.8) в (4.13) и умножив все члены на Ь/с, имеем

д w Ь д2 w 1 — Требуя совпадения с уравнением (4.13), получим два соотношения

А = 1, ьс"-1 = 1.

а2

Выразим параметры а и с через Ь:

а = лГь, с = ¿1-". Подставляя эти величины в (4.8), находим инвариантное преобразование

? = Ь, }у = Ь~:й>, (4-14)

содержащее свободный параметр С = Ь.

Перейдем теперь к определению инвариантов. Из первых двух соотношений (4.14) исключаем Ь:

х=(1)1/2Х => хГ1'2 = И~~1'2 => 11=хГ1'2. (4.15)

Из второго и последнего соотношений (4.14) исключаем Ь:

/1\ТГ - -!- - 1 1

и- = И- => \vtn-1 = И?п-1 => /2 = \vin-l . (4.16)

Соотношение (4.6) с учетом (4.15)—(4.16) определяет вид автомодельное решения

и> = Л^Ф(г), 1 = хГ1/2. (4.17)

Подставив (4.17) в (4.13), после несложных преобразований приходим к обыкновенному дифференциальному уравнению

—!— Ф - -гФС = аФ" + М>". 1-й 2

Замечание: Преобразование (4.8) часто удобнее записывать в виде

г = Сх х = Скх, w = CnW (С > 0),

а затем определять значения постоянных к и т, при которых сохраняется вид исходного уравнения (С — свободный параметр). Соответствующие таким преобразованиям точные решения называются автомодельными [6-8].

Примеры некоторых автомодельных решений, указанные в [9-11], приведены в таблице 2 (используется сокращение ОДУ — обыкновенное дифференциальное уравнение).

Таблица 2

Некоторые нелинейные уравнения математической физики, которые допускают автомодельные решения

Уравнение Его название Вид решения Итоговое ОДУ

£=&[/(><] Нестационарное уравнение теплопроводности п> = г = хг1'2 [Ди>)м/]' + \Ы = 0

Уравнение Бюргерса и- = г1!2 и(г), г = хг1'2 аи" + Ъии'+ + \ги' + ¿и = 0

г/ду>\д2у> дг ~}Удх> дх2 Уравнение фильтрации и- = ¿/^(г), г = хг1'2 2/(и')и" + 1и' -и = 0

д2ж _ д г Волновое уравнение п> = г = х^ (г2и/)' =

д2м> , д2м> _ п Уравнение теплопроводности с источником ■м = хАи(г), г = у/х (1 + 12)и"-Щ^ги'+ + Ш^и-аи" = 0

д2\\> д2\\> _ п 1 Ти 1 <-, о — V дх1 ду ду1 Уравнение околозвукового газового потока м = х~зк~2и(£), ь г = ту, к — любое к+1 и и + к+11 и - Ъкги' + 3(3к + 2)и = 0

Уравнение Кортевега -де Фриза = г2/3и(г), г = хг1/3 аи'" + Ъии'+ + ¿ги' + §и = 0

д\у д2\у _ ду дхду _ д2\у _ 1'.' дх ду2 ду3 Уравнение пограничного слоя ж = х1+1и(г), г = хку, X — любое (2Х + 1)(и')2- - (X + \)ии" = аи'"

Аналогичный подход с успехом может использоваться для построения точных решений систем уравнений с частными производными.

Пример 3: Рассмотрим систему уравнений стационарного ламинарного гидродинамического пограничного слоя на плоской пластине:

ди ди д2и дх ду

Сделаем в (4.18) растяжение независимых и зависимых переменных по правилу

х = Сх, у = Ску, и = С" и, V = С V. (4.19)

Умножив полученные уравнения на подходящие постоянные множители, имеем

-д^+Сп-т-к+1^= с-т-2к+1уд2 Й

д_х д у д (4.20)

ди +сп-т-к+1^_=0 У '

дх д у

Потребуем, чтобы вид уравнений преобразованной системы (4.20) совпал с видом уравнений исходной системы (4.18). Это условие дает два линейных алгебраических уравнения: n - m - k + 1 = 0, -2k - m + 1 = 0. Разрешив их относительно m и n, получим

m = 1 - 2k, n = -k, (4.21)

где показатель k может быть выбран произвольно. Подставим (4.21) в (4.19):

х = Cx, y = Ôy, и = Ô~2ku, v = Ckv.

Находим отсюда инварианты путем исключения параметра C:

У = )ку => ух~к = У*к = h= %,\ v х'

и = (^]~2ки => u£k~l = ulk~l = 1г\ V х'к

V = V => vi = vl = I3. Решение ищется в виде I2 = U(I1), I3 = V(I{), т.е.

их2k-1 = U(Z), vxk = V(Z)

или

u(x, y) = x1~2kU (Z), v(x, y) = x~kV(Ç), Z = yx~k, (4.22)

где k — произвольная постоянная. Подставив (4.22) в исходную систему (4.18), получим систему обыкновенных дифференциальных уравнений для функций U = U(Z),

V = V (Z):

U[(1 - 2k)U - kZUZ] + VUZ = vUZZ, (1 - 2k)U - kZUZ + VZ = 0.

Пример 4: Рассмотрим уравнение Клейна-Гордона с нелинейностью экспоненциального вида

d2w 32w 1

= + <4 23»

Ищем инвариантное преобразование в виде комбинации преобразований сдвига и масштабирования:

x = ax, t = bt, w = w + c. (4.24)

Подставив (4.24) в (4.23) и умножив на b2, получим

д2 w b2 д2 w ,9 w , w

Чтобы это уравнение совпало с (4.23), надо положить

i, = 1, ъ2еи = 1. a2

Отсюда имеем

2

а = Ь, с = - - 1п й X

и преобразование (4.24) принимает вид

х = Ъх, ? = ¿Г, м>=м>--\п.Ь, (4-25)

л

где Ь > 0 — произвольная постоянная.

Исключая из соотношений (4.24) параметр Ь, находим инварианты

хх 2 — 2 — ,

/1 = - = —, /2 = и-+ - = и-+-1ги. (4.26)

г г л л

Соотношение (4.6) с учетом (4.26) определяет вид автомодельного решения

2х = Ф(г) - -1т, г = ~. (4.27)

л г

Подставляя (4.27) в (4.23), получим обыкновенное дифференциальное уравнение

2

Замечание: Точно таким же образом ищется точное решение стационарного уравнения теории горения с экспоненциальным тепловыделением

д2™ д2™ 1 + =

дх2 ду2

Пример 5: Рассмотрим нелинейное уравнение нестационарной теплопроводности

д™ д , п д™ ч

Ищем инвариантное преобразование в виде комбинации преобразований сдвига и масштабирования:

х = ах, г = г + Ь, ™ = с™. (4.29)

Подставив (4.29) в (4.28) и поделив на с, получим

дм с" д/-рдм\ дг а2 д х д х'

Требование совпадения с уравнением (4.28) дает одно условие

сп

— = 1 => а = с"/2, а2

где с — произвольная постоянная. При этом другой параметр Ь также остается произвольным, т.е. полученное инвариантное преобразование

х = сп/г х, г = х+ Ь, ™ = с™ (4.30)

является двухпараметрическим.

Поскольку для применимости используемого метода требуется однопарамет-рическое инвариантное преобразование положим

Ь = /(с), (4.31)

где вид функции / будет определяться далее. Исключая параметр с из двух последних соотношений (4.30) с учетом (4.31), имеем

Это равенство будет определять инвариант 1(г, = 1(г, если выполняется соотношение

Д™ \ —

■=) = ф(м0-ф(м0, (4.33)

которое можно рассматривать как функциональное уравнение. Нетрудно проверить, что уравнение (4.33) имеет решение

/(с) = А 1п с, (4.34)

где А — произвольная постоянная. Подставив (4.31) в (4.30) с учетом (4.34), получим

х = сп/2 X, г = г + А 1п с, ™ = с№. (4.35)

При А ф 0 выразим из второго соотношения (4.35) параметр с через г и подставим в оставшиеся соотношения. В результате находим инварианты

/, = хека = х«К /2 = ыет = где к = - — .

2А

Используя теперь формулу (4.6), определяем вид точного решения уравнения (4.28):

™ = е~тФ(1), I = хекпг. (4.36)

Подставив (4.36) в (4.28), получим обыкновенное дифференциальное уравнение для функции Ф(г):

-2кФ + кп1Ф[ = а(н>Х);.

Замечание: Преобразование (4.30) при Ь = с сохраняет вид уравнения (4.28), но не имеет инвариантов. Отсюда следует, что не каждое преобразование, сохраняющее вид уравнения, имеет инварианты.

Краткие выводы

Описанный метод исследования нелинейных уравнений с частными производными является естественным обобщением метода подобия [6,8] и использует идеи метода группового анализа, но значительно проще последнего. Содержание разд. 4.2 (основанного на результатах [7,10-14]) на начальном этапе может излагаться как общий эвристический подход, эффективность которого проверяется на конкретных примерах. Для успешного использования данного метода надо уметь решать лишь простейшие алгебраические уравнения и уметь дифференцировать (и не требуется понимать и использовать сложную терминологию группового анализа дифференциальных уравнений).

Таблица 3

Инвариантные решения, которые могут быть получены путем использования комбинаций преобразований сдвига и масштабирования, сохраняющих вид уравнений (С — произвольная постоянная, С > 0)

№ Инвариантные преобразования Вид инвариантных решений Примеры уравнений

1 1 = 1+ Ск, х = х+ СХ уу = и (г), г = кх-Х1 !г = т^т

2 1 = а, х = Скх, № = Ст,Н> ™ = 1 = хгк Уравнения см. в таблице 2

3 1 = 1 + 1п С, х = Скх, № = Ст,Н> ™ = 1 = хе~кг (к = ^тп, т — любое)

4 1 = а, х = х + к\пС, № = Ст,Н> ™ = 1 = х-к\п1 ~ ~ аТх{™ —) (т = -1/п, к — любое)

5 1 = а, х=С$х, ж = ж + а 1п С w = и{г) + а\т, г = хг13 ды _ д /'.'.' \ д1 ~ дх\с дх> (а = 2(3 - 1, (3 — любое)

6 1 = а, х = х + (31п С, ж = ж + а 1п С = и(1)+а\т, г = х-|31ги /• д2У! \2 д2У! д2У! Л \dxdt' дх2 дг2 ~ и (а, (3 — любые)

7 1 = 1 + С, х = х + СХ, ж = м> + Ск ц; = 11(1) + Ы, 1 = Х-Х1 /у- ду>\д2у> дг ~ J V дх > дх2 (к, X — любые)

8 1 = 1 + 1п С, х = х+кЫС, ж = Ст ж и' = етги(1), 1 = х-Ы /• д2\ч \2 д2\ч д2\ч _ Удхдг> дх2 дг2 ~ и [к, т — любые)

Общие выводы и замечания

Описанный метод с успехом может использоваться для решения самых различных типов математических уравнений. В случае обыкновенных дифференциальных уравнений и уравнений математической физики данный метод намного проще для понимания, чем классический метод группового анализа [3-5,7,10-14], и может трактоваться как его частный случай. При использовании описанного простого метода надо уметь решать лишь самые простые алгебраические уравнения (и системы) и уметь дифференцировать. Важно отметить, что эффективность используемого в статье метода для обыкновенных дифференциальных уравнений практически не уступает классическому методу группового анализа [4,5], а для нелинейных уравнений математической физики позволяет находить все наиболее распространенные решения (которые по самым грубым оценкам составляют не менее половины решений, которые можно найти с помощью метода группового анализа; эту оценку можно получить путем исследования представленных

в [9-11] результатов). Указанные обстоятельства говорят о целесообразности введения данного метода в стандартные курсы лекций по обыкновенным дифференциальным уравнениям и уравнениям математической физики, которые читаются студентам, специализирующимся в области прикладной математики, физики и механики (этот метод с успехом можно включать также в соответствующие специальные курсы, читаемые в некоторых технических и педагогических вузах). Кроме того, этот простой метод полезно излагать на предварительном этапе (в качестве введения) до чтения основного материала в существенно более сложных для восприятия спецкурсах по методам группового анализа дифференциальных уравнений.

Изложенный в данной статье материал представляет собой замкнутый цикл лекций, который является первой частью специального курса по нелинейным уравнениям математической физики и механики, читаемого студентам кафедры "Прикладная математика" Московского государственного технического университета имени Н.Э.Баумана.

Литература

[1] Болтянский, В.Г. Симметрия в алгебре / В.Г.Болтянский, Н.Я. Вилен-кин. 2-е изд. - М.: Наука, 2002. - 240 с.

[2] Кудряшов, Н.А. Симметрия алгебраических и дифференциальных уравнений / Н.А. Кудряшов // Соросовский образовательный журнал. -1998. - №9. - С. 104-110.

[3] Ибрагимов, Н.Х. Азбука группового анализа / Н.Х.Ибрагимов. - М.: Знание, 1989. - 48 с.

[4] Ибрагимов, Н.Х. Опыт группового анализа обыкновенных дифференциальных уравнений / Н.Х.Ибрагимов. - М.: Знание, 1991. - 48 с.

[5] Polyanin, A.D. Handbook of Exact Solutions for Ordinary Differential Equations, 2nd Edition / A.D. Polyanin, V.F. Zaitsev. - Boca Raton: Chapman & Hall/CRC Press, 2003. - 816 pp.

[6] Баренблатт, Г.И. Подобие, автомодельность, промежуточная асимптотика / Г.И. Баренблатт. - Л.: Гидрометеоиздат, 1978. - 218 с.

[7] Полянин, А.Д. Методы решения нелинейных уравнений математической физики и механики / А.Д. Полянин, В.Ф.Зайцев, А.И. Журов. -М.: Физматлит, 2005. - 256 с.

[8] Седов, Л.И. Методы подобия и размерности в механике / Л.И. Седов. - М.: Наука, 1972.

[9] Полянин, А.Д. Справочник по нелинейным уравнениям математической физики: Точные решения / А.Д. Полянин, В.Ф.Зайцев. - М.: Физматлит. 2002. - 432 с.

[10] Ibragimov, N.H. CRC Handbook of Lie Group to Differential Equations, Vol. 1. Symmetries, Exact Solutions and Conservation Laws / N.H. Ibragimov (Editor). - Boca Raton: CRC Press, 1994. - 432 pp.

[11] Polyanin, A.D. Handbook of Nonlinear Partial Differential Equations / A.D. Polyanin, V.F. Zaitsev. - Boca Raton: Chapman & Hall/CRC Press, 2004. - 840 pp.

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

[13] Дородницын, В.А. Симметрия в решениях уравнений математической физики / В.А.Дородницын, Г.Г. Еленин. - М.: Знание, 1984. - 64 с.

[14] Polyanin, A.D. Handbook of Mathematics for Engineers and Scientists / A.D. Polyanin, A.V. Manzhirov. Boca Raton: Chapman & Hall/CRC Press, 2007. - 1544 pp.

Поступила в редакцию 16/ VIII/2008;

в окончательном варианте — 16/VIII/2008.

ELEMENTARY THEORY OF USING INVARIANTS FOR SOLVING MATHEMATICAL EQUATIONS

© 2008 A.D. Polyanin2

The paper describes a general scheme for studying mathematical equations. The scheme is based on the use of invariants and permits one to simplify algebraic equations and systems, reduce the order of (or integrate) ordinary differential equations, and also obtain exact solutions of nonlinear equations of mathematical physics. The invariants are constructed by a search for transformations preserving the form of the equations in question. (The notions and complicated techniques of group analysis are not used here.) Numerous examples of solution of specific algebraic and differential equations are given. It is important to note that the application of the simplest scaling and translation transformations permits one to give a unified description of a substantially larger number of solvable ordinary differential equations (or equations admitting order reduction) than in classical and special textbooks. To use the method, one should only be able to solve the simplest algebraic equations and compute derivatives. This shows that it is advisable to include the method not only in standard lecture courses on ordinary differential equations and equations of mathematical physics for students specializing in applied mathematics, physics, and mechanics but also in special courses read in some technical and pedagogical higher education institutions.

Keywords and phrases: invariant, algebraic equation, differential equation,

group analysis.

Paper received 16/VIII/2008.

Paper accepted 16/VIII/2008.

2Polyanin Andrey Dmitrievich ([email protected]), Ishlinsky Institute for Problems

in Mechanics of the Russian Academy of Sciences, Moscow, 119526, Russia.