Просмотр кода страницы Групповой анализ дифференциальных уравнений

'''Групповой анализ дифференциальных уравнений''' — раздел математики, изучающий свойства симметрии [[Дифференциальное уравнение|дифференциальных уравнений]] относительно различных преобразований зависимых и независимых переменных. Включает в себя методы и прикладные аспекты [[Дифференциальная геометрия|дифференциальной геометрии]], теории [[Группа Ли|групп]] и [[Алгебра Ли|алгебр Ли]], [[Вариационное исчисление|вариационного исчисления]] и является, в свою очередь, эффективным инструментом исследования в теории [[Обыкновенное дифференциальное уравнение|ОДУ]], [[Дифференциальное уравнение в частных производных|ДУЧП]] и [[Математическая физика|математической физике]].

== Мотивировка ==
Если дифференциальное уравнение после некоторой замены переменных переходит само в себя (с точностью до тождественных преобразований), то эта замена переводит любое решение уравнения снова в решение, вообще говоря, не совпадающее с исходным. Все такие замены образуют [[Группа (математика)|группу]], называемую [[Группы симметрии|группой симметрии]] дифференциального уравнения, или группой, допускаемой дифференциальным уравнением. Таким образом, знание группы симметрии и некоторых частных решений позволяет строить семейства решений, получаемые из исходных применением всех преобразований группы. Кроме того, если некоторое решение уравнения [[Инвариант (математика)|инвариантно]] относительно группы (или некоторой её [[Подгруппа|подгруппы]]), данный факт накладывает определённые условия на его вид, что позволяет ожидать упрощения исходного уравнения при его ограничении на такие инвариантные решения (в частности, уменьшения числа независимых переменных). Эти рассуждения приводят к задаче об общих методах нахождения допускаемой группы данного дифференциального уравнения. С другой стороны, по заданной группе преобразований в принципе может быть построено множество дифференциальных уравнений, допускающих её в качестве своей группы симметрии, что особенно актуально для фундаментальных разделов [[Теоретическая физика|теоретической физики]].

Хорошо развитые методы [[Теория групп|теории групп]] и дифференциальной геометрии позволяют придать изложенным соображениям строгие формулировки и конструктивно решить ряд сопутствующих задач, а также существенно расширяют арсенал средств для исследования качественного поведения решений дифференциальных уравнений, численного интегрирования и т.&nbsp;п.

== Определения ==
Пусть <math>x = (x^1, x^2, ... x^n) \in X = R^n</math> и <math>u = (u^1, u^2, ..., u^m) \in U = R^m</math> обозначают совокупности соответственно независимых и зависимых переменных некоторой системы дифференциальных уравнений порядка <math>s</math>

:{{EF|:|<math>F(x, u, \underset{1}{u}, ..., \underset{s}{u}) = 0,\qquad F=(F^1, F^2, ..., F^m),</math>|style=|ref=1|center=}}

а <math>\underset{k}{u} \in \underset{k}{U}</math> — совокупность всевозможных производных <math>u</math> по <math>x</math> порядка <math>k</math>. Система уравнений ({{eqref|1|}}) определяет в пространстве <math>\underset{s}{Z} = X\times U\times \underset{1}{U}\times ... \times \underset{s}{U}</math> некоторое подмногообразие <math>M</math>.

Пусть группа Ли <math>G</math> [[Действие группы|действует]] в пространстве <math>\underset{0}{Z}\equiv Z = X\times U</math> независимых и зависимых переменных преобразованиями

:{{EF|:|<math>x' = \phi(x, u, g),\quad u' = \psi(x, u, g), \quad g \in G.</math>|style=|ref=2|center=}}

Посредством пересчёта производных к преобразованным переменным преобразования ({{eqref|2|}}) единственным образом продолжаются на всё пространство <math>\underset{s}{Z}</math>:

:<math>{u'}^\alpha_{a_1...a_k} \equiv \frac{\partial u'^\alpha}{\partial x'^{a_1}...\partial x'^{a_k}} = \psi^\alpha_{a_1...a_k}(x, u, \underset{1}{u}, ..., \underset{k}{u}, g),\quad k=1, ..., s.</math>

