Anharmonischer Oszillator

Der anharmonische Oszillator i​st ein schwingungsfähiges physikalisches System, b​ei dem d​ie Rückstellkraft n​icht proportional z​ur Auslenkung a​us der Ruhelage ist. Das h​at zur Folge, d​ass die f​reie Schwingung n​icht streng sinusförmig verläuft. Mechanische Beispiele s​ind etwa Pendel (Anharmonizität bemerkbar b​ei größerer Auslenkung), transversal schwingende Saiten, Kippschwingungen (das Kippeln e​ines aufrecht stehenden Gegenstands), Hüpfen e​ines Balls a​uf einer ebenen Fläche.

Bei genauer Betrachtung s​ind praktisch a​lle realen schwingungsfähigen Systeme anharmonisch. Die meisten nähern s​ich aber e​inem harmonischen Oszillator an, j​e kleiner d​ie Auslenkungen a​us der Ruhelage sind, w​eil dann d​ie Näherung e​iner linearen Rückstellkraft i​mmer besser zutrifft (für d​ie Mechanik s​iehe Hookesches Gesetz, Mathematisches Pendel). Bei solchen anharmonischen Oszillatoren verlaufen kleine f​reie Schwingungen näherungsweise sinusförmig u​nd mit e​iner bestimmten Eigenfrequenz, d​er Grundfrequenz d​es Oszillators.

Beim anharmonischen Oszillator treten i​m Vergleich z​um harmonischen Oszillator grundsätzlich n​eue Phänomene auf:

• Die Abweichung von der Sinusform bedeutet, dass die Schwingung auch Oberschwingungen (akustisch: Obertöne) der Grundfrequenz enthält. Bei harmonischer Anregung mit der Grundfrequenz klingen sie auch nach dem Einschwingvorgang nicht ab, sondern schwingen zusammen mit der Grundschwingung weiter.
• Bei unsymmetrischem Kraftgesetz verschiebt sich der Mittelpunkt der Schwingung gegenüber der Ruhelage. Dies tritt z. B. bei den Schwingungen zwischen den Atomen der festen Körper auf und ist die Ursache von deren thermischer Ausdehnung.
• Bei periodischer Anregung sind in der entstehenden erzwungenen Schwingung nicht nur, wie beim harmonischen Oszillator, die Anregungs- und die Eigenfrequenz enthalten, sondern zusätzlich deren Differenzfrequenz und Summenfrequenz sowie alle anderen ganzzahligen Kombinationen davon. Technische Anwendung findet dies z. B. in der nichtlinearen Optik bei der Frequenzverdopplung bei Laserlicht.
• Bei äußerer Anregung mit einer periodischen Kraft kann der anharmonische Oszillator auch mit einer chaotischen Bewegung reagieren, wenn die Anfangsbedingungen entsprechend gewählt sind oder seine Parameter außerhalb bestimmter Grenzen liegen (z. B. bei zu geringer Dämpfung). Siehe dazu: Doppelpendel

Bewegungsgleichung

Illustration der anharmonischen Kraftgesetze F(x) in Beispiel A (rot) und B (grün) für den Fall k=1 und l=0,2, im Vergleich zum harmonischen Kraftgesetz F(x)=-kx (blau). Außerdem sind die zugehörigen Potentiale V(x) mit F(x)=-dV/dx im grauen Kasten dargestellt.

Trajektorien des harmonischen Oszillators und der zwei diskutierten anharmonischen Oszillatoren im Vergleich für die Parameter m=k=1, l=0,2, c=0 und verschiedene Anfangsbedingungen ${\displaystyle x(0)}$ und ${\displaystyle m\cdot {\dot {x}}(0)=p(0)}$. Die rechte Spalte zeigt Phasenraumportraits der drei Oszillatoren.

Es gibt eine Vielzahl unterschiedlicher anharmonischer Oszillatoren und entsprechend viele Bewegungsgleichungen. Außerdem gibt es Unterschiede, je nachdem ob das Problem aus der Mechanik, Akustik oder Elektrik stammt. Ihnen gemeinsam ist, dass ihre Rückstellkraft nicht wie beim harmonischen Oszillator gemäß ${\displaystyle F(x)=-kx}$ nur linear von der Auslenkung ${\displaystyle x}$ abhängt, sondern auch von höheren Potenzen von ${\displaystyle x}$ ("nichtlineares Kraftgesetz"). In der nebenstehenden Abbildung sind die einfachsten Beispiele zusammen mit dem genäherten linearen Kraftgesetz (blaue Gerade) dargestellt:

A) Ein einfaches asymmetrisches Kraftgesetz (rote Kurve im Bild rechts) lautet ${\displaystyle F(x)=-kx-lx^{2}}$.

Eine Umformung ${\displaystyle F(x)=-kx(1+{\tfrac {l}{k}}x)}$ lässt erkennen, dass der nichtlineare Term die Rückstellkraft zur einen Seite hin zunehmend schwächt und zur anderen Seite hin stärkt, sofern die Auslenkungen nicht zu groß werden (${\displaystyle |{\tfrac {l}{k}}x|<1}$). Wie man an den Trajektorien ${\displaystyle x(t)}$ in der zweiten Abbildung sieht, führt dies zu einer im Vergleich zum harmonischen Oszillator stärkeren Auslenkung des Oszillators in die Richtung der abgeschwächten Rückstellkraft und zu einer verlangsamten Periode. Ein solches asymmetrisches Kraftgesetz kann etwa als zweite Näherung (nach dem linearen Kraftgesetz) an das unten erwähnte Morse-Potential in der Molekülphysik genutzt werden.

B) Ein einfaches symmetrisches Kraftgesetz (grüne Kurve im Bild rechts) lautet ${\displaystyle F(x)=-kx-lx^{3}}$.

Je nach Vorzeichen von ${\displaystyle l}$ wächst die Rückstellkraft, im Vergleich zum harmonischen Oszillator, zunehmend stärker (${\displaystyle l>0}$) oder schwächer (${\displaystyle l<0}$) an. Im ersten Fall erfolgen die Schwingungen bei gleicher Amplitude schneller, im zweiten Fall langsamer.

Zum Beispiel ist dies Kraftgesetz für ${\displaystyle l=-k/6}$ eine Näherung für das Schwere- oder mathematische Pendel, dessen Schwingungsperiode mit zunehmender Amplitude zunimmt. Die Größe ${\displaystyle x}$ steht in diesem Fall für den dimensionslosen Auslenkungswinkel des Pendels. Eine solche Näherung ergibt sich, wenn man die Rückstellkraft ${\displaystyle F(x)=-k\cdot \sin x}$ des mathematischen Pendels in einer Taylorreihe bis zur dritten Ordnung in ${\displaystyle x}$ entwickelt:

${\displaystyle F(x)=-k\cdot \sin x\approx \underbrace {-k\cdot x+{\frac {k}{6}}x^{3}} _{\text{1. und 3. Ordnung}}-{\frac {k}{120}}x^{5}+\ldots }$

Die Bewegungsgleichungen mit dem Dämpfungsterm ${\displaystyle c{\dot {x}}}$ ergeben sich dann aus dem Newton'schen Gesetz

${\displaystyle m{\ddot {x}}=-c{\dot {x}}+F(x)}$

wobei ${\displaystyle {\dot {x}}}$ die erste und ${\displaystyle {\ddot {x}}}$ die zweite Ableitung der Funktion ${\displaystyle x(t)}$ nach der Zeit ${\displaystyle t}$ bezeichnet, zu:

Fall A: ${\displaystyle m{\ddot {x}}+c{\dot {x}}+kx+lx^{2}=0}$

Fall B: ${\displaystyle m{\ddot {x}}+c{\dot {x}}+kx+lx^{3}=0}$

Für ${\displaystyle l=0}$ erhält man wieder in beiden Fällen das lineare Kraftgesetz sowie die Differentialgleichung des gedämpften harmonischen Oszillators. Auch für kleine Auslenkungen ist die Lösung des anharmonischen Oszillators wieder nahezu harmonisch, wie die ersten Trajektorien in der Abbildung rechts zeigen.

