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

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

При этом существенное значение будут иметь рассуждения пп. 27 , 28 об инвариантности условий стационарности; с механической точки

зрения наиболее важный случай будет тот, когда полная энергия $H$ имеет стационарное значение в абсолютном смысле или в зависимости от инвариантных соотношений или каких-либо интегралов. Решения, к которым мы таким образом приходим, мы будем называть, следуя Раусу, стационарными.
52. Статические решения. Чтобы начать с простого, но не лишенного, однако, интереса случая, возьмем снова каноническую систему, характеристическая функция которой не зависит от $t$. В этом случае существует интеграл $H=$ const, и, согласно следствию п. 27 , соответствующее условие стационарности $\delta H=0$ позволяет написать $2 n$ инвариантных соотношений
\[
\frac{\partial H}{\partial p_{h}}=0, \frac{\partial H}{\partial q_{h}}=0 \quad(h=1,2, \ldots, n),
\]

которые приводят каноническую систему к виду $\dot{p}_{h}=\dot{q}_{h}=0$ и показывают, что решения, при которых удовлетворяется условие $\delta H=0$, все являются как раз такими, при которых отдельные $p$ и $q$ сохраняют постоянные значения. Следовательно, мы имеем дело со статическим решением в узком смысле п. 17 гл. VI; так как число $2 n$ уравнений (104) как раз равно числу постоянных $\tilde{p}^{0}, q^{0}$, то эти уравнения, за исключением случаев несовместности и неопределенности, пригодны для определения искомых решений.

На любое из этих решений $\sigma$ распространяется замечание, вытекающее из теоремы Дирихле для динамического случая, а именно, что возможно указать чисто качественное условие устойчивости, т.е. условие, выражаемое посредством одних только соотношений неравенства. Действительно, таким является в силу уравнений (104) условие, что $H$ имеет для решения $\bar{\sigma}$ действительный максимум или минимум (см. п. 7 и гл. VII, пп. 5-6, 17); замечание о лагранжевых системах с кинетическим потенциалом, не зависящим от времени, в конце упомянутого п. 17, гл. VI, таким образом, будет вполне оправдано, так как, как это непосредственно следует из п. 1 той же самой главы, всякая такая лагранжева система определяет каноническую систему с характеристической функцией, не зависящей от $t$, и обратно.
53. СЛУЧАЙ НЕСКОЛЬКИХ иНВАРИАнТНЫХ СОотнОшЕНИЙ, НАХОДящихСя в инволюции. Переходя после этого к более общим предположениям, докажем, что если для указанной канонической системы порядка $2 n$ с характеристической функцией, не зависящей от времени, известны $m<n$ инвариантных соотношений, находящихся в инволюции и разрешимых относительно $m$ переменных $p$, то можно определить $\infty^{m}$ частных решений данной системы посредством интегрирования приведённой системы дифференциальных уравнений порядка $m$.

Действительно, пусть известные инвариантные соотношения даны в разрешенном виде
\[
p_{r}=\varphi_{r}\left(p_{m+1}, \ldots, p_{n} ; q_{1}, \ldots, q_{n}\right) \quad(r=1,2, \ldots, m) ;
\]

как мы знаем из п. 30 , если наложим на функцию $H$ условие стационарности, имея в виду эти $m$ соотношений (т. е. если мы исключим из $H$ величины $p_{1}, p_{2}, \ldots, p_{m}$ посредством равенств (105) и, обозначив через $\bar{H}$ приведенную таким образом функцию, положим $\delta \bar{H}=0$ ), то явные уравнения, эквивалентные этому условию стационарности, сведутся к $2(n-m)$ уравнениям
\[
\frac{\partial \bar{H}}{\partial p_{m+j}}=0, \quad \frac{\partial \bar{H}}{\partial q_{m+j}}=0 \quad(j=1,2, \ldots, n-m),
\]

так как из них непосредственно следует, что остальные производные $\partial \bar{H} / \partial q_{r}(r=1,2, \ldots, m)$ обращаются в нуль.

Если, обращаясь к условиям, образующим нормальный случай, предположим, что уравнения (106) разрешимы относительно $p_{m+j}$, $q_{m+j}(j=1,2, \ldots, n-m)$, то эти уравнения и равенства (105) позволяют полностью выразить все $p$ и $q_{m+1}, \ldots, q_{n}$ через $q_{1}, \ldots, q_{m}$.

