Хори [55] и Кэмел [57] независимо друг от друга развили методы теории возмущений для общих негамильтоновых систем, обобщив соответствующий метод, пригодный только для гамльтоновых систем дифференциальных уравнений. Ясно, что такое обобщение не является таким уж необходимым, ибо, как уже говорилось, любую систему можно свести к гамильтоновой увеличением ее порядка вдвое и взедением котангепциального пространства Дирака. Увеличение вдвое числа рассматриваемых уравнений окупается тем, что необходимо найти только две функции: новый гамильтониан и генератор преобразования. При
непосредственном подходе приходится иметь дело со столькими неизвестными функциями, сколько есть переменных; действительно, при использовании результатов § 7 главы I это становится ясно сразу же. Здесь мы опишем такой способ наиболее близким к работе Кэмела [57] об̈разом.

Рассмотрим систему $n$ дифференциальных уравнений первого порядка
\[
\dot{x}=f(x, \varepsilon)
\]

и предположим, что $\boldsymbol{f}(\boldsymbol{x} \varepsilon)$ – вещественная аналитическая по каждой из $n+1$ переменных $x_{1}, \ldots, x_{n}$, $\varepsilon$ вектор-функция в некоторой области $\Omega\left\{x \in D \subset R^{n},|\varepsilon|<\varepsilon_{0}\right\}$. Правая часть уравнения (2.9.1) может быть разложена при достаточно малых $\varepsilon$ в сходящийся степенной ряд
\[
\dot{\boldsymbol{x}}=\sum_{k \geqslant 0} \frac{\varepsilon^{k}}{k !} f^{(k)}(x),
\]

где
\[
f^{(k)}(x)=\left.\frac{\partial^{k} f}{\partial \varepsilon^{k}}\right|_{\varepsilon=0} .
\]

Очевидно, функции $\boldsymbol{f}^{(k)}(\boldsymbol{x})$ в области $D$ будут вещественными аналитическими функциями. В конечном счете это требование аналитичности может быть ослаблено, если на каждом шаге вычисления приближений восцользоваться олерацией сглаживания, однако для общего понимания метода это не имеет значения. Мы не будем рассматривать неавтономных систем, хотя и есть такие случаи, когда $t$ не может рассматриваться в качестве дополнительной $x$-переменной. Таковы, например, случаи, когда исследуется асимптотическое поведение, устойчивость или периодические решения.

Если уравнения (2.9.1) или (2.9.2) нельзя проинтегрировать в общем виде, то можно искать такое преобразование к новым $n$ переменным $\xi$, скажем,
\[
x=\boldsymbol{x}(\xi, \varepsilon),
\]

что дифференциальные уравнения относительно $\xi$
\[
\dot{\xi}=\boldsymbol{g}(\xi, \varepsilon),
\]

получающиеся из (2.9.3) и (2.9.1), поддаются более простому цзучению. Ясно, что сформулированная так проблема слишком обща, чтобы определить, какими свойствами должны обладать эти преобразования переменных. Один из способов убедиться в этом, разумеется, состоит в предположении о том, что при $\varepsilon=0$ общее решение уравнения (2.9.1) известно, т. е. уравнение
\[
\dot{\boldsymbol{y}}=\boldsymbol{f}(\boldsymbol{y}, 0)=\boldsymbol{f}^{(0)}(\boldsymbol{y})
\]

интетрируемо. Затем можно заинтересоваться вопросом о существовании преобразования (2.9.3), которое переводит уравнение (2.9.1) в уравнение $(2.9 .5)$, точнее, в уравнение
\[
\dot{\xi}=f^{(0)}(\xi) .
\]

Так как при $\varepsilon=0$ преобразование (2.9.3), очевидно, будет тождественным, то мы приходим к задаче нахождения преобразования, близкого к тождественному, т. е.
\[
\boldsymbol{x}=\boldsymbol{\xi}+\varepsilon h(\boldsymbol{\xi}, \varepsilon),
\]