Durch die Nichtlinearität der Differentialgleichungen wird das Superpositionsprinzip außer Kraft gesetzt. Das bedeutet, dass nicht mehr jedes Vielfache ${\displaystyle \alpha \cdot x(t)}$ einer Lösung ${\displaystyle x(t)}$ auch eine Lösung der Differentialgleichung ist und allgemeiner, dass mit zwei Lösungen ${\displaystyle x_{1}(t),\;x_{2}(t)}$ nicht jede Linearkombination ${\displaystyle x(t)=\alpha _{1}\cdot x_{1}(t)+\alpha _{2}\cdot x_{2}(t)}$ auch eine Lösung ist (wobei ${\displaystyle \alpha ,\,\alpha _{1},\,\alpha _{2}}$ beliebige feste Zahlen sind). Die Lösung der Bewegungsgleichung ist meist ein elliptisches Integral und daher in geschlossener Form mit elementaren Funktionen nicht darstellbar. Dieser Artikel konzentriert sich auf periodische Bewegungen des anharmonischen Oszillators. Dabei wird die Dämpfung stellenweise vernachlässigt, d. h. ${\displaystyle c=0}$ gesetzt. Nur so ergeben sich, falls keine äußere Kraft einwirkt, periodische Bewegungen im strengen Sinne. Zu den chaotischen Bewegungsformen siehe Chaosforschung und dort angegebene Fachliteratur.

Näherungslösung

Im Fall einer schwachen anharmonischen Störung, d. h. ${\displaystyle |kx|\gg |lx^{2}|}$ bzw. ${\displaystyle |kx|\gg |lx^{3}|}$, kann man die Lösung durch Störungsrechnung erhalten. Dazu gibt man ${\displaystyle x}$ in Form einer Potenzreihe eines Störparameters ${\displaystyle \epsilon =-l/m}$ an:

${\displaystyle {\begin{aligned}x(t)&=x_{0}(t)&+&\;\epsilon \cdot x_{1}(t)&+&\;\epsilon ^{2}\cdot x_{2}(t)&+&\ldots \\\omega &=\omega _{0}&+&\;\epsilon \cdot \omega _{1}&+&\;\epsilon ^{2}\cdot \omega _{2}&+&\ldots \end{aligned}}}$

Dabei ist das erste Glied die (den Anfangsbedingungen angepasste) Lösung für den harmonischen Fall ${\displaystyle l=0}$, z. B. ${\displaystyle x_{0}(t)=A\cos \omega _{0}t}$, wenn der Oszillator zur Zeit ${\displaystyle t=0}$ bei der Auslenkung ${\displaystyle x=A}$ mit Geschwindigkeit ${\displaystyle {\dot {x}}=0}$ freigegeben wird und die Grundfrequenz ${\displaystyle \omega _{0}={\sqrt {k/m}}}$ hat.

Nach Einsetzen von ${\displaystyle x(t)}$ in die Bewegungsgleichung, wobei ${\displaystyle l}$ durch ${\displaystyle \epsilon }$ ausgedrückt wird, ergibt sich eine Potenzreihe in ${\displaystyle \epsilon }$, deren Koeffizienten sämtlich gleich Null zu setzen sind. So erhält man Differentialgleichungen für die einzelnen Näherungsfunktionen ${\displaystyle x_{i}(t)}$, die rekursiv zu lösen sind. Konkret hat die Differentialgleichung für ${\displaystyle x_{i+1}(t)}$ die Form der Bewegungsgleichung für einen harmonischen Oszillator, der von einer externen Kraft, welche durch die vorherigen ${\displaystyle x_{i}}$ gegeben ist, zu erzwungenen Schwingungen angeregt wird.

Im Fall A f​olgt im ersten Schritt

${\displaystyle x(t)=A\cos \omega _{0}t+\epsilon {\frac {A^{2}}{6\omega _{0}^{2}}}(3-4\cos \omega _{0}t+\cos 2\omega _{0}t)+\epsilon ^{2}(\ldots )+\ldots }$

Hier tritt bereits die verdoppelte Grundfrequenz auf. Der mathematische Grund dafür lässt sich zurückverfolgen zum Auftreten des quadratischen Terms ${\displaystyle x^{2}(t)}$ in der Bewegungsgleichung, der durch die trigonometrische Identität ${\displaystyle 2\cos ^{2}\phi =\cos(2\phi )+1}$ linearisiert wird. Bei den weiteren Näherungen ergeben sich Summanden mit entsprechend höheren Vielfachen der Grundfrequenz, insgesamt also ein ganzes Spektrum von Oberschwingungen.

Im Fall B ergibt sich im ersten Schritt eine instabile Lösung, denn ${\displaystyle x_{1}(t)}$ enthält einen zu ${\displaystyle t\cos \omega _{0}t}$ proportionalen Term. Dieser lässt sich jedoch eliminieren, wenn ${\displaystyle \omega _{1}=-{\tfrac {3}{8}}{\tfrac {A^{2}}{\omega _{0}}}}$ gesetzt wird. So ergibt sich schon in 1. Näherung eine Abhängigkeit der Oszillationsfrequenz ${\displaystyle \omega }$ von der Amplitude ${\displaystyle A}$ (wie im Fall A erst im 2. Schritt).

