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

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

Для определения приближенных решений, соответствующих установившемуся режиму в колебательной системе (стационарным колебаниям), необходимо приравнять правую часть уравнения (5.1) нулю, так как при стационарном режиме амплитуда постоянна и, следовательно, производная от нее равна нулю. Из полученного алгебраического уравнения находим стационарные значения амплитуды. Однако для построения приближенных решений, соответствующих непосредственно стационарным колебаниям, можно указать более простой способ, чем изложенный выше. Рассмотрим сначала уравнение консервативной колебательной системы (2.1), которое можно записать в виде
$$\frac{d^{2}x}{d{t}^2}+{\omega }^{2}x=sf(x).$$

Согласно результатам § 2 стационарное решение этого уравнения во втором приближении имеет вид:
$$\begin{array}{c}x=a\cos (\omega t+\phi )-\frac{\epsilon }{{\omega }^2}\sum _{\begin{array}{c}n=0\\ neq1\end{array}}^{\mathrm{\infty }}\frac{f_n(a)\cos n(\omega t+\phi )}{n^2-1},\\ {\omega }_{\mathrm{II}}^2(a)={\omega }^2-\epsilon \frac{f_1(a)}{a}+s^2\dots ,\end{array}\}$$

где $f_n(a)(n=0,1,2,\dots )$ — коэффициенты Фурье в разложении
$$\begin{array}{rl}f(a\cos \psi )& =\sum _{n=0}^{\mathrm{\infty }}{f}_n(a)\cos n\psi ,\\ f_n(a)& =\frac{1}{\pi }{\int }_0^{2\pi }f(a\cos \psi )\cos n\psi d\psi ,\end{array}$$
$a$ и $p$-постоянные интегрирования, определяющиеся начальными значениями.

Исходя из выражений (6.2), естественно для получения высших приближений, соответствующих стационарному режиму, воспользоваться следующим приемом.

Представим репение уравнения (6.1) в виде $x=z(\overline{\omega }t+\phi )$, где $z(\overline{\omega t}+\phi )$ — периодическая функция $\overline{\omega t}+\phi$ с периодом $2\pi$.

Заметим, что $x=z(\overline{\omega }t+\phi )$ будет удовлетворять уравнению только тогда, когда $z(\overrightarrow{\omega }t+\phi )$ удовлетворяет уравнению
$$\overline{{\omega }^2}\frac{d^{2}z}{d{\psi }^2}+{\omega }^{2}z=\epsilon f(z).$$

Решение уравнения (6.4) $z=z(\psi ),\psi =\overline{\omega }t+\phi$, а также выражение для частоты колебания $\overline{\omega }$ естественно искать в виде разложений
$$\begin{array}{rl}z(\psi )& =\sum _{n=0}^{\mathrm{\infty }}{\mathrm{e}}^n{z}_n(\psi ),\\ {\overline{\omega }}^2& =\sum _{n=0}^{\mathrm{\infty }}{\mathrm{s}}^n{\alpha }_n,\end{array}$$

коэффициенты которых определим, подставив (6.5) и (6.6) в (6.4) и приравняв коэффициенты при одинаковых степенях малого параметра $\epsilon$, причем потребуем, чтобы $z_n(\psi )$ были периодическими функциями $\psi$ с периодом $2\pi$.
Произведя подстановку, иолучаем следующие уравнения:
$$\begin{array}{l}{\alpha }_0\frac{d^2{z}_0}{d{\psi }^2}+{\omega }^2{z}_0=0,\\ {\alpha }_0\frac{d^2{z}_1}{d{\psi }^2}+{\omega }^2{z}_1=f\left(z_0\right)-{\alpha }_1\frac{d^2{z}_0}{d{\psi }^2},\\ {\alpha }_0\frac{d^2{z}_2}{d{\psi }^2}+{\omega }^2{z}_2=f^{\prime }\left(z_0\right)z_1-{\alpha }_2\frac{d^2{z}_0}{d{\psi }^2}-{\alpha }_1\frac{d^2{z}_1}{d{\psi }^2},\\ {\alpha }_0\frac{d^2{z}_3}{d{\psi }^2}+{\omega }^2{z}_3=f^{\prime }\left(z_0\right)z_2+\frac{1}{2}{f}^{\prime \prime }\left(z_0\right)z_1^2-{\alpha }_3\frac{d^2{z}_0}{d{\psi }^2}-{\alpha }_2\frac{d^2{z}_1}{d{\psi }^2}-{\alpha }_1\frac{d^2{z}_2}{d{\psi }^2},\\ \text{. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . }\end{array}$$