а функция $h(\xi, \varepsilon)$ предполагается аналитической в некоторой области по каждой из $n+1$ переменных $\xi$, $\varepsilon$, причем эта область содержит точку $\varepsilon=0$. Очевидно, что при достаточно малых $\varepsilon$ вблизи $\varepsilon=0$ преобразование будет обратимым. Можно записать
\[
\boldsymbol{x}=\boldsymbol{\xi}+\sum_{k \geqslant 1} \frac{\varepsilon^{k}}{k !} \boldsymbol{E}_{k}(\xi),
\]

и цреобразованная система дифференциальных уравнений в общем случае будет иметь вид
\[
\dot{\xi}=\varphi(\xi, \varepsilon)=\sum_{k \geqslant 0} \frac{\varepsilon^{k}}{k !} \varphi^{(k)}(\xi)
\]

где
\[
\varphi^{(k)}(\xi)=\left.\frac{\partial^{k} \varphi}{\partial \varepsilon^{k}}\right|_{\varepsilon=0} .
\]

Теперь нашей задачей является получение по данному преобразованию (2.9.8) функций $\varphi^{(k)}(\xi)$ в $(2.9 .9)$ из функций $f^{(k)}(\boldsymbol{x})$ в (2.9.2). Очевидно, это можно сделать разными способами, однако, если необходимы члены больпого порядка и их систематическое исследование, то можно рекомендовать рекуррентный алгоритм, аналогичный алгоритму § 7 главы I. Дифференцирование соотношения (2.9.8) по времени $t$ дает
\[
\dot{x}=\dot{\xi}+\sum_{k \geqslant 1} \frac{\varepsilon^{k}}{k !} \frac{\partial E_{k}}{\partial \xi} \dot{\xi},
\]

и, подставляя сюда (2.9.2) и (2.9.9), находим
\[
\sum_{k \geqslant 0} \frac{\varepsilon^{k}}{k !} f^{(k)}(x)=\sum_{k \geqslant 0} \frac{\varepsilon^{k}}{k !} \varphi^{(k)}(\xi)+\sum_{k \geqslant 1} \frac{\varepsilon^{k}}{k !} \frac{\partial E_{k}}{\partial \xi} \sum_{j \geqslant 0} \frac{\varepsilon^{j}}{j !} \varphi^{(j)}(\xi) .
\]

Из соотношения (1.7.2) теперь мы видим, что
\[
f(\boldsymbol{x}(\xi, \varepsilon), \varepsilon)=\sum_{n=0}^{\infty} \frac{\varepsilon^{n}}{n !} f_{n}(\xi),
\]

и в нашем распоряжении есть рекуррентные соотношения для определения функций $f_{n}(\xi)$, например, уравнения (1.7.4) пли (1.7.15), или заключительные соотношения из § 7 главы I. Из (2.9.10) теперь следует, что
\[
f_{n}(\xi)=\boldsymbol{\varphi}^{(n)}(\xi)+\sum_{m=1}^{n} C_{n}^{m} \frac{\partial \boldsymbol{E}_{m}(\xi)}{\partial \xi} \boldsymbol{\varphi}^{(n-m)}(\xi) .
\]

Если теперь рассмотреть (1.7.22), то получим
\[
E_{n}(\xi)=-T_{n}(\xi)-\sum_{m=1}^{n-1} C_{n-1}^{m-1} T_{m}(\xi) \frac{\partial E_{n-m}(\xi)}{\partial \xi},
\]

или в обозначениях (1.7.9)
\[
\boldsymbol{E}_{n}(\xi)=-\boldsymbol{T}_{n}(\xi)-\sum_{m=1}^{n-1} C_{n-1}^{m-1} L_{m} \boldsymbol{E}_{n-m}(\xi) .
\]

Преобразование, обратное к (2.9.8), запишем в виде
\[
\xi=\boldsymbol{x}+\sum_{k \geqslant 1} \frac{\varepsilon^{k}}{k !} \boldsymbol{X}^{(k)}(x)
\]

