范德珀尔振子
定义
范德珀尔振子是非线性动力学理论中的一个重要发现，作为一个非保守非线性振荡器，其基本模型为一个二阶非线性方程： [math]\displaystyle{ \ddot{x} - \mu (1 - x^2)\dot{x} + x = 0, \quad \mu \gt 0 }[/math] 作为非线性动力学系统的一个典型案例，其动力学特性由阻尼系数μ的大小决定，在阻尼系数不同时，系统通过霍普夫分岔从平凡平衡点通过振动，以位移依赖的非线性阻尼作为能量调节手段，在平衡点附近负阻尼发散，远离平衡点时正阻尼收敛，逐渐收敛于唯一极限环，从而实现近似正弦振荡或自激振荡。是解释霍普夫分岔产生的极限环现象的典型模型。
历史
19世纪末20世纪初，动力学理论与电子电气技术持续发展，如1880年Gérard-Lescuyer的串联电机、1901年Duddell的音乐电弧、1907年de Forest发明真空三极管，都展示了“自激振荡”现象，为后续研究提供了实验基础，而面对自激振动现象，传统线性振动理论无法解释极限环及分岔现象，传统摄动分析方法面临级数发散问题，同时缺失解释不同自激振动现象的统一理论体系。在此背景下，时任飞利浦公司的电气工程师的荷兰物理学家巴尔塔萨·范德珀尔与1920年采用非线性动力学分析方法，引入非线性阻尼概念，在研究真空电子管振荡电路时提出了这一二阶非线性模型[8]，该模型最初用于描述早期无线电设备中三极管电路的自持振荡现象，其电路是一个由电感器、电容器和非线性电阻器组成的非线性电路，是首个被深入研究的“自激振荡”非线性模型，也是极限环的典型代表之一。范德珀尔在研究该系统时开创性地提出了驰豫振荡（张弛振荡，relaxation oscillation）这一概念，并系统阐述了此现象；同时，也开创了“非线性震荡”的系统性研究方法，范德珀尔在研究此系统中引入了平均法，他与Van der Mark 于1927年在《Nature》上发表了关于“不规则噪声”的报告，后被公认为混沌现象最早的实验观察之一；此外，该发现还间接催生了后续混沌理论的诞生。
主要内容
数学模型
范德珀尔振子（Van der Pol oscillator）是研究极限环非线性系统的经典实例，作为一个耗散性非线性震荡系统，其数学模型为一个二阶非线性齐次方程，：
[math]\displaystyle{ \ddot{x} - \mu (1 - x^2)\dot{x} + x = 0, \quad \mu \gt 0 }[/math]
其中x为位置坐标，作为时间t的函数，μ是标量参数，决定阻尼的非线性和强度。
也可写为状态空间方程组的形式：
\begin{cases} \ddot{x} - \mu(1 - x^2)\dot{x} + x = 0, \\ \ddot{y} + \mu(1 - x^2)\dot{y} + y = 0. \end{cases}
其中： （1）x是位移（或电压等物理量）。 （2）x ̇是速度，x ̈是加速度 （3）μ是非线性阻尼参数（μ≥0），控制着非线性效应的强弱，振子动力学特性由Mu单独决定。
我们可通过如下步骤来探究：考虑如图所示的RLC电路：它由电感L、电容器C、电阻R组成，我们假设电压源是电池（E（t）=E常数）。当电池开关关闭时，根据基尔霍夫电压定律：L I' + R I + \frac{1}{C} Q = E，电路中开始流动电流I(t)，其中Q作为电容器中的电荷，其一阶导数等于I,如果我们对等式两边求导，我们将得到式：LI"+RI'+\frac{1}{C}I=0,这是一个具有恒定系数的二阶线性常微分方程，它表示一个阻尼谐波振荡器。 然而，范德珀尔考虑的电路显示的是有源元件（真空管阵列-半导体），而不是无源电阻。通过这种方式，当电流低时，半导体就像在系统中泵送能量一样，当电流过高时，半导体会抑制系统的能量（不像电阻器只是消散能量）。能量注入和能量吸收之间的相互作用导致电压和电流的周期性振荡。半导体的作用由函数I^2-α所刻画，其中α是电流的阈值水平，电流的方程变成下式：LI'+(I^2-α)I+\frac{1}{C}Q=E,为了简单起见，我们假设没有外部电压源。通过取方程两侧的导数，我们确实得到了范德珀尔振子的方程：LI"+3I'(I^2-\frac{α}{3})+\frac{1}{C}I=0 。自引入以来，范德波尔方程一直被用作物理学、电子学、生物学、神经学、社会学和经济学中振荡过程的基本模型。该模型很快被推广为受迫范德珀尔振子系统。
物理意义与动力学行为
范德珀尔振子在非线性阻尼参数不同时，存在以下三种情况：
（1）当 [math]\displaystyle{\quad \mu \gt 0 }[/math] 时：方程退化为简谐振荡器[math]\displaystyle{ \ddot{x} + x = 0, \quad \mu \gt 0 }[/math]，系统做等幅正弦振荡，没有能量耗散也没有能量注入。
（2）当[math]\displaystyle{1 > mu >> 0, \quad \mu \gt 0 }[/math]时：非线性较弱，系统表现出接近正弦波的振荡，但振幅会缓慢调整，最终稳定在一个极限环上。
（3）当[math]\displaystyle{mu >> 1, \quad \mu \gt 0 }[/math]时：非线性很强，系统表现出弛豫振荡（Relaxation Oscillation）。波形不再是正弦波，而是呈现“缓慢充电、快速放电”的锯齿状或脉冲状特征，最终也将收敛于一个极限环。这在生物神经元的动作电位模型中非常常见。
延伸
受迫范德珀尔振子
当方程右边不为0的时候，即系统受到外部激励（External Excitation）时，则成为受迫范德珀尔振子（Van der Pol-Duffing oscillator）： [math]\displaystyle{ \ddot{x} - \mu (1 - x^2)\dot{x} + x = Fcos(omegat)}[/math] 其中F是正弦驱动函数的振幅，ω是其角频率。 此系统经历一系列分岔（Bifurcation）过程，这个系统从稳定的极限环走向混沌。
霍普夫分叉
当μ从小于零移动到大于零时，原点的螺旋汇变成了螺旋源，一个极限周期骤然出现，半径为2。这是因为过渡不是通用：当 μ= 0时，微分方程都变成线性，原点变成圆形节点。
在霍普夫分叉中，极限环应当具有正比于[math]\displaystyle{ ε^(1/2) }[/math]的振幅，现令 ε=μ， 我们可能会尝试使用变量的变化[math]\displaystyle{ u\gt μ^(1/2)x }[/math]将范德珀尔振子转换为霍普夫分叉,得到：
[math]\displaystyle{ \ddot{u} +u+ u^2*dot{u} - ε *dot{u}= 0}[/math]
此式即为一个霍普夫分叉的模型。
量子物理领域的延伸——范德珀尔振子的哈密尔顿量
人们还可以为范德波尔振荡器编写一个与时间无关的哈密顿形式，方法是使用辅助二阶非线性微分方程将其增强为四维自主动力系统，如下所示： \ddot{x} - \mu(1 - x^2)\dot{x} + x = 0 \ddot{y} + \mu(1 - x^2)\dot{y} + y = 0 请注意，由于x和y变量的时间演变之间的单向耦合，原始范德波尔振荡器的动力学不受影响。可以通过这个方程组的哈密顿函数H证明： \[H(x, y, p_x, p_y) = p_x p_y + xy - \mu (1 - x^2)y p_y,\] 其中\[p_x=dot{y}+\mu(1-x^2)y]和\[p_y=dot{x}]（即为状态空间方程组中的二式）分别是对应于x和y的共轭矩。 原则上，这可能会导致范德波尔振荡器的量子化，这样的哈密顿量也将具有时间依赖参数的极限周期系统的几何相位与相应的哈密顿量系统的汉奈角连接起来[23]。
关于量子振荡器 量子范德波尔振荡器是经典范德波尔振荡器的量子力学版本，使用林德布拉德方程来研究其量子动力学和量子同步。 请注意，上述具有辅助二阶方程的哈密顿方法产生无界相空间轨迹，因此不能用于量子化范德波尔振荡器。在弱非线性极限（即μ→0）中，范德波尔振荡器简化为Stuart-Landau方程。Stuart-Landau方程实际上描述了弱非线性极限中的一整类极限周期振荡器。经典Stuart-Landau方程的形式要简单得多，也许不足为奇的是，可以用Lindblad方程量化，Lindblad方程也比van der Pol振荡器的Lindblad方程更简单。量子Stuart-Landau模型在量子同步研究中发挥了重要作用[25][26]（它通常被称为范德波尔振荡器，尽管它不能与范德波尔振荡器唯一关联）。经典Stuart-Landau模型（μ→0）与更一般的极限周期振荡器（任意μ）之间的关系也在相应的量子模型中进行了数值证明。 然而，上述借助辅助变量构建哈密顿量的方法虽在数学上具有巧思，却止步于量子化的物理门槛——其产生的无界相空间轨迹无法构造物理上允许的束缚态，因而不能作为量子化的有效出发点。这恰恰揭示了经典非保守系统（如范德波尔振子）量子化的本质困难：系统因非线性阻尼与外界环境存在持续的能量交换，必须将其视为一个开放量子系统。 要跨越这一鸿沟，理论物理学家转向了Lindblad主方程形式，通过引入增益和耗散两个独立的退相干通道，成功构建了量子范德波尔振子的标准模型。这一理论突破为此后的实验物理铺平了道路。近年来，随着离子阱、超导电路等量子模拟平台的飞速发展，该理论模型终于得以在实验室中精确“搭建”。 2025年，中国科学技术大学林毅恒教授团队联合深圳大学及德国斯图加特大学等机构的研究人员，正是基于上述开放量子系统的理论框架，首次在离子阱实验中实现了量子范德波尔振子。相关成果以“Experimental realization and synchronization of a quantum van der Pol oscillator”为题，于2025年10月10日发表在国际学术期刊《Science Advances》上。该研究演示了量子极限环的形成、量子同步及其增强等非经典行为。
应用与实例
电子电气领域
范德珀尔振子在电子电气领域有着广泛而根基性的应用，将LRC振荡电路中的震荡过程和电阻分别与振子震荡过程和非线性阻尼相类比，即可发现其非线性动力学性质，其核心的极限环和非线性阻尼特性直接对应于含负阻元件的振荡电路行为。例如，在使用真空管的早期电路中，当系统在极限环附近受到外部周期信号驱动时，振荡频率会被“拖入”驱动信号的频率，两者形成同步——这一现象被称为频率牵引或注入锁定。 范德珀尔方程最初即源于对含非线性负阻元件（真空管、隧道二极管等）的LC振荡电路的分析。对电路应用基尔霍夫定律，可推导出完全相同的数学形式。方程中的非线性阻尼项精确描述了负阻元件的电流-电压特性——小信号时呈现负电阻（注入能量），大信号时呈现正电阻（耗散能量），这正是电子管放大电路能够产生稳定正弦波振荡的物理机制。 范德珀尔模型已被证明是分析现代RC振荡器的通用振荡器模型，适用于弱非线性和强非线性两种 regimes。在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方程可以在很好的近似程度上被约化为一个受外部电流驱动的广义范德珀尔振荡器，该电流作为外部常力驱动动作电位的产生。FitzHugh和Nagumo认识到范德珀尔电振荡器与振荡神经元之间的相似性。该模型至今仍是计算神经科学中的标准模型。
催生混沌理论的诞生
范德珀尔振子及受迫范德珀尔振子呈现的分叉性质，为最初被人类观测到的混沌现象之一，叩开了混沌世界的大门，后续基于此系统的相关研究，使混沌理论进一步发展，如1945年英国数学家卡特赖特（M.Cartwright）和李特伍德（J.Littlewood）在对受迫范德珀尔振子的研究中，在理论层面首次证明受迫范德珀尔方程存在奇异解，预示了混沌的可能性。随后，日本学者上田睆亮（Ueda Yoshisuke）于1961年首次在模拟计算机上观察到混沌吸引子，随着后续研究的深入，混沌理论的面纱也逐渐被科学家们所揭晓。