Определяя из первых $N+1$ уравнений системы (6.7) функции $z_0,z_1,z_2,\dots ,z_N$, а также величины ${\alpha }_0,{\alpha }_1,{\alpha }_2,\dots ,{\alpha }_N$, можем составить выражение
$$x=\sum _{n=0}^N{\epsilon }^n{z}_n(\overline{\omega }t+\phi ),$$

где
$${\overline{\omega }}^2=\sum _{n=0}^N{\epsilon }^n{a}_n,$$

которое будет удовлетворять уравнению (6.1) с точностью до величин порядка малости ${\epsilon }^{N+1}$ и, следовательно, может рассматриваться как $N+1$-е приближение решения уравнения (6.1), соответствующее стационарным колебаниям. Определение функций $z_n(\overline{\omega }t+\phi )$ и величин ${\alpha }_n(n=0,1,2,\dots )$ из уравнений (6.7) может быть произведено, вообще говоря, неоднозначно. Для того чтобы эти величины были определены однозначно, необходимо наложить некоторые дополнительные условия.

Потребуем, чтобы $z_n(\psi )(n=1,2,3,\dots )$ не имели в своем составе основной гармоники аргумента $\psi$. Из первого уравнения (6.7) находим:
$$z_0(\psi )=a\cos \psi ,{\alpha }_0={\omega }^2.$$

Подставляя (6.9) во второе уравнение системы (6.7), имеем:
$${\omega }^2(\frac{d^2{z}_1}{d{\psi }^2}+z_1)=f(a\cos \psi )+{\alpha }_{1}a\cos \psi$$

или, учитывая (6.3),
$${\omega }^2(\frac{d^2{z}_1}{d{\psi }^2}+z_1)=\sum _{\begin{array}{c}n=0\\ neq1\end{array}}^{\mathrm{\infty }}{f}_n(a)\cos n\psi +({\alpha }_{1}a+f_1(a))\cos \psi .$$

Принимая во внимание требование периодичности функции $z_1(\psi )$, приравниваем в правой части уравнения (6.11) нулю коэффициент при первой гармонике аргумента $\psi$. В результате получаем уравнение для определения ${\alpha }_1$ :
$$a_{1}a+f_1(a)=0,$$

откуда находим:
$${\alpha }_1=\frac{-f_1(a)}{a}.$$

Подставляя найденное значение ${\alpha }_1$ в правую часть уравнения (6.11), имеем:
$$\frac{d^2{z}_1}{d{\psi }^2}+z_1=\frac{1}{{\omega }^2}\sum _{\begin{array}{c}n=0\\ neq1\end{array}}^{\mathrm{\infty }}{f}_n(a)\cos n\psi$$

откуда находим для $z_1(\psi )$ следующее выражение:
$$z_1(\psi )=\frac{1}{{\omega }^2}\sum _{\begin{array}{c}n=0\\ neq1\end{array}}^{\mathrm{\infty }}\frac{f_n(a)\cos n\psi }{1-n^2}.$$

При этом как выражение (6.13), так и выражение для ${\alpha }_1$ совнадают с выражениями, найденными согласно общему методу.

Продолжая изложенный процесс, можем последовательно определить все функции $z_1,z_2,z_3,\dots$ и величины ${\alpha }_1,{\alpha }_2,{\alpha }_3,\dots$ до сколь угодно большого значения индекса и тем самым построить приближенные решения, удовлетворяющие рассматриваемому уравнению (6.1) с точностью до любой степени є.

В качестве примера определим стационарное решение в третьем улучшенном приближении (с точностью до величин порядка малости ${\epsilon }^3$ включительно) для колебательной системы, описываемой уравнением вида
$$\frac{d^{2}x}{d{t}^2}+x+\epsilon x^3=0.$$

Для определения функций $z_0,z_1,z_2,z_3$ и величин ${\alpha }_0,{\alpha }_1,{\alpha }_2,{\alpha }_3$ на основании (6.7) получаем уравнения:
$$\begin{array}{l}{\alpha }_0\frac{d^2{z}_0}{d{\psi }^2}+z_0=0\\ {\alpha }_0\frac{d^2{z}_1}{d{\psi }^2}+z_1=-z^3-{\alpha }_1\frac{d^2{z}_0}{d{\psi }^2}\\ {\alpha }_0\frac{d^2{z}_2}{d{\psi }^2}+z_2=-3z_0^2{z}_1-{\alpha }_2\frac{d^2{z}_0}{d{\psi }^2}-{\alpha }_1\frac{d^2{z}_1}{d{\psi }^2},\\ {\alpha }_0\frac{d^2{z}_3}{d{\psi }^2}+z_3=-3z_0^2{z}_2-3z_1{z}_0-{\alpha }_3\frac{d^2{z}_0}{d{\psi }^2}-{\alpha }_2\frac{d^2{z}_1}{d{\psi }^2}-{\alpha }_1\frac{d^2{z}_2}{d{\psi }^2}.\end{array}$$