Anwendungen

Real i​n technischen Geräten verbaute Federn weisen i​m Allgemeinen, zuweilen a​uch durch d​ie Konstruktion beabsichtigt, n​ur in gewissen Grenzen e​ine lineare Beziehung zwischen Rückstellkraft u​nd Auslenkung auf. Die Dynamik e​ines Systems m​it solchen Federn f​olgt dann d​en nichtlinearen Bewegungsgleichungen, w​ie sie o​ben eingeführt wurden.

Das Morse-Potential (blau) im Vergleich zum quadratischen Potential des harmonischen Oszillators (grün). Eingezeichnet sind auch die Energiestufen, die beim Harmonischen Oszillator äquidistant (${\displaystyle \hbar \omega }$) sind, beim Morsepotential hingegen mit zunehmender Energie immer weniger Abstand haben. (${\displaystyle D_{e}}$ ist die Bindungsenergie, gerechnet vom Minimum des Potentials. ${\displaystyle D_{0}}$ ist die tatsächlich zur Flucht aus der Potentialmulde benötigte Dissoziationsenergie, wegen der Nullpunktenergie (${\displaystyle \nu =0}$) kleiner als ${\displaystyle D_{e}}$.)

Wichtige Anwendungen für anharmonische Schwingungen finden s​ich etwa i​n der Molekülphysik b​ei der Schwingung zweiatomiger Moleküle, o​der in d​er Festkörperphysik b​ei wärmebedingten Schwingungen v​on Atomen. Die Anharmonizität bildet s​ich durch d​ie unterschiedlichen Effekte b​ei Annäherung (elektrostatische Abstoßung, teilweise d​urch die Elektronen d​er Atomhülle abgeschirmt, a​ber durch d​as Pauli-Prinzip verstärkt) u​nd Entfernung (Rückstellkraft d​urch die kovalente Bindung d​er Atome) d​er Atomrümpfe aus. Wie i​n der Abbildung rechts gezeigt, können solche Schwingungen e​twa in e​inem Morse-Potential berechnet werden.

Erzwungene anharmonische Schwingungen

Bei der Bewegung unter dem Einfluss einer zeitabhängigen äußeren Kraft ${\displaystyle F(t)}$ unterscheiden sich anharmonischer und harmonischer Oszillator grundsätzlich voneinander. Beispielsweise kann, auch nach Beendigung eines Einschwingvorgangs, der anharmonische Oszillator mit anderen Frequenzen schwingen als in der erregenden Kraft vertreten. Es kann auch bei langsamer Variation der Erregerfrequenz zu sprunghafter Änderung der Amplitude kommen. Bei gekoppelten anharmonischen Oszillatoren kann eine minimale, kaum messbare Veränderung der Frequenz oder Amplitude zu starken Veränderungen des Verhaltens aller gekoppelten Elemente führen. Außerdem treten in solchen gekoppelten Systemen weitere Effekte wie das Frequenzziehen und die Frequenzmitnahme auf. Diese Phänomene sind auch praktisch von großem Interesse, da reale Oszillatoren sich nur solange annähernd harmonisch verhalten, wie bestimmte Grenzen für Auslenkung und/oder Frequenz eingehalten werden.

Diese Phänomene lassen s​ich ganz allgemein a​uf die nicht-lineare Form d​er Bewegungsgleichung zurückführen, w​omit das Superpositionsprinzip h​ier nicht m​ehr angewendet werden kann. Einige d​er Konsequenzen:

• Es gibt kein allgemein gültiges Lösungsverfahren wie bei den erzwungenen Schwingungen des harmonischen Oszillators. Man muss die Bewegungsgleichung numerisch integrieren oder analytische Lösungen verwenden, die nur unter Näherungsannahmen zu gewinnen sind.
• Der Einschwingvorgang (bei periodischer Anregung) ist nicht durch die Überlagerung der stationären Schwingung mit einer freien Schwingung gegeben.
• Die eingeschwungene stationäre Schwingung ist nicht immer unabhängig von den Anfangsbedingungen.
• Nach Fouriertransformation bleiben die Bewegungsgleichungen für verschiedene Frequenzen miteinander gekoppelt.