так что
\[
\boldsymbol{X}^{(n)}(x)=\boldsymbol{T}_{n}(x)-\sum_{m=1}^{n-1} C_{n-1}^{m} L_{m} \boldsymbol{X}_{m, n-m}(x),
\]

где использованы введенные в (1.7.20) и (1.7.21) обозначения, т. e.
\[
\begin{array}{l}
\boldsymbol{X}_{p, q}(\boldsymbol{x})=-\sum_{m=1}^{p} C_{p-1}^{m-1} L_{m} \boldsymbol{X}_{p-m, q}(\boldsymbol{x}), \\
\boldsymbol{X}_{0, q}(\boldsymbol{x})=\boldsymbol{X}^{(q)}(\boldsymbol{x}) .
\end{array}
\]

Наконец, находим
\[
\varphi^{(n)}(\xi)=f^{(n)}(\xi)+\sum_{j=1}^{n} C_{n}^{j}\left[f_{j, n-j}(\xi)-\frac{\partial E_{j}}{\partial \xi} \varphi^{(n-j)}(\xi)\right],
\]

что и является искомым рекуррентным соотношением. Очевидно, уравнение (2.9.16) содержит коэффициенты $T_{n}$, определяющие отображение $(2.9 .8)$, т. е. коэффициенты $X^{(n+1)}(x)$ разложения
\[
\frac{\partial \xi}{\partial \varepsilon}=\sum_{n \geqslant 0} \frac{\varepsilon^{n}}{n !} \boldsymbol{X}^{(n+1)}(x),
\]

получаемого из (2.9.13), и коэффициенты разложения
\[
\frac{\partial \xi}{\partial \varepsilon}=\boldsymbol{T}(\xi, \varepsilon)=\sum_{n \geqslant 0} \frac{\varepsilon^{n}}{n !} \boldsymbol{T}_{n+1}(\xi)
\]

так же, как и в (1.7.6), (1.7.7) и (1.7.8). На каждом шаге приближения функции $T_{n}(\xi)$ в зависимости от выдвинутых требований должны быть выбраны соответствующим образом. Эти неизвестные функции $T_{n}(\xi)$ можно выделить в уравнении (2.9.16) стедующим образом:
\[
\frac{\partial T_{n}}{\partial \xi} \varphi^{(0)}(\xi)-\frac{\partial \varphi^{(0)}(\xi)}{\partial \xi} T_{n}(\xi)=\varphi^{(n)}-f^{(n)}-G_{n}(\xi),
\]

где функции $\boldsymbol{G}_{n}(\xi)$ зависят от всех предыдущих приближений. Кәмел получил для них такие формулы:
\[
\begin{array}{l}
\boldsymbol{G}_{n}(\xi)= \\
\quad=\frac{\partial E_{n}^{*}(\xi)}{\partial \xi} \boldsymbol{\varphi}^{(0)}(\xi)-f_{n, 0}^{*}(\xi)+\sum_{m=1}^{n-1} C_{n}^{m}\left[\frac{\partial E_{m}}{\partial \xi} \varphi^{(n-m)}(\xi)-f_{m, n-m}^{r}(\xi)\right],
\end{array}
\]

где
\[
\begin{aligned}
\boldsymbol{E}_{n}^{*}(\boldsymbol{\xi}) & =\boldsymbol{E}_{n}(\boldsymbol{\xi}) \quad \text { при } \quad \boldsymbol{T}_{n}=0, \\
\boldsymbol{f}_{n, 0}^{*}(\boldsymbol{\xi}) & =\boldsymbol{f}_{n, 0}(\boldsymbol{\xi}) \quad \text { при } \quad \boldsymbol{T}_{n}=0, \\
f_{p, q}(\boldsymbol{\xi}) & =-\sum_{m=1}^{p} C_{p-1}^{m-1} L_{m} \boldsymbol{f}_{p-m, q}(\boldsymbol{\xi}), \\
\boldsymbol{f}_{0, q}(\boldsymbol{\xi}) & =\boldsymbol{\varphi}^{(q)}(\boldsymbol{\xi}) .
\end{aligned}
\]