Из первого уравнения находим:
$$z_0(\psi )=a\cos \psi ,{\alpha }_0=1.$$

После этого второе уравнение можно записать в виде
$$\frac{d^2{z}_1}{d{\psi }^2}+z_1=-a^3{\cos }^3\psi +{\alpha }_{1}a\cos \psi ,$$

или
$$\frac{d^2{z}_1}{d{\psi }^2}+z_1=-\frac{a^3}{4}\cos 3\psi +({\alpha }_{1}a-\frac{3}{4}{a}^3)\cos \psi ,$$

откуда имеем:
$$\begin{array}{rl}{\alpha }_1& =\frac{3}{4}{a}^2,\\ z_1(\psi )& =\frac{a^3}{32}\cos 3\psi .\end{array}\}$$

Подставляя (6.16) и (6.19) в третье уравнение системы (6.15), получаем:
$$\frac{d^2{z}_2}{d{\psi }^2}+z_2=-\frac{3a^5}{32}{\cos }^2\psi \cos 3\psi +\frac{27}{128}{a}^5\cos 3\psi +{\alpha }_{2}a\cos \psi ,$$

или
$$\frac{d^2{z}_2}{d{\psi }^2}+z_2=\frac{21}{128}{a}^5\cos 3\psi -\frac{3a^5}{128}\cos 5\psi +({\alpha }_{2}a-\frac{3a^5}{128})\cos \psi ,$$

откуда находим:
$$\begin{array}{rl}{\alpha }_2& =\frac{3a^4}{128},\\ z_2(\psi )& =-\frac{21}{1024}{a}^5\cos 3\psi +\frac{a^5}{1024}\cos 5\varphi .\end{array}\}$$

После этого последнее уравнение системы (6.15) может быть записано в виде:
$$\begin{array}{l}\frac{d^2{z}_3}{d{\psi }^2}+z_3=-3a^2{\cos }^2\psi (-\frac{21}{1024}{a}^5\cos 3\psi +\frac{a^5}{1024}\cos 5\psi )-\\ -3\frac{a^6}{1024}{\cos }^{2}3\varphi a\cos \psi +{\alpha }_{3}a\cos \psi +\frac{3a^4}{128}\frac{9a^3}{32}\cos 3\psi -\\ -\frac{3}{4}{a}^2(\frac{21}{1024}{a}^{5}9\cos 3\psi -\frac{a^5}{1024}25\cos 5\psi ),\end{array}$$

или
$$\begin{array}{rl}\frac{d^2{z}_3}{d{\psi }^2}+z_3=-\frac{1059a^7}{2048}\cos 3\psi +& \frac{177}{2048}{a}^7\cos 5\psi -\\ & -\frac{3}{2048}{a}^7\cos 7\psi +({\alpha }_{3}a+\frac{57}{4096}{a}^7)\cos \psi ,\end{array}$$

из которого находим:
$$\begin{array}{rl}{\alpha }_3& =-\frac{57}{4096}{a}^6,\\ z_3(\psi )& =\frac{1059}{2048\cdot 8}{a}^7\cos 3\psi -\frac{177}{2048\cdot 24}{a}^7\cos 5\psi +\frac{3}{2048\cdot 48}{a}^7\cos 7\psi .\end{array}$$

Таким образом, принимая во внимание (6.8), (6.19), (6.22) и (6.25), получаем приближенное стационарное решение уравнения (6.14) с точностью до величин порядка малости ${}^3$ включительно в виде
$$\begin{array}{l}x=a\cos (\overline{\omega }t+\phi )+\epsilon \frac{a^3}{32}[1-\epsilon \frac{21}{32}{a}^2+{\epsilon }^2\frac{1059}{512}{a}^4]\cos 3(\overline{\omega }t+\phi )+\\ +{\epsilon }^2\frac{a^5}{1024}[1-\epsilon \frac{177}{48}{a}^2]\cos 5(\overline{\omega }t+\phi )+{\epsilon }^3\frac{3a^7}{2048\cdot 48}\cos 7(\overline{\omega }t+\phi ),\end{array}$$

