范德珀尔振子

范德珀尔振子

概述

范德珀尔振子是非线性动力学理论中的一个重要发现。作为一个非保守非线性振荡器，其基本模型为一个二阶非线性方程： [math]\displaystyle{ \ddot{x} - \mu (1 - x^2)\dot{x} + x = 0, \quad \mu > 0 }[/math] 作为非线性动力学系统的一个典型案例，其动力学特性由位移相关的非线性阻尼系数决定。通过中间的非线性阻尼项 [math]\displaystyle{ -\mu (1 - x^2)\dot{x} }[/math]决定不同的动力学特性，该模型将线性理论、分岔理论、极限环、快慢动力学乃至混沌的边缘整合了起来，实现了由传统线性振荡理论向非线性动力学，乃至于混沌理论的飞跃（当方程右侧不为0时，模型变为受迫范德珀尔振子）。 在阻尼系数不同时，系统通过霍普夫分岔从平凡平衡点产生振荡。其能量调节机制完全取决于位移 x 的大小：当位移较小（|x| < 1）时，总阻尼项为负，系统注入能量（负阻尼发散）；当位移较大（|x| > 1）时，总阻尼项为正，系统耗散能量（正阻尼收敛）[1]。 这种机制最终使系统收敛于唯一极限环，从而实现近似正弦振荡或自激振荡。它是解释霍普夫分岔产生极限环现象的典型模型。

历史

19世纪末20世纪初，动力学理论与电子电气技术持续发展，如1880年Gérard-Lescuyer的串联电机[2]、1901年Duddell的音乐电弧[3]、1907年de Forest发明真空三极管[4]，都展示了“自激振荡”现象。面对自激振动现象，传统线性振动理论无法解释极限环及分岔现象，传统摄动分析方法面临级数发散问题，同时缺失解释不同自激振动现象的统一理论体系。 在此背景下，时任飞利浦公司电气工程师的荷兰物理学家巴尔塔萨·范德珀尔于1920年采用非线性动力学分析方法，引入非线性阻尼概念，在研究真空电子管振荡电路时提出了这一二阶非线性模型[5]。该模型最初用于描述早期无线电设备中三极管电路的自持振荡现象，是首个被深入研究的“自激振荡”非线性模型。 范德珀尔在研究该系统时开创性地提出了弛豫振荡（张弛振荡，Relaxation Oscillation）这一概念，并系统阐述了此现象；同时，也开创了“非线性振荡”的系统性研究方法。他与Van der Mark于1927年在《Nature》上发表了关于“不规则噪声”的报告[6]，后被公认为混沌现象最早的实验观察之一。此外，该发现还间接催生了后续混沌理论的诞生。

主要内容

数学模型

阻尼系数不同时的范德珀尔振子动力学相图，可见当μ非0时，随着振子不断振荡，最终都会收敛于一个极限环，且当μ显著大于时，系统出现了节律性的弛豫震荡的现象。

范德珀尔振子（Van der Pol oscillator）是研究极限环非线性系统的经典实例。作为一个耗散性非线性振荡系统，其数学模型为一个二阶非线性齐次方程：

[math]\displaystyle{ \ddot{x} - \mu (1 - x^2)\dot{x} + x = 0, \quad \mu > 0 }[/math]

其中 x 为位置坐标（作为时间 t 的函数），\mu 是标量参数，决定阻尼的非线性和强度。 也可写为二维状态空间方程组的形式：

[math]\displaystyle{ \begin{cases} \dot{x} = y, \\ \dot{y} = \mu(1 - x^2)y - x. \end{cases} }[/math]

其中： （1）x 是位移（或电压等物理量）。 （2）y = \dot{x} 是速度，\dot{y} = \ddot{x} 是加速度。 （3）\mu 是非线性阻尼参数（\mu \geq 0），控制着非线性效应的强弱。

物理意义与动力学行为

范德珀尔振子在非线性阻尼参数不同时，存在以下三种情况：

（1）当 [math]\displaystyle{ \mu = 0 }[/math] 时：方程退化为简谐振荡器 [math]\displaystyle{ \ddot{x} + x = 0 }[/math]，系统做等幅正弦振荡，没有能量耗散也没有能量注入。