Полное описание этого метода было дано в работах Кәмела [57] и Хенрарда [51] и в работе Хори [55]. Кәмел показал, что подход, основанный на преобразованиях Ли, содержит в себе как частные случаи важные методы разложения по двум переменным и подбора асимнтотических решений, развитые в работе Кеворкяна [58]. Этот вопрос здесь не рассматривается, так как он весьма подробно изложен в работе Коула (см. [18.5]). Может быть, стоило бы отметить, что алгоритм Депри для преобразований Ли, которые генерируются функцией малого параметра, примененный к функциям Гамильтона, зависящим от этого же параметра, также может быть здесь получен, если с самого начала забыть про наличие этого параметра, что и было показано выше и следует из работы Мерсмана [73]. Аналогичным образом описанный выше метод может быть упрощен введением операторов и функций, не являющихся функциями параметра, с последующим разложением в степенные ряды по $\varepsilon$ всех полученных результатов. Такой подход был осуществлен Хори [55], и здесь мы ограничимся кратким изложением его результатов.

Рассмотрим набор из $n$ переменных $\xi_{1}, \ldots, \xi_{n}$ и операторы $T_{k}(\xi)$, и пусть
\[
D_{\xi}=\sum_{k=1}^{n} T_{k}(\xi) \frac{\partial}{\partial \xi_{k}} .
\]

Рассмотрим отображение
\[
x_{j}=\xi_{j}+\sum_{p \geqslant 1} \frac{1}{p !} D^{p-1} T_{j}(\xi),
\]

где
\[
D_{\xi}^{0}=1, \quad D_{\xi}^{1}=D_{\xi}, \quad D_{\xi}^{p}=D_{\xi} D_{\xi}^{p-1} .
\]

Это отображение аналогично отображению (2.8.1) и соответствующим определениям. В частности, функция $T_{k}(\xi)$ играет роль функции $\partial S / \partial \xi_{k}$. Рассмотрим также отображение вещественной аналитической функции $f(\boldsymbol{x})$, зависящей от $n$ переменных $x_{1}, \ldots$ $\ldots, x_{n}$, в $\xi$-пространство по формулам
\[
f(x)=f(\xi)+\sum_{p \geqslant 1} \frac{1}{p !} D_{\xi}^{v} f(\xi) .
\]

В действительности формулы (2.9.20) являются следствием формул (2.9.21). Определим обратное отображение
\[
\begin{aligned}
T_{k}^{-1}(x) & =\left.T_{k}(\xi)\right|_{\xi=x}, \\
D_{x} & =\sum_{k=1}^{n} T_{k}^{-1}(x) \frac{\partial}{\partial x_{k}},
\end{aligned}
\]

так что, обращая отображение (2.9.21), имеем
\[
f(\xi)=f(x)+\sum_{p \geqslant 1} \frac{(-1)^{p}}{p !} D_{x}^{p} f(x),
\]

что является непосредственным обобщением преобразования Ли. Все выписанные соотношения в действительности в том или ином виде содержатся в предыдущем изложении (где есть зависимость от $\varepsilon$ ), а их доказательство получается сразу же.
Уравнение
\[
\dot{x}_{k}=\dot{f}_{k}(x)
\]

после применения преобразования (2.9.20), генерируемого функцией $T_{k}$, переходит в уравнение
\[
\dot{\xi}_{k}=\varphi_{k}(\xi) .
\]

Используя (2.9.24), получаем преобразование, обратное к (2.9.20), в виде
\[
\xi_{j}=x_{j}+\sum_{p \geqslant 1} \frac{(-1)^{p}}{p !} D_{x}^{p-1} T_{j}^{-1}(x) .
\]