где $a$ и $\phi$-произвольные постоянные, а частота колебаний $\overline{\omega }$ определяется выражением
$$\overline{{\omega }^2}=1+\frac{3}{4}\approx a^2+\frac{3}{128}{\epsilon }^2{a}^4-\frac{57}{4096}{\epsilon }^3{a}^6.$$

Перейдем теперь к построению приближенных решений для стационарных колебаний в неконсервативных системах. Для этого рассмотрим уравнение вида (1.1)
$$\frac{d^{2}x}{d{t}^2}+{\omega }^{2}x=\epsilon f(x,\frac{dx}{dt}).$$

Стационарное решение (улучшенное первое приближение) этого уравнения согласно § 2 может быть записано в виде
$$x=a\cos \psi +\frac{g_0(a)}{{\omega }^2}-\frac{1}{{\omega }^2}\sum _{n=2}^{\mathrm{\infty }}\frac{g_n(a)\cos n\psi +h_n(a)\sin n\psi }{n^2-1},$$

где $g_n(a)$ и $h_n(a)(n=0,2,3,\dots )$ определяются согласно формулам (1.29), $\psi =\overline{\omega }(a)t+\phi$, а $a$ и $\overline{\omega }(a)$ должны быть найдены из следующих выражений:
$$\begin{array}{c}{A}_1(a)=0,\\ \overline{\omega }(a)=\omega +\epsilon B_1(a),\end{array}\}$$
$A_1(a)$ и $B_1(a)$ определяются согласно формулам (1.27).
Для случая консервативных колебательных систем, как мы уже убедились, $A_1(a)$ тождественно равно нулю, и поэтому выражение для приближенного стационарного решения (6.29) зависит от двух произвольных постоянных $a$ и $\phi$.

Рассмотрим теперь случай, когда $A_1(a)$ не обращается тождественно в нуль ни в каком интервале значений $a$. Предположим также, что функция $A_1(a)$ имеет лишь простые корни. В этом случае каждому корню $A_1(a)$ соответствует некоторое стационарное состояние, причем выражение (6.29) для этого стационарного состояния зависит лишь от одной шроизвольной постоянной, а именно от $\phi$.

Приступая к изложению формальной методики построения стационарных решений уравнения (1.1), воспользуемся в основном способом, примененным выше для консервативных колебательных систем.

Решение уравнения (6.28), соответствующее стационарным колебаниям, представим в виде
$$x=z(\overline{\omega }t+\phi )\text{, }$$

где $\phi$ — произвольная постоянная, $\overline{\omega }$ — частота колебания, $z(\psi )$ — периодическая функция $\psi$ с периодом $2\pi$.
Как и выше, функция $z(\psi )$ должна удовлетворять уравнению
$${\overline{\omega }}^2\frac{d^{2}z}{d{\psi }^2}+{\omega }^{2}z=\epsilon f(z,\overline{\omega }\frac{dz}{d\psi }).$$

Будем искать функцию $z(\psi )$ и частоту колебаний $\overline{\omega }$ в виде разложений
$$\begin{array}{rl}z(\psi )& =\sum _{n=0}^{\mathrm{\infty }}{\epsilon }^n{z}_n(\psi ),\\ \overline{\omega }& =\sum _{n=0}^{\mathrm{\infty }}{\epsilon }^n{\omega }_n,\end{array}$$

где $z_n(\psi )$ — периодические функции $\psi$ с периодом $2\pi$.
Подставляя значения $z(\varphi )$ и $\overline{\omega }$ (6.33) в уравнение (6.32) и приравнивая коэффициенты при одинаковых степенях ぇ, получим следующую систему уравнений:
$$\begin{array}{l}{\omega }_0^2\frac{d^2{z}_0}{d{\psi }^2}+{\omega }^2{z}_0=0,\\ {\omega }_0^2\frac{d^2{z}_1}{d{\psi }^2}+{\omega }^2{z}_1=f(z_0,{\omega }_0\frac{d{z}_0}{d\psi })-2{\omega }_0{\omega }_1\frac{d^2{z}_0}{d{\psi }^2},\\ {\omega }_0^2\frac{d^2{z}_2}{d{\psi }^2}+{\omega }^2{z}_2=f_z^{\prime }(z_0,{\omega }_0\frac{d{z}_0}{d\psi })z_1+f_{z^{\prime }}^{\prime }(z_0,{\omega }_0\frac{d{z}_0}{d\psi }){\omega }_0\frac{d{z}_1}{d\psi }+\\ +f_{z^{\prime }}^{\prime }(z_0,{\omega }_0\frac{d{z}_0}{d\psi }){\omega }_1\frac{d{z}_0}{d\psi }-2{\omega }_0{\omega }_2\frac{d^2{z}_0}{d{\psi }^2}-2{\omega }_0{\omega }_1\frac{d^2{z}_1}{d{\psi }^2}-{\omega }_1^2\frac{d^2{z}_0}{d{\psi }^2},\end{array}$$

