Oscilador de Van der Pol

Sistema dinámico oscilante con amortiguamiento no lineal

En el estudio de sistemas dinámicos , el oscilador de van der Pol (llamado así por el físico holandés Balthasar van der Pol ) es un sistema oscilante no conservativo con amortiguamiento no lineal . Evoluciona en el tiempo según la ecuación diferencial de segundo orden donde x es la coordenada de posición —que es una función del tiempo t— y μ es un parámetro escalar que indica la no linealidad y la fuerza del amortiguamiento. $${\displaystyle {d^{2}x \over dt^{2}}-\mu (1-x^{2}){dx \over dt}+x=0,}$$

Diagrama de fase del oscilador de van der Pol , con μ que varía de 0,1 a 3,0. Las líneas verdes son las líneas nulas de x .

El mismo gráfico de fase del oscilador, pero con la transformada de Liénard .

El oscilador de Van der Pol simulado con Brain Dynamics Toolbox ^[1]

Evolución del ciclo límite en el plano de fases . El ciclo límite comienza como un círculo y, al variar μ , se vuelve cada vez más agudo. Ejemplo de oscilador de relajación .

Historia

El oscilador de Van der Pol fue propuesto originalmente por el ingeniero eléctrico y físico holandés Balthasar van der Pol mientras trabajaba en Philips . ^[2] Van der Pol encontró oscilaciones estables, ^[3] que posteriormente llamó oscilaciones de relajación ^[4] y que ahora se conocen como un tipo de ciclo límite , en circuitos eléctricos que emplean tubos de vacío . Cuando estos circuitos se activan cerca del ciclo límite , se arrastran , es decir, la señal de activación arrastra la corriente junto con ella. Van der Pol y su colega, van der Mark, informaron en la edición de septiembre de 1927 de Nature que a ciertas frecuencias de activación se escuchaba un ruido irregular , ^[5] que luego se descubrió que era el resultado del caos determinista . ^[6]

La ecuación de Van der Pol tiene una larga historia de uso tanto en las ciencias físicas como en las biológicas . Por ejemplo, en biología, Fitzhugh ^[7] y Nagumo ^[8] extendieron la ecuación en un campo plano como modelo para los potenciales de acción de las neuronas . La ecuación también se ha utilizado en sismología para modelar las dos placas en una falla geológica [ ^9] y en estudios de fonación para modelar los osciladores de las cuerdas vocales derecha e izquierda ^[10] .

Forma bidimensional

El teorema de Liénard se puede utilizar para demostrar que el sistema tiene un ciclo límite. Aplicando la transformación de Liénard , donde el punto indica la derivada temporal, el oscilador de Van der Pol se puede escribir en su forma bidimensional: ^[11] ${\displaystyle y=x-x^{3}/3-{\dot {x}}/\mu }$

${\displaystyle {\dot {x}}=\mu \left(x-{\tfrac {1}{3}}x^{3}-y\right)}$

${\displaystyle {\dot {y}}={\frac {1}{\mu }}x}$.

Otra forma comúnmente utilizada basada en la transformación conduce a: ${\displaystyle y={\dot {x}}}$

${\displaystyle {\dot {x}}=y}$

${\displaystyle {\dot {y}}=\mu (1-x^{2})y-x}$.

Resultados para el oscilador no forzado

Oscilación de relajación en el oscilador de Van der Pol sin forzamiento externo. El parámetro de amortiguamiento no lineal es igual a μ = 5 .

^[12]

• Cuando μ = 0 , es decir, no hay función de amortiguamiento, la ecuación se convierte en Esta es una forma del oscilador armónico simple y siempre hay conservación de energía . $${\displaystyle {\frac {d^{2}x}{dt^{2}}}+x=0.}$$