Группа <math>G</math> называется '''группой симметрии системы ({{eqref|1|}})''', если многообразие <math>M</math> является инвариантным многообразием <math>s</math>-го продолжения действия ({{eqref|2|}}), то есть действия ({{eqref|2|}}), продолженного на производные до порядка <math>s</math> включительно. Действие каждой [[Однопараметрическая группа|однопараметрической подгруппы]] <math>e^{t A}</math>, <math>A\in\mathcal{G} = \mathrm{Lie}\,G</math> (см. {{iw|Экспоненциальное отображение (теория Ли)|экспоненциальное отображение|en|Exponential map (Lie theory)}}) группы <math>G</math> в пространстве <math>X\times U</math> [[Генераторы группы|генерируется]] векторным полем (здесь и ниже подразумевается [[Соглашение Эйнштейна|правило суммирования Эйнштейна]])

:{{EF|:|<math>Y_A = \xi^a(x, u)\frac{\partial}{\partial x^a} + \eta^\alpha(x, u)\frac{\partial}{\partial u^\alpha},\quad \xi^a = \left.\frac{d\phi^a(x, u, e^{t A})}{d t}\right|_{t = 0},\quad\eta^\alpha = \left.\frac{d\psi^\alpha(x, u, e^{t A})}{d t}\right|_{t = 0}.</math>|style=|ref=3|center=}}

Соответствующий генератор действия подгруппы <math>e^{t A}</math>, продолженного на пространство <math>\underset{s}{Z}</math>,

:{{EF|:|<math>\underset{s\ \ \ }{Y_A} = \xi^a(x, u)\frac{\partial}{\partial x^a} + \eta^\alpha(x, u)\frac{\partial}{\partial u^\alpha} + \sum\limits_{|I| = 1}^s\theta^\alpha_I(x, u, \underset{1}{u}, ..., \underset{|I|}{u})\frac{\partial}{\partial u^\alpha_I},\quad \theta^\alpha_I = \left.\frac{d\psi^\alpha_I(x, u, \underset{1}{u}, ..., \underset{|I|}{u}, e^{t A})}{d t}\right|_{t = 0},</math>|style=|ref=4|center=}}

где <math>I</math> — [[мультииндекс]], называется '''<math>s</math>-м продолжением''' генератора <math>Y_A</math>. По аналогии путём формального прибавления к ряду ({{eqref|4|}}) неограниченного числа слагаемых с производными высших порядков вводится понятие '''бесконечного продолжения''' <math>\underset{\infty\ \ \ }{Y_A}</math>. При этом вопрос о сходимости данного ряда не возникает, так как на практике всегда приходится иметь дело с функциями, зависящими от производных конечного порядка.

== Коэффициенты продолженных генераторов ==
Явный вид коэффициентов продолженного генератора находится дифференцированием ограничений

:<math>d{u'}^\alpha = {u'}^\alpha_a d{x'}^a,\quad d{u'}^\alpha_a = {u'}^\alpha_{ab} d{x'}^b</math>

и т. д., накладываемых на координаты в пространстве <math>\underset{s}{Z}</math>, по параметру преобразования <math>t</math> при <math>t = 0</math>. Например, для нахождения коэффициентов при <math>\partial/\partial u^\alpha_a</math> рассмотрим соотношения