Решаем первое из этих уравнений, полагая
$$z_0=a\cos \psi ,{\omega }_0=\omega ,$$

где $a$-некоторая, пока не определенная постоянная.
Подставив значения $z_0(\varphi )$ и ${\omega }_0(6.35)$ в правую часть второго уравнения системы (6.34) и учитывая (1.16), имеем
$${\omega }^2(\frac{d^2{z}_1}{d{\psi }^2}+z_1)=\sum _{n=0}^{\mathrm{\infty }}[g_n(a)\cos n\psi +h_n(a)\sin n\psi ]+2\omega {\omega }_{1}a\cos \dot{\psi }.$$

Чтобы это уравнение имело относительно $z_1(\psi )$ периодическое решение (т. е. чтобы в выражении для $z_1(\psi )$ не появилось секулярных членов), необходимо приравнять нулю коэффициенты при основной гармонике, входящей в правую часть (6.36).
Приравнивая, получаем уравнения
$$\begin{array}{rl}{h}_n(a)& =A_1(a)=0,\\ {\omega }_1& =-\frac{g_1(a)}{2a\omega },\end{array}\}$$

определяющие $a$ и ю ${}_1$. После этого уравнение (6.36) примет вид
$${\omega }^2(\frac{d^2{z}_1}{d{\psi }^2}+z_1)=g_0(a)+\sum _{n=2}^{\mathrm{\infty }}[g_n(a)\cos n\dot{\psi }+h_n(a)\sin n\psi ].$$

Решая өго, находим:
$$z_1(\psi )=a_1\cos \psi +\frac{g_0(a)}{{\omega }^2}+\frac{1}{{\omega }^2}\sum _{n=2}^{\mathrm{\infty }}\frac{g_n(a)\cos n\psi +h_n(a)\sin n\psi }{1-n^2},$$

где $a_1$ — произвольная пэстоянная, которая должна быть определена из условия периодичности $z_2(\varphi )$.

Покажем, как определяется $a_1$. Для этого выражение для $z_1(\psi )$ представим в виде
$$z_1(\psi )=a_1\cos \psi +u(a,\psi ),$$

где
$$u(a,\psi )=\frac{g_0(a)}{{\omega }^2}+\frac{1}{{\omega }^2}\sum _{n=2}^{\mathrm{\infty }}\frac{g_n(a)\cos n\psi +h_n(a)\sin n\psi }{1-n^2}$$

является известной периодической функцией $\psi$.
Подставляя значение $z_1(\psi )(6.40)$ в третье уравнение системы (6.34), получим:
$$\begin{array}{l}{\omega }^2(\frac{d^2{z}_2}{d{\psi }^2}+z_2)=f_z^{\prime }(a\cos \psi ,-a\omega \sin \psi )a_1\cos \psi -\\ -f_{z^{\prime }}^{\prime }(a\cos \psi ,-a\omega \sin \psi )a_1\omega \sin \psi +2\omega {\omega }_1{a}_1\cos \psi +\\ +2\omega {\omega }_{2}a\cos \psi +v(\psi ),\end{array}$$

где
$$\begin{array}{l}v(\psi )=f_z^{\prime }(a\cos \psi ,-a\omega \sin \psi )u+f_{z^{\prime }}^{\prime }(a\cos \psi ,-a\omega \sin \psi )\omega \frac{du}{d\psi }+\\ +f_{z^{\prime }}^{\prime }(a\cos \psi ,-a\omega \sin \psi ){\omega }_{1}a\sin \psi +{\omega }_1^{2}a\cos \psi .\end{array}$$

Напишем разложение для $v(\psi )$ в ряд Фурье:
$$v(\psi )=\sum _{n=0}^{\mathrm{\infty }}[v_n(a)\cos n\psi +w_n(a)\sin n\psi ].$$

