Bestimmung der Stabilität aus den Eigenwerten zweier DE zweiter Ordnung

Ich analysiere die Stabilität eines doppelten physikalischen Pendels und habe das dimensionslose Gleichungssystem bestimmt

$$\begin{align} \label{eq:dimensionless} \begin{aligned} & \ddot{\phi} \cos{(\theta - \phi)} + \dot{\phi}^2 \sin{(\theta - \phi)} + \alpha \ddot{\theta} + \beta \sin{\theta} = 0 \end{aligned} \\[5pt] \begin{aligned} & \ddot{\theta} \cos{(\theta - \phi)} - \dot{\theta}^2 \sin{(\theta - \phi)} + \gamma \ddot{\phi} + \sin{\phi} = 0 \end{aligned} \end{align}$$ Wo $\alpha$, $\beta$Und $\gamma$definiert sind $$\begin{align} \label{eq:dim_parameters} \alpha &= \frac{\big ( I_1 + (m_1 + m_2) L_{1}^2 \big )}{m_2 L_{1} L_{2}} \\[5pt] \beta &= \frac{L_1 (m_1 + m_2)}{m_2 L_{2}} \\[5pt] \gamma &= \frac{(I_2 + m_2 L_{2}^2)}{m_2 L_{1} L_{2}} \, . \end{align}$$ Einstellung $\ddot{\phi} = \ddot{\theta} = \dot{\phi} = \dot{\theta} = 0$, finden wir, dass die stationären Punkte des Systems bei sind $(0, 0), \, (0, \pi), \, (\pi, 0), \, (\pi, \pi)$.

An diesen Punkten möchte ich die Stabilität des Systems finden. So habe ich versucht, sie zu finden, bin mir aber etwas unsicher, ob meine Arbeit korrekt ist:

Ich habe das Gleichungssystem linearisiert

$$\begin{align} \begin{aligned} & \ddot{\phi} + \alpha \ddot{\theta} + \beta \theta = 0 \\[5pt] & \ddot{\theta} + \phi + \gamma \ddot{\phi} = 0 \, . \end{aligned} \end{align}$$ Umstellen und Ersetzen der Gleichungen ineinander $$\begin{align} \begin{aligned} \nonumber & \ddot{\theta} = \frac{\beta \gamma}{1 - \alpha \gamma} \theta - \frac{1}{1 - \alpha \gamma} \phi \end{aligned} \\[5pt] \begin{aligned} \nonumber & \ddot{\phi} = \frac{\alpha}{1 - \alpha \gamma} \phi - \frac{\beta}{1 - \alpha \gamma} \theta \, . \end{aligned} \end{align}$$ Dies als Matrixgleichung ausdrücken $$\begin{align} \begin{aligned} \label{eq:matrix} \frac{d}{dt} \begin{bmatrix} \theta \\ \dot{\theta} \\ \phi \\ \dot{\phi} \end{bmatrix} &= \begin{bmatrix} 0 & 1 & 0 & 0 \\ \frac{\beta \gamma}{1 - \alpha \gamma} & 0 & -\frac{1}{1 - \alpha \gamma} & 0 \\ 0 & 0 & 0 & 1 \\ -\frac{\beta}{1 - \alpha \gamma} & 0 & \frac{\alpha}{1 - \alpha \gamma} & 0 \end{bmatrix} \begin{bmatrix} \theta \\ \dot{\theta} \\ \phi \\ \dot{\phi} \end{bmatrix} \, . \end{aligned} \end{align}$$ Bestimmung der charakteristischen Gleichung für die Eigenwerte der Matrix $$\begin{align} \begin{aligned} \nonumber & \text{det}(A - \lambda I) = \begin{vmatrix} -\lambda & 1 & 0 & 0 \\ \frac{\beta \gamma}{1 - \alpha \gamma} & -\lambda & -\frac{1}{1 - \alpha \gamma} & 0 \\ 0 & 0 & -\lambda & 1 \\ -\frac{\beta}{1 - \alpha \gamma} & 0 & \frac{\alpha}{1 - \alpha \gamma} & -\lambda \end{vmatrix} = 0 \\[5pt] & \implies \lambda^4 - \frac{\alpha + \beta \gamma}{1 - \alpha \gamma} \lambda^2 + \frac{\alpha \beta \gamma - \beta}{(1 - \alpha \gamma)^2} = 0 \\[5pt] & \implies \lambda^4 - \frac{\alpha + \beta \gamma}{1 - \alpha \gamma} \lambda^2 - \frac{\beta}{1 - \alpha \gamma} = 0 \end{aligned} \end{align}$$ Das Auflösen nach den Eigenwerten unter Verwendung der quadratischen Formel ergibt $$\begin{align} \begin{aligned} \lambda^2 &= \frac{(\alpha + \beta \gamma) \pm \sqrt{(-(\alpha + \beta \gamma))^2 - 4 (1 - \alpha \gamma) (-\beta)}}{2 (1 - \alpha \gamma)} \\[5pt] &= \frac{(\alpha + \beta \gamma) \pm \sqrt{\alpha^2 - 2 \alpha \beta \gamma + \beta^2 \gamma^2 + 4 \beta}}{2 (1 - \alpha \gamma)} \\[5pt] &= \frac{(\alpha + \beta \gamma) \pm \sqrt{(\alpha - \beta \gamma)^2 + 4 \beta}}{2 (1 - \alpha \gamma)} \end{aligned} \end{align}$$

Ich bin mir nicht sicher, ob das funktioniert. Die Stabilität an diesen Punkten sollte sicherlich nicht von den Parametern abhängen. Ich fühle, dass$(0, 0)$sollte immer stabil sein.

Habe ich den Jacobi nicht richtig gefunden? Oder soll ich vielleicht die Tatsache ausnutzen, dass die Parameter positiv sind? Jede Hilfe oder Klärung wäre sehr willkommen.