Рассмотрим еще колебательную систему, описываемую уравнением вида
$$\frac{d^{2}x}{d{t}^2}+{\omega }^{2}x=\epsilon f(x)\frac{dx}{dt},$$

которое также является частным случаем уравнения (1.1).
Заметим, что ранее рассмотренное уравнение (3.1) может быть приведено к виду (4.1).
В самом деле, полагая
$$\frac{dx}{dt}=y$$

и дифференцируя уравнение (3.1), получаем:
$$m\frac{d^{2}y}{d{t}^2}+ky=\epsilon F^{\prime }(y)\frac{dy}{dt}.$$

Сопоставляя уравнение (4.1) с (1.1), имеем:
$$f(x,\frac{dx}{dt})=f(x)\frac{dx}{dt}$$

поэтому, для того чтобы воспользоваться формулами (1.21)-(1.28), необходимо разложить в ряд Фу́рье выражение
$$f(a\cos \psi )a\omega \sin \psi .$$

Для упрощения этой операции рассмотрим функцию
$$F^{\ast }(x)={\int }_0^{x}f(x)dx$$

и разложение в ряд Фурье:
$$F^{\ast }(a\cos \psi )=\sum _{n=0}^{\mathrm{\infty }}{F}_n^{\ast }(a)\cos n\psi .$$

Дифференцируя (4.3) по ф, на основании (4.2) получим:
$$f(a\cos \psi )a\omega \sin \psi =\sum _{n=0}^{\mathrm{\infty }}\omega n{F}_n^{\ast }(a)\sin n\psi .$$

Сопоставляя (4.4) с (1.16) п (1.17), находим:
$$A_1(a)=\frac{1}{2}{F}_1^{\ast }(a),B_1(a)=0,$$

откуда в первом приближении имеем:
$$x=a\cos \psi ,$$

где $a$ и $ф$ должны удовлетворять уравнениям
$$\begin{array}{l}\frac{da}{dt}=F_2^{\epsilon }{F}_1^{\ast }(a)\\ \frac{d\psi }{dt}=\omega .\end{array}\}$$

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

Исходя из выражений (3.2), (3.6), (3.7), (3.8) и (4.4) и учитывая, что
$$\begin{array}{l}\frac{1}{\pi }{\int }_0^{\pi }f(a^{-}\cos \psi )\cos \psi \sin \psi d\psi =0,\\ \frac{1}{\pi }{\int }_0^{\pi }f(a\cos \psi ){\cos }^2\psi d\psi =\frac{d{F}_1^{\ast }(a)}{da},\end{array}$$

можем написать:
$$x=a\cos \psi +\frac{\epsilon }{\omega }\sum _{n=2}^{\mathrm{\infty }}\frac{n{F}_n^{\ast }(a)\sin n\psi }{n^2-1},$$

где $a$ и $\psi$ определяются уравнениями
$$\begin{array}{l}\frac{da}{dt}=\frac{\epsilon }{2}{F}_1^{\ast }(a)\\ \frac{d\psi }{dt}=\omega +{\epsilon }^2{B}_2(a)\end{array}\}$$

а $B_2(a)$ имеет следующий вид:
$$B_2(a)=-\frac{1}{8a\omega }{F}_1^{\ast }(a)\frac{d{F}_1^{\ast }(a)}{da}-\frac{1}{2\omega a^2}\sum _{n=2}^{\mathrm{\infty }}\frac{n^2{F}_n^{\ast 2}(a)}{n^2-1},$$

Сопоставляя полученные приближенные решения с решениями уравнения (3.1), найденными в предыдущем параграфе, убеждаемся в их полной идентичности.

Таким образом, система, описываемая уравнением (4.1), также является квазиизохронной.
В качестве примера рассмотрим уравнение Ван-дер-Поля:
$$\frac{d^{2}x}{d{t}^2}-\epsilon (1-x^2)\frac{dx}{dt}+x=0.$$

Сопоставляя (4.10) и (4.1), имеем:
$$f(x)=1-x^2,$$

и потому
$$F^{\ast }(x)=x-\frac{x^3}{3},$$

после чего находим разложение (4.3) для нашего случая:
$$F^{\ast }(a\cos \psi )=a(1-\frac{a^2}{4})\cos \psi -\frac{a^3}{12}\cos 3\psi ,$$

