Oscilador de Van der Pol

Sistema dinámico oscilante con amortiguación 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 conservador con amortiguación no lineal . Evoluciona en el tiempo según la ecuación diferencial de segundo orden.

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

xcoordenadafuncióntμescalar que indica la no linealidad

Gráfico de fase del oscilador de van der Pol , con μ variando de 0,1 a 3,0. Las líneas verdes son las x - clinclinas nulas .

El mismo gráfico de fase del oscilador, pero con 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. Un 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 conducen cerca del ciclo límite , quedan arrastrados , es decir, la señal de activación arrastra la corriente consigo. Van der Pol y su colega, van der Mark, informaron en la edición de septiembre de 1927 de Nature que en ciertas frecuencias del motor se escuchaba un ruido irregular , ^[5] que luego se descubrió que era el resultado de un 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] ampliaron 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 del tiempo, 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 forzado externo. El parámetro de amortiguación no lineal es igual a μ = 5 .

^[12]

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

• 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, tal solución existe 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 en serie
  $${\displaystyle T={\frac {2\pi }{1-\mu ^{2}/16+17\mu ^{4}/3072+O(\mu ^{6})}}.}$$
  Consulte el método de Poincaré-Lindstedt para una derivación de orden 2. Consulte el capítulo 10 de ^[14] para una derivación de hasta orden 3, y ^[15] para una derivación numérica de hasta orden 164.
• Para μ grandes , 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
  $${\displaystyle {\dot {x}}=\mu \left(x-{\frac {1}{3}}x^{3}-y\right)\!,\quad {\dot {y}}={\frac {1}{\mu }}x.}$$
  De esta forma, el oscilador completa un ciclo de la siguiente manera:
  • Ascendiendo lentamente por la rama derecha de la curva cúbica desde (2, –2/3) hasta (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:
  $${\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 }$$
  Los órdenes superiores del período del ciclo son
  $${\displaystyle T=(3-2\ln 2)\mu +3\alpha \mu ^{-1/3}-{\frac {23}{\mu ^{-1}}}\ln \mu +O(\ln \mu ^{-1}).}$$
  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 tipográfico de 3 α a 2 α . )
• La amplitud del ciclo es ^[18]
  $${\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 μ se mueve de menos de cero a más de cero, el sumidero espiral en el origen se convierte en una fuente espiral y aparece un ciclo límite "de la nada" con radio 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 un tamaño, podemos intentar convertir esto en una bifurcación de Hopf utilizando el cambio de variables que da ${\displaystyle \propto \varepsilon ^{1/2},}$ ${\displaystyle u=\varepsilon ^{1/2}x,}$

$${\displaystyle {\ddot {u}}+u+u^{2}{\dot {u}}-\varepsilon {\dot {u}}=0}$$

^[19]

Hamiltoniano para oscilador de Van der Pol

Las condiciones iniciales elegidas al azar 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 aumentá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.}$

Tenga en cuenta que la dinámica del oscilador de Van der Pol original no se ve afectada debido al acoplamiento unidireccional entre las evoluciones temporales de las variables xey . Se puede demostrar que una H hamiltoniana para este sistema de ecuaciones es ^[20]

