Формула интегрирования по частям для интегралов от функций одной переменной известна всем со школьной скамьи. В вузах она, как правило, изучается в координатном виде в\mathbb{R}^mв курсах математического анализа и УРЧП. В этой статье мы приведем инвариантную геометрическую версию этой формулы и отметим некоторые конструкции теории динамических систем, с которыми она естественным образом связана.

Напомним некоторые определения.

Возьмём гладкое многообразие M с локальными координатами x=(x^1,\ldots,x^m) и зададим на нём векторное поле v(x)=(v^1,\ldots, v^m)(x) (все объекты далее мы считаем гладкими).

Будем считать, что для t \in (-\varepsilon, \varepsilon) определено семейство диффеоморфизмов g^t: M \to M, удовлетворяющее следующей задаче Коши:

\frac{d}{dt}g^t(x) = v(g^t(x)), \quad g^0(x) = x.

Это семейство называется фазовым потоком поля v.

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

g^t(x) = (g^t_1(x), \dots, g^t_m(x)).

Производной Ли вдоль поля v от дифференциальной k-формы

\omega = \sum_{i_1 < \dots < i_k} \omega_{i_1 \dots i_k}(x) dx^{i_1} \wedge \dots \wedge dx^{i_k}

называется k-форма

L_v \omega =  \left. \frac{d}{dt} \right|_{t=0} \sum_{i_1 < \dots < i_k} \omega_{i_1 \dots i_k}(g^t(x)) d(g^t_{i_1}(x)) \wedge \dots \wedge d(g^t_{i_k}(x)).

Кроме линейности, оператор L_v обладает ещё рядом замечательных свойств, которые мы будем восстанавливать по мере надобности, а сейчас отметим лишь одно: на 0 -форму, т.е. функцию f:M\to\mathbb R , производная Ли действует так:

L_vf=v^i\frac{\partial f}{\partial x^i}.

По повторяющимся индексам ведётся суммирование.

Итак, для функций производная Ли — это производная в силу системы ОДУ с векторным полем v.

Нам ещё понадобится оператор внутреннего дифференцирования i_v. Этот оператор действует на k-формы,k>0. Он по определению линеен и задаётся на базисных k-формах следующим образом:

i_v(dx^{i_1} \wedge \dots \wedge dx^{i_k}) = \sum_{s=1}^k (-1)^{s+1} v^{i_s} dx^{i_1} \wedge \dots \wedge \widehat{dx^{i_s}} \wedge \dots \wedge dx^{i_k}.

Символ \widehat{\dots} означает, что соответствующий сомножитель пропущен.

Основная теорема.

Теорема. Предположим, что на M задана m-форма

\psi=\rho(x)dx^1\wedge\ldots\wedge dx^m,

причём

L_v\psi=0.\qquad (1)

(Об этом интересном предположении поговорим подробнее ниже.)

Пусть Q\subset M — компактное m-мерное подмногообразие с гладкой границей \partial Q, которое получается в результате замыкания открытого множества.

Тогда для любых гладких функций f,g:M\to \mathbb{R} верна формула:

\int_Q (L_vf)g\psi=\int_{\partial Q}fg i_v\psi-\int_{Q}f(L_v g)\psi.\qquad(2)

Это и есть инвариантная дифференциально-геометрическая версия формулы интегрирования по частям.

На самом деле в этой теореме достаточно, чтобы все объекты f,g,\psi,v были определены во внутренности Q и вели себя регулярно во всём Q вплоть до границы включительно.

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

L_v\psi=\frac{\partial (\rho v^i)}{\partial x^i}dx^1\wedge\ldots\wedge dx^m.

Чаще всего в учебниках по анализу и УРЧП формула (2) встречается для случая:

M=\mathbb{R}^m,\quad v=\frac{\partial}{\partial x^j},\quad \psi=dx^1\wedge\ldots\wedge dx^m.

Тогда она приобретает вид:

\int_Q \frac{\partial f}{\partial x^j}gdx^1\wedge\ldots\wedge dx^m=\int_{\partial Q}fg n_jdS-\int_{Q}\frac{\partial g}{\partial x^j}fdx^1\wedge\ldots\wedge dx^m,

где n_j — компонента вектора внешней единичной нормали к поверхности \partial Q, а dS — элемент площади на этой поверхности.

Равенство (1) имеет следующий динамический и геометрический смысл. Предположим, что \rho \ge 0 на M. Тогда на M определено понятие объёма (меры) множества:

\mu(D)=\int_D\psi.

Если равенство (1) выполнено, то поток g^t сохраняет объём любого (измеримого по Лебегу) множества:

\mu(D)=\mu(g^t(D)),\quad \forall t\in\mathbb{R}.

Отметим, что если многообразие Q не имеет границы (\partial Q=\emptyset), то в силу равенства (2) оператор L_v оказывается кососимметричным на пространстве гладких функций в смысле L^2_\mu(Q)-скалярного произведения.

Это весьма неслучайный факт. Предположим дополнительно, что g^t(Q)\subset Qи рассмотрим семействоU^t линейных операторов на пространствеL^2_\mu(Q):

(U^tf)(x):=f(g^t(x)).

В силу отмеченного выше свойства фазового потока сохранять объем, каждый из операторовU^tоказывается изометрией:

\|U^t\|=1.

Если функцияfгладкая, то

\frac{d}{dt}U^tf=L_v f,\quad U^0=\mathrm{id}.

В этом смысле мы можем написать

U^t=e^{tL_v}.

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

Доказательство теоремы.

Оператор L_v подчиняется правилу Лейбница:

L_v(fg\psi)=(L_vf)g\psi+f(L_vg)\psi+fgL_v\psi.

Последнее слагаемое справа равно нулю в силу (1).

Воспользуемся магической формулой Картана (Cartan’s magic formula):

L_v(fg\psi)=di_v(fg\psi)+i_vd(fg\psi).

Последнее слагаемое в ней равно нулю как внешний дифференциал от m-формы на m-мерном многообразии.

Для завершения доказательства нам остаётся проинтегрировать получившееся равенство

di_v(fg\psi)=(L_vf)g\psi+f(L_vg)\psi

по Q и воспользоваться формулой Стокса:

\int_Q di_v(fg\psi)=\int_{\partial Q}i_v(fg\psi).

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