согласно которому получим:
$$\begin{array}{l}{F}_1^{\ast }(a)=a(1-\frac{a^2}{4}),\\ F_3^{\ast }(a)=-\frac{a^3}{12},\\ F_n^{\ast }(a)=0,\text{ если }neq1,neq3.\end{array}\}$$

Таким образом, учитывая (4.6), в первом приближении имеем:
$$x=a\cos \psi \text{, }$$

где $a$ и $\psi$ должны быть определены из системы уравнений:
$$\begin{array}{l}\frac{da}{dt}=\frac{\epsilon a}{2}(1-\frac{a^2}{4}),\\ \frac{d\psi }{dt}=1.\end{array}\}$$

Итак, в первом приближении получаем гармоническое колебание, обладающее постоянной частотой $\omega =1$, амплитуда которого изменяется в соответствии с первым дифференциальным уравнөнием системы (4.13). Чтобы найти в явном виде закон зависимости амплитуды колебания от времени, необходимо решить это уравнение. Умножая обе его части на $a$, имеем:
$$\frac{d{a}^2}{dt}=\epsilon (1-\frac{a^2}{4})a^2,$$

откуда
$$\frac{d{a}^2}{(1-\frac{a^2}{4})a^2}=\epsilon dt,$$

или
$$\frac{d{a}^2}{4-a^2}+\frac{d{a}^2}{a^2}=\epsilon dt,$$

что дает
$$\ln \frac{a^2}{4-a^2}=\ln \frac{a_0^2}{4-a_0^2}+\epsilon t,$$

где $a_0$ — начальное значение амплитуды.
Из (4.15) окончательно находим:
$$a=\frac{a_0{e}^{\frac{1}{2}\epsilon t}}{\sqrt{1+\frac{1}{4}{a}_0^2(e^{8t}-1)}}.$$

Подставив (4.16) в (4.12), имеем выражение для первого приближения в явном виде:
$$x=\frac{a_0{e}^{\frac{1}{2}\epsilon t}}{\sqrt{1+\frac{1}{4}{a}_0^2(e^{et}-1)}}\cos (\omega t+\theta ).$$

Как видно из (4.17), если начальное значение амплитуды $a_0$ равно нулю, то амплитуда останется равной нулю для любого $t$, и мы получим $x=0$, т. е. тривиальное решение уравнения Ван-дер-Поля. Это тривиальное решение, очевидно, соответствует статическому режиму, т. е. отсутствию колебаний в системе.

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

Из (4.17) также замечаем, что если $a_0=2$, то $a=2$ для пюбых $t⩾0$. Это решение соответствует стационарному (установивпемуся) динамическьму режиму:
$$x=2\cos (t+\theta ).$$

В отличие от статического динамический режим обладает сильной устойчивостью, заключающейся в том, что каково бы ни было значение $a_{0}eq0$, малое или большое, все равно $a(t)\to 2$ при $t\to \mathrm{\infty }$.

Иначе говоря, любое колебание при увеличении $t$ приближается к стационарному колебанию (4.18).

Заметим, что только в первом приближении можно представить стационарный режим (4.18) как гармоническое колебание с частотой $\omega =1$ и амплитудой, равной 2 . В действительности же стационарный режим не гармонический.

Перейдем теперь к построению второго приближения. Согласно (4.7), (4.8) и (4.11) находим:
$$x=a\cos \psi -\frac{\epsilon a^3}{32}\sin 3\psi ,$$

где $a$ и $\psi$ должны быть определены из уравнений
$$\begin{array}{l}\frac{da}{dt}=\frac{\epsilon a}{2}(1-\frac{a^2}{4}),\\ \frac{d\psi }{dt}=1-{ϵ}^2(\frac{1}{8}-\frac{a^2}{8}+\frac{7a^4}{256}).\end{array}\}$$

Для стационарных колебаний во втором приближении получим:
$$x=2\cos (\omega t+\theta )-\frac{8}{4}\sin 3(\omega t+\theta ),$$

причем
$$\omega =1-\frac{{\epsilon }^2}{16}\cdot$$

На рассмотренном нами простом примере колебательной самовозбуждающейся системы, описываемой уравнением Ван-дер-Поля, замечаем коренное отличие этой системы от колебательных консервативных систем, описываемых уравнением вида (2.1).

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

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