Так как из (2.9.25) следует, что
\[
\frac{d}{d t}=\sum_{k=1}^{n} f_{k} \frac{\partial}{\partial x_{k}}
\]

то для произвольной функцип имеем
\[
\frac{d}{d t} F(x)=\dot{F}(x)=\sum_{k=1}^{n} \frac{\partial F}{\partial x_{k}} \dot{x}_{k}=\left(\sum_{k=1}^{n} f_{h} \frac{\partial}{\partial x_{k}}\right) F .
\]

Вычисление функции $\varphi_{k}(\xi)$ в (2.9.26) проведем следующим образом. Дифференцирование соотношения (2.9.27) дает
\[
\dot{\xi}_{j}=\dot{x}_{j}+\sum_{p \geqslant 1} \frac{(-1)^{p}}{p !} \frac{d}{d t}\left\{D_{x}^{p-1} T_{j}^{-1}(x)\right\},
\]

и, подставляя сюда (2.9.25) и (2.9.23), ваходим
\[
\dot{\xi}_{j}=f_{j}(x)+\sum_{p \geqslant 1} \frac{(-1)^{p}}{p !} \sum_{k=1}^{n} f_{k}(x) \frac{\partial}{\partial x_{k}}\left(D_{\boldsymbol{x}}^{p-1} T_{j}^{-1}(x)\right)
\]

или, используя (2.9.21) и (2.9.20),
\[
\begin{array}{l}
\dot{\xi}_{j}=f_{j}(\xi)+\sum_{p \geqslant 1} \frac{1}{p !} D_{\xi}^{p} f_{j}(\xi)+\sum_{p \geqslant 1} \frac{(-1)^{p}}{p !} \sum_{k=1}^{n} f_{k}(\xi) \frac{\partial}{\partial \xi_{k}}\left(D_{\xi}^{p-1} T_{j}(\xi)\right)+ \\
+\sum_{p \geqslant 1} \frac{(-1)^{p}}{p !} \sum_{q \geqslant 1} \frac{1}{q !} D_{\xi}^{q}\left\{\sum_{k=1}^{n} f_{k}(\xi) \frac{\hat{o}}{\partial \xi_{k}}\left(D^{p-1} T_{j}(\xi)\right)\right\}=\varphi_{j}(\xi) . \quad \text { (2.9.29) }
\end{array}
\]

Теперь рассмотрим ряды
\[
\begin{array}{l}
f_{j}=f_{j}^{(0)}+f_{j}^{(1)}+\ldots, \\
\varphi_{j}=\varphi_{j}^{(0)}+\varphi_{j}^{(1)}+\ldots \\
T_{j}=T_{j}^{(0)}+T_{j}^{(1)}+\ldots
\end{array}
\]

и будем искать такие операторы $T_{j}$, чтобы функции $\varphi_{j}$ имели желаемый вид. Очевидно, предшолагаем, что уравнения
\[
\dot{y}_{k}=f_{k}^{(0)}(y)
\]

имеют известное общее решение. Разложение (2.9.30) функций $f_{i}$ не обязательно подразумевает разложение в степенные ряды по некоторому малому параметру и не обязательно имеет бесконечное число членов. Действительно, в обычном случае для возмущений интегрируемой системы (2.9.31) имеем $f_{j}^{(k)}=0$ при $k \geqslant 2$, т. e.
\[
f_{j}=f_{j}^{(0)}+f_{j}^{(1)} .
\]