:<math>d{u'}^\alpha = \left(\frac{\partial {u'}^\alpha}{\partial x^a} + \frac{\partial {u'}^\alpha}{\partial u^\beta} u^\beta_a\right) dx^a = {u'}^\alpha_a d{x'}^a = {u'}^\alpha_b\left(\frac{\partial {x'}^b}{\partial x^a} + \frac{\partial {x'}^b}{\partial u^\beta} u^\beta_a\right) dx^a.</math>

Приравнивая коэффициенты при <math>dx^a</math> и дифференцируя их по <math>t</math> при <math>t = 0</math>, с учётом выражений ({{eqref|3|}} — {{eqref|4|}}) имеем

:<math>\frac{\partial \eta^\alpha}{\partial x^a} + \frac{\partial \eta^\alpha}{\partial u^\beta} u ^\beta_a = \theta^\alpha_a + u^\alpha_b\left(\frac{\partial \xi^b}{\partial x^a} + \frac{\partial \xi^b}{\partial u^\beta} u^\beta_a\right),</math>

откуда

:<math>\theta^\alpha_a = D_a\eta^\alpha - u^\alpha_b D_a\xi^b,</math>

где введено обозначение

:<math>D_a = \frac{\partial}{\partial x^a} + u^\beta_a \frac{\partial}{\partial u^\beta} + u^\beta_{ab} \frac{\partial}{\partial u^\beta_b} + ...</math>

для оператора полной производной по координате <math>x^a</math>. Аналогичным образом могут быть найдены общие рекуррентное и явное выражения для коэффициентов произвольного порядка:

:{{EF|:|<math>\theta^\alpha_{aI} = D_a\theta^\alpha_I - u^\alpha_{bI} D_a\xi^b,\quad \theta^\alpha_I = D_I\left(\eta^\alpha - u^\alpha_b \xi^b\right) + \xi^b u^\alpha_{bI}.</math>|style=|ref=5|center=}}

Инфинитезимальным критерием инвариантности системы ({{eqref|1|}}) является условие

:<math>\left.\underset{s\ \ \ }{Y_A} F\right|_{F = 0} = 0,</math>

которое должно выполняться для любого элемента <math>A</math> из окрестности нуля алгебры Ли <math>\cal G</math>. Так как данное условие содержит не только переменные <math>x</math> и <math>u</math>, от которых зависят коэффициенты генератора <math>Y</math>, но и производные, вообще говоря, до порядка <math>s</math> включительно, которые в данном случае фигурируют как независимые переменные, при любых значениях которых условие обязано выполняться, то оно распадается на систему, как правило, переопределённых линейных дифференциальных уравнений на коэффициенты <math>\xi^a</math>, <math>\eta^\alpha</math>. Решив эту систему, можно в принципе восстановить (локальное) действие группы <math>G</math> в пространстве <math>X\times U</math>, а затем и в <math>\underset{s}{Z}</math>.

== Дифференциальные инварианты ==
'''Дифференциальным инвариантом порядка <math>s</math>''' группы <math>G</math> называется [[Дифференцируемая функция|дифференцируемая]] функция на <math>\underset{s}{Z}</math>, зависящая от производных до порядка <math>s</math> включительно, инвариантная относительно <math>s</math>-го продолжения действия этой группы. Дифференциальные инварианты порядка <math>s</math> удовлетворяют системе линейных уравнений первого порядка
:<math>\underset{s\ \ \ }{Y_i} J = 0,\qquad i = 1, ..., \dim G,</math>
где <math>Y_i</math> — базис генераторов группы <math>G</math> на <math>Z</math>. Из общей теории таких систем следует, что произвольный инвариант может быть выражен через некоторый минимальный набор <math>N - r</math> функционально независимых инвариантов, где <math>N = n + m\,C_s^{n + s}</math> — число независимых переменных и <math>r</math> — число независимых уравнений в системе, равное максимальному рангу матрицы её коэффициентов.

Значительная часть приложений группового анализа основаны на следующей теореме.

{{quotation|'''Теорема.''' Пусть <math>J_1</math>, <math>J_2</math>, ..., <math>J_{N - r}</math> — полный набор функционально независимых дифференциальных инвариантов порядка до <math>s</math> включительно группы <math>G</math>. Тогда система уравнений ({{eqref|1|}}) допускает <math>G</math> в качестве своей группы симметрий тогда и только тогда, когда она может быть записана в виде
:<math>\tilde F(J_1, J_2, ..., J_{N - r}) = 0</math>
с некоторыми новыми функциями <math>\tilde F</math>.}}

Таким образом, знание дифференциальных инвариантов позволяет найти общий вид уравнений, инвариантных относительно заданной группы, а анализ структуры алгебры Ли группы симметрий даёт возможность выбрать замену переменных, приводящую заданное уравнение к возможно более простому виду, например, допускающему понижение порядка (см. раздел «[[Групповой анализ дифференциальных уравнений#Приложения|Приложения]]»).

== Инвариантное дифференцирование ==
'''Оператором инвариантного дифференцирования''' группы <math>G</math> называется дифференциальный оператор, который при действии на дифференциальный инвариант этой группы даёт дифференциальный инвариант более высокого порядка. Из определения следует, что оператор  <math>\delta</math> является оператором инвариантного дифференцирования группы  <math>G</math> тогда и только тогда, когда он коммутирует с любым генератором <math>\underset{\infty}{Y}</math> продолженного действия этой группы:

:<math>\left[\delta,\,\underset{\infty}{Y}\right] = 0.</math>

Для любой группы  <math>G</math> преобразований пространства <math>Z = X\times U</math> существуют <math>n = \dim X</math> операторов инвариантного дифференцирования первого порядка, линейно независимых над [[Поле (алгебра)|полем]] инвариантов данной группы. Эти инварианты имеют вид  <math>\delta=\lambda^a D_a</math> и с учётом ({{eqref|5|}}) удовлетворяют системе уравнений

:<math>\underset{\kappa\ \ }{Y_i}\lambda^a = \lambda^b D_b\xi^a_i,\qquad i = 1, ..., \dim G.</math>

Число <math>\kappa</math> является наименьшим порядком продолжения группы <math>G</math>, ранг которого максимален, то есть равен <math>\dim G</math>. Поле дифференциальных инвариантов имеет конечный набор образующих в том смысле, что произвольный дифференциальный инвариант может быть получен конечным числом действий, включающих функциональные операции и применение операторов инвариантного дифференцирования первого порядка, из базиса дифференциальных инвариантов порядка <math>\kappa + 1</math>.

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

Ниже приведены основные результаты группового анализа применительно к ОДУ.

=== Понижение порядка ===
Если обыкновенное дифференциальное уравнение
:<math>F(x, u', ..., u^{(s)}) = 0</math>
допускает однопараметрическую группу симметрий с генератором
:{{EF|:|<math>Y = \xi(x, u)\frac{\partial}{\partial x} + \eta(x, u)\frac{\partial}{\partial u},</math>|style=|ref=5|center=}}
то переходом к переменным, {{iw|Теорема о выпрямлении векторного поля|выпрямляющим|en|Straightening theorem for vector fields}} векторное поле ({{eqref|5|}}), его порядок может быть понижен на единицу. В частности, уравнение первого порядка, разрешённое относительно производной, при этом условии интегрируется в квадратурах.

Последнее утверждение может быть сформулировано альтернативным образом в терминах интегрирующего множителя.

=== Интегрирующий множитель ===
Обыкновенное дифференциальное уравнение
:<math>P(x, u)\,dx + Q(x, u)\,du = 0</math>
допускает однопараметрическую группу симметрий с генератором ({{eqref|6|}}) тогда и только тогда, когда функция
:<math>\mu = \frac{1}{P\,\xi + Q\,\eta}</math>
является для этого уравнения [[Простейшие дифференциальные уравнения первого порядка#Интегрирующий множитель|интегрирующим множителем]].

=== Теорема Ли ===
Приведённые выше результаты обобщает следующая теорема.

{{quotation|'''Теорема (С. Ли).''' Если система <math>s</math> обыкновенных дифференциальных уравнений первого порядка допускает [[Разрешимая группа|разрешимую]] <math>r</math>-параметрическую группу симметрий, действующую регулярно и имеющую <math>r</math>-мерные [[Действие группы#Орбиты|орбиты]], то её общее решение можно получить из общего решения редуцированной системы <math>s - r</math> уравнений первого порядка. В частности, если <math>r = s</math>, то система интегрируется в квадратурах.}}

Ввиду соответствия между уравнениями порядка <math>s</math> и системами <math>s</math> уравнений первого порядка аналогичная теорема справедлива также для одного уравнения порядка <math>s</math>.

== Литература ==
* {{книга|автор=Л. В. Овсянников|заглавие=Групповой анализ дифференциальных уравнений|ссылка=|место={{М}}|издательство=Наука. Гл. ред. физ.-мат. лит.|год=1978|страниц=400|страницы={{{страницы|}}}|isbn=}}
* {{книга|автор=П. Олвер|заглавие=Приложения групп Ли к дифференциальным уравнениям. Пер. с англ.|ссылка=|место={{М}}|издательство=Мир|год=1989|страниц=639|страницы={{{страницы|}}}|isbn=5-03-001178-1}}
* {{книга|автор=Н. Х. Ибрагимов|заглавие=Группы преобразований в математической физике|ссылка=|место={{М}}|издательство=Наука. Гл. ред. физ.-мат. лит.|год=1983|страниц=280|страницы={{{страницы|}}}|isbn=}}

== Примечания ==
{{примечания}}

[[Категория:Дифференциальные уравнения]]
[[Категория:Теория групп]]
[[Категория:Математическая физика]]