# Van der Pol-vergelijking

In de dynamica beschrijft de Van der Pol-vergelijking een trilling met niet-lineaire demping, waarbij energie verloren gaat, die dus niet conservatief is. De oscillator wordt Van der Pol-oscillator genoemd en wordt beschreven door de tweede-ordedifferentiaalvergelijking
    ${\frac {\mathrm {d} ^{2}x}{\mathrm {d} t^{2}}}-\mu (1-x^{2}){\frac {\mathrm {d} x}{\mathrm {d} t}}+x=A\sin(\omega t)$
waarin de variabele $x$ de plaatscoördinaat is en $A\sin(\omega t)$ de voorgeschreven excitatie. Beide zijn een functie van de tijd $t$. De scalairen $A$, de amplitude, $\omega$ de hoeksnelheid en $\mu$ als maat voor de niet-lineaire demping, zijn de parameters van de differentiaalvergelijking.
De Van der Pol-vergelijking vindt haar oorsprong in de elektrodynamica, en is genoemd naar haar bedenker, de Nederlandse elektronicapionier Balthasar van der Pol (1889–1959). Van der Pol ontdekte dat het gedrag van een met een triode opgewekte elektrische trilling afhing van niet-lineaire termen in de karakteristiek van de triode. De triode was één van de eerste elektronische componenten. Hij stelde de differentiaalvergelijking op, waarmee dit verschijnsel, dat ook in veel andere niet-lineaire systemen optreedt, adequaat kon worden beschreven.
Met de stelling van Liénard kan bewezen worden, dat er een zogenaamde limietcyclus bestaat voor de niet-geëxciteerde Van der Pol-oscillator $A=0$, aldus resulterende als voorbeeld voor een Liénard-systeem. Uit de differentiaalvergelijking voor de oscillator blijkt dat er in de limietcyclus twee grensgevallen zijn: vrijwel harmonische sinusvormige trillingen $0<\mu \ll 1$, en relaxatietrillingen $\mu \gg 1$. Deze laatste, sterk niet-lineaire trillingsvorm dankt haar naam aan Van der Pol: dit vanwege de overeenkomst van deze oscillatievorm, gedurende een halve periode, met die van een (ont)ladende RC-schakeling, met relaxatietijd $\tau =RC$.
De geëxciteerde Van der Pol-oscillator wordt veel gebruikt in de niet-lineaire dynamica, als een voorbeeld van een natuurkundig systeem dat aanleiding kan geven tot chaotisch gedrag, alsmede tot andere niet-lineaire trillingsverschijnselen zoals bifurcaties en fasevergrendeling, phase locking.

## Betekenis

De Van der Pol-vergelijking is voortgevloeid uit onderzoek dat werd verricht aan het NatLab, en is een voorbeeld van het feit dat toegepast onderzoek — indien de onderzoeker voldoende vrijheid wordt gelaten — ook tot resultaten kan leiden die in meer algemene zin bruikbaar zijn. Zo heeft deze niet-lineaire differentiaalvergelijking vele nieuwe toepassingsgebieden gevonden buiten de elektrodynamica, zoals in de plasmafysica. Deze eponieme vergelijking is een begrip in de theorie over natuurkundige trillingsverschijnselen en in de chaostheorie.

## Gedrag van de niet-geëxciteerde vergelijking

* $\mu =0$: er is geen demping, en de vergelijking wordt die van de lineaire harmonische oscillator:
       ${\mathrm {d} ^{2}x \over \mathrm {d} t^{2}}+x=0$
    de afwezigheid van demping in deze lineaire vergelijking resulteert in behoud van energie.
* $\mu >0$: het systeem is gedempt voor grote uitwijkingen $|x|$, maar instabiel voor kleine $|x|$. Voor grote waarden van de tijd $t$ nadert het systeem naar een relaxatietrilling, in welke limietcyclus er weer energiebehoud is.
Voor kleine demping $\mu \ll 1$ zijn de oscillaties bijna sinusvormig. Voor grotere waarden van $\mu$, rond 1, worden de trillingen laagfrequenter en hoekiger: minder sinusvormig, maar met een steiler pulsfront gevolgd door een gedempte daaltijd. Bij grote waarden $\mu \gg 1$ krijgt men de zo kenmerkende, afgeronde blokgolfvorm van deze relaxatietrillingen.

## Gedrag bij excitatie