Кроме того, имеем:
$f_z^{\prime }(a\cos \psi ,-a\omega \sin \psi )\cos \psi -f_{z^{\prime }}^{\prime }(a\cos \psi ,-a\omega \sin \psi )\omega \sin \psi =$
$$=\sum _{n=0}^{\mathrm{\infty }}[g_n^{\prime }(a)\cos n\psi +h_n^{\prime }(a)\sin n\psi ]$$

и поэтому уравнение (6.41) можно эаписать в виде
$$\begin{array}{l}{\omega }^2(\frac{d^2{z}_2}{d{\psi }^2}+z_2)=a_1\sum _{n=0}^{\mathrm{\infty }}[g_n^{\prime }(a)\cos n\psi +h_n^{\prime }(a)\sin n\psi ]+\\ +2\omega ({\omega }_1{a}_1+{\omega }_{2}a)\cos \psi +\sum _{n=0}^{\mathrm{\infty }}[v_n(a)\cos n\psi +w_n(a)\sin n\psi ].\end{array}$$

Для того чтобы $z_2(\psi )$ являлось периодической функцией, необходимо в правой части уравнения (6.45) приравнять нулю коэффидиенты шри $\sin \psi$ и $\cos \psi$. Получаем:
$$\begin{array}{rl}{a}_1{h}_n^{\prime }(a)& =-w_1(a),\\ a_1{g}_1^{\prime }(a)+2\omega ({\omega }_1{a}_1+{\omega }_{2}a)& =-v_1(a).\end{array}\}$$

Первое из этих уравнений определяет $a_1$ :
$$a_1=-\frac{w_1(a)}{h_1^{\prime }(a)},$$

причем $h_1^{\prime }(a)=A_1^{\prime }(a)eq0$, так как мы предположили, что функция $A_1(a)$ имеет только простые корни.
Из второго уравнения находим ${\omega }_2$ :
$${\omega }_2=-\frac{1}{a}[\frac{v_1(a)+a_1{g}_1^{\prime }(a)}{2\omega }+{\omega }_1{a}_1].$$

После этого уравнение (6.45) можем записать в виде
$$\begin{array}{r}{w}^2(\frac{d^2{z}_2}{d{\psi }^2}+z_2)=(v_0(a)+a_1{g}_0^{\prime }(a))+\sum _{n=2}^{\mathrm{\infty }}\{(v_n(a)+a_1{g}_n^{\prime }(a))\cos n\psi +\\ +(w_n(a)+a_1{h}_n^{\prime }(a))\sin n\psi \}.\end{array}$$

Решением этого уравнения будет:
$$\begin{array}{c}{z}_2(\psi )=a_2\cos \psi +\frac{1}{{\omega }^2}(v_0(a)+a_1{g}_0^{\prime }(a))+\frac{1}{{\omega }^2}\sum _{n=2}^{\mathrm{\infty }}\{(v_n(a)+a_1{g}_n^{\prime }(a))\cos n\psi +\\ +(w_n(a)+a_1{h}_n^{\prime }(a))\sin n\psi \}\frac{1}{1-n^2},\end{array}$$

где $a_2$-неопределенная постоянная, которую определяем из условия периодичности функции $z_3(\varphi )$ и т. д. Продолжая изложенный процесс, можем построить приближенные решения для стациопарного режима с любой наперед заданной степенью точности. Например, во втором приближении согласно (6.31), (6.39) и (6.37) имеем:
$$\begin{array}{l}x=(a+\epsilon a_1)\cos (\overline{\omega }t+\phi )+\frac{\epsilon }{{\omega }^2}{g}_0(a)+\\ +\frac{\epsilon }{{\omega }^2}\sum _{n=2}^{\mathrm{\infty }}\frac{g_n(a)\cos n(\overline{\omega }t+\phi )+h_n(a)\sin n(\overline{\omega }t+\phi )}{1-n^2},\end{array}$$

где
$$\overline{\omega }=\omega +\epsilon \frac{g_1(a)}{2\omega a},$$

а амплитуда должна быть определена из уравнения
$$h_n(a)=0\text{. }$$

Сравнивая полученное выражение для $x$ (6.51) с ранее найденным выражением для стационарного режима (6.29), замечаем, что единственное отличие этих двух репений состоит в том, что амплитуда первой гармоники в (6.29) равна $a$, где $a$-корень первого уравнения (6.30), а в формуле (6.51) амплитуда первой гармоники равна $a+\epsilon a_1$. Однако такое расхождение полностыю совпадает с нашим замечапием по этому поводу, приведенным в 1 .