• Cuando μ > 0 , todas las condiciones iniciales convergen a un ciclo límite globalmente único. Cerca del origen, el sistema es inestable y, lejos del origen, el sistema está amortiguado. ${\displaystyle x={\tfrac {dx}{dt}}=0,}$
• El oscilador de Van der Pol no tiene una solución analítica exacta. ^[13] Sin embargo, sí existe una solución para el ciclo límite si f ( x ) en la ecuación de Lienard es una función constante por partes.
• El período en μ pequeño tiene expansión serial . Véase el método de Poincaré–Lindstedt para una derivación hasta el orden 2. Véase el capítulo 10 de ^[14] para una derivación hasta el orden 3, y ^[15] para una derivación numérica hasta el orden 164. $${\displaystyle T={\frac {2\pi }{1-\mu ^{2}/16+17\mu ^{4}/3072+O(\mu ^{6})}}.}$$
• Para valores grandes de μ , el comportamiento del oscilador tiene un ciclo de acumulación lenta y liberación rápida (un ciclo de acumulación de tensión y liberación de tensión, por lo tanto, una oscilación de relajación ). Esto se ve más fácilmente en la forma En esta forma, el oscilador completa un ciclo de la siguiente manera: $${\displaystyle {\dot {x}}=\mu \left(x-{\frac {1}{3}}x^{3}-y\right)\!,\quad {\dot {y}}={\frac {1}{\mu }}x.}$$
  • Ascendiendo lentamente la rama derecha de la curva cúbica desde (2, –2/3) a (1, 2/3) . ${\displaystyle y=x-{\tfrac {x^{3}}{3}},}$
  • Moviéndose rápidamente hacia la rama izquierda de la curva cúbica, de (1, 2/3) a (–2, 2/3) .
  • Repita los dos pasos en la rama izquierda.
• El término principal en el período del ciclo se debe al lento ascenso y descenso, que se puede calcular como: Los órdenes superiores del período del ciclo son donde α ≈ 2.338 es la raíz más pequeña de Ai(– α ) = 0 , donde Ai es la función de Airy . (Sección 9.7 ^[16] ) ( ^[17] contiene una derivación, pero tiene un error de imprenta de 3 α a 2 α .) Esto fue derivado por Anatoly Dorodnitsyn . ^{[18] [19]} $${\displaystyle T=2\int dt=2\int \mu {\frac {dy}{x}}=2\mu \int _{2}^{1}{\frac {dy}{dx}}{\frac {dx}{x}}=(3-2\ln 2)\mu }$$ $${\displaystyle T=(3-2\ln 2)\mu +3\alpha \mu ^{-1/3}-{\frac {22}{9\mu }}\ln \mu +{\frac {0.0087\cdots }{\mu }}+O(\mu ^{-4/3}).}$$
• La amplitud del ciclo es ^[20] $${\displaystyle 2+{\frac {\alpha }{3}}\mu ^{-4/3}-{\frac {16}{27}}\mu ^{-2}\ln \mu +O(\mu ^{-2})}$$

Bifurcación de Hopf

A medida que μ pasa de menos de cero a más de cero, el sumidero espiral en el origen se convierte en una fuente espiral y aparece "de la nada" un ciclo límite con un radio de dos. Esto se debe a que la transición no es genérica: cuando ε = 0 , tanto la ecuación diferencial se vuelve lineal como el origen se convierte en un nodo circular.

Sabiendo que en una bifurcación de Hopf , el ciclo límite debe tener tamaño , podemos intentar convertir esto en una bifurcación de Hopf usando el cambio de variables que da Esto de hecho es una bifurcación de Hopf. ^[21] ${\displaystyle \propto \varepsilon ^{1/2},}$ ${\displaystyle u=\varepsilon ^{1/2}x,}$ $${\displaystyle {\ddot {u}}+u+u^{2}{\dot {u}}-\varepsilon {\dot {u}}=0}$$

Hamiltoniano para oscilador de Van der Pol

Las condiciones iniciales elegidas aleatoriamente son atraídas hacia una órbita estable.

También se puede escribir un formalismo hamiltoniano independiente del tiempo para el oscilador de Van der Pol ampliándolo a un sistema dinámico autónomo de cuatro dimensiones utilizando una ecuación diferencial no lineal auxiliar de segundo orden de la siguiente manera:

${\displaystyle {\ddot {x}}-\mu (1-x^{2}){\dot {x}}+x=0,}$

${\displaystyle {\ddot {y}}+\mu (1-x^{2}){\dot {y}}+y=0.}$

Obsérvese que la dinámica del oscilador Van der Pol original no se ve afectada debido al acoplamiento unidireccional entre las evoluciones temporales de las variables x e y . Se puede demostrar que un hamiltoniano H para este sistema de ecuaciones es ^[22]

