更改

=概述及历史背景=

=概述及历史背景=

范德珀尔振子（Van der Pol oscillator）是非线性动力学中最经典、最具代表性的自激振荡模型之一。它通过一个简明的二阶非线性微分方程，特别是其中关键的振幅相关非线性阻尼项，就能够展示线性振子所没有的稳定极限环、振幅不同而产生的分岔及其能量流向、自激振荡以及随非线性增强而出现的松弛震荡（亦称弛豫振荡）等核心非线性现象，因此常被视为理解极限环理论和非线性振荡的"原型系统"。

范德珀尔振子（Van der Pol oscillator）是非线性动力学中最经典、最具代表性的自激振荡模型之一。它通过一个简明的二阶非线性[[微分方程]]，特别是其中关键的振幅相关非线性阻尼项，就能够展示线性振子所没有的稳定[[极限环]]、振幅不同而产生的分岔及其能量流向、自激振荡以及随非线性增强而出现的[[松弛震荡]]（亦称弛豫振荡）等核心非线性现象，因此常被视为理解极限环理论和非线性振荡的"原型系统"。

这一模型的诞生有着深刻的时代背景。19世纪末20世纪初，随着动力学理论与电子电气技术的持续发展，1880年Gérard-Lescuyer的串联电机<ref name="Gerard-Lescuyer_Sur_un_paradoxe_electrodynamique" />、1901年Duddell的音乐电弧<ref name="Nature_Musical_Arcs" />以及1907年de Forest发明的真空三极管<ref name="De_Forest_The_Audion" />等装置，都展示了"自激振荡"现象。然而，面对这些自激振动现象，传统线性振动理论无法解释极限环及分岔现象，传统摄动分析方法也面临级数发散问题，始终缺失解释不同自激振动现象的统一理论体系。

这一模型的诞生有着深刻的时代背景。19世纪末20世纪初，随着动力学理论与电子电气技术的持续发展，1880年Gérard-Lescuyer的串联电机<ref name="Gerard-Lescuyer_Sur_un_paradoxe_electrodynamique" />、1901年Duddell的音乐电弧<ref name="Nature_Musical_Arcs" />以及1907年de Forest发明的真空三极管<ref name="De_Forest_The_Audion" />等装置，都展示了"自激振荡"现象。然而，面对这些自激振动现象，传统线性振动理论无法解释极限环及分岔现象，传统摄动分析方法也面临级数发散问题，始终缺失解释不同自激振动现象的统一理论体系。

在此背景下，时任飞利浦公司电气工程师的荷兰物理学家巴尔塔萨·范德珀尔（Balthasar van der Pol）于1920年采用非线性动力学分析方法，引入非线性阻尼概念，提出了这一二阶非线性模型<ref name="Van_Der_Pol_A_theory_of_the_amplitude_of_free_and_forced_triode_vibrations" />。该模型最初用于描述早期无线电设备中三极管电路的自持振荡现象，是首个被深入研究的"自激振荡"非线性模型。范德珀尔在研究该系统时开创性地提出了松弛震荡（张弛振荡，Relaxation Oscillation）这一概念，并系统阐述了此现象，同时开创了"非线性振荡"的系统性研究方法。他与Van der Mark于1927年在《Nature》上发表了关于"不规则噪声"的报告<ref name="Van_Der_Pol_Van_Der_Mark_Frequency_demultiplication" />，后被公认为混沌现象最早的实验观察之一，间接催生了后续混沌理论的诞生。

在此背景下，时任飞利浦公司电气工程师的荷兰物理学家巴尔塔萨·范德珀尔（Balthasar van der Pol）于1920年采用非线性动力学分析方法，引入非线性阻尼概念，提出了这一二阶非线性模型<ref name="Van_Der_Pol_A_theory_of_the_amplitude_of_free_and_forced_triode_vibrations" />。该模型最初用于描述早期无线电设备中三极管电路的自持振荡现象，是首个被深入研究的"自激振荡"非线性模型。范德珀尔在研究该系统时开创性地提出了松弛震荡（张弛振荡，Relaxation Oscillation）这一概念，并系统阐述了此现象，同时开创了"非线性振荡"的系统性研究方法。他与Van der Mark于1927年在《Nature》上发表了关于"不规则噪声"的报告<ref name="Van_Der_Pol_Van_Der_Mark_Frequency_demultiplication" />，后被公认为[[混沌]]现象最早的实验观察之一，间接催生了后续混沌理论的诞生。

