133
个编辑
更改
→范德珀尔振子
=概述及历史背景=
范德珀尔振子（Van der Pol oscillator）是非线性动力学中最经典、最具代表性的自激振荡模型之一。它通过一个简明的二阶非线性微分方程，就能够展示线性振子所没有的稳定极限环、自激振荡以及随非线性增强而出现的松弛震荡（亦称弛豫振荡）等核心非线性现象，因此常被视为理解极限环理论和非线性振荡的“原型系统”。
这一模型的诞生有着深刻的时代背景。19世纪末20世纪初，随着动力学理论与电子电气技术的持续发展，1880年Gérard-Lescuyer的串联电机[1]、1901年Duddell的音乐电弧[2]以及1907年de Forest发明的真空三极管[3]等装置，都展示了“自激振荡”现象。然而，面对这些自激振动现象，传统线性振动理论无法解释极限环及分岔现象，传统摄动分析方法也面临级数发散问题，始终缺失解释不同自激振动现象的统一理论体系。
[math]\displaystyle{ \ddot{x} - \mu (1 - x^2)\dot{x} + x = 0, \quad \mu > 0 }[/math]
在此背景下，时任飞利浦公司电气工程师的荷兰物理学家巴尔塔萨·范德珀尔（Balthasar van der Pol）于1920年采用非线性动力学分析方法，引入非线性阻尼概念，在研究真空电子管振荡电路时提出了这一二阶非线性模型[4]。该模型最初用于描述早期无线电设备中三极管电路的自持振荡现象，是首个被深入研究的“自激振荡”非线性模型。范德珀尔在研究该系统时开创性地提出了松弛震荡（张弛振荡，Relaxation Oscillation）这一概念，并系统阐述了此现象，同时开创了“非线性振荡”的系统性研究方法。他与Van der Mark于1927年在《Nature》上发表了关于“不规则噪声”的报告[5]，后被公认为混沌现象最早的实验观察之一，间接催生了后续混沌理论的诞生。
后来人们发现，这种机制并不局限于电子电路，而广泛出现在生物节律、神经活动、心脏起搏、机械振动等系统中，因此范德珀尔振子也成为从具体物理问题走向一般非线性动力学理论的一个标志性模型。
作为一个非保守非线性振荡器，范德珀尔振子的基本模型为一个二阶非线性微分方程：
x ̈-μ(1-x^2)x ̇+x=0, μ>0
其中 μ 为非线性阻尼系数，中间的非线性阻尼项 -μ(1-x^2)x ̇ 决定了系统的全部动力学特性。当方程右侧不为0时，模型变为受迫范德珀尔振子，可进一步展示分岔、混沌等更复杂的非线性现象。
该模型将线性理论、分岔理论、极限环、快慢动力学乃至混沌的边缘整合了起来，实现了由传统线性振荡理论向非线性动力学、乃至于混沌理论的飞跃。
=主要内容=
==数学模型==
[[文件:62ed93b61faa264fce7d907959fdbe8f.png|缩略图|阻尼系数不同时的范德珀尔振子动力学相图，可见当μ非0时，随着振子不断振荡，最终都会收敛于一个极限环，且当μ显著大于时，系统出现了节律性的弛豫震荡的现象。]]
[[文件:62ed93b61faa264fce7d907959fdbe8f.png|缩略图|阻尼系数不同时的范德珀尔振子动力学相图，可见当μ非0时，随着振子不断振荡，最终都会收敛于一个极限环，且当μ显著大于时，系统出现了节律性的弛豫震荡的现象。]]
范德珀尔振子（Van der Pol oscillator）是研究[[极限环]][[非线性系统]]的经典实例。作为一个耗散性非线性振荡系统，其数学模型为一个二阶非线性齐次方程：
范德珀尔振子（Van der Pol oscillator）是研究[[极限环]][[非线性系统]]的经典实例。作为一个耗散性非线性振荡系统，其数学模型为一个二阶非线性齐次方程：
（3）\mu 是非线性阻尼参数（\mu \geq 0），控制着非线性效应的强弱。
（3）\mu 是非线性阻尼参数（\mu \geq 0），控制着非线性效应的强弱。
==物理意义与动力学行为==
===阻尼参数不同时的振动性质===
范德珀尔振子在非线性阻尼参数不同时，存在以下三种情况：
范德珀尔振子在非线性阻尼参数不同时，存在以下三种情况：
（3）当 [math]\displaystyle{ \mu \gg 1 }[/math] 时：非线性很强，系统表现出松弛震荡（Relaxation Oscillation）。波形不再是正弦波，而是呈现“缓慢充电、快速放电”的锯齿状或脉冲状特征，最终也将收敛于一个极限环。这在生物[[神经元]]的[[动作电位]]模型中非常常见。
（3）当 [math]\displaystyle{ \mu \gg 1 }[/math] 时：非线性很强，系统表现出松弛震荡（Relaxation Oscillation）。波形不再是正弦波，而是呈现“缓慢充电、快速放电”的锯齿状或脉冲状特征，最终也将收敛于一个极限环。这在生物[[神经元]]的[[动作电位]]模型中非常常见。
===位移不同时的能量传递===
在不同阻尼系数时，系统通过[[霍普夫分岔]]从平凡平衡点产生振荡。其能量调节机制完全取决于位移 x 的大小：当位移较小（|x| < 1）时，总阻尼项为负，系统注入能量（[[负阻尼]]发散）；当位移较大（|x| > 1）时，总阻尼项为正，系统[[耗散]]能量（正阻尼收敛）[1]。 这种机制最终使系统收敛于唯一极限环，从而实现近似正弦振荡或自激振荡。它是解释霍普夫分岔产生极限环现象的典型模型。
在不同阻尼系数时，系统通过[[霍普夫分岔]]从平凡平衡点产生振荡。其能量调节机制完全取决于位移 x 的大小：当位移较小（|x| < 1）时，总阻尼项为负，系统注入能量（[[负阻尼]]发散）；当位移较大（|x| > 1）时，总阻尼项为正，系统[[耗散]]能量（正阻尼收敛）[6]。 这种机制最终使系统收敛于唯一极限环，从而实现近似正弦振荡或自激振荡。它是解释霍普夫分岔产生极限环现象的典型模型。
[[文件:阻尼系数不同时的范德珀尔振子位移随时间变化图.png|缩略图|阻尼系数不同的范德珀尔振子位移随时间变化图]]
[[文件:阻尼系数不同时的范德珀尔振子位移随时间变化图.png|缩略图|阻尼系数不同的范德珀尔振子位移随时间变化图]]
[[文件:阻尼系数不同的范德珀尔振子机械能随时间变化图.png|缩略图|阻尼系数不同的范德珀尔振子机械能随时间变化图]]
[[文件:阻尼系数不同的范德珀尔振子机械能随时间变化图.png|缩略图|阻尼系数不同的范德珀尔振子机械能随时间变化图]]
==模型的推导过程==
我们可通过如下步骤来探究[7]：考虑如图所示的RLC电路：它由电感L、电容器C、电阻R组成，我们假设电压源是电池（[math]\displaystyle{ E(t)=E }[/math]常数）。当电池开关关闭时，根据基尔霍夫电压定律：
我们可通过如下步骤来探究[7]：考虑如图所示的RLC电路：它由电感L、电容器C、电阻R组成，我们假设电压源是电池（[math]\displaystyle{ E(t)=E }[/math]常数）。当电池开关关闭时，根据基尔霍夫电压定律：
自引入以来，范德波尔方程一直被用作物理学、电子学、生物学、神经学、社会学和经济学中振荡过程的基本模型。该模型很快被推广为受迫范德珀尔振子系统。
自引入以来，范德波尔方程一直被用作物理学、电子学、生物学、神经学、社会学和经济学中振荡过程的基本模型。该模型很快被推广为受迫范德珀尔振子系统。
=延伸=
==受迫范德珀尔振子==
当方程右边不为0的时候，即系统受到[[外部激励]]（External Excitation）时，则成为受迫范德珀尔振子（Forced Van der Pol oscillator）：
当方程右边不为0的时候，即系统受到[[外部激励]]（External Excitation）时，则成为受迫范德珀尔振子（Forced Van der Pol oscillator）：
[math]\displaystyle{ \ddot{x} - \mu (1 - x^2)\dot{x} + x = F \cos(\omega t) }[/math]
[math]\displaystyle{ \ddot{x} - \mu (1 - x^2)\dot{x} + x = F \cos(\omega t) }[/math]
（注：受迫VDP振子因引入时间维度而升维至三维及以上，打破了[[庞加莱-本迪克松定理]]的限制，从而触及混沌理论的门槛。日本学者上田睆亮于1961年在此系统中首次观察到了[[混沌吸引子]][8]。）
（注：受迫VDP振子因引入时间维度而升维至三维及以上，打破了[[庞加莱-本迪克松定理]]的限制，从而触及混沌理论的门槛。日本学者上田睆亮于1961年在此系统中首次观察到了[[混沌吸引子]][8]。）
==霍普夫分叉==
当\mu从小于零移动到大于零时，原点的螺旋汇变成了螺旋源，一个极限周期骤然出现，半径为2。这是因为过渡不是通用的：当 [math]\displaystyle{ \mu = 0 }[/math] 时，微分方程变成线性，原点变成[[中心]]（Center）（而非节点）。
当\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{ \varepsilon^{1/2} }[/math] 的振幅，现令 [math]\displaystyle{ \varepsilon = \mu }[/math]，我们可能会尝试使用变量的变化 [math]\displaystyle{ u = \mu^{1/2} x }[/math] 将范德珀尔振子转换为霍普夫分叉，得到：
此式即为一个霍普夫分叉的模型[9]。
此式即为一个霍普夫分叉的模型[9]。
==量子物理领域——范德珀尔振子的[[哈密顿量]]==
人们还可以为范德波尔振荡器编写一个与时间无关的哈密顿形式，方法是使用辅助二阶非线性微分方程将其增强为四维自主[[动力系统]]，如下所示：
人们还可以为范德波尔振荡器编写一个与时间无关的哈密顿形式，方法是使用辅助二阶非线性微分方程将其增强为四维自主[[动力系统]]，如下所示：
其中 [math]\displaystyle{ p_x = \dot{y} + \mu(1 - x^2)y }[/math] 和 [math]\displaystyle{ p_y = \dot{x} }[/math] 分别是对应于x和y的共轭矩。原则上，这可能会导致范德波尔振荡器的量子化，这样的哈密顿量也将具有时间依赖参数的极限周期系统的[[几何相位]]与相应的哈密顿量系统的[[汉奈角]]连接起来[11]。
其中 [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方程量化。
量子范德波尔振荡器是经典范德波尔振荡器的量子力学版本，使用[[林德布拉德方程]]来研究其量子动力学和[[量子同步]][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）之间的关系也在相应的量子模型中进行了数值证明。
量子Stuart-Landau模型在量子[[同步]]研究中发挥了重要作用[13][14]。经典Stuart-Landau模型（[math]\displaystyle{ \mu \to 0 }[/math]）与更一般的极限周期振荡器（任意\mu）之间的关系也在相应的量子模型中进行了数值证明。
2025年，中国科学技术大学林毅恒教授团队联合深圳大学及德国斯图加特大学等机构的研究人员，基于上述开放量子系统的理论框架，首次在离子阱实验中实现了量子范德波尔振子。相关成果于2025年10月10日发表在国际学术期刊《Science Advances》上，演示了量子极限环的形成、量子同步及其增强等非经典行为[15]。
2025年，中国科学技术大学林毅恒教授团队联合深圳大学及德国斯图加特大学等机构的研究人员，基于上述开放量子系统的理论框架，首次在离子阱实验中实现了量子范德波尔振子。相关成果于2025年10月10日发表在国际学术期刊《Science Advances》上，演示了量子极限环的形成、量子同步及其增强等非经典行为[15]。
=应用与实例=
==电子电气领域==
范德珀尔振子在电子电气领域有着广泛而根基性的应用。其核心的极限环和非线性阻尼特性直接对应于含负阻元件的振荡电路行为。例如，在使用真空管的早期电路中，当系统在极限环附近受到外部周期信号驱动时，振荡频率会被“拖入”驱动信号的频率，两者形成同步——这一现象被称为[[频率牵引]]或[[注入锁定]]。
范德珀尔振子在电子电气领域有着广泛而根基性的应用。其核心的极限环和非线性阻尼特性直接对应于含负阻元件的振荡电路行为。例如，在使用真空管的早期电路中，当系统在极限环附近受到外部周期信号驱动时，振荡频率会被“拖入”驱动信号的频率，两者形成同步——这一现象被称为[[频率牵引]]或[[注入锁定]]。
范德珀尔方程最初即源于对含非线性负阻元件（真空管、隧道二极管等）的LC振荡电路的分析。方程中的非线性阻尼项精确描述了负阻元件的电流-电压特性——小信号时呈现负电阻（注入能量），大信号时呈现正电阻（耗散能量），这正是电子管放大电路能够产生稳定正弦波振荡的物理机制。范德珀尔模型已被证明是分析现代RC振荡器的通用振荡器模型。在CMOS纳米技术中，范德珀尔参数与多种电路实现之间的对应关系也已得到系统阐述。注入锁定的范德珀尔振荡器可以用简单的现代半导体电路实现。
范德珀尔方程最初即源于对含非线性负阻元件（真空管、隧道二极管等）的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”的论文，建立了心脏的电子电路模型，将心跳视为松弛震荡[16]。
范德珀尔振子的极限环概念精准对应于心脏[[窦房结]]（SA node）和[[房室结]]（AV node）等天然起搏器的自发节律放电。范德珀尔和van der Mark早在1928年即于《Philosophical Magazine》发表了题为“The heartbeat considered as a relaxation oscillation, and an electrical model of the heart”的论文，建立了心脏的电子电路模型，将心跳视为松弛震荡[16]。
此后，改良的范德珀尔振荡器被广泛用于模拟心脏起搏器的动作电位时序。Grudziński与Żebrowski在2004年设计了修正的范德波尔振荡器模型，通过在经典范德波尔振荡器的相空间中增加一个鞍点和一个节点，并引入非对称阻尼项，使得起搏节点的放电频率可在不改变不应期长度的条件下独立调节[17]。
此后，改良的范德珀尔振荡器被广泛用于模拟心脏起搏器的动作电位时序。Grudziński与Żebrowski在2004年设计了修正的范德波尔振荡器模型，通过在经典范德波尔振荡器的相空间中增加一个鞍点和一个节点，并引入非对称阻尼项，使得起搏节点的放电频率可在不改变不应期长度的条件下独立调节[17]。
在神经科学中，FitzHugh于1961年将范德珀尔方程与[[Hodgkin-Huxley模型]]结合，提出了[[FitzHugh-Nagumo模型]]（亦称[[Bonhoeffer-van der Pol模型]]）[18]。该模型的相平面中存在一个稳定的极限环，对应于神经元的周期性脉冲发放。研究表明，FitzHugh-Nagumo方程可以在很好的近似程度上被约化为一个受外部电流驱动的广义范德珀尔振荡器。
在神经科学中，FitzHugh于1961年将范德珀尔方程与[[Hodgkin-Huxley模型]]结合，提出了[[FitzHugh-Nagumo模型]]（亦称[[Bonhoeffer-van der Pol模型]]）[18]。该模型的相平面中存在一个稳定的极限环，对应于神经元的周期性脉冲发放。研究表明，FitzHugh-Nagumo方程可以在很好的近似程度上被约化为一个受外部电流驱动的广义范德珀尔振荡器。
==催生混沌理论的诞生==
范德珀尔振子及受迫范德珀尔振子呈现的分叉性质，为最初被人类观测到的混沌现象之一，叩开了混沌世界的大门。1945年英国数学家卡特赖特（M.Cartwright）和李特伍德（J.Littlewood）在对受迫范德珀尔振子的研究中，在理论层面首次证明受迫范德珀尔方程存在奇异解，预示了混沌的可能性[19]。随后，日本学者上田睆亮（Ueda Yoshisuke）于1961年首次在模拟计算机上观察到混沌[[吸引子]][8]，随着后续研究的深入，混沌理论的面纱也逐渐被科学家们所揭晓。
范德珀尔振子及受迫范德珀尔振子呈现的分叉性质，为最初被人类观测到的混沌现象之一，叩开了混沌世界的大门。1945年英国数学家卡特赖特（M.Cartwright）和李特伍德（J.Littlewood）在对受迫范德珀尔振子的研究中，在理论层面首次证明受迫范德珀尔方程存在奇异解，预示了混沌的可能性[19]。随后，日本学者上田睆亮（Ueda Yoshisuke）于1961年首次在模拟计算机上观察到混沌[[吸引子]][8]，随着后续研究的深入，混沌理论的面纱也逐渐被科学家们所揭晓。
=参考文献=
1. Strogatz, Steven (2019). "Example 7.3.1". Nonlinear dynamics and chaos: with applications to
1.Gérard-Lescuyer, J.-M.-A. (1880). Sur un paradoxe électrodynamique.
2.Gérard-Lescuyer, J.-M.-A. (1880). Sur un paradoxe électrodynamique.
2.Nature (1901). Musical Arcs. Nature, 63(1640), 542–543.
3.Nature (1901). Musical Arcs. Nature, 63(1640), 542–543.
3.De Forest, L. (1906). The Audion: A New Receiver for Wireless Telegraphy. Transactions of the American Institute of Electrical Engineers.
4.De Forest, L. (1906). The Audion: A New Receiver for Wireless Telegraphy. Transactions of the American Institute of Electrical Engineers.
4.VAN DER POL B. A theory of the amplitude of free and forced triode vibrations[J]. Radio Review, 1920, 1: 701-710.
5.VAN DER POL B. A theory of the amplitude of free and forced triode vibrations[J]. Radio Review, 1920, 1: 701-710.
5.VAN DER POL B, VAN DER MARK J. Frequency demultiplication[J]. Nature, 1927, 120(3019): 363-364. DOI:10.1038/120363a0.
6.VAN DER POL B, VAN DER MARK J. Frequency demultiplication[J]. Nature, 1927, 120(3019): 363-364. DOI:10.1038/120363a0.
6. 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).
7.https://www.phys.uconn.edu/~rozman/Courses/P2200_13F/downloads/vanderpol/vanderpol-oscillator-draft.pdf
7.https://www.phys.uconn.edu/~rozman/Courses/P2200_13F/downloads/vanderpol/vanderpol-oscillator-draft.pdf
19.CARTWRIGHT M L, LITTLEWOOD J E. On non-linear differential equations of the second order: I. The equation \ddot{y} - k(1 - y^2)\dot{y} + y = b\lambda k \cos(\lambda t + a), k large[J]. Journal of the London Mathematical Society, 1945, s1-20(3): 180-189. DOI:10.1112/jlms/s1-20.3.180.
19.CARTWRIGHT M L, LITTLEWOOD J E. On non-linear differential equations of the second order: I. The equation \ddot{y} - k(1 - y^2)\dot{y} + y = b\lambda k \cos(\lambda t + a), k large[J]. Journal of the London Mathematical Society, 1945, s1-20(3): 180-189. DOI:10.1112/jlms/s1-20.3.180.
注：10—15源因子维基百科，“范德珀尔振子”词条。
==署名==
==署名==