${\displaystyle H(x,y,p_{x},p_{y})=p_{x}p_{y}+xy-\mu (1-x^{2})yp_{y},}$

donde y son los momentos conjugados correspondientes a x e y , respectivamente. Esto puede, en principio, conducir a la cuantificación del oscilador de Van der Pol. Un hamiltoniano de este tipo también conecta ^[23] la fase geométrica del sistema de ciclo límite que tiene parámetros dependientes del tiempo con el ángulo de Hannay del sistema hamiltoniano correspondiente. ${\displaystyle p_{x}={\dot {y}}+\mu (1-x^{2})y}$ ${\displaystyle p_{y}={\dot {x}}}$

Oscilador cuántico

El oscilador cuántico de van der Pol, que es la versión mecánico cuántica del oscilador clásico de van der Pol, se ha propuesto utilizando una ecuación de Lindblad para estudiar su dinámica cuántica y sincronización cuántica. ^[24] Nótese que el enfoque hamiltoniano anterior con una ecuación auxiliar de segundo orden produce trayectorias ilimitadas en el espacio de fases y, por lo tanto, no se puede utilizar para cuantificar el oscilador de van der Pol. En el límite de no linealidad débil (es decir, μ→ 0), el oscilador de van der Pol se reduce a la ecuación de Stuart-Landau . De hecho, la ecuación de Stuart-Landau describe una clase completa de osciladores de ciclo límite en el límite débilmente no lineal. La forma de la ecuación clásica de Stuart-Landau es mucho más simple y, tal vez no sea sorprendente, se puede cuantificar mediante una ecuación de Lindblad que también es más simple que la ecuación de Lindblad para el oscilador de van der Pol. El modelo cuántico de Stuart-Landau ha desempeñado un papel importante en el estudio de la sincronización cuántica ^{[25] [26]} (donde a menudo se lo ha denominado oscilador de van der Pol, aunque no se lo puede asociar únicamente con el oscilador de van der Pol). La relación entre el modelo clásico de Stuart-Landau ( μ→ 0) y osciladores de ciclo límite más generales ( μ arbitrario ) también se ha demostrado numéricamente en los modelos cuánticos correspondientes. ^[24]

Oscilador de Van der Pol forzado

Comportamiento caótico en el oscilador de Van der Pol con forzamiento sinusoidal. El parámetro de amortiguamiento no lineal es igual a μ = 8,53 , mientras que el forzamiento tiene amplitud A = 1,2 y frecuencia angular ω = 2π/10 .

El oscilador de Van der Pol forzado, o impulsado, toma la función "original" y agrega una función impulsora A sin( ωt ) para dar una ecuación diferencial de la forma:

${\displaystyle {d^{2}x \over dt^{2}}-\mu (1-x^{2}){dx \over dt}+x-A\sin(\omega t)=0,}$

donde A es la amplitud , o desplazamiento , de la función de onda y ω es su velocidad angular .

Cultura popular

Circuito eléctrico que involucra un triodo , lo que resulta en un oscilador de Van der Pol forzado. ^[27] El circuito contiene: un triodo, una resistencia R , un capacitor C , un conjunto de inductores acoplados con autoinductancia L e inductancia mutua M. En el circuito RLC en serie hay una corriente i , y hacia el ánodo ("placa") del triodo una corriente i _a , mientras que hay un voltaje u _g en la rejilla de control del triodo. El oscilador de Van der Pol es forzado por una fuente de voltaje de CA E _s .

El autor James Gleick describió un oscilador de Van der Pol de tubo de vacío en su libro de 1987 Chaos: Making a New Science . ^[28] Según un artículo del New York Times ^{, [29]} Gleick recibió un oscilador electrónico moderno de Van der Pol de un lector en 1988.

Véase también

• Mary Cartwright , matemática británica, una de las primeras en estudiar la teoría del caos determinista, particularmente aplicada a este oscilador. ^[30]

Referencias

