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

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

дифференциальных уравнений
\[
\dot{x}=\Lambda x+\ldots
\]

Будем говорить, что набор (комплексных) собственных значений $\lambda_{1}, \ldots, \lambda_{n}$ матрицы $\Lambda$ является резонансным, если $\sum m_{i} \lambda_{i}=0$ при некоторых натуральных значениях $m_{i}$. Отметим, что при исследовании системы (8.1) (например, в теории нормальных форм) обычно используется более слабое условие резонағсности: $\sum m_{i} \lambda_{i}=0$ при некоторых целых $m_{i} \geqslant 0$ и $\sum\left|m_{i}\right|
eq 0$.

Предложение 10. Если набор $\lambda_{1}, \ldots, \lambda_{n}$ нерезонансный, то в малой окрестности точки $x=0$ уразнение (8.1) не имеет интегрального инварианта с аналитической плотностью.

Условие нерезонансности выполнено, например, в случае, когда $\operatorname{Re} \lambda_{i} \geqslant 0(\leqslant 0)$ и $\sum \operatorname{Re} \lambda_{i}>0(<0)$.

ДокаЗаТЕЛЬСТВО. Разложим плотность $M(x)$ в сходящийся ряд по однородным формам:
\[
M=M_{s}+M_{s+1}+\ldots, \quad s \geqslant 0 .
\]

Ясно, что функция $M_{s}$ является плотностью интегрального инварианта для линейной системы $\dot{x}=\Lambda x$. Можно сразу считать, что $\Lambda$ приведена к канонической жордановой форме. Расположим мономы формы $M_{s}$ в некотором лексикографическом порядке:
\[
M_{s}=\sum_{\substack{m_{i} \geqslant 0 \\ m_{1}+\ldots+m_{n}=s}} a_{m_{1} \ldots m_{n}} x_{1}^{m_{1}} \ldots x_{n}^{m_{n}} .
\]

Очевидно, что $\operatorname{div} M_{s}(\Lambda x)$ является некоторой формой той же степени. Приравнивая нулю ее коэффициенты, получим линейную однородную систему уравнений относительно $a_{m_{1} \ldots m_{n}}$. Определитель этой системы равен произведению
\[
\prod_{\substack{m_{i} \geqslant 0 \\ m_{1}+\ldots+m_{n}=s}}\left[\left(m_{1}+1\right) \lambda_{1}+\ldots+\left(m_{n}+1\right) \lambda_{n}\right] .
\]

Согласно предположению, это произведение отлично от нуля. Следовательно, все $a_{m_{1} \ldots m_{n}}=0$. Что и требовалось.

ЗАМЕЧАНИЕ. При более сильном условии отсутствия резонансных соотношений в традиционном смысле уравнения (8.1) не имеют первых интегралов, аналитических в окрестности точки $x=0$.

В качестве примера рассмотрим задачу о постоянных вращениях выпуклого твердого тела с аналитической выпуклой границей на горизонтальной абсолютно шероховатой плоскости (см. [4]). Движение такого тела описывается системой шести дифференциальных уравнений, имеющих интеграл энергии и геометрический интеграл. В частном случае, когда одна из главных центральных осей инерции тела ортогональна его поверхности, мы имеем однопараметрическое семейство стационарных вращений вокруг вертикально расположенной оси инерции. Стационарным движением отвечєют особые точки уравнений движения. Характеристическое уравнение имеет следующий вид:
\[
\lambda^{2}\left(a_{4} \lambda^{4}+a_{3} \lambda^{3}+a_{2} \lambda^{2}+a_{1} \lambda^{1}+a_{0}\right)=0 .
\]

Коэффициенты $a_{s}$ довольно сложно выражаются через многочисленные параметры задачи; фактически, они могут принимать произвольные значения. Наличие двойного нулевого корня связано с существованием двух независимых интегралов: в точках, отвечающих перманентным вращениям, дифференциалы интеграла энергии и геометрического интеграла в общем случае независимы. Фиксируя уровни первых интегралов, мы будем иметь дифференциальные уравнения на четырехмерных многообразиях, не имеющие в общем случае инвариантной меры с аналитической плотностью. Следовательно, в окрестности стационарных движений исходные уравнения также не имеют инвариантной меры.

Рассмотрим теперь задачу об инвариантной мере для систем дифференциальных уравнений, близких к интегрируемым системам, удовлетворяющим условиям теоремы 1. В окрестности инвариантных торов невозмущенной интегрируемой системы в качестве независимых переменных естественно принять постоянные первых интегралов $I_{1}, \ldots, I_{n-2}$ и угловые координаты $x, y \bmod 2 \pi$ на самих инвариантных торах. В этих переменных возмущенная система будет иметь следующий вид:
\[
\begin{array}{c}
\dot{I}_{s}=\varepsilon f_{s}(I, x, y)+\ldots, \quad s=1, \ldots, n-2, \\
\dot{x}=\frac{\lambda(I)}{\Phi(I, x, y)}+\varepsilon X(I, x, y)+\ldots, \quad \dot{y}=\frac{\mu(I)}{\Phi(I, x, y)}+\varepsilon Y(I, x, y)+\ldots
\end{array}
\]

Пусть функции, входящие в правые части этих дифференциальных уравнений, являются аналитическими в прямом произведении $D \times T^{2}$, где $D-$ некоторая область в $\mathbb{R}^{n-2}=\left\{I_{1}, \ldots, I_{n-2}\right\}, T^{2}=\{x, y$ $\bmod 2 \pi\} ; \varepsilon-$ малый параметр. Естественно поставить задачу о существовании у системы (8.2) инвариантной меры, плотность которой аналитична по переменным $I, x, y, 2 \pi$-периодична по $x, y$ и аналитически зависит от параметра $\varepsilon$ :
\[
M=M_{0}+\varepsilon M_{1}+\varepsilon^{2} M_{2}+\ldots
\]

