Der Duffing-Oszillator ist ein nichtlinearer Oszillator. Er kann als Erweiterung des harmonischen Oszillators, dessen Potential das lineare hookesche Gesetz zu Grunde liegt, um eine kubische Rückstellkraft betrachtet werden. Sein Verhalten wird durch folgende Differentialgleichung mit den zeitlichen Ableitungen von x beschrieben:
- $ {\ddot {x}}+\delta {\dot {x}}+\alpha x+\beta x^{3}=\gamma \cos(\omega _{0}t) $
$ \delta $ ist die Dämpfung, $ \gamma ,\omega _{0} $ sind die Amplitude und Frequenz der Anregung, $ \alpha ,\beta $ sind systemspezifische Parameter, welche die nichtlineare, rücktreibende Kraft charakterisieren.
Duffing-Oszillator ohne Anregung
Die Zustandsraumdarstellung des ungetriebenen Duffing-Oszillators $ {\ddot {x}}+\delta {\dot {x}}+\alpha x+\beta x^{3}=0 $ ist
- $ {\begin{bmatrix}{\dot {x}}_{1}\\{\dot {x}}_{2}\\\end{bmatrix}}={\begin{bmatrix}x_{2}\\-\delta x_{2}-\alpha x_{1}-\beta x_{1}^{3}\\\end{bmatrix}} $
Für den stationären Fall gilt
- $ {\begin{bmatrix}0\\0\\\end{bmatrix}}={\begin{bmatrix}x_{2}\\-\delta x_{2}-\alpha x_{1}-\beta x_{1}^{3}\\\end{bmatrix}} $
und damit
- $ x_{2}=0\, $ und $ \alpha x_{1}+\beta x_{1}^{3}=0 $.
Die Gleichung liefert für $ x_{1}\; $ drei stationäre Lösungen
- $ x_{1_{0}}=0,x_{1_{1,2}}=\pm {\sqrt {-{\frac {\alpha }{\beta }}}} $
Diese sind nur dann reell, wenn $ \alpha \cdot \beta <0 $ ist. Zur Beurteilung, welche dieser stationären Lösungen stabil sind wird das Differentialgleichungssystem um diese Punkte linearisiert. Die Jacobi-Matrix des Systems
- $ {\textbf {J}}={\begin{bmatrix}0&1\\-\alpha -3\beta x_{1}^{2}&-\delta \\\end{bmatrix}} $
hat für $ x_{1_{0}}\; $ die Eigenwerte
- $ \lambda _{0}={\frac {-\delta \pm {\sqrt {\delta ^{2}-4\alpha }}}{2}} $
und für $ x_{1_{1,2}}\; $ die Eigenwerte
- $ \lambda _{1}={\frac {-\delta \pm {\sqrt {\delta ^{2}+8\alpha }}}{2}} $.
Die Bedingung $ \alpha \cdot \beta <0 $ liefert zwei Fälle.
Fall 1: $ \alpha >0\; $ und $ \beta <0\; $
- $ \lambda _{0}\; $ hat negative Realteile, d.h. dieser Punkt ist stabil.
- $ \lambda _{1}\; $ hat einen positiven Realteil, d.h. diese Punkte sind instabil.
Fall 2: $ \alpha <0\; $ und $ \beta >0\; $
- $ \lambda _{0}\; $ hat positive Realteile, d.h. dieser Punkt ist instabil.
- $ \lambda _{1}\; $ hat negative Realteile, d.h. diese Punkte sind stabil.
Die Differenzialgleichung
- $ {\ddot {x}}+\delta {\dot {x}}-ax+bx^{3}=0 $
mit $ \delta >0,a>0,b>0\; $ beschreibt den stabilen Duffing-Oszillator.
