I am solving a physical problem with a known periodic solution. When I simulate the behaviour of the system numerically, with full blown differential equations, I get stable, but rather complicated looking oscillations about a certain point in parameter space. I constructed a simplified set of differential equations by perturbing the stable position by a small change and then obtained its equations of motion.
I have 2 coordinates $\epsilon$ and $\beta$, 1 characteristic angular velocity $\Omega$ and a small constant $q$ ($0 \lt q \ll 1$). The set looks as follows:
$\ddot{\epsilon} = 2 \Omega \dot{\beta} + 3 \Omega^2 \epsilon - \frac{9}{4} q \Omega^2 \beta$
$\ddot{\beta} = -2 \Omega \dot{\epsilon} -\frac{9}{4} q \Omega^2 \epsilon + \frac{9}{4} q \Omega^2 \beta$
I was wondering whether it can be shown that it allows for periodic solutions (or stable) and what their form could be. If it only were for the second derivatives and the actual variables, I could try to diagonalise the system, but I am confused by the first derivatives.
Thanks a lot.
SSF
A system of the form $$ \mathbf{\ddot{x}}+\mathbf{A}\mathbf{\dot{x}}+\mathbf{B}\mathbf{x}=0\,, $$ where $\mathbf{x}$ is a vector and $\mathbf{A}, \mathbf{B}$ are matrices, can be solved by substituting $\mathbf{x}=\mathbf{x_0} e^{\lambda t}$ with constant $\mathbf{x_0}$. Then you are solving $$ \left(\lambda^2 \mathbf{I}+\lambda \mathbf{A}+\mathbf{B}\right)\mathbf{x_0}=0\,. $$ In two dimensions you get four linearly independent solutions. You find the $\lambda$s by solving $$ \det\left(\lambda^2 \mathbf{I}+\lambda \mathbf{A}+\mathbf{B}\right)=0\,, $$ which in your case is a quartic. The system is stable if $\Re{\lambda}\leq0$ for every root.