De Van der Pol-vergelijking kan bij sinusvormige forcering, $A\neq 0$ en $\omega \neq 0$, aanleiding geven tot chaotisch gedrag in de oplossing. Dit werd reeds door Van der Pol en zijn collega Van der Mark opgemerkt, en gerapporteerd in hun artikel in Nature in 1927:
   … often an irregular noise is heard in the telephone receivers before the frequency jumps, however this is a subsidiary phenomenon … (… vaak wordt voordat de frequentie verspringt een onregelmatig geluid in de telefoonhoorns gehoord, maar dit is echter een bijkomend fenomeen …)
Van der Pol en Van der Mark gebruikten telefoonhoorns om de signalen van hun triode-opstelling hoorbaar te maken.
Het gedrag van de geforceerde Van der Pol-vergelijking wordt bepaald door de parameters $A$ amplitude, $\omega$ hoeksnelheid en $\mu$, de niet-lineaire dempingsparameter, alsmede door de beginvoorwaarden. Voor verschillende waarden van de parameters wordt verschillend gedrag gevonden:
* Voor kleine amplitude, $A\ll 1$, als de excitatiefrequentie $\omega$ ongelijk is aan en geen is veelvoud van de trillingsfrequentie van de oscillator in vrije trilling — zodat geen resonantie optreedt — domineert de vrije trilling en geeft de excitatie slechts aanleiding tot kleine verstoringen op de limietcyclus.
* Voor grotere amplitudes is gecompliceerd gedrag waarneembaar, waaronder bifurcaties en chaos, in verschillende gebieden van de $A-\omega -\mu$ ruimte. Fasevergrendeling, phase locking treedt op, waarbij de oscillator trilt met een frequentie $m\omega /n$, met $m$ en $n$ positieve gehele getallen. Bij $m=1$ en $n=1$ spreekt men van subharmonische oscillaties, en voor $n=1$ en $m>1$ is sprake van superharmonischen, of boventonen. Maar, bij de Van der Pol-oscillator kunnen veel verschillende combinaties van $m$ en $n$ optreden in verschillende delen van de parameterruimte.
* Voor grote waarden van de amplitude, $A\gg 1$, trilt het systeem met de forceringsfrequentie $\omega$.

## Afleiding van de vergelijking voor een triode-oscillator

Hiernaast wordt een triode-schakeling weergegeven voor het realiseren van een Van der Pol-oscillator. De met wederzijdse inductie $M$ gekoppelde spoelen $L$, waarvan de primaire vanuit de sinusvormige bron $E_{s}$ met hoekfrequentie $\omega$ wordt geëxciteerd, resulteren in een elektrische stroom $i$ in de LCR-kring. De secundaire keten, met anodestroom $i_{a}$ van de triode, wordt gevoed vanaf de pluspool, terwijl de spanning op het triodestuurrooster bestaat uit een (negatieve) bias, waarop een (meegekoppelde) stuurroosteroscillatie $u_{g}$ vanuit de LCR-kring gesuperponeerd is.
Daarop kan met de spanningswet van Kirchhoff en Van der Pol's modellering van het triode-gedrag, de oscillator worden beschreven. Voor de LCR-kring geldt:
    $L{\frac {\mathrm {d} i}{\mathrm {d} t}}+R\,i+u_{g}-M{\frac {\mathrm {d} i_{a}}{\mathrm {d} t}}=P_{0}\sin(\omega t)$,
evenals de condensatorstroom:
    $i=C{\frac {\mathrm {d} u_{g}}{\mathrm {d} t}}$
Terwijl voor de triode geldt, volgens de betrekking van Barkhausen met steilheid $S$:
    ${\frac {\mathrm {d} i_{a}}{\mathrm {d} t}}=S{\frac {\mathrm {d} u_{g}}{\mathrm {d} t}}+{\frac {1}{R_{i}}}{\frac {\mathrm {d} u_{a}}{\mathrm {d} t}}$
en inwendige weerstand $R_{i}$ alsmede anodespanning $u_{a}$. Voor zijn triode-oscillatorschakeling vindt Van der Pol voor de anodestroom $i_{a}$ een niet-lineaire vorm:
    $i_{a}=S\,u_{g}\left(1-{\frac {u_{g}^{2}}{3K^{2}}}\right)$