后来人们发现，这种机制并不局限于电子电路，而广泛出现在生物节律、神经活动、心脏起搏、机械振动等系统中，因此范德珀尔振子也成为从具体物理问题走向一般非线性动力学理论的一个标志性模型。

后来人们发现，这种机制并不局限于电子电路，而广泛出现在生物节律、神经活动、心脏起搏、机械振动等系统中，因此范德珀尔振子也成为从具体物理问题走向一般非线性动力学理论的一个标志性模型。

作为一个非保守非线性振荡器，范德珀尔振子的基本模型为一个二阶非线性微分方程：

作为一个非保守非线性振荡器，范德珀尔振子的基本模型为一个二阶非线性微分方程：

其中  μ 为非线性阻尼系数，中间的非线性阻尼项  -μ(1-x^2)ẋ决定了系统的全部动力学特性。当方程右侧不为0时，模型变为受迫范德珀尔振子，可进一步展示分岔、混沌等更复杂的非线性现象。

其中  μ 为非线性阻尼系数，中间的非线性阻尼项  -μ(1-x^2)ẋ决定了系统的全部动力学特性。当方程右侧不为0时，模型变为受迫范德珀尔振子，可进一步展示分岔、混沌等更复杂的非线性现象。

=主要内容=

=主要内容=

==数学模型==

==数学模型==

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

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

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

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

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

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

===位移不同时的能量传递===

===位移不同时的能量传递===

在不同阻尼系数时，系统通过[[霍普夫分岔]]从平凡平衡点产生振荡。其能量调节机制完全取决于位移 x 的大小：当位移较小（|x| < 1）时，总阻尼项为负，系统注入能量（[[负阻尼]]发散）；当位移较大（|x| > 1）时，总阻尼项为正，系统[[耗散]]能量（正阻尼收敛）<ref name="Strogatz_Example_7.3.1" />。 这种机制最终使系统收敛于唯一极限环，从而实现近似正弦振荡或自激振荡。它是解释霍普夫分岔产生极限环现象的典型模型。

在不同阻尼系数时，系统通过[[霍普夫分岔]]从平凡平衡点产生振荡。其能量调节机制完全取决于位移 x 的大小：当位移较小（|x| < 1）时，总阻尼项为负，系统注入能量（[[负阻尼]]发散）；当位移较大（|x| > 1）时，总阻尼项为正，系统[[耗散]]能量（正阻尼收敛）<ref name="Strogatz_Example_7.3.1" />。 这种机制最终使系统收敛于唯一极限环，从而实现近似正弦振荡或自激振荡。它是解释[[霍普夫分岔]]产生极限环现象的典型模型。

[[文件:阻尼系数不同时的范德珀尔振子位移随时间变化图.png|缩略图|阻尼系数不同的范德珀尔振子位移随时间变化图]]

[[文件:阻尼系数不同时的范德珀尔振子位移随时间变化图.png|缩略图|阻尼系数不同的范德珀尔振子位移随时间变化图]]

[[文件:阻尼系数不同的范德珀尔振子机械能随时间变化图.png|缩略图|阻尼系数不同的范德珀尔振子机械能随时间变化图]]

[[文件:阻尼系数不同的范德珀尔振子机械能随时间变化图.png|缩略图|阻尼系数不同的范德珀尔振子机械能随时间变化图]]

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