Amplitudensprünge

Als Beispiel sei das symmetrische Kraftgesetz (der obige Fall B) untersucht: Die Bewegungsgleichung lautet (nach Division durch ${\displaystyle m}$, mit ${\displaystyle \gamma ={\frac {c}{2m}}}$ sowie ${\displaystyle f(t)={\frac {F(t)}{m}}}$):

${\displaystyle {\ddot {x}}+2\gamma {\dot {x}}+\omega _{0}^{2}x+\epsilon x^{3}=f(t)}$

Ausgehend v​on der Annahme e​iner stationären harmonischen Schwingung

${\displaystyle x(t)=A\cdot \sin \Omega t}$

ergibt s​ich daraus für d​ie verursachende Kraft

(**) ${\displaystyle \;f(t)=(\omega _{0}^{2}-\Omega ^{2})A\sin \Omega t\;+\;2\gamma \Omega A\cdot \cos \Omega t\;+\;\epsilon A^{3}\left({\tfrac {3}{4}}\sin \Omega t-{\tfrac {1}{4}}\sin(3\Omega t)\right)}$.

Der mit ${\displaystyle 3\Omega }$ oszillierende Anteil der Kraft rührt von der Umformung ${\displaystyle \sin ^{3}\phi \equiv {\tfrac {3}{4}}\sin \phi -{\tfrac {1}{4}}\sin(3\phi )}$ her. Dieser Anteil wird im Weiteren vernachlässigt. Die Kraft lässt sich dann näherungsweise zu

${\displaystyle f(t)=f_{0}\cdot \sin(\Omega t+\varphi )}$

zusammenfassen, so dass hier eine mit der Frequenz ${\displaystyle \Omega }$ harmonisch schwingende Kraft eine harmonische Schwingung derselben Frequenz erzeugt. Dabei ist ${\displaystyle \varphi }$ eine Phasenverschiebung und die Amplitude der Kraft durch

${\displaystyle f_{0}^{2}=\left((\omega _{0}^{2}-\Omega ^{2})A\;+\;{\tfrac {3}{4}}\epsilon A^{3}\right)^{2}\;+\;4\,\gamma ^{2}\Omega ^{2}A^{2}}$

gegeben. Diese Gleichung lässt sich zwar nicht in der für die Resonanzkurve üblichen Form ${\displaystyle A(\Omega )}$ umstellen. Jedoch kann man sie nach ${\displaystyle \Omega ^{2}}$ auflösen und erhält den Zusammenhang zwischen Erregerfrequenz ${\displaystyle \Omega }$ und stationärer Schwingungsamplitude ${\displaystyle A}$ (für die durch ${\displaystyle f_{0}}$ gegebene Kraftamplitude) in der Form:

${\displaystyle (\Omega _{1/2})^{2}=\omega _{0}^{2}-2\gamma ^{2}\;+\;{\tfrac {3}{4}}\epsilon A^{2}\;\pm \;{\sqrt {{\frac {f_{0}^{2}}{A^{2}}}+\;4\,\gamma ^{2}\left(\gamma ^{2}\omega _{0}^{2}-{\tfrac {3}{4}}\epsilon A^{2}\right)}}}$.

Die beiden Lösungen ${\displaystyle \Omega _{1/2}}$, die hier durch die Lösung einer quadratischen Gleichung entstehen, drücken aus, dass im Allgemeinen zwei Erregerfrequenzen zu gleich großer Amplitude der stationären Schwingung führen, wie es auch schon beim harmonischen Oszillator links und rechts der Resonanzspitze der Fall ist. Für den harmonischen Fall, ${\displaystyle \epsilon =0}$, stimmt diese Formel mit der Resonanzkurve der harmonischen erzwungenen Schwingungen überein, die ihr Amplitudenmaximum bei ${\displaystyle \Omega ^{2}=\omega _{0}^{2}-2\gamma ^{2}}$ hat und zu beiden Seiten hin symmetrisch abfällt. Neu ist beim anharmonischen Oszillator, dass die Resonanzfrequenz sich mit steigender Amplitude verschiebt (Term ${\displaystyle {\tfrac {3}{4}}\epsilon A^{2}}$ in der Formel). Dadurch kann sich im Schaubild der Resonanzkurve ${\displaystyle A(\Omega )}$ die ganze Resonanzspitze derartig krümmen, dass sie in bestimmten Frequenzbereichen eine S-förmige Gestalt annimmt, also trotz gleicher Kraftamplitude ${\displaystyle f_{0}}$ und Erregerfrequenz ${\displaystyle \Omega }$ bis zu drei verschiedene mögliche Werte für die stationäre Amplitude ${\displaystyle A}$ anzeigt. Wird bei langsamer, stetiger Variation der Erregerfrequenz solch ein Bereich erreicht, springt die Amplitude von einem Ast der Resonanzkurve auf einen anderen: Die Schwingung "kippt".