Невозмущенная задача имеет инвариантную меру с плотностью $M_{0}$. В соответствии с известным принципом усреднения, усредним правые части системы (8.2) по мере $d m=\Phi d x \wedge d y$. В результате получим замкнутую систему уравнений для изменения медленных переменных $I$ в области
\[
\begin{array}{c}
\dot{I}_{s}=\varepsilon F_{s}(I), \quad 1 \leqslant s \leqslant n-2, \\
F_{s}=\frac{1}{\Lambda} \int_{0}^{2 \pi} \int_{0}^{2 \pi} f_{s} \Phi d x d y, \quad \Lambda=\int_{0}^{2 \pi} \int_{0}^{2 \pi} \Phi d x d y .
\end{array}
\]

Предложение 11. Предположим, что $m \lambda(I)+n \mu(i)
ot \equiv 0$ в области $D$ при всех целых $m$, $n$ не равных одновременно нулю. Если усредненная система (8.4) не имеет инвариантной меры с аналитической плотностью, то исходная система (8.2) также не имеет инвариантной меры с плотностью (8.4).

Система (8.4) проще системы (8.2); достаточное условие отсутствия инвариантной меры у (8.4) дает предложение 10.

ДОКАЗАТЕЛЬСТВО ПРЕДЛОЖЕНИЯ 11. Коэффициенты $M_{0}$ и $M_{1}$ разложения (8.4) удовлетворяют уравнениям
\[
\begin{array}{c}
\lambda \frac{\partial}{\partial x} \frac{M_{0}}{\Phi}+\mu \frac{\partial}{\partial y} \frac{M_{0}}{\Phi}=0 \\
\sum_{s} \frac{\partial}{\partial I_{s}}\left(M_{0} f_{s}\right)+\frac{\partial}{\partial x} M_{0} X+\frac{\partial}{\partial y} M_{0} Y=-\left(\lambda \frac{\partial}{\partial x} \frac{M_{1}}{\Phi}+\mu \frac{\partial}{\partial y} \frac{M_{1}}{\Phi}\right) .
\end{array}
\]

Поскольку $\lambda / \mu$ иррационально почти при всех $I \in D$, то из уравнения (8.5) вытекает, что $M_{0}=\Gamma(I) \Phi$. Подставляя это соотношение в уравнение (8.6) и усредняя затем по переменным $x, y$, получим следующее уравнение
\[
\sum_{s} \frac{\partial}{\partial I_{s}} \Gamma F_{s}=0 .
\]

Следовательно, функция $Г-$ плотность интегрального инварианта для системы (8.4). Осталось показать, что $\Gamma
ot \equiv 0$. Если это не так, то $M_{0}=0$. Но тогда функция $M_{1}+\varepsilon M_{2}+\ldots$ является плотностью инвариантной меры для системы (8.2). Если $M_{1} \equiv 0$, то эту операцию можно повторить еще раз. Предложение доказано.

ЗАМЕЧАНИЕ. Можно показать, что (в условиях предложения 11), если усредненная система (8.4) не имеет аналитического первого интеграла в области $D$, то исходная система (8.2) не имеет интеграла в виде ряда $g_{0}+\varepsilon g_{1}+\ldots$ с аналитическими в $D \times T^{2}$ коэффициентами $g_{s}$.

Рассмотрим более подробно частный случай, когда $n=3$. Индекс $s$ можно опустить. Пусть $F(I)
ot \equiv 0$. Если в некоторой точке интервала $D$ функция $F(I)=0$, то уравнение (8.4) не имеет, очевидно, инвариантной меры. Так мы будем считать, что $F(I)
eq 0$ в $D$. Рассмотрим следующие разложения Фурье:
\[
\begin{aligned}
\frac{X \Phi}{F} & =\sum X_{m n}(I) \exp i(m x+n y), \\
\frac{Y \Phi}{F} & =\sum Y_{m n}(I) \exp i(m x+n y), \\
\frac{f \Phi}{F} & =\sum f_{m n}(I) \exp i(m x+n y) .
\end{aligned}
\]

Резонансным множеством $\Delta$ назовєм множество точек $I \in D$, для которых
\[
\sum_{|m|+|n|
eq 0}\left|\frac{a_{m n}}{m \lambda+n \mu}\right|^{2}=\infty, \quad a_{m n}=\frac{d f_{m n}}{d I}+i\left(m X_{m n}+n Y_{m n}\right) .
\]

Предложение 12.
1) Предположим, что $\lambda(I) / \mu(I)
ot \equiv$ const;
2) пересечение $\Delta \cap D$ не пусто.

Тогда уравнения (8.2) не имеют интегрального инварианта с плотностью (8.3).

Действительно, из соотношения (8.7) получим, что $\Gamma=c / F$, где $c=$ $=$ const. Пусть
\[
\frac{M_{1}}{\Phi}=\sum b_{m n}(I) \exp i(m x+n y) .
\]

Уравнение (8.6) даст нам серию соотношений
\[
-(m \lambda+n \mu) b_{m n}=c a_{m n} .
\]

Пусть $I \in \Delta$. Тогда из условия
\[
\sum\left|b_{m n}\right|^{2}<\infty
\]

получим, что $c=0$. Что и требовалось.
Автор благодарен В.Ф.Журавлеву, прочитавшему рукопись и сделавшему ряд замечаний.