[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的共轭矩。原则上，这可能会导致范德波尔振荡器的量子化，这样的哈密顿量也将具有时间依赖参数的极限周期系统的[[几何相位]]与相应的哈密顿量系统的[[汉奈角]]连接起来<ref name="Chattopadhyay_Finding_the_Hannay_angle" />。

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

===关于量子振荡器===

===关于量子振荡器===

量子范德波尔振荡器是经典范德波尔振荡器的量子力学版本，使用[[林德布拉德方程]]来研究其量子动力学和[[量子同步]]<ref name="Chia_Relaxation_oscillations_and_frequency_entrainment_in_quantum_mechanics" />。请注意，上述具有辅助二阶方程的哈密顿方法产生无界相空间轨迹，因此无法构造物理上允许的束缚态，不能直接用于量子化范德波尔振荡器。在弱非线性极限（即 [math]\displaystyle{ \mu \to 0 }[/math]）中，范德波尔振荡器简化为[[Stuart-Landau方程]]。Stuart-Landau方程实际上描述了弱非线性极限中的一整类极限周期振荡器。经典Stuart-Landau方程的形式要简单得多，可以用Lindblad方程量化。

量子范德波尔振荡器是经典范德波尔振荡器的量子力学版本，使用[[林德布拉德方程]]来研究其量子动力学和[[量子同步]]<ref name="Chia_Relaxation_oscillations_and_frequency_entrainment_in_quantum_mechanics" />。请注意，上述具有辅助二阶方程的哈密顿方法产生无界[[相空间]]轨迹，因此无法构造物理上允许的束缚态，不能直接用于量子化范德波尔振荡器。在弱非线性极限（即 [math]\displaystyle{ \mu \to 0 }[/math]）中，范德波尔振荡器简化为[[Stuart-Landau方程]]。[[Stuart-Landau方程]]实际上描述了弱非线性极限中的一整类极限周期振荡器。经典Stuart-Landau方程的形式要简单得多，可以用[[Lindblad方程]]量化。

量子Stuart-Landau模型在量子[[同步]]研究中发挥了重要作用<ref name="Walter_Quantum_Synchronization_of_a_Driven_Self-Sustained_Oscillator" /><ref name="Lee_Quantum_Synchronization_of_Quantum_van_der_Pol_Oscillators_with_Trapped_Ions" />。经典Stuart-Landau模型（[math]\displaystyle{ \mu \to 0 }[/math]）与更一般的极限周期振荡器（任意\mu）之间的关系也在相应的量子模型中进行了数值证明。

量子Stuart-Landau模型在量子[[同步]]研究中发挥了重要作用<ref name="Walter_Quantum_Synchronization_of_a_Driven_Self-Sustained_Oscillator" /><ref name="Lee_Quantum_Synchronization_of_Quantum_van_der_Pol_Oscillators_with_Trapped_Ions" />。经典Stuart-Landau模型（[math]\displaystyle{ \mu \to 0 }[/math]）与更一般的极限周期振荡器（任意\mu）之间的关系也在相应的量子模型中进行了数值证明。

然而，上述借助辅助变量构建哈密顿量的方法虽在数学上具有巧思，却止步于量子化的物理门槛——其产生的无界相空间轨迹无法构造物理上允许的束缚态。这恰恰揭示了经典非保守系统（如范德波尔振子）量子化的本质困难：系统因非线性阻尼与外界环境存在持续的能量交换，必须将其视为一个[[开放量子系统]]。

然而，上述借助辅助变量构建哈密顿量的方法虽在数学上具有巧思，却止步于量子化的物理门槛——其产生的无界相空间轨迹无法构造物理上允许的束缚态。这恰恰揭示了经典非保守系统（如范德波尔振子）量子化的本质困难：系统因非线性阻尼与外界环境存在持续的能量交换，必须将其视为一个[[开放量子系统]]。

要跨越这一鸿沟，理论物理学家转向了Lindblad主方程形式，通过引入增益和耗散两个独立的退相干通道，成功构建了量子范德珀尔振子的标准模型。近年来，随着离子阱、超导电路等量子模拟平台的飞速发展，该理论模型终于得以在实验室中精确"搭建"<ref name="Chia_Relaxation_oscillations_and_frequency_entrainment_in_quantum_mechanics" />。

要跨越这一鸿沟，理论物理学家转向了Lindblad主方程形式，通过引入增益和耗散两个独立的退相干通道，成功构建了量子范德珀尔振子的标准模型。近年来，随着[[离子阱]]、[[超导电路]]等量子模拟平台的飞速发展，该理论模型终于得以在实验室中精确"搭建"<ref name="Chia_Relaxation_oscillations_and_frequency_entrainment_in_quantum_mechanics" />。

2025年，中国科学技术大学林毅恒教授团队联合深圳大学及德国斯图加特大学等机构的研究人员，基于上述开放量子系统的理论框架，首次在离子阱实验中实现了量子范德珀尔振子。相关成果于2025年10月10日发表在国际学术期刊《Science Advances》上，演示了量子极限环的形成、量子同步及其增强等非经典行为<ref name="CAS_cn_2025_van_der_Pol" />。

2025年，中国科学技术大学林毅恒教授团队联合深圳大学及德国斯图加特大学等机构的研究人员，基于上述[[开放量子系统]]的理论框架，首次在离子阱实验中实现了量子范德珀尔振子。相关成果于2025年10月10日发表在国际学术期刊《Science Advances》上，演示了[[量子极限环]]的形成、[[量子同步]]及其增强等非经典行为<ref name="CAS_cn_2025_van_der_Pol" />。

=应用与实例=

=应用与实例=

==电子电气领域==

==电子电气领域==

范德珀尔振子在电子电气领域有着广泛而根基性的应用。其核心的极限环和非线性阻尼特性直接对应于含负阻元件的振荡电路行为。例如，在使用真空管的早期电路中，当系统在极限环附近受到外部周期信号驱动时，振荡频率会被"拖入"驱动信号的频率，两者形成同步——这一现象被称为[[频率牵引]]或[[注入锁定]]。

范德珀尔振子在电子电气领域有着广泛而根基性的应用。其核心的极限环和非线性阻尼特性直接对应于含[[负阻]]元件的振荡电路行为。例如，在使用真空管的早期电路中，当系统在极限环附近受到外部周期信号驱动时，振荡频率会被"拖入"驱动信号的频率，两者形成[[同步]]——这一现象被称为[[频率牵引]]或[[注入锁定]]。

范德珀尔方程最初即源于对含非线性负阻元件（真空管、隧道二极管等）的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"的论文，建立了心脏的电子电路模型，将心跳视为松弛震荡<ref name="Van_Der_Pol_Van_Der_Mark_The_heartbeat_considered_as_a_relaxation_oscillation" />。

范德珀尔振子的极限环概念精准对应于心脏[[窦房结]]（SA node）和[[房室结]]（AV node）等天然起搏器的自发节律放电。范德珀尔和van der Mark早在1928年即于《Philosophical Magazine》发表了题为"The heartbeat considered as a relaxation oscillation, and an electrical model of the heart"的论文，建立了心脏的电子电路模型，将心跳视为松弛震荡<ref name="Van_Der_Pol_Van_Der_Mark_The_heartbeat_considered_as_a_relaxation_oscillation" />。

此后，改良的范德珀尔振荡器被广泛用于模拟心脏起搏器的动作电位时序。Grudziński与Żebrowski在2004年设计了修正的范德波尔振荡器模型，通过在经典范德波尔振荡器的相空间中增加一个鞍点和一个节点，并引入非对称阻尼项，使得起搏节点的放电频率可在不改变不应期长度的条件下独立调节<ref name="Grudzinski_Zebrowski_Modeling_cardiac_pacemakers_with_relaxation_oscillators" />。

此后，改良的范德珀尔振荡器被广泛用于模拟心脏起搏器的[[动作电位]]时序。Grudziński与Żebrowski在2004年设计了修正的范德波尔振荡器模型，通过在经典范德波尔振荡器的相空间中增加一个鞍点和一个节点，并引入非对称阻尼项，使得起搏节点的放电频率可在不改变[[不应期]]长度的条件下独立调节<ref name="Grudzinski_Zebrowski_Modeling_cardiac_pacemakers_with_relaxation_oscillators" />。

在神经科学中，FitzHugh于1961年将范德珀尔方程与[[Hodgkin-Huxley模型]]结合，提出了[[FitzHugh-Nagumo模型]]（亦称[[Bonhoeffer-van der Pol模型]]）<ref name="Fitzhugh_Impulses_and_physiological_states_in_theoretical_models_of_nerve_membrane" />。该模型的相平面中存在一个稳定的极限环，对应于神经元的周期性脉冲发放。研究表明，FitzHugh-Nagumo方程可以在很好的近似程度上被约化为一个受外部电流驱动的广义范德珀尔振荡器。

在神经科学中，FitzHugh于1961年将范德珀尔方程与[[Hodgkin-Huxley模型]]结合，提出了[[FitzHugh-Nagumo模型]]（亦称[[Bonhoeffer-van der Pol模型]]）<ref name="Fitzhugh_Impulses_and_physiological_states_in_theoretical_models_of_nerve_membrane" />。该模型的相平面中存在一个稳定的极限环，对应于[[神经元]]的周期性脉冲发放。研究表明，FitzHugh-Nagumo方程可以在很好的近似程度上被约化为一个受外部电流驱动的广义范德珀尔振荡器。

==催生混沌理论的诞生==

==催生混沌理论的诞生==