Past exam of the differential equations and applications course of the University of Oxford / 2008 / A / 1

Set $y=\dot{x}$. The phase-plane system is $\dot{x}=y$, $\dot{y}=-9x+ϵ(1-4x^2)y$. The origin is its only equilibrium, and its Jacobian matrix is $\left(\begin{array}{cc}0& 1\\ -9& ϵ\end{array}\right)$. ThusFor the stated small positive $ϵ$, the eigenvalues have positive real part and nonzero imaginary part: the origin is an unstable focus.

Use the Poincaré-Lindstedt method. Introduce $\theta =\Omega t$ and $x=x_0(\theta )+ϵx_1(\theta )+\cdots$. The leading equation gives ${\omega }_0=3$ and $x_0=A_0\cos \theta$, after fixing the arbitrary phase. At order $ϵ$,Here $4{\cos }^2\theta \sin \theta =\sin \theta +\sin 3\theta$. Periodicity requires the two resonant coefficients to vanish. With $A_0>0$, this givesTherefore the leading limit cycle is $x(t)=\cos (3t)+O(ϵ)$, with frequency $3+O({ϵ}^2)$. For example the nonresonant first correction contains $-ϵ\sin (9t)/24$; its homogeneous sine/cosine pieces depend on phase and amplitude conventions. The absence of a first-order frequency shift is a solvability condition, not an assumption.

Use the method of multiple scales with independent fast and slow variables $t$ and $\tau =ϵt$, and write $x_0=R(\tau )\cos \Theta$, $\Theta =3t+\Phi (\tau )$. At the next order, $x_{1,tt}+9x_1=-2x_{0,t\tau }-(4x_0^2-1)x_{0,t}$. Its resonant terms areSuppressing secular terms yieldsThe leading initial data are $\boxed{R_0(0)=4,{\Phi }_0(0)=0}$, up to adding an integer multiple of $2\pi$ to the phase. Thus $x(t)=R_0(ϵt)\cos [3t+{\Phi }_0(ϵt)]+O(ϵ)$. The amplitude decreases from four towards the stable value one, consistent with the selected limit cycle; the zero-amplitude state is unstable. Order-$ϵ$ corrections enforce the exact initial derivative after differentiating the slow amplitude term.