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

УДК 517.9

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

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

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

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

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

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

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

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

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

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

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

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

X + а X + 1 = 0

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

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

4

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

х4 + ах1 + 1 = 0.

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

у'1х — = 0 (1.2)

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

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

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

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

такое же уравнение

Х х— у = °.

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

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

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

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

I = у — х = X — X.

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

^2 = у х = Ъ 13 = ухх = ~Ххх ....

Пример 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 + а2 х2 + а0 = 0, (2.1)

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

Замена

х = — х (2.2)

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

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

х2 = X.

(2.3)

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

а2пгп + а2п-2 гп-1 + «2п-4 гп-2 + ••• + а4г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 + сХ + Ьх + а = 0.

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

а(х2 н—+ Ь(х н—) + с = 0.

х2 х

(2.7)

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

х х2 х

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

аі2 + Ьі + с - 2а = 0.

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

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

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

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

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

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

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

Нетрудно проверить, что преобразование

Р2и+\( х) = 0,

Р2и+\( х) = (х + 1)Q2n( х),

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

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

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

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

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

а0 х2п + а1 х2п 1 + ••• + ап-1 хп+1 + апхп+

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

х=

(2.9)

х

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

г = и— = х + — . (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Л2)

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

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

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

(2.13)

получим

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

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

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

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

22

х2+у2 = а’ (214)

х3 + у3 = Ь. (2.14)

Переходя к переменным (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. Преобразования, сохраняющие вид уравнения, и их инварианты

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

р(х> у у'х,---,у(х1)) = 0 (з.1)

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

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

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

р (х, у, ту,..., Уу?) = 0. (3.3)

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

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

у = 8(х) (3.4)

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

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

у = 8( х). (3.5)

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

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

у2х - Ух = о (3.6)

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

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

У =

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

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

у = Аех + Ь, А = в~а,

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

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

1( х У, у'х) =1 (х, у

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

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

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

На первом этапе ищется преобразование

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

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

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

II = 1\ (х, у), 12 = 12(х, у, у'х)■ (3.9)

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

и = 12, z = 1\, и = и(г).

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

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

z = II, z = z( х).

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

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

Процедура определения инвариантов

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

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

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

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

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

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

, А2~, ,, А2~г/ (п) А2—(п) /пил

Ух = ~гУъ УXX = ~аУ~х~х Ух - (3-п)

А 1 А 1 А 1

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

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

Мх - В1 \

А )■

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

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

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

ух = р(х, у'х). (3.12)

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

х ух.

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

и = у'х, z = х, и = и^).

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

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

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

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

у> ух.

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

и = у'х, z = у, и = и(£).

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

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

у"=^(х,^) (3.14)

не меняется при растяжении зависимой переменной: у ау. При этом две ком-

бинации из трех величин х, у, у'х остаются остаются неизменными:

Ух

х, —.

у

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

ух / ч

и = —, 1 = X, и = и(1).

у

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

и'х = У^-(Ц2 = У-^-и\

х у у у

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

и'х = р(х, и) - и2.

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

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

уу”х - (у'х)2 = ky3eXx■ (3.15)

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

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

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

у!х- (Ух)2 = ae}■bkXie}■х

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

аекЬ = 1 => Ь = --\па. (3-17)

Л

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

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

у = у(Л(х-х> у^^ = у<Л~х.

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

Р = уeЛx. (3.18)

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

ух = аЪ

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

у' у

— = 4? = /2. (3.19)

у у

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

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

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

и'х=У^-(Ц = У^-и\ (3.21)

у Ку* у

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

и'х = = 0^*У + екХУхЖ = + е1ХУу)иг = + ги)“г (З-22)

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

у" о у" т

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

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

(Л + и)и'1 = k■

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

У^ = ур(^у£\

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

и2 + ^ + zu)u'z = р^, и).

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

Ух \а2х + |32;у + У2>

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

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

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

сдх + Р1У + У1 _ а(щ х + в у) + сцЬ + Р1С + 71

а2х + (Згу + 72 а(а2 х + у) + оф + (32с + у2' '

/ = 7«і£+Ріу + уі \ , .

-1 \а2х+|32у + у2/

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

а1Ь + в с + у1 = ауь а2Ь + &с + 42 = ау2. (3.26)

Эти соотношения можно рассматривать как систему двух линейных алгебраических уравнений первого порядка для определения коэффициентов Ь и с. При

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

г = а1Х + + , где г = г(х). (3.27)

а2 х + в2у + 42

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

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

(а2Р1 - аф2)х + (31у2 - |32у1 , _ а2г - а!

(Ргг-РО2 ^ ' т + (32г - РГ

В таблице 1 представлены некоторые обыкновенные дифференциальные уравнения второго порядка, допускающие понижение порядка путем использования простейших инвариантных преобразований. Для уравнений первого порядка (когда функция р(и, V, м>) не зависит от третьего аргумента) уравнения, указанные в таблице 1, решаются путем перехода от у к новой зависимой переменной z = Ь(х,у), где ^ —первый инвариант.

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

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

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

Таблица 1

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

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

Р(х,/х,ух'х) = 0 У=У + а 1\ - X Ь =УХ

р(у,У„У~) = о х = х + а \=У 12=У'Х

F(ax + |3у + у, у'х, ухх) = 0 х = х + а|3, у = у-аа 1\ = ах + (Зу + у ь = ух

У у" F їх,—,— =0 У у У = ау 1\ - X Ь = ~ У

F(y,Лy',x2y") = 0 х- ах Ь= У Ь = хух

Риу,ук,Щ=0 х = х — - 1п а, У = ау 1у = екху 1г = ^ У

Р(екху,екхУх,екхухх) = 0 - 1, х = х - - та, к-У = аЗ 1у = екХу 12 = екху'х

Р(хеку, ху'х, хгу”х) = 0 х = ах, У=У~ї1па 1у = хекУ 12 = ХУ'х

F(x/^y,x/^+1y',x/^+2y") = 0 х - ах, у - аку 11= Ху 12 = хк+1у'х

/ , г2!;" \ р хпут,ху\ Ухх =0 \ У х = атх, у а ' у 1 = х”/1 1 _^Ух 2 у

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

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

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

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

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

= 0 (42)

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

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

F( х,(, —=, —, ) = 0. (4.3)

' дх <9г дИ? дхдс дР- '

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

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

I(х, г, w) = 1( х,г, w).

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

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

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

На первом этапе ищется преобразование

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

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

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

II = II(х, у, w), 12 = 12(х, у, w). (4.5)

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

менных и точное решение уравнения (4.2) ищется в виде

!2 = ф(!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)

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

г = г + Ск, х = х + Ск, w = w, (4.10)

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

х - хг - г 1 . -

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

к к

Вторым инвариантом здесь является ^ = 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 д^ \

а? *•«?•**•*•■■■) = “■ <412>

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

F(Ф, кФ', -кФ', к2Ф'', -ккФ'', к2Ф'', ...) = 0.

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

дw д1^

-=<,^ + Р,,Л (4.18)

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

д w b д2 w 1 —

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

А = 1, ъсп~х = 1.

а2

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

а = л/b, с = bi-n.

Подставляя эти величины в (4.8), находим инвариантное преобразование

х= V& х, t = bt, w = bi-n w, (4-14)

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

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

х=ф1/2х => => h=xrll2. (4.15)

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

/1\тг- -4- -1 1

w = \~j "w => wtn~1 = vrfn-i => /г = wtn-i. (4.16)

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

w = ?1-»Ф(г), z = xt~1^2. (4-17)

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

—Ф - -гФ' = аФ" + ЬФп.

1 -п 2 zz

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

t = Ct, x = CkX, w = CnW (C > 0),

а затем определять значения постоянных к и m, при которых сохраняется

вид исходного уравнения (C — свободный параметр). Соответствующие та-

ким преобразованиям точные решения называются автомодельными [6-8].

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

Таблица 2

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

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

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

Уравнение Бюргерса Ы = г1/2и(г), г = хг1'2 аи" + Ьии'+ + 1 г г/ + \и = 0

д\\> _ г( д\\> \ 32уу дг ~ •* V дх / дх2 Уравнение фильтрации и- = ^2и(г), г = хг4'2 2/(и')и" + їй' - и = 0

д2и> _ д г £ґЛА,\дпі дг1 ~ дхУ^ дх і Волновое уравнение И1 = г = х^ (г2м/)' =

+ II а 3 Уравнение теплопроводности с источником ж = хіі и(г), г =у/х (1 + 12)и" - 2^іи' + + (і_„)2и аи - 0

о и с? К + Ой|°й Уравнение околозвукового газового потока ’и> = х~зк~2и(і), 1 = Ху, к — любое МИМ + £+1г “ - 5£гг/ + 3(3£ + 2)и = 0

дю _ „дъу/ . илл,дю -дї-аи+Ьц;-дї Уравнение Кортевега -де Фриза и1 = Г2//Зи(г), г = хг1/3 аи’” + Ъии'+ + |гг/ + |и = 0

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

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

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

ди ди д2и

<«»>

Э~Х*Э1~

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

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

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

+ с„_т_,+1-£и = с_т_2к+ 1у£и дх ду дХ (4.20)

ди п_т_к+1дУ_ у >

дх ду

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

получим

т = 1 — 2к, п = —к, (4.21)

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

х = С X, у = С^у, и = С1~2ки, V = Ск X.

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

у = (*)*у ух~к = у~хк = /!=<;;

4 X7

и = (-) и => и£к~1 = иУ-1 = /2;

4 Х'к

V = (|) V => у^=у1 = /з.

Решение ищется в виде 12 = и(Т\), 1з = V(/1), т.е.

их2к—1 = и&), ухк = V(Z)

или

и(х, у) = х1~2ки (£), у(х, у) = x~kV (О, £ = ух~к, (4.22)

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

V = V©:

и [(1 — 2к)и — к^Щ ] + VU^Q = уи",

(1 — 2к)и — к^Щ + VI. = 0.

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

д2^ cP^w т

^ = + ^ <4 23»

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

х = ах, Г = ЙХ w = w + с. (4.24)

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

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

4 = 1, й Vе = 1.

а2

Отсюда имеем

2

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

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

х = Ьх, Г = Ь, м>=м>-—\п.Ь, (4-25)

Л

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

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

хх 2 — 2 — ,

1\ = - = —, І2 = и1 + - 1п Г = и1 + - 1п Г. (4.26)

її Л Л

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

2х н-= Ф(г) - - 1пГ, г =-. (4.27)

Л ї

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

2

(г2Ф^+-=аФ"+(3/ф.

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

д2н д2н т

и УУ _ 1н,

Эх2 ду2 1 '

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

дн д , п дн \

а = &)' (4'28)

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

х = ах, Г = I + Ь, н = сн. (4.29)

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

д н сп д ,-п д н

дГ =

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

сп

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

а2

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

х = сп/г х, Г = Г + Ь, н = сН (4.30)

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

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

Ь = /(с),

(4.31)

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

' = /(’,)■ (4-32>

Это равенство будет определять инвариант I(t, w) = I(t, w), если выполняется соотношение

4w \ —

■=) = ф(н’)-ф(н’), (4.33)

w

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

f(c) = A ln с, (4.34)

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

х = cn/2 х, t = t + A ln с, w = cw. (4.35)

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

/| = хекм = х£ш, /2 = we2kt = wikt, где к = ——.

’ Л 2А

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

w = в~2к,Ф{1), z = xeknt. (4.36)

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

для функции Ф^):

-2кФ + knzФZ = a(wnwZ)Z.

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

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

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

Таблица 3

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

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

1 1 = 1+ Ск, х= х+ СХ ж = и (г), г = кх-Х1 £ = £[/(><]

2 Г = Сг, х = Ск х, ж = Стй> и- = Ги(1), 1 = хгк У равнения см. в таблице 2

3 1 = 1 + 1п С, х = Ск х, ™ = Стч; \у = етЧ7(1), 1 = хе~кг дм _ д (..п ЗиЛ ТЙ ~ д!) (к = |тп, т — любое)

4 г = Сг, х = х + к\пС, ж = Стй> ™ = 1ти{£), 1 = х-к\п1 д\\> _ „ д (,,,п <9иЛ :ш' ~ Тх) (т = -1/и, к — любое)

5 г = Сг, х = С& х, ж = ж + а 1п С и1 = и {г) + а 1п Г, г = хГ$ д\\> _ д (р\ч <9и> \ дг ~ дхУ^ дх) (а = 2(3 - 1, (3 — любое)

6 г = Сг, х = х + (31п С, ж = ж + а 1п С 1л> = и(1) + а1п?, г = х-|31ги (дгу/ \2 дгу/ дгу/ Г) У'дхдг> дх2 812 " и (а, (3 — любые)

7 ? = ? + С, х = х + СХ, ж = и* + Ск и1 = и (г) + Ы, г = х-Х1 д\у _ /?/ дм \ д2\у дг ~ ^ дх' дх2 (к, X — любые)

8 1 = 1 + 1п С, х = х+к\пС, ж = Стй> ы = етги(1), 1 = х-Ы /• д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.