（2）当 [math]\displaystyle{ 0 < \mu \ll 1 }[/math] 时：非线性较弱，系统表现出接近正弦波的振荡，但振幅会缓慢调整，最终稳定在一个极限环上。

（3）当 [math]\displaystyle{ \mu \gg 1 }[/math] 时：非线性很强，系统表现出弛豫振荡（Relaxation Oscillation）。波形不再是正弦波，而是呈现“缓慢充电、快速放电”的锯齿状或脉冲状特征，最终也将收敛于一个极限环。这在生物神经元的动作电位模型中非常常见。

阻尼系数不同的范德珀尔振子位移随时间变化图

阻尼系数不同的范德珀尔振子机械能随时间变化图

模型的推导过程

我们可通过如下步骤来探究[7]：考虑如图所示的RLC电路：它由电感L、电容器C、电阻R组成，我们假设电压源是电池（[math]\displaystyle{ E(t)=E }[/math]常数）。当电池开关关闭时，根据基尔霍夫电压定律：

[math]\displaystyle{ L I' + R I + \frac{1}{C} Q = E }[/math]

电路中开始流动电流 [math]\displaystyle{ I(t) }[/math]，其中Q作为电容器中的电荷，其一阶导数等于I。如果我们对等式两边求导，我们将得到式：

[math]\displaystyle{ L I + R I' + \frac{1}{C} I = 0 }[/math]

这是一个具有恒定系数的二阶线性常微分方程，它表示一个阻尼谐波振荡器。 然而，范德珀尔考虑的电路显示的是有源元件（真空管阵列-半导体），而不是无源电阻。通过这种方式，当电流低时，半导体就像在系统中泵送能量一样，当电流过高时，半导体会抑制系统的能量。能量注入和能量吸收之间的相互作用导致电压和电流的周期性振荡。半导体的作用由函数 [math]\displaystyle{ I^2 - \alpha }[/math] 所刻画，其中 \alpha 是电流的阈值水平，电流的方程变成下式：

[math]\displaystyle{ L I' + (I^2 - \alpha)I + \frac{1}{C} Q = E }[/math]

为了简单起见，我们假设没有外部电压源。通过取方程两侧的导数，我们确实得到了范德珀尔振子的方程：

[math]\displaystyle{ L I + 3 I' \left(I^2 - \frac{\alpha}{3}\right) + \frac{1}{C} I = 0 }[/math]

自引入以来，范德波尔方程一直被用作物理学、电子学、生物学、神经学、社会学和经济学中振荡过程的基本模型。该模型很快被推广为受迫范德珀尔振子系统。

延伸

受迫范德珀尔振子

当方程右边不为0的时候，即系统受到外部激励（External Excitation）时，则成为受迫范德珀尔振子（Forced Van der Pol oscillator）： [math]\displaystyle{ \ddot{x} - \mu (1 - x^2)\dot{x} + x = F \cos(\omega t) }[/math] 其中F是正弦驱动函数的振幅，\omega是其角频率。此系统经历一系列分岔（Bifurcation）过程，这个系统从稳定的极限环走向混沌。 （注：受迫VDP振子因引入时间维度而升维至三维及以上，打破了庞加莱-本迪克松定理的限制，从而触及混沌理论的门槛。日本学者上田睆亮于1961年在此系统中首次观察到了混沌吸引子[8]。）

霍普夫分叉

当\mu从小于零移动到大于零时，原点的螺旋汇变成了螺旋源，一个极限周期骤然出现，半径为2。这是因为过渡不是通用的：当 [math]\displaystyle{ \mu = 0 }[/math] 时，微分方程变成线性，原点变成中心（Center）（而非节点）。 在霍普夫分叉中，极限环应当具有正比于 [math]\displaystyle{ \varepsilon^{1/2} }[/math] 的振幅，现令 [math]\displaystyle{ \varepsilon = \mu }[/math]，我们可能会尝试使用变量的变化 [math]\displaystyle{ u = \mu^{1/2} x }[/math] 将范德珀尔振子转换为霍普夫分叉，得到：

[math]\displaystyle{ \ddot{u} + u + u^2 \dot{u} - \varepsilon \dot{u} = 0 }[/math]

此式即为一个霍普夫分叉的模型[9]。

量子物理领域——范德珀尔振子的哈密顿量

人们还可以为范德波尔振荡器编写一个与时间无关的哈密顿形式，方法是使用辅助二阶非线性微分方程将其增强为四维自主动力系统，如下所示：

[math]\displaystyle{ \ddot{x} - \mu(1 - x^2)\dot{x} + x = 0 }[/math] [math]\displaystyle{ \ddot{y} + \mu(1 - x^2)\dot{y} + y = 0 }[/math]

请注意，由于x和y变量的时间演变之间的单向耦合，原始范德波尔振荡器的动力学不受影响。可以通过这个方程组的哈密顿函数H证明[10]：

[math]\displaystyle{ H(x, y, p_x, p_y) = p_x p_y + xy - \mu (1 - x^2)y p_y }[/math]

其中 [math]\displaystyle{ p_x = \dot{y} + \mu(1 - x^2)y }[/math] 和 [math]\displaystyle{ p_y = \dot{x} }[/math] 分别是对应于x和y的共轭矩。原则上，这可能会导致范德波尔振荡器的量子化，这样的哈密顿量也将具有时间依赖参数的极限周期系统的几何相位与相应的哈密顿量系统的汉奈角连接起来[11]。

关于量子振荡器

量子范德波尔振荡器是经典范德波尔振荡器的量子力学版本，使用林德布拉德方程来研究其量子动力学和量子同步[12]。请注意，上述具有辅助二阶方程的哈密顿方法产生无界相空间轨迹，因此无法构造物理上允许的束缚态，不能直接用于量子化范德波尔振荡器。在弱非线性极限（即 [math]\displaystyle{ \mu \to 0 }[/math]）中，范德波尔振荡器简化为Stuart-Landau方程。Stuart-Landau方程实际上描述了弱非线性极限中的一整类极限周期振荡器。经典Stuart-Landau方程的形式要简单得多，可以用Lindblad方程量化。 量子Stuart-Landau模型在量子同步研究中发挥了重要作用[13][14]。经典Stuart-Landau模型（[math]\displaystyle{ \mu \to 0 }[/math]）与更一般的极限周期振荡器（任意\mu）之间的关系也在相应的量子模型中进行了数值证明。 然而，上述借助辅助变量构建哈密顿量的方法虽在数学上具有巧思，却止步于量子化的物理门槛——其产生的无界相空间轨迹无法构造物理上允许的束缚态。这恰恰揭示了经典非保守系统（如范德波尔振子）量子化的本质困难：系统因非线性阻尼与外界环境存在持续的能量交换，必须将其视为一个开放量子系统。 要跨越这一鸿沟，理论物理学家转向了Lindblad主方程形式，通过引入增益和耗散两个独立的退相干通道，成功构建了量子范德波尔振子的标准模型。近年来，随着离子阱、超导电路等量子模拟平台的飞速发展，该理论模型终于得以在实验室中精确“搭建”[12]。 2025年，中国科学技术大学林毅恒教授团队联合深圳大学及德国斯图加特大学等机构的研究人员，基于上述开放量子系统的理论框架，首次在离子阱实验中实现了量子范德波尔振子。相关成果于2025年10月10日发表在国际学术期刊《Science Advances》上，演示了量子极限环的形成、量子同步及其增强等非经典行为[15]。

应用与实例

电子电气领域

范德珀尔振子在电子电气领域有着广泛而根基性的应用。其核心的极限环和非线性阻尼特性直接对应于含负阻元件的振荡电路行为。例如，在使用真空管的早期电路中，当系统在极限环附近受到外部周期信号驱动时，振荡频率会被“拖入”驱动信号的频率，两者形成同步——这一现象被称为频率牵引或注入锁定。 范德珀尔方程最初即源于对含非线性负阻元件（真空管、隧道二极管等）的LC振荡电路的分析。方程中的非线性阻尼项精确描述了负阻元件的电流-电压特性——小信号时呈现负电阻（注入能量），大信号时呈现正电阻（耗散能量），这正是电子管放大电路能够产生稳定正弦波振荡的物理机制。范德珀尔模型已被证明是分析现代RC振荡器的通用振荡器模型。在CMOS纳米技术中，范德珀尔参数与多种电路实现之间的对应关系也已得到系统阐述。注入锁定的范德珀尔振荡器可以用简单的现代半导体电路实现。

生物医学领域

范德珀尔振子的极限环概念精准对应于心脏窦房结（SA node）和房室结（AV node）等天然起搏器的自发节律放电。范德珀尔和van der Mark早在1928年即于《Philosophical Magazine》发表了题为“The heartbeat considered as a relaxation oscillation, and an electrical model of the heart”的论文，建立了心脏的电子电路模型，将心跳视为弛豫振荡。 此后，改良的范德珀尔振荡器被广泛用于模拟心脏起搏器的动作电位时序。Gois等在2004年设计了修正的范德珀尔振荡器模型，通过在经典范德珀尔振荡器的相空间中增加一个鞍点和一个节点，并引入非对称阻尼项，使得起搏节点的放电频率可在不改变不应期长度的条件下独立调节。 在神经科学中，FitzHugh于1961年将范德珀尔方程与Hodgkin-Huxley模型结合，提出了FitzHugh-Nagumo模型（亦称Bonhoeffer-van der Pol模型）。该模型的相平面中存在一个稳定的极限环，对应于神经元的周期性脉冲发放。研究表明，FitzHugh-Nagumo方程可以在很好的近似程度上被约化为一个受外部电流驱动的广义范德珀尔振荡器。

催生混沌理论的诞生

范德珀尔振子及受迫范德珀尔振子呈现的分叉性质，为最初被人类观测到的混沌现象之一，叩开了混沌世界的大门。1945年英国数学家卡特赖特（M.Cartwright）和李特伍德（J.Littlewood）在对受迫范德珀尔振子的研究中，在理论层面首次证明受迫范德珀尔方程存在奇异解，预示了混沌的可能性。随后，日本学者上田睆亮（Ueda Yoshisuke）于1961年首次在模拟计算机上观察到混沌吸引子，随着后续研究的深入，混沌理论的面纱也逐渐被科学家们所揭晓。

参考文献

1. Strogatz, Steven (2019). "Example 7.3.1". Nonlinear dynamics and chaos: with applications to physics, biology, chemistry, and engineering (2nd ed.). Boca Raton. ISBN 978-0-367-09206-1. OCLC 1112373147 (https://search.worldcat.org/oclc/1112373147). 2.Gérard-Lescuyer, J.-M.-A. (1880). Sur un paradoxe électrodynamique. 3.Nature (1901). Musical Arcs. Nature, 63(1640), 542–543. 4.De Forest, L. (1906). The Audion: A New Receiver for Wireless Telegraphy. Transactions of the American Institute of Electrical Engineers. 5.VAN DER POL B. A theory of the amplitude of free and forced triode vibrations[J]. Radio Review, 1920, 1: 701-710. 6.VAN DER POL B, VAN DER MARK J. Frequency demultiplication[J]. Nature, 1927, 120(3019): 363-364. DOI:10.1038/120363a0. 7.https://www.phys.uconn.edu/~rozman/Courses/P2200_13F/downloads/vanderpol/vanderpol-oscillator-draft.pdf 8.C. Hayashi, H. Shibayama, and Y. Ueda. "Quasi-periodic Oscillations in Self-oscillatory Systems with External Force." (In Japanese), IECE Technical Report, Nonlinear Theory, 16 December 1961. 9. Strogatz, Steven (2019). "Example 8.4.1". Nonlinear dynamics and chaos: with applications to physics, biology, chemistry, and engineering (2nd ed.). Boca Raton. ISBN 978-0-367-09206-1. OCLC 1112373147 (https://search.worldcat.org/oclc/1112373147). 10. Shah, Tirth; Chattopadhyay, Rohitashwa; Vaidya, Kedar; Chakraborty, Sagar (2015). "Conservative perturbation theory for nonconservative systems" (https://journals.aps.org/pre/ab stract/10.1103/PhysRevE.92.062927). Physical Review E. 92 (6) 062927. arXiv:1512.06758 (ht tps://arxiv.org/abs/1512.06758). Bibcode:2015PhRvE..92f2927S (https://ui.adsabs.harvard.ed u/abs/2015PhRvE..92f2927S). doi:10.1103/physreve.92.062927 (https://doi.org/10.1103%2Fph ysreve.92.062927). PMID 26764794 (https://pubmed.ncbi.nlm.nih.gov/26764794). S2CID 14930486 (https://api.semanticscholar.org/CorpusID:14930486). 11. Chattopadhyay, Rohitashwa; Shah, Tirth; Chakraborty, Sagar(2018). "Finding the Hannay angle in dissipative oscillatory systems via conservative perturbation theory" (https://journals.a ps.org/pre/abstract/10.1103/PhysRevE.97.062209). Physical Review E. 97 (6)062209.arXiv:1610.05218(https://arxiv.org/abs/1610.05218). Bibcode:2018PhRvE..97f2209C (https://ui.adsabs.harvard.edu/abs/2018PhRvE..97f2209C). doi:10.1103/PhysRevE.97.062209 (https://d oi.org/10.1103%2FPhysRevE.97.062209). PMID 30011548 (https://pubmed.ncbi.nlm.nih.gov/3 0011548). S2CID 51635019 (https://api.semanticscholar.org/CorpusID:51635019). 12. Chia, A.; Kwek, L. C.; Noh, C. (2020-10-16). "Relaxation oscillations and frequency entrainment in quantum mechanics". Physical Review E. 102 (4) 042213. arXiv:1711.07376 (htt ps://arxiv.org/abs/1711.07376). Bibcode:2020PhRvE.102d2213C (https://ui.adsabs.harvard.ed u/abs/2020PhRvE.102d2213C). doi:10.1103/physreve.102.042213 (https://doi.org/10.1103%2F physreve.102.042213). ISSN 2470-0045 (https://search.worldcat.org/issn/2470-0045). PMID 33212685 (https://pubmed.ncbi.nlm.nih.gov/33212685). S2CID 224801468 (https://api.se manticscholar.org/CorpusID:224801468). 13. Walter, Stefan; Nunnenkamp, Andreas; Bruder, Christoph (2014-03-06). "Quantum Synchronization of a Driven Self-Sustained Oscillator". Physical Review Letters. 112 (9) 094102. arXiv:1307.7044 (https://arxiv.org/abs/1307.7044). Bibcode:2014PhRvL.112i4102W (h ttps://ui.adsabs.harvard.edu/abs/2014PhRvL.112i4102W). doi:10.1103/physrevlett.112.094102 (https://doi.org/10.1103%2Fphysrevlett.112.094102). ISSN 0031-9007 (https://search.worldcat. org/issn/0031-9007). PMID 24655255 (https://pubmed.ncbi.nlm.nih.gov/24655255). S2CID 7950471 (https://api.semanticscholar.org/CorpusID:7950471). 14. Lee, Tony E.; Sadeghpour, H. R. (2013-12-04). "Quantum Synchronization of Quantum van der Pol Oscillators with Trapped Ions". Physical Review Letters. 111 (23) 234101. arXiv:1306.6359 (https://arxiv.org/abs/1306.6359). Bibcode:2013PhRvL.111w4101L (https://ui.adsabs.harvard.e du/abs/2013PhRvL.111w4101L). doi:10.1103/physrevlett.111.234101 (https://doi.org/10.1103% 2Fphysrevlett.111.234101). ISSN 0031-9007 (https://search.worldcat.org/issn/0031-9007). PMID 24476274 (https://pubmed.ncbi.nlm.nih.gov/24476274). S2CID 33622111 (https://api.se manticscholar.org/CorpusID:33622111). 15.