Подставляя ряды (2.9.30) в (2.9.29) и приравнивая члены одинакового порядка, можно получить рекуррентный алгоритм вычисления неизвестных функций $\varphi_{j}^{(k)}$ и $T_{j}^{(h)}$. В этом отношении явное использование параметра є для представления членов разного порядка является весьма удобным, хотя и не обязательным. Это означает, что приравнивание членов одинакового порядка лучше заменить приравниванпем коэффициентов при одинаковых степенях $\varepsilon$, если положить $f_{j}^{(k)}=O\left(\varepsilon^{k}\right), \varphi_{j}^{(k)}=O\left(\varepsilon^{k}\right), T_{j}^{(k)}=$ $=O\left(\varepsilon^{k}\right)$,
Первые несколько приближений имеют вид
\[
\begin{array}{c}
f_{j}^{(0)}=\varphi_{j}^{(0)}, \\
\sum_{k=1}^{n}\left\{-f_{k}^{(0)} \frac{\partial T_{j}^{(1)}}{\partial \xi_{k}}+T_{k}^{(1)} \frac{\partial f_{j}^{(0)}}{\partial \xi_{k}}\right\}+f_{j}^{(1)}=\varphi_{j}^{(1)}, \\
\sum_{k=1}^{n}\left\{-f_{k}^{(0)} \frac{\partial T_{j}^{(2)}}{\partial \xi_{k}}+T_{k}^{(2)} \frac{\partial f_{j}^{(0)}}{\partial \xi_{k}}\right\}+\frac{1}{2 !} \sum_{k=1}^{n} T_{k}^{(1)} \frac{\partial}{\partial \xi_{k}}\left(\varphi_{j}^{(1)}+f_{j}^{(1)}\right)- \\
-\frac{1}{2 !} \sum_{k=1}^{n} \frac{\partial T_{j}^{(1)}}{\partial \xi_{k}}\left(\varphi_{k}^{(1)}+f_{k}^{(1)}\right)+f_{j}^{(2)}=\varphi_{j}^{(2)},
\end{array}
\]

а в общем случае
\[
\sum_{k=1}^{n}\left\{-f_{k}^{(0)} \frac{\partial T_{i}^{(p)}}{\partial \xi_{k}}+T_{k}^{(p)} \frac{\partial f_{j}^{(0)}}{\partial \xi_{k}}\right\}+\cdots+f_{j}^{(p)}=\varphi_{j}^{(p)},
\]

что, по суцеству, эквивалентно выписанным выше соотношениям (2.9.17). Здесь все функции зависят от переменных $\xi$. Теперь мы введем важное понятие дополнительной системы, определив
\[
\frac{d \xi_{j}}{d \tau}=f_{j}^{(0)}(\xi)
\]

с общим решением
\[
\xi_{j}=\xi_{j}(\tau)
\]

так что
\[
\sum_{k=1}^{n}\left\{-f_{k}^{(0)} \frac{\partial T_{j}^{(p)}}{\partial \xi_{k}}+T_{k}^{(p)} \frac{\partial f_{j}^{(0)}}{\partial \xi_{k}}\right\}=-\frac{d T_{j}^{(p)}}{d \tau}+\sum_{k=1}^{n} T_{k}^{(p)} \frac{\partial f_{j}^{(0)}}{\partial \xi_{k}} .
\]

Общее уравнение (2.9.32) на каждом шаге вычисления приближения приводится к линейной системе дифференциальных уравнений относительно $T_{j}^{(p)}(\xi)$
\[
-\frac{d T_{j}^{(p)}}{d \tau}+\sum_{k=1}^{n} \frac{\partial f_{j}^{(0)}(\xi(\tau))}{\partial \xi_{k}} T_{k}^{(p)}(\tau)+F_{j}^{(p)}(\tau)=\varphi_{j}^{(p)},
\]

где все $\xi$ заменены решениями (2.9.34) дополнительной системы. Ясно, что уравнение (2.9.35) является непосредственным обобщением уравнения (2.8.9). Заметим, что как и в обычных методах усреднения, функции $\varphi_{i}^{(p)}$ надо выбирать так, чтобы в $T_{j}^{(p)}(\tau)$ отсутствовали секулярные члены, т. е. надо, чтобы предел
\[
\lim _{\tau \rightarrow \infty} T_{j}^{(p)}(\tau)
\]