Постоянное же значение амплитуда будет сохранять только в том случае, когда оба упомянутые количества энергии точно уравновешивают друг друга:

Построим теперь приближенные решения для уравнения Ван-дерПоля, воспользовавшись принципом усреднения.

Для этого необходимо уравнение (4.10) привести к стандартной форме. Это легко сделать, если вместо неизвестной функции $x$ ввести две новые функции $a$ и $\theta$ посредством следующих формул замены переменных:
$$\begin{array}{c}x=a\cos (t+\theta ),\\ \frac{dx}{dt}=-a\sin (t+\theta ).\end{array}$$

Дифференцируя (4,22) и сравнивая с (4.23), получаем:
$$\frac{da}{dt}\cos (t+\theta )-a\frac{d\theta }{dt}\sin (t+\theta )=0.$$

Дифференцируя (4,23) и принимая во внимание (4.22), (4.10), имеем:
$$\frac{da}{dt}\sin (t+\theta )+a\frac{d\theta }{dt}\cos (t+\theta )=\epsilon [1-a^2{\cos }^2(t+\theta )]a\sin (t+\theta ).$$

Разрешив соотношения (4.24), (4.25) относительно производных, приходим к системе двух уравнений в стандартной форме:
$$\begin{array}{l}\frac{da}{dt}=\epsilon [1-a^2{\cos }^2(t+\theta )]a{\sin }^2(t+\theta ),\\ \frac{d\theta }{dt}=\epsilon [1-a^2{\cos }^2(t+\theta )]\sin (t+\theta )\cos (t+\theta ),\end{array}\}$$

или
$$\begin{array}{l}\frac{da}{dt}=s\{\frac{a}{2}(1-\frac{a^2}{4})-\frac{a}{2}\cos 2(t+\theta )+\frac{a^3}{8}\cos 4(t+\theta )\},\\ \frac{d\theta }{dt}=\epsilon \{\frac{1}{2}(1-\frac{a^2}{2})\sin 2(t+\theta )-\frac{a^2}{8}\sin 4(t+\theta )\}.\end{array}\}$$

Применяя принцип усреднения, получаем в первом приближении
$$a=a_1,\theta ={\theta }_1,$$

причем
$$\frac{d{a}_1}{dt}=\frac{s{a}_1}{2}(1-\frac{a_1^2}{4}),\frac{d{\theta }_1}{dt}=0,$$

так как
$$\underset{t}{M}\{\cos 2(t+\theta )\}=\underset{t}{M}\{\sin 2(t+\theta )\}=\underset{t}{M}\{\cos 4(t+\theta )\}=\underset{t}{M}\{\sin 4(t+\theta )\}=0.$$

Как видно, уравнения первого приближения (4.28) совпадают с полученными выше уравнениями (4.13).
Улучшенное первое приближение, очевидно, будет:
$$\begin{array}{l}a=a_1-\frac{\epsilon a_1}{4}\sin 2(t+{\theta }_1)+\frac{\epsilon a_1^3}{32}\sin 4(t+{\theta }_1),\\ \theta ={\theta }_1-\frac{s}{4}(1-\frac{a_1^2}{2})\cos 2(t+{\theta }_1)+\frac{\epsilon a_1^2}{32}\cos 4(t+{\theta }_1).\end{array}\}$$

Для стационарного режима, как и выше, имеем:
$$a(t)\to 2\text{ при }t\to \mathrm{\infty },$$

и, следовательно, для установившегося колебательного режима при $a_1=2$ формулы (4.29) улучшенного первого приближения дают:
$$\begin{array}{l}a=2-\frac{s}{2}\sin 2(t+{\theta }_1)+\frac{\epsilon }{4}\sin 4(t+{\theta }_1),\\ \theta ={\theta }_1+\frac{s}{4}\cos 2(t+{\theta }_1)+\frac{\epsilon }{8}\cos 4(t+{\theta }_1).\end{array}\}$$

Подставляя эти значения в (4.22), получим:
$$\begin{array}{rl}x=& [2-\frac{\epsilon }{2}\sin 2(t+{\theta }_1)+\\ & +\frac{\epsilon }{4}\sin 4(t+{\theta }_1)]\cos (t+{\theta }_1+\frac{\epsilon }{4}\cos 2(t+{\theta }_1)+\frac{\epsilon }{8}\cos 4(t+{\theta }_1)),\end{array}$$