С другой стороны, если первоначальную каноническую систему представим себе разбитой на две частичные системы
\[
\begin{aligned}
\dot{p}_{h} & =-\frac{\partial H}{\partial q_{h}} & & (h=1,2, \ldots, n), \\
\dot{q}_{m+j} & =\frac{\partial H}{\partial p_{m+j}} & & (j=1,2, \ldots, n-m), \\
\dot{q}_{r} & =\frac{\partial H}{\partial p_{r}} & & (r=1,2, \ldots, m),
\end{aligned}
\]

то уравнения ( $\left.5^{\prime}\right)$ будут обязательно удовлетворяться только что указанными $2(n-m)$ разрешенными уравнениями; а в силу тех же явных выражений $p_{h}$ и $q_{m+j}$ равенства (5\”) сводятся к системе $m$ уравнений относительно $q_{1}, \ldots, q_{m}$, т. е. к системе порядка $m$. Интегрирование этой системы введет $m$ произвольных постоянных; соответствующие уравнения интегралов, присоединенные к равенствам (105) и (106), как раз дадут искомый класс $\infty_{\varepsilon}^{m}$ решений для заданной канонической системы. Один частный случай, заслуживающий упоминания, в котором $m=n$, получен на основании замечания п. 40 , так как форма уравнений (71), позволяет видеть, что речь идет о системе, находящейся в инволюции.

Естественно, что общность решений, к которой мы приходим таким способом; возрастает, когда какое-нибудь из известных инвариантных соотношений (105) является истинным интегралом и поэтому содержит произвольную постоянную; поэтому если при приведении мы прямо используем $m$ интегралов, находящихся попарно в инволюции, то придем к классу $\infty^{2 m}$ частных решений,

54. Движение по Раусу. Результат предыдущего пункта находит важное применение, когда каноническая система имеет $m$ игнорируемых координат $q_{1}, q_{2}, \ldots, q_{m}$. Как мы уже знаем (п. 42), в этом случае существуют $m$ интегралов обобщенных количеств движения
\[
p_{r}-p_{r}^{0}=0 \quad(r=1,2, \ldots, m) ;
\]

соответствующая приведенная характеристическая функция $\vec{H}$ является не чем иным, как первоначальной функцией $H$, в которой вместо каждого $p_{r}$ подставлено соответствующее постоянное значение $p_{r}^{0}$. Если теперь сопоставим условия стационарности (106) с каноническими уравнениями
\[
\dot{p}_{m+j}=-\frac{\partial H}{\partial q_{m+j}}, \quad \dot{q}_{m+j}=\frac{\partial H}{\partial p_{m-j}} \quad(j=1,2, \ldots, n-m),
\]

то увидим, что искомые стационарные решения будут меростатическими (см. гл. VI, п. 24); точнее, каждое такое решение состоит из решения уравнений (5а), статического относительно $p_{m+j}, q_{m+j}$, к которому надо присоединить постоянные значения также и для $p_{\boldsymbol{r}}$, а $q_{r}$ на основании остальных уравнений
\[
\dot{q}_{r}=\frac{\partial H}{\partial p_{r}} \quad(r=1,2, \ldots, m)
\]

будут линейными функциями от $t$.
Этот тип семейств $\infty^{2 m}$ стационарных решений был изучен Раусом в частном предположении динамического случая; и поэтому движения, которые определяются этими решениями, называются движениями по Payсу.
55. Заслуживает внимания то обстоятельство, что с теоретической точки зрения рассмотренный в п. 53 случай оказывается только более общим случаем Рауса, разобранным в предыдущем пункте. Действительно, как это доказывается в теории преобразований прикосновения, инвариантные соотношения (105), находящиеся в инволюции, можно всегда привести надлежащим (вполне) каноническим преобразованием переменных $p, q$ к простейшему виду
\[
p_{r}=0 \quad(r=1,2, \ldots, m) ;
\]

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

Это замечание приводит к почти непосредственному доказательству результата п. 30. Достаточно заметить, что, когда выполнено преобразование, приводящее равенства (105) к виду (105′), инвариантность этих последних соотношений в силу канонических уравнений
\[
\dot{p}_{r}=-\frac{\partial H}{\partial q_{r}} \quad(r=1,2, \ldots, m)
\]

влечет за собой то, что $m$ производных $\partial H / \partial q_{r}$ должны тождественно обращаться в нуль в силу равенств ( $105^{\prime}$ ), откуда, полагая преобразование выполненным в обратном порядке, мы прямо увидим, что условия стационарности характеристической функции $H$, соответственно равенствам (105), выражаются только посредством $2(n-m)$ соотношений (106).