Zodat, voor de gehele schakeling, de volgende differentiaalvergelijking resulteert voor de spanningsoscillatie $u_{g}$ op het stuurrooster:
    $LC{\frac {\mathrm {d} ^{2}u_{g}}{\mathrm {d} t^{2}}}-\left[(MS-RC)-MS{\frac {u_{g}^{2}}{K^{2}}}\right]{\frac {\mathrm {d} u_{g}}{\mathrm {d} t}}+u_{g}=P_{0}\sin(\omega t)$,
hetgeen — na normalisatie, en mits $MS>RC$ — resulteert in voornoemde Van der Pol-vergelijking.

## Relatie tot de Rayleigh-vergelijking

Lord Rayleigh introduceerde het volgende niet-lineaire model voor een systeem dat instabiel is bij lage snelheid $\mathrm {d} y/\mathrm {d} t$ en gedempt bij grotere snelheid:
    ${\frac {\mathrm {d} ^{2}y}{\mathrm {d} t^{2}}}-\mu \left[1-\left({\frac {\mathrm {d} y}{\mathrm {d} t}}\right)^{2}\right]{\frac {\mathrm {d} y}{\mathrm {d} t}}+y=0$
De Rayleigh-vergelijking is direct om te zetten in de Van der Pol-vergelijking via de transformatie
    ${\frac {\mathrm {d} y}{\mathrm {d} t}}={\frac {1}{\sqrt {3}}}\,x$
Lord Rayleigh onderzocht het optreden van een limietcyclus voor kleine en positieve waarden van de dempingsparameter $\mu \ll 1$, als zwak niet-lineaire oplossing van de naar hem genoemde differentiaalvergelijking. De Van der Pol-vergelijking wordt soms — vanwege deze relatie tot de eerdere Rayleigh-vergelijking — ook wel de Rayleigh–van der Pol-vergelijking genoemd.

### voetnoten

* ↑ (en) B vd Pol. A theory of the amplitude of free and forced triode vibrations, 1920. Een theorie voor de amplitude van vrije en gedwongen triode trillingen, in Radio Review 1, blz 701–710, 754–762
* ↑ (en) B. van der Pol & J. van der Mark (1928): "The heartbeat considered as a relaxation oscillation, and an electrical model of the heart". De hartslag beschouwd als relaxatie-oscillatie en een elektrisch model van het hart. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science Series 7, Suppl. 6, pp. 763–775, 1928.
* 1 2 (en) J.M.T. Thompson & H.B. Stewart (2002): Nonlinear dynamics and chaos. John Wiley and Sons, 2e editie, ISBN 0471876844, pp. 2 & 80–82.
* ↑ (en) B. van der Pol & J. van der Mark (1927): "Frequency demultiplication", Nature 120, DOI:10.1038/120363a0, blz 363–364.
* ↑ (en) H.J. Pain (2005): The physics of vibrations and waves. John Wiley and Sons, 6$^{de}$ druk, ISBN 047001296X , pp. 467–468.
* ↑ (en) U. Parlitz & W. Lauterborn (1987): "Period-doubling cascades and devil's staircases of the driven van der Pol oscillator", Physical Review A 36(3), DOI:10.1103/PhysRevA.36.1428, pp. 1428–1434.
* ↑ (en) R. Mettin, U. Parlitz & W. Lauterborn (1993) "Bifurcation structure of the driven van der Pol oscillator". International Journal of Bifurcation and Chaos 3(6), DOI:10.1142/S0218127493001203, pp. 1529–1555.
* 1 2 3 K Tomita. Periodically forced nonlinear oscillators, 1986. in Chaos, blz 213-214 ISBN 0719018110
* 1 2 (en) Lord Rayleigh (John William Strutt) (1896): The Theory of Sound. Cambridge:Cambridge University Press, 2e editie. Zie §68a.
* ↑ (en) P.S. Landa (2001): Regular and chaotic oscillations. Springer, ISBN 3540410015, pp. 89–92.

### websites

* Uitgebreide Engelse behandeling van de Oscillator-werking op Scholarpedia
* Instituut voor Sociale Geschiedenis – H.B.G. Casimir, Pol, Balthasar van der (1889–1959), in Biografisch Woordenboek van Nederland
* (en) notities betreffende Van der Pol-oscillator, op de website van Illinois State University
* (en) Interactieve demonstraties betreffende de Van der Pol-oscillator — Java applets op de website van het California Institute of Technology
* (en) R. Grimshaw. Nonlinear ordinary differential equations, 1993. blz 153–163
* (en) F. Verhulst. Nonlinear differential equations and dynamical systems, 1996. ISBN 3540609342, blz 15–16, 168–170 en 255–259