Subharmonische Anregung

Im vorigen Abschnitt wurde ein Beitrag zur Kraft, der mit der Frequenz ${\displaystyle 3\Omega }$ oszilliert, in der Gleichung (**) einfach weggelassen. Das ist nicht immer gerechtfertigt, denn dieser Beitrag kann bei bestimmten Bedingungen auch die Hauptrolle spielen. Wenn gilt, dass

${\displaystyle \gamma =0}$   sowie   ${\displaystyle (\omega _{0}^{2}-\Omega ^{2})\;+\;{\frac {3}{4}}\epsilon A^{2}=0}$,

dann verschwinden aus Gleichung (**) alle mit ${\displaystyle \Omega }$ periodischen Terme. Es bleibt:

${\displaystyle x(t)=A\cdot \sin \Omega t}$

ist e​ine Lösung z​ur externen Kraft

${\displaystyle f(t)=\epsilon A^{3}{\frac {1}{4}}\sin(3\Omega t)}$.

Beispiel: Das mathematische Pendel mit Grundfrequenz ${\displaystyle \omega _{0}={\sqrt {g/L}}}$ (${\displaystyle g}$ Erdbeschleunigung, ${\displaystyle L}$ Pendellänge), das Kraftgesetz genähert durch Wahl des Parameters ${\displaystyle \epsilon =-\omega _{0}^{2}/6}$, werde angetrieben durch eine externe Kraft

${\displaystyle F(t)=-{\frac {1}{24}}m\omega _{0}^{2}\cdot \sin(3\Omega t)\quad {\text{mit}}\quad \Omega ={\sqrt {7/8}}\omega _{0}}$.

Dann schwingt es mit Amplitude ${\displaystyle A=1\approx 57^{\circ }}$ und einer dreifach untersetzten Frequenz wie ${\displaystyle x(t)=A\sin(\Omega t)}$. Um dies Verhalten beobachten zu können, muss man allerdings entweder die richtigen Anfangsbedingungen treffen oder das Abklingen von zusätzlichen Eigenschwingungen abwarten, was wegen der Annahme einer vernachlässigbar geringen Dämpfung sehr lange dauern kann.

Intermodulation

Mit Intermodulation wird das Phänomen bezeichnet, dass der Oszillator bei Anregung mit zwei Frequenzen ${\displaystyle \Omega _{1}\neq \Omega _{2}}$ mit einer Schwingung antwortet, in der auch Kombinationsfrequenzen ${\displaystyle \Omega =|n\Omega _{1}\;\pm \;m\Omega _{2}|}$ (${\displaystyle n}$ und ${\displaystyle m}$ ganzzahlig) vertreten sind. Während des Einschwingvorgangs, der aber wegen der Dämpfung gewöhnlich nur kurz dauert, sind solche Frequenzvielfache und Kombinationsfrequenzen auch in Bezug auf die Grundfrequenz ${\displaystyle \omega _{0}}$ vorhanden. In der Akustik können sie als hörbare Töne auftreten, die ihre Ursache darin haben, dass das Trommelfell oder eine Lautsprechermembran über diejenige Auslenkung hinaus erregt wird, bis zu der das lineare Kraftgesetz anwendbar ist.

Literatur

Friedhelm Kuypers: Klassische Mechanik. 8. Auflage. Wiley-VCH, Weinheim 2007, ISBN 978-3-527-40721-7.

Dieter Guicking: Schwingungen (Kap. 5 Nichtlineare u​nd parametrische Systeme). Wiley-VCH, Weinheim 2016, ISBN 978-3-658-14136-3.

Siehe auch

• Double-Well Potential
• Duffing-Oszillator
• Van-der-Pol-System
• Methode der harmonischen Balance
• Harmonischer Oszillator
• Harmonischer Oszillator (Quantenmechanik)
• Störungstheorie (Quantenmechanik)