или, пренебрегая членами второго порядка малости, после элементарных преобразований получаем улучшенное приближение
$$x=2\cos (t+{\theta }_1)-\frac{s}{4}\sin 3(t+{\theta }_1),$$

которое совпадает с выражением улучшенного приближения, найденным нами ранее.

Прежде чем перейти к специальному рассмотрению стационарных амплитуд и их устойчивости, рассмотрим еще один классический пример автоколебательной системы — колебания часового маятника, возбуждаемого импульсами.
В этом случае имеем уравнение
$$J\frac{d^{2}x}{d{t}^2}+\{\lambda \frac{dx}{dt}-I\frac{\frac{dx}{dt}+\left|\frac{dx}{dt}\right|}{2}\delta (x-x_0)\}+kx=0,$$

где $x_0$-значение угла отклонения маятника, при котором на маятник посылается импульс $I,\delta (x)$ — «несобственная функция», определяемая соотнопениями:
$$\begin{array}{rl}{\int }_{-0}^{+0}\delta (x)dx& =1,\\ \delta (x)=0& \text{ при }xeq0.\end{array}\}$$

Чтобы привести уравнение (4.33) к виду (1.1), полагаем
$$\begin{array}{c}{\omega }^2=\frac{k}{J},\\ \epsilon f(x,\frac{dx}{dt})=\frac{\lambda }{J}\frac{dx}{dt}-\frac{I}{J}\frac{\frac{dx}{dt}+\left|\frac{dx}{dt}\right|}{2}\delta (x-x_0).\end{array}$$

На основании (4.34) и (4.35) можем написать:
$$\begin{array}{c}\epsilon {\int }_0^{2\pi }f(a\cos \psi ,-a\omega \sin \psi )\sin \psi d\psi =\frac{1}{J}{\int }_0^{2\pi }\lambda a\omega {\sin }^2\psi d\psi -\\ -\frac{Ia\omega }{J}{\int }_0^{2\pi }\delta (a\cos \psi -x_0)\frac{\sin \psi +|\sin \psi |}{2}\sin \psi d\psi =\\ =\frac{\lambda a\omega \pi }{J}-\frac{Ia\omega }{J}{\int }_{-\frac{\pi }{2}}^{\frac{\pi }{2}}\delta (a\cos \psi -x_0){\sin }^2\psi d\psi .\end{array}$$

Пусть фа является корнем уравнения
$$a\cos \psi =x_0,$$

лежащим между 0 и $\frac{\pi }{2}$.
Тогда для $a⩾x_0$ имеем:
$$\begin{array}{rl}{\int }_{-\frac{\pi }{2}}^{\frac{\pi }{2}}\delta (a\cos \psi -x_0){\sin }^2\psi d\psi & ={\int }_{{\psi }_a^{-0}}^{{\psi }_a+0}\delta (a\cos \psi -x_0){\sin }^2\psi d\psi =\\ & =\frac{1}{a}{\int }_{x_0-0}^{x_0+0}\delta (x-x_0)\sin {\psi }_{a}dx=\frac{\sin {\psi }_a}{a}=\frac{1}{a}\sqrt{1-\frac{x_0^2}{a^2}}.\end{array}$$

Если же $a<x_0$, то, очевидно,
$${\int }_{-\frac{\pi }{2}}^{\frac{\pi }{2}}\delta (a\cos \psi -x_0){\sin }^2\psi d\psi =0.$$

Таким образом, на основании (4.36), (1.24) и (1.27) получаем уравнения для определения мгновенной амплитуды в первом ґпрближении
$$\begin{array}{r}\frac{da}{dt}=-\frac{\lambda }{2J}a,\text{ если }a<x_0;\\ \frac{da}{dt}=-\frac{\lambda }{2J}a+\frac{I}{2\pi J}\sqrt{1-\frac{x_0^2}{a^2}},\text{ если }a⩾x_0,\end{array}\}$$

и совершенно аналогично уравнение для полной фазы
$$\frac{d\psi }{dt}=\omega (a),$$