был конечным. В простейшем случае функция $f_{j}^{(0)}$ линейно зависит от переменных $\xi$, так что система уравнений (2.9.35) явлнется линейной неоднородной системой с постоянными коәффициентами для приближения любого порядка. Если имеет место не такой случай, а, например, функции $\partial f_{j}^{(0)} /\left.\partial \xi_{k}\right|_{\xi=\xi(\tau)}$ являются периодическими или условно-периодическими функциями $\tau$, то интегрирование уравнений (2.9.36) представляет собой нетривиальную задачу. Следовательно, желательно произвести такое разложение функций $f_{j}(\xi)$, чтобы все $f_{j}^{(0)}(\xi)$ были линейными.

Уравнение Ван дер Поля. В качестве примера рассмотрим уравнение
\[
\ddot{x}+\varepsilon\left(1-x^{2}\right) \dot{x}+x=0,
\]

которое можно записать так
\[
\dot{x}_{1}=x_{2}, \quad \dot{x}_{2}=-x_{1}-\varepsilon\left(1-x_{1}^{2}\right) x_{2} .
\]

Здесь мы имеем
\[
\begin{array}{ll}
f_{1}^{(0)}=x_{2}, \quad f_{2}^{(0)}=-x_{1}, & f_{1}^{(1)}=f_{1}^{(2)}=\ldots=f_{1}^{(p)}=\ldots=0, \\
f_{2}^{(1)}=-\varepsilon\left(1-x_{1}^{2}\right) x_{2}, & f_{2}^{(2)}=f_{2}^{(3)}=\ldots=f_{2}^{(p)}=\ldots=0 .
\end{array}
\]

Дополнительная система
\[
\frac{d \xi_{j}}{d \tau}=\cdot f_{j}^{(0)}=\varphi_{j}^{(0)}
\]

имеет решение вида
\[
\xi_{1}=\alpha \cos (\tau+\beta), \quad \xi_{2}=-\alpha \sin (\tau+\beta),
\]

где $\alpha, \beta$-скалярные постоянные. Уравнения первого приближения принимают вид
\[
\begin{array}{r}
-\frac{d T_{1}^{(1)}}{d \tau}+T_{2}^{(1)}=\varphi_{1}^{(1)}, \\
-\frac{d T_{2}^{(1)}}{d \tau}-T_{1}^{(1)}+\varepsilon\left[1-\alpha^{2} \cos ^{2}(\tau+\beta)\right] \alpha \sin (\tau+\beta)=\varphi_{2}^{(1)},
\end{array}
\]

или
\[
\begin{array}{l}
\frac{d_{2} T_{1}^{(1)}}{d \tau^{2}}+T_{1}^{(1)}=\varepsilon\left[\left(1-\frac{\alpha^{2}}{4}\right) \alpha \sin (\tau+\beta)-\frac{\alpha^{3}}{4} \sin (3 \tau+3 \beta)\right]- \\
-\frac{d \varphi_{1}^{(1)}}{d \tau}-\varphi_{2}^{(1)} .
\end{array}
\]

Для того чтобы не появилось секулярных членов, в уравнении для $T_{1}^{(1)}$ должен отсутствовать член $\sin (\tau+\beta)$. Один из возможных способов. выбора произвольных функций описывается формулами
\[
\begin{array}{l}
\frac{d^{2} T_{1}^{(1)}}{d \tau^{2}}+T_{1}^{(1)}=-\varepsilon \frac{\alpha^{3}}{4} \sin (3 \tau+3 \beta), \\
\frac{d \varphi_{1}^{(1)}}{d \tau}+\varphi_{2}^{(1)}=\varepsilon\left(1-\frac{\alpha^{2}}{4}\right) \alpha \sin (\tau+\beta), \\
\varphi_{2}^{(1)}=\frac{d \varphi_{1}^{(1)}}{d \tau}=\frac{1}{2} \varepsilon\left(1-\frac{\alpha^{2}}{4}\right) \alpha \sin (\tau+\beta),
\end{array}
\]