${\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. En principio, esto puede conducir a la cuantificación del oscilador de Van der Pol. Un hamiltoniano de este tipo también conecta ^[21] 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ánica 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. ^[22] Tenga en cuenta que el enfoque hamiltoniano anterior con una ecuación auxiliar de segundo orden produce trayectorias de espacio de fase ilimitadas y, por lo tanto, no se puede utilizar para cuantificar el oscilador de van der Pol. En el límite de la 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 toda una clase 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, puede cuantificarse 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 ^{[23] [24]} (donde a menudo se le ha llamado oscilador de van der Pol, aunque no puede asociarse ú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. ^[22]

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 , dando como resultado un oscilador de Van der Pol forzado. ^[25] El circuito contiene: un triodo, una resistencia R , un condensador 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 triodo ("placa") 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 CA E _s .

El autor James Gleick describió un oscilador de Van der Pol de tubo de vacío en su libro de 1987 Caos: Haciendo una nueva ciencia . ^[26] Según un artículo del New York Times ^{, [27]} Gleick recibió un oscilador electrónico moderno de Van der Pol de un lector en 1988.

Ver 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. ^[28]

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). "Baltasar Van Der Pol". Revista de la Sociedad Matemática de Londres . Wiley. T1-35 (3): 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 del triodo libre y forzada", Radio Review (más tarde Wireless World) 1 701–710 (1920)
4. van der Pol, Balth. (1926). "Sobre las "oscilaciones de relajación"". Revista filosófica y revista científica de Londres, Edimburgo y Dublín . Informa UK Limited. 2 (11): 978–992. doi :10.1080/14786442608564127. ISSN  1941-5982.
5. VAN DER POL, BALTH; VAN DER MARK, J. (1927). "Desmultiplicación de frecuencias". Naturaleza . Springer Science y Business Media LLC. 120 (3019): 363–364. Código Bib :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 membrana nerviosa". Revista Biofísica . Elsevier BV. 1 (6): 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 pulso activa que simula el axón nervioso". Actas del IRE . Instituto de Ingenieros Eléctricos y Electrónicos (IEEE). 50 (10): 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 Bib : 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ág. 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, Dakota del Norte; Vakakis, AF (1 de septiembre de 2003). "Sobre la falta de soluciones analíticas del oscilador de Van der Pol". ZAMM . Wiley. 83 (9): 611–615. Código Bib : 2003ZaMM...83..611P. doi :10.1002/zamm.200310040. ISSN  0044-2267. S2CID  120504403.
14. Verhulst, Fernando (1996). Ecuaciones diferenciales no lineales y sistemas dinámicos. Texto universitario. Berlín, Heidelberg: Springer Berlín 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 en series de potencias para la frecuencia y período del ciclo límite de la ecuación de Van Der Pol". Revista SIAM de Matemática Aplicada . 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 la perturbación. Steven A. Orszag. Nueva York, NY. ISBN 978-1-4757-3069-2. OCLC  851704808.{{cite book}}: Mantenimiento CS1: falta el editor de la ubicación ( enlace )
17. Grimshaw, R. (1993). Ecuaciones diferenciales ordinarias no lineales. Boca Ratón, Florida: CRC Press. págs. 161-163. ISBN 0-8493-8607-1. OCLC  28275539.
18. 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.
19. 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 Ratón. ISBN 978-0-367-09206-1. OCLC  1112373147.{{cite book}}: Mantenimiento CS1: falta el editor de la ubicación ( enlace )
20. Shah, Tirth; Chattopadhyay, Rohitashwa; Vaidya, Kedar; Chakraborty, Sagar (2015). "Teoría de la perturbación conservadora para sistemas no conservadores". Revisión física E. 92 (6): 062927. arXiv : 1512.06758 . Código bibliográfico : 2015PhRvE..92f2927S. doi : 10.1103/physreve.92.062927. PMID  26764794. S2CID  14930486.
21. Chattopadhyay, Rohitashwa; Shah, Tirth; Chakraborty, Sagar (2018). "Encontrar el ángulo de Hannay en sistemas oscilatorios disipativos mediante la teoría de la perturbación conservadora". Revisión física E. 97 (6): 062209. arXiv : 1610.05218 . Código Bib : 2018PhRvE..97f2209C. doi : 10.1103/PhysRevE.97.062209. PMID  30011548. S2CID  51635019.
22. Chia, A.; Kwek, LC; Noh, C. (16 de octubre de 2020). "Oscilaciones de relajación y arrastre de frecuencia en mecánica cuántica". Revisión física E. 102 (4): 042213. arXiv : 1711.07376 . Código bibliográfico : 2020PhRvE.102d2213C. doi :10.1103/physreve.102.042213. ISSN  2470-0045. PMID  33212685. S2CID  224801468.
23. Walter, Stefan; Nunnenkamp, ​​Andreas; Bruder, Christoph (6 de marzo de 2014). "Sincronización cuántica de un oscilador autosostenido impulsado". Cartas de revisión física . 112 (9): 094102. arXiv : 1307.7044 . Código Bib : 2014PhRvL.112i4102W. doi : 10.1103/physrevlett.112.094102. ISSN  0031-9007. PMID  24655255. S2CID  7950471.
24. Lee, Tony E.; Sadeghpour, recursos humanos (4 de diciembre de 2013). "Sincronización cuántica de osciladores cuánticos de van der Pol con iones atrapados". Cartas de revisión física . 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.
25. K. Tomita (1986): "Osciladores no lineales forzados periódicamente". En: Caos , Ed. Arun V. Holden. Prensa de la Universidad de Manchester, ISBN 0719018110 , págs. 
26. Gleick, James (1987). Caos: haciendo una nueva ciencia . Nueva York: Penguin Books. págs. 41–43. ISBN 0-14-009250-1.
27. Colman, David (11 de julio de 2011). "No hay silencio sin ruido". New York Times . Consultado el 11 de julio de 2011 .
28. 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 . Wiley. T1-20 (3): 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 la Wayback Machine.