1. Heitmann, S., Breakspear, M (2017–2022) Caja de herramientas de dinámica cerebral. bdtoolbox.org doi.org/10.5281/zenodo.5625923
2. Cartwright, ML (1960). "Balthazar Van Der Pol". Revista de la Sociedad Matemática de Londres . s1-35 (3). Wiley: 367–376. doi :10.1112/jlms/s1-35.3.367. ISSN  0024-6107.
3. B. van der Pol: "Una teoría de la amplitud de las vibraciones de triodo libre y forzado", Radio Review (más tarde Wireless World) 1 701–710 (1920)
4. van der Pol, Balth. (1926). "Sobre las oscilaciones de relajación"". Revista y diario filosófico de Londres, Edimburgo y Dublín . 2 (11). Informa UK Limited: 978–992. doi :10.1080/14786442608564127. ISSN  1941-5982.
5. VAN DER POL, BALTH; VAN DER MARK, J. (1927). "Desmultiplicación de frecuencias". Nature . 120 (3019). Springer Science and Business Media LLC: 363–364. Código Bibliográfico :1927Natur.120..363V. doi :10.1038/120363a0. ISSN  0028-0836. S2CID  186244992.
6. Kanamaru, T., "Oscilador de Van der Pol", Scholarpedia , 2 (1), 2202, (2007).
7. FitzHugh, Richard (1961). "Impulsos y estados fisiológicos en modelos teóricos de la membrana nerviosa". Revista biofísica . 1 (6). Elsevier BV: 445–466. Código Bibliográfico :1961BpJ.....1..445F. doi :10.1016/s0006-3495(61)86902-6. ISSN  0006-3495. PMC 1366333 . PMID  19431309. 
8. Nagumo, J.; Arimoto, S.; Yoshizawa, S. (1962). "Una línea de transmisión de pulsos activos que simula el axón nervioso". Actas del IRE . 50 (10). Instituto de Ingenieros Eléctricos y Electrónicos (IEEE): 2061–2070. doi :10.1109/jrproc.1962.288235. ISSN  0096-8390. S2CID  51648050.
9. Cartwright, Julyan HE ; Eguíluz, Víctor M.; Hernández-García, Emilio; Piro, Oreste (1999). "Dinámica de medios elásticos excitables". Revista Internacional de Bifurcación y Caos . 09 (11): 2197–2202. arXiv : chao-dyn/9905035 . Código Bibliográfico :1999IJBC....9.2197C. doi :10.1142/s0218127499001620. ISSN  0218-1274. S2CID  9120223.
10. Lucero, Jorge C.; Schoentgen, Jean (2013). Modelado de asimetrías de cuerdas vocales con osciladores de van der Pol acoplados. Actas de reuniones sobre acústica. Vol. 19. p. 060165. doi :10.1121/1.4798467. ISSN  1939-800X.
11. Kaplan, D. y Glass, L., Comprensión de la dinámica no lineal , Springer, 240–244, (1995).
12. Grimshaw, R., Ecuaciones diferenciales ordinarias no lineales , CRC Press , 153–163, (1993), ISBN 0-8493-8607-1 . 
13. Panayotounakos, DE; Panayotounakou, ND; Vakakis, AF (1 de septiembre de 2003). "Sobre la falta de soluciones analíticas del oscilador de Van der Pol". ZAMM . 83 (9). Wiley: 611–615. Bibcode :2003ZaMM...83..611P. doi :10.1002/zamm.200310040. ISSN  0044-2267. S2CID  120504403.
14. Verhulst, Ferdinand (1996). Ecuaciones diferenciales no lineales y sistemas dinámicos. Universitext. Berlín, Heidelberg: Springer Berlin Heidelberg. doi :10.1007/978-3-642-61453-8. ISBN 978-3-540-60934-6.
15. Andersen, CM; Geer, James F. (junio de 1982). "Expansiones de series de potencias para la frecuencia y el período del ciclo límite de la ecuación de Van Der Pol". Revista SIAM de Matemáticas Aplicadas . 42 (3): 678–693. doi :10.1137/0142047. ISSN  0036-1399.
16. Bender, Carl M. (1999). Métodos matemáticos avanzados para científicos e ingenieros I: métodos asintóticos y teoría de perturbaciones. Steven A. Orszag. Nueva York, NY. ISBN 978-1-4757-3069-2.OCLC 851704808  .{{cite book}}: Mantenimiento de CS1: falta la ubicación del editor ( enlace )
17. Grimshaw, R. (1993). Ecuaciones diferenciales ordinarias no lineales. Boca Raton, Fla.: CRC Press. pp. 161–163. ISBN 0-8493-8607-1.OCLC 28275539  .
18. A, Dorodnicyn A. (1947). "Solución asintótica de la ecuación de van der Pol". Prikl. Mat. Mech. (en ruso). 11 : 313–328.
19. Dorodnicyn, AA Solución asintótica de la ecuación de Van Der Pol. Biblioteca de investigación del NIST. Oficina Nacional de Normas.
20. Zonneveld, JA (1966). "Soluciones periódicas de la ecuación de Van der Pol". Indagationes Mathematicae (Actas) . 69 : 620–622. doi : 10.1016/s1385-7258(69)50068-x . ISSN  1385-7258.
21. Strogatz, Steven (2019). "Ejemplo 8.4.1". Dinámica no lineal y caos: con aplicaciones a la física, la biología, la química y la ingeniería (2.ª ed.). Boca Raton. ISBN 978-0-367-09206-1.OCLC 1112373147  .{{cite book}}: Mantenimiento de CS1: falta la ubicación del editor ( enlace )
22. Shah, Tirth; Chattopadhyay, Rohitashwa; Vaidya, Kedar; Chakraborty, Sagar (2015). "Teoría de perturbación conservativa para sistemas no conservativos". Physical Review E . 92 (6): 062927. arXiv : 1512.06758 . Bibcode :2015PhRvE..92f2927S. doi :10.1103/physreve.92.062927. PMID  26764794. S2CID  14930486.
23. Chattopadhyay, Rohitashwa; Shah, Tirth; Chakraborty, Sagar (2018). "Encontrar el ángulo de Hannay en sistemas oscilatorios disipativos mediante la teoría de perturbación conservativa". Physical Review E . 97 (6): 062209. arXiv : 1610.05218 . Bibcode :2018PhRvE..97f2209C. doi :10.1103/PhysRevE.97.062209. PMID  30011548. S2CID  51635019.
24. Chia, A.; Kwek, LC; Noh, C. (16 de octubre de 2020). "Oscilaciones de relajación y arrastre de frecuencia en mecánica cuántica". Physical Review E . 102 (4): 042213. arXiv : 1711.07376 . Bibcode :2020PhRvE.102d2213C. doi :10.1103/physreve.102.042213. ISSN  2470-0045. PMID  33212685. S2CID  224801468.
25. Walter, Stefan; Nunnenkamp, ​​Andreas; Bruder, Christoph (6 de marzo de 2014). "Sincronización cuántica de un oscilador autosostenido impulsado". Physical Review Letters . 112 (9): 094102. arXiv : 1307.7044 . Código Bibliográfico :2014PhRvL.112i4102W. doi :10.1103/physrevlett.112.094102. ISSN  0031-9007. PMID  24655255. S2CID  7950471.
26. Lee, Tony E.; Sadeghpour, HR (4 de diciembre de 2013). "Sincronización cuántica de osciladores de van der Pol cuánticos con iones atrapados". Physical Review Letters . 111 (23): 234101. arXiv : 1306.6359 . Código Bibliográfico :2013PhRvL.111w4101L. doi :10.1103/physrevlett.111.234101. ISSN  0031-9007. PMID  24476274. S2CID  33622111.
27. K. Tomita (1986): "Osciladores no lineales forzados periódicamente". En: Chaos , Ed. Arun V. Holden. Manchester University Press, ISBN 0719018110 , págs. 213-214. 
28. Gleick, James (1987). Caos: la creación de una nueva ciencia . Nueva York: Penguin Books. pp. 41–43. ISBN 0-14-009250-1.
29. Colman, David (11 de julio de 2011). "No hay silencio sin ruido". New York Times . Consultado el 11 de julio de 2011 .
30. Cartwright, ML ; Littlewood, JE (1945). "Sobre ecuaciones diferenciales no lineales de segundo orden: I. la ecuación y¨ − k (1-y 2 )y˙ + y = b λk cos(λl + α), k Grande". Revista de la Sociedad Matemática de Londres . s1-20 (3). Wiley: 180–189. doi :10.1112/jlms/s1-20.3.180. ISSN  0024-6107.

Enlaces externos

• "Ecuación de Van der Pol", Enciclopedia de Matemáticas , EMS Press , 2001 [1994]
• Oscilador de Van der Pol en Scholarpedia
• Demostraciones interactivas del oscilador Van Der Pol Archivado el 14 de febrero de 2017 en Wayback Machine