Это доказательство оказывается, несомненно, более простым, чем доказательство II. 30 ; нужно, однако, заметить, что оно опирается на теорию преобразований прикосновения, которую мы здесь не затрагивали. Во всяком случае, даже если мы отвлечемся от этого несущественного обстоятельства, теоретическая возможность сведения $m$ инвариантных соотношений, находящихся в инволюции (105), к соотношениям (105′) не лишает интереса рассуждения, которые мы развили в пп. 28,30 , относя систему к совершенно общим координатам.

Существенная цель нашего исследования состояла в определении некоторых классов решений простыми средствами, или, по крайней мере, более простыми, чем полное интегрирование заданной системы дифференциальных уравнений, каковым является интегрирование приведенной системы порядка $m<2 n$. Далее, определение преобразования прикосновения, пригодного для приведения случая п. 53 к случаю Рауса, вообще говоря, требует операций порядка более высокого, чем $m$, так что его нельзя рассматривать как полезное орудие для вычисления, хотя совершенно законно и даже удобно пользоваться им как средством для доказательства.
56. ЗАМЕЧАНИЯ О ДЕЙСТВИТЕЛЬНОМ ПОСТРОЕНИИ СТАЦИОНАРНЫХ РЕШЕНИЙ. Чтобы выразить условия стационарности функции $f_{0}$ соответственно некоторому числу $m$ соотношений
\[
f_{r}=0 \quad(r=1,2, \ldots, m),
\]

инвариантных относительно заданной системы дифференциальных уравнений, $n$-го порядка (96), нет нужды строго следовать способу, указанному в наших теоретических выкладках в п. 30 , т. е. прежде всего разрешать систему (107) относительно $n$ аргументов и после исключения этих аргументов из $f_{0}$ приравнивать нулю виртуальную вариацию приведенной таким образом функции $\overline{f_{0}}$. Наоборот, как известно из анализа, можно избежать предварительного решения уравнения (107), прибегая к классическому методу множителей.

Этот метод состоит, как известно, в присоединении к уравнениям (107) символического уравнения
\[
\delta f_{0}+\sum_{r=1}^{n} \mu_{r} \delta f_{r}=0,
\]

где $\delta f_{0}, \delta f_{r}$ обозначают виртуальные вариации, а множители $\mu$ должны рассматриваться как вспомогательные неизвестные. Исключая эти

множители $\mu$ из $n$ уравнении, объединенных в уравнении (108), и присоединяя уравнения, которые таким образом получатся, к уравнениям (107), мы придем к системе, эквивалентной системе, состоящей из уравнений (107) и из уравнения $\delta f_{0}=0$.

Другими словами, можно сказать, что уравнения (107), (108), по существу, содержат две группы уравнений $(A)$ и (B), из которых уравнения ( $A$ ) составляют систему, инвариантную относительно заданной системы дифференциальных уравнений с одними только $x$, а уравнения ( $B$ ) дают определение множителей $\mu$ в функция от $x$. Отсюда следует, что если продифференцируем по $t$ систему уравнений (107), (108) или эквивалентную ей систему ( $A$ ), ( $B$ ) и примем, конечно, во внимание систему (36), то частичная система ( $A$ ) в силу своего инвариантного характера не дает места никакому новому соотношению, тогда как система ( $B$ ) приведет к такому же числу уравнений ( $B^{\prime}$ ), которые определят производные от множителей $\mu$ в функциях от $x$. Таким образом, можно также сказать, что система (107), $(108)$, как и эквивалентная ей система $(A),(B)$, инвариантна относительно расширенной системы дифференциальных уравнений с $n+m$ неизвестными $\boldsymbol{x}$ и $\mu$, которая получается путем присоединения к уравнениям (36) уравнений ( $\left.B^{\prime}\right)$, так как, дифференцируя уравнения $(A)$ и (B) и принимая во внимание уравнения (36), $\left(B^{\prime}\right)$, мы не придем к какому-нибудь соотношению между $x, \mu, t$, отличному от уравнений (107), (108).

Это и есть результат, от которого следует исходить при фактическом приведении системы дифференциальных уравнений и идти к определению стационарных решений, о которых идет речь.

Необходимо признать, что, с теоретической точки зрения, способ множителей как способ, преобразующий первоначальную задачу, сводящуюся к системе дифференциальных уравнений с $n$ неизвестными, в аналогичный вопрос, связанный с системой уравнений с $n+m$ неизвестными, не представляет преимуществ по сравнению с первоначальным способом. Однако вместе с указанным выше преимуществом, заключающимся в том, что его применение позволяет избежать предварительного решения уравнени́й (107), он соединяет еще достоинство особенной алгоритмическон ясности, которая, как мы увидим, будет ценна в механических приложениях, так как допускает прямое и изящное истолкование природы движения.