тде
$$\dot{\omega }(a)=\{\begin{array}{lll}\omega ,& \text{ если }& a<x_0;\\ \omega -\frac{I{x}_0}{2\pi J{a}^2},& \text{ если }& a>x_0.\end{array}$$

Из уравнения (4.38) следует, очевидно, что при достаточно малом начальном значении амплитуды $a_0$, например при
$$a⩽x_0,$$

амплитуда $a(t)$ будет убывать:
$$a(t)\to 0\text{ при }t\to \mathrm{\infty },$$

и колебания с течением времени затухнут.
Данная колебательная система не имеет, следовательно, самовозбуждения, и статический режим (равновесие)
$$a=0$$

является устойчивым.
Однако мы можем поставить вопрос о существовании стационарных (установившихся) динамических режимов, соответствующих постоянным, не равным нулю значениям амплитуды. Очевидно, что такие значения амплитуды $а$ должны удовлетворять уравнению
$$-\frac{\lambda }{2J}a+\frac{I}{2\pi J}\sqrt{1-\frac{x_0^2}{a^2}}=0,$$

так как для стационарных амплитуд $\frac{da}{dt}=0$.
Таким образом, вопрос о существовании стационарных динамических режимов связан с вопросом о существовании вещественных положительных корней у уравнения стационарных амплитуд (4.39). Рассмотрим это уравнение.
Переписывая его в форме
$$1-\frac{x_0^2}{a^2}={\left(\frac{\pi \lambda }{I}\right)}^2{a}^2,$$

получим биквадратное уравнение для $a$ :
$$a^4-{\left(\frac{I}{\pi \lambda }\right)}^2{a}^2+{\left(\frac{I}{\pi \lambda }\right)}^2{x}_0^2=0,$$

из которого находим два корня:
$$\begin{array}{l}{a}_1^2=\frac{1}{2}{\left(\frac{I}{\pi \lambda }\right)}^2-\sqrt{\frac{1}{4}{\left(\frac{I}{\pi \lambda }\right)}^4-x_0^2{\left(\frac{I}{\pi \lambda }\right)}^2};\\ a_2^2=\frac{1}{2}{\left(\frac{I}{\pi \lambda }\right)}^2+\sqrt{\frac{1}{4}{\left(\frac{I}{\pi \lambda }\right)}^4-x_0^2{\left(\frac{I}{\pi \lambda }\right)}^2}.\end{array}$$

Корни эти будут вещественны, если
$$x_0<\frac{1}{2}\frac{I}{\pi \lambda }.$$

Если
$$x_0>\frac{1}{2}\frac{I}{\pi \lambda },$$

то оба корня (4.40) и (4.41) будут комплексными, и поэтому в случае (4.43) стационарные динамические режимы невозможны. В этом случае единственным возможным стационарным режимом может быть положение равновесия $a=0$.

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

При выполнении условия (4.42), как это следует из (4.38), процесе изменения амплитуды будет происходить следующим образом:
$$\begin{array}{rlrl}\text{ если }& a(0)<a_1,& \text{ то }a(t)& \to 0;\\ & & \to \mathrm{\infty }\\ \text{ если }& a(0)>a_1,& \text{ то }a(t)& \to a_2,\\ t& \to \mathrm{\infty }\end{array}$$

причем если $a(0)>a_2$, то амплитуда монотонно убывает до значения $a_{\epsilon }$, а если $a_1<a(0)<a_2$, то амплитуда монотонно возрастает, стремясь к тому же значению $a_2$.

Итак, если параметры колебательной системы удовлетворяют условию (4.42), в системе возможно устойчивое стационарное колебание с постоянной амплитудой:
$$a_2=\sqrt{\frac{1}{2}{\left(\frac{I}{\pi \lambda }\right)}^2+\sqrt{\frac{1}{4}{\left(\frac{I}{\pi \lambda }\right)}^4-x_0^2{\left(\frac{I}{\pi \lambda }\right)}^2}},$$

но для того, чтобы возбудить это колебание, необходимо, чтобы начальная амплитуда колебаний была больше $a_1$, т. е. необходим достаточно сильный начальный толчок, который создал бы отклонение маятника, превышающее $a_1$.
Например, если начальные условия при $t=0$ будут
$$x=0,\frac{dx}{dt}=A,$$

то мы видим, что
$${\psi }_0=\frac{3\pi }{2},A=a_0\omega$$

и, следовательно, условие возбуждения колебаний будет
$$A>a_1{}^{\omega }.$$