так что
\[
\begin{array}{l}
\varphi_{2}^{(1)}=-\varepsilon\left[1-\frac{1}{4}\left(\xi_{1}^{2}+\xi_{2}^{2}\right)\right] \xi_{2}, \\
\varphi_{1}^{(1)}=-\varepsilon\left[1-\frac{1}{4}\left(\xi_{1}^{2}+\xi_{2}^{2}\right)\right] \xi_{1}
\end{array}
\]

и, следовательно,
\[
\begin{array}{l}
T_{1}^{(1)}=\frac{\varepsilon \alpha^{3}}{32} \sin (3 \tau+3 \beta)=\frac{\varepsilon}{32} \xi_{2}\left(\xi_{2}^{2}-3 \xi_{1}^{2}\right) \\
T_{2}^{(1)}=\frac{\varepsilon}{2} \xi_{1}\left(-1+\frac{7}{16} \xi_{1}^{2}-\frac{5}{16} \xi_{2}^{2}\right) .
\end{array}
\]

Таким образом, уравнения в новых переменных в первом приближении имеют вид
\[
\begin{array}{l}
\frac{d \xi_{1}}{d t}=\xi_{2}-\varepsilon\left[1-\frac{1}{4}\left(\xi_{1}^{2}+\xi_{2}^{2}\right)\right] \xi_{1} \\
\frac{d \xi_{2}}{d t}=-\xi_{1}-\varepsilon\left[1-\frac{1}{4}\left(\xi_{1}^{2}+\xi_{2}^{2}\right)\right] \xi_{2} .
\end{array}
\]

Легко проверить, что уравнение для $\xi_{2}$ получается из уравнения для $\xi_{1}$ при замене $\xi_{2} \rightarrow-\xi_{1}, \xi_{1} \rightarrow \xi_{2}$, что объясняется сделанным выбором $\varphi_{2}^{(1)}=d \varphi_{1}^{(1)} / d \tau$. Если положить
\[
u^{2}=\xi_{1}^{2}+\xi_{2}^{2},
\]

то найдем, что
\[
\frac{d u^{2}}{d t}=-2 \varepsilon\left(1-\frac{u^{2}}{4}\right) u^{2}
\]

и, следовательно,
\[
u^{2}=\xi_{1}^{2}+\xi_{2}^{2}=\frac{4 k e^{-\varepsilon t}}{k e^{-\varepsilon t} \pm 1},
\]

а знак + или – надо выбирать в зависимости от знака постоянной $k$, т. е. в зависимости от начальных условий так, чтобы величина $u^{2}$ была положительной.
При $\varepsilon>0$ мы получаем асимптотическое поведение:
\[
u^{2} \rightarrow 0 \quad \text { при } \quad t \rightarrow \infty,
\]

что описывает хорошо известное демпфированное движение по направлению к фокусу. Если $\varepsilon<0$, то $u^{2} \rightarrow 4$ при $t \rightarrow \infty$ и мы имеем предельный цикл в уравнении Ван дер Поля. То, что первое приближение (по $\varepsilon$ ) дает возможность получить полную информацию об асимптотическом поведении системы, объясняется тем, что в любом приближении уравнения для $\xi_{1}, \xi_{2}$ имеют тот же вид, что и уравнения первого приближения, т. е.
\[
\begin{array}{c}
\dot{\xi}_{1}=\left[1+\varepsilon^{2} f_{2}\left(u^{2}\right)+\varepsilon^{4} f_{4}\left(u^{2}\right)+\ldots\right] \xi_{2}- \\
-\varepsilon\left[1-\frac{u^{2}}{4}+\varepsilon^{2} g_{3}\left(u^{2}\right)+\varepsilon^{4} g_{5}\left(u^{2}\right)+\ldots\right] \xi_{1}, \\
\dot{\xi}_{2}=\dot{\xi}_{1} \quad\left(\xi_{2} \rightarrow-\xi_{1}, \quad \xi_{1} \rightarrow \xi_{2}\right),
\end{array}
\]

так что ошисанные выше асимптотические свойства сохраняются.