“范德珀尔振子”的版本间的差异
| 第4行： | 第4行： | ||
范德珀尔振子（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" /> | + | 在此背景下，时任飞利浦公司电气工程师的荷兰物理学家巴尔塔萨·范德珀尔（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" />，后被公认为[[混沌]]现象最早的实验观察之一，间接催生了后续[[混沌理论]]的诞生。 |
后来人们发现，这种机制并不局限于电子电路，而广泛出现在生物节律、神经活动、心脏起搏、机械振动等系统中，因此范德珀尔振子也成为从具体物理问题走向一般非线性动力学理论的一个标志性模型。 | 后来人们发现，这种机制并不局限于电子电路，而广泛出现在生物节律、神经活动、心脏起搏、机械振动等系统中，因此范德珀尔振子也成为从具体物理问题走向一般非线性动力学理论的一个标志性模型。 | ||
作为一个[[非保守非线性振荡器]]，范德珀尔振子的基本模型为一个二阶非线性[[微分方程]]： | 作为一个[[非保守非线性振荡器]]，范德珀尔振子的基本模型为一个二阶非线性[[微分方程]]： | ||
| + | |||
[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] | ||
| + | |||
其中 [math]\displaystyle{ \mu }[/math] 为非线性阻尼系数，中间的非线性阻尼项 [math]\displaystyle{ -\mu(1-x^2)\dot{x} }[/math] 决定了系统的全部动力学特性，该模型将线性理论、[[分岔理论]]、[[极限环]]、[[快慢动力学]]乃至[[混沌的边缘]]整合了起来，实现了由传统线性振荡理论向非线性动力学、乃至于[[混沌]]理论的飞跃。 | 其中 [math]\displaystyle{ \mu }[/math] 为非线性阻尼系数，中间的非线性阻尼项 [math]\displaystyle{ -\mu(1-x^2)\dot{x} }[/math] 决定了系统的全部动力学特性，该模型将线性理论、[[分岔理论]]、[[极限环]]、[[快慢动力学]]乃至[[混沌的边缘]]整合了起来，实现了由传统线性振荡理论向非线性动力学、乃至于[[混沌]]理论的飞跃。 | ||
=主要内容= | =主要内容= | ||
==数学模型== | ==数学模型== | ||
范德珀尔振子（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] | ||
| + | |||
其中 [math]\displaystyle{ x }[/math] 为位置坐标（作为时间 [math]\displaystyle{ t }[/math] 的函数），[math]\displaystyle{ \mu }[/math] 是标量参数，二者决定阻尼的非线性和强度。 | 其中 [math]\displaystyle{ x }[/math] 为位置坐标（作为时间 [math]\displaystyle{ t }[/math] 的函数），[math]\displaystyle{ \mu }[/math] 是标量参数，二者决定阻尼的非线性和强度。 | ||
也可写为二维状态空间方程组的形式： | 也可写为二维状态空间方程组的形式： | ||
| + | |||
[math]\displaystyle{ \dot{x}=y,\quad \dot{y}=\mu(1-x^2)y-x }[/math] | [math]\displaystyle{ \dot{x}=y,\quad \dot{y}=\mu(1-x^2)y-x }[/math] | ||
| + | |||
其中： | 其中： | ||
| + | |||
（1）[math]\displaystyle{ x }[/math] 是位移（或电压等物理量）。 | （1）[math]\displaystyle{ x }[/math] 是位移（或电压等物理量）。 | ||
| + | |||
（2）[math]\displaystyle{ y=\dot{x} }[/math] 是速度，[math]\displaystyle{ \dot{y}=\ddot{x} }[/math] 是加速度。 | （2）[math]\displaystyle{ y=\dot{x} }[/math] 是速度，[math]\displaystyle{ \dot{y}=\ddot{x} }[/math] 是加速度。 | ||
| + | |||
（3）[math]\displaystyle{ \mu }[/math] 是非线性阻尼参数（[math]\displaystyle{ \mu\ge 0 }[/math]），控制着非线性效应的强弱。 | （3）[math]\displaystyle{ \mu }[/math] 是非线性阻尼参数（[math]\displaystyle{ \mu\ge 0 }[/math]），控制着非线性效应的强弱。 | ||
| + | |||
当方程左侧为0的时候，模型为自由范德珀尔振子，不受外界噪声干扰；当方程右侧不为0时，模型变为[[受迫范德珀尔振子]]，可进一步展示[[分岔]]、[[混沌]]等更复杂的非线性现象。 | 当方程左侧为0的时候，模型为自由范德珀尔振子，不受外界噪声干扰；当方程右侧不为0时，模型变为[[受迫范德珀尔振子]]，可进一步展示[[分岔]]、[[混沌]]等更复杂的非线性现象。 | ||
==物理意义与动力学行为== | ==物理意义与动力学行为== | ||
| 第26行： | 第36行： | ||
[[文件:阻尼系数不同的范德珀尔振子相图.png|缩略图|阻尼系数不同的范德珀尔振子振动相图]] | [[文件:阻尼系数不同的范德珀尔振子相图.png|缩略图|阻尼系数不同的范德珀尔振子振动相图]] | ||
[[文件:不同起振点的范德珀尔振子.png|缩略图|同一阻尼系数，不同起振点的范德珀尔振子收敛于同一极限环]] | [[文件:不同起振点的范德珀尔振子.png|缩略图|同一阻尼系数，不同起振点的范德珀尔振子收敛于同一极限环]] | ||
| − | 范德珀尔振子作为一个非线性振荡系统，其动力学特性与阻尼系数 [math]\displaystyle{ \mu }[/math] 的符号密切相关。将方程线性化于原点附近（即忽略 [math]\displaystyle{ x^2 }[/math] 项），可得线性近似方程 [math]\displaystyle{ \ddot{x}-\mu\dot{x}+x=0 }[/math] | + | 范德珀尔振子作为一个非线性振荡系统，其动力学特性与阻尼系数 [math]\displaystyle{ \mu }[/math] 的符号密切相关。将方程线性化于原点附近（即忽略 [math]\displaystyle{ x^2 }[/math] 项），可得线性近似方程 |
| + | [math]\displaystyle{ \ddot{x}-\mu\dot{x}+x=0 }[/math] | ||
| + | 其特征方程为 [math]\displaystyle{ \lambda^2-\mu\lambda+1=0 }[/math]，特征根为 [math]\displaystyle{ \lambda=(\mu\pm\sqrt{\mu^2-4})/2 }[/math]。由此可以看出，原点处的稳定性完全由 [math]\displaystyle{ \mu }[/math] 的符号决定： | ||
| + | |||
当 [math]\displaystyle{ \mu<0 }[/math] 时，特征根的实部为负，原点是一个稳定的焦点（螺旋汇），[[相空间]]中所有[[轨迹]]均螺旋收敛于原点，系统最终静止，不存在极限环。此时若存在极限环，则必为不稳定的（排斥子），任何微小的扰动都会使系统远离该极限环。 | 当 [math]\displaystyle{ \mu<0 }[/math] 时，特征根的实部为负，原点是一个稳定的焦点（螺旋汇），[[相空间]]中所有[[轨迹]]均螺旋收敛于原点，系统最终静止，不存在极限环。此时若存在极限环，则必为不稳定的（排斥子），任何微小的扰动都会使系统远离该极限环。 | ||
| + | |||
当 [math]\displaystyle{ \mu=0 }[/math] 时，特征根为纯虚数 [math]\displaystyle{ \pm i }[/math]，原点退化为[[线性中心]]，系统退化为保守的简谐振荡器 [math]\displaystyle{ \ddot{x}+x=0 }[/math]，相轨迹为一族同心圆，既不发散也不收敛。[math]\displaystyle{ \mu=0 }[/math] 因此成为[[系统动力学]]性质发生突变的[[临界点]]——[[霍普夫分岔]]点。 | 当 [math]\displaystyle{ \mu=0 }[/math] 时，特征根为纯虚数 [math]\displaystyle{ \pm i }[/math]，原点退化为[[线性中心]]，系统退化为保守的简谐振荡器 [math]\displaystyle{ \ddot{x}+x=0 }[/math]，相轨迹为一族同心圆，既不发散也不收敛。[math]\displaystyle{ \mu=0 }[/math] 因此成为[[系统动力学]]性质发生突变的[[临界点]]——[[霍普夫分岔]]点。 | ||
| + | |||
当 [math]\displaystyle{ \mu>0 }[/math] 时，特征根的实部为正，原点变为不稳定的焦点（螺旋源），任何从原点附近出发的轨迹都会被”推开”，螺旋向外运动。然而，由于非线性阻尼项的存在，轨迹并不会无限发散，而是最终被”捕获”到一个稳定的闭合轨道——即稳定的极限环上。此时系统呈现自激振荡，振幅由非线性机制自动锁定。 | 当 [math]\displaystyle{ \mu>0 }[/math] 时，特征根的实部为正，原点变为不稳定的焦点（螺旋源），任何从原点附近出发的轨迹都会被”推开”，螺旋向外运动。然而，由于非线性阻尼项的存在，轨迹并不会无限发散，而是最终被”捕获”到一个稳定的闭合轨道——即稳定的极限环上。此时系统呈现自激振荡，振幅由非线性机制自动锁定。 | ||
| + | |||
上述 [math]\displaystyle{ \mu }[/math] 穿过零点时原点稳定性翻转、极限环随之产生或消失的现象，正是[[超临界霍普夫分岔]]的典型图景。值得注意的是，[math]\displaystyle{ \mu>0 }[/math] 与 [math]\displaystyle{ \mu<0 }[/math] 两种情形在时间反演下互为镜像：将 [math]\displaystyle{ t\to -t }[/math] 代入方程，[math]\displaystyle{ \mu }[/math] 的符号恰好翻转，因此一者为另一者的逆过程，体现了系统在时空意义上的对称性。从相图上看，当 [math]\displaystyle{ \mu>0 }[/math] 时，轨迹随机从平面内各点螺旋逼近极限环的位置（类似逆时针书写”6”字形，从“尾巴”出发，最后在下面的小圈处反复“画圈”），最终被极限环”捕获”并形成稳定循环；当 [math]\displaystyle{ \mu<0 }[/math] 时则相反，轨迹从极限环螺旋向内或向外随机收敛于平面内各点（类似顺时针书写“6”字形，从下面“小圈”左上角和“尾巴”衔接的地方出发，回到这里，再画出一条长长的“尾巴”）。分岔在此扮演了”入口”的角色——没有分岔，原点便是稳定的，所有轨迹终将被吸入口径为零的平衡点，极限环与后续的自持振荡便无从谈起。 | 上述 [math]\displaystyle{ \mu }[/math] 穿过零点时原点稳定性翻转、极限环随之产生或消失的现象，正是[[超临界霍普夫分岔]]的典型图景。值得注意的是，[math]\displaystyle{ \mu>0 }[/math] 与 [math]\displaystyle{ \mu<0 }[/math] 两种情形在时间反演下互为镜像：将 [math]\displaystyle{ t\to -t }[/math] 代入方程，[math]\displaystyle{ \mu }[/math] 的符号恰好翻转，因此一者为另一者的逆过程，体现了系统在时空意义上的对称性。从相图上看，当 [math]\displaystyle{ \mu>0 }[/math] 时，轨迹随机从平面内各点螺旋逼近极限环的位置（类似逆时针书写”6”字形，从“尾巴”出发，最后在下面的小圈处反复“画圈”），最终被极限环”捕获”并形成稳定循环；当 [math]\displaystyle{ \mu<0 }[/math] 时则相反，轨迹从极限环螺旋向内或向外随机收敛于平面内各点（类似顺时针书写“6”字形，从下面“小圈”左上角和“尾巴”衔接的地方出发，回到这里，再画出一条长长的“尾巴”）。分岔在此扮演了”入口”的角色——没有分岔，原点便是稳定的，所有轨迹终将被吸入口径为零的平衡点，极限环与后续的自持振荡便无从谈起。 | ||
===阻尼参数大小不同时的振动非线性性质=== | ===阻尼参数大小不同时的振动非线性性质=== | ||
当阻尼系数 [math]\displaystyle{ \mu\ge 0 }[/math] 时，系统呈现稳定或趋于稳定的振动。系统的动力学特性随 [math]\displaystyle{ \mu }[/math] 的大小不同而呈现三种典型情形，分别对应模型非线性的不同强度，且随着 [math]\displaystyle{ \mu }[/math] 的增加，极限环的非线性动力学特征也在不断增强： | 当阻尼系数 [math]\displaystyle{ \mu\ge 0 }[/math] 时，系统呈现稳定或趋于稳定的振动。系统的动力学特性随 [math]\displaystyle{ \mu }[/math] 的大小不同而呈现三种典型情形，分别对应模型非线性的不同强度，且随着 [math]\displaystyle{ \mu }[/math] 的增加，极限环的非线性动力学特征也在不断增强： | ||
| + | |||
（1）当 [math]\displaystyle{ \mu=0 }[/math] 时：方程退化为简谐振荡器 [math]\displaystyle{ \ddot{x}+x=0 }[/math]，系统做等幅正弦振荡。作为一个保守系统，没有能量[[耗散]]也没有能量注入。 | （1）当 [math]\displaystyle{ \mu=0 }[/math] 时：方程退化为简谐振荡器 [math]\displaystyle{ \ddot{x}+x=0 }[/math]，系统做等幅正弦振荡。作为一个保守系统，没有能量[[耗散]]也没有能量注入。 | ||
| + | |||
（2）当 [math]\displaystyle{ 0<\mu\ll 1 }[/math] 时：非线性较弱，系统表现出接近正弦波的振荡，但振幅会缓慢调整，最终稳定在一个极限环上。此时极限环的形状接近圆形，振荡频率接近线性固有频率 [math]\displaystyle{ \omega\approx 1 }[/math]，收敛速度正比于 [math]\displaystyle{ \mu }[/math]。这一 regime 可用弱非线性摄动方法（如[[平均法]]、多尺度法）进行解析近似。 | （2）当 [math]\displaystyle{ 0<\mu\ll 1 }[/math] 时：非线性较弱，系统表现出接近正弦波的振荡，但振幅会缓慢调整，最终稳定在一个极限环上。此时极限环的形状接近圆形，振荡频率接近线性固有频率 [math]\displaystyle{ \omega\approx 1 }[/math]，收敛速度正比于 [math]\displaystyle{ \mu }[/math]。这一 regime 可用弱非线性摄动方法（如[[平均法]]、多尺度法）进行解析近似。 | ||
| + | |||
（3）当 [math]\displaystyle{ \mu\gg 1 }[/math] 时：非线性很强，系统表现出[[松弛振荡]]（Relaxation Oscillation）。波形不再是正弦波，而是呈现”缓慢充电、快速放电”的锯齿状或脉冲状特征：系统在 [math]\displaystyle{ |x|<1 }[/math] 区域内缓慢积累能量（准静态充电），一旦越过 [math]\displaystyle{ |x|=1 }[/math] 的阈值，非线性阻尼迅速转为正阻尼，系统在极短时间内释放能量（快速放电），随后又进入下一轮充电-放电循环。最终系统也将收敛于一个极限环，但该极限环在相空间中呈现明显的”圆角矩形”形状，左右两侧近似垂直的跳跃段对应于速度的瞬时剧变。这在生物[[神经元]]的[[动作电位]]模型中非常常见。 | （3）当 [math]\displaystyle{ \mu\gg 1 }[/math] 时：非线性很强，系统表现出[[松弛振荡]]（Relaxation Oscillation）。波形不再是正弦波，而是呈现”缓慢充电、快速放电”的锯齿状或脉冲状特征：系统在 [math]\displaystyle{ |x|<1 }[/math] 区域内缓慢积累能量（准静态充电），一旦越过 [math]\displaystyle{ |x|=1 }[/math] 的阈值，非线性阻尼迅速转为正阻尼，系统在极短时间内释放能量（快速放电），随后又进入下一轮充电-放电循环。最终系统也将收敛于一个极限环，但该极限环在相空间中呈现明显的”圆角矩形”形状，左右两侧近似垂直的跳跃段对应于速度的瞬时剧变。这在生物[[神经元]]的[[动作电位]]模型中非常常见。 | ||
| + | |||
[[文件:阻尼系数不同时的范德珀尔振子位移随时间变化图.png|缩略图|阻尼系数不同的范德珀尔振子位移随时间变化图]] | [[文件:阻尼系数不同时的范德珀尔振子位移随时间变化图.png|缩略图|阻尼系数不同的范德珀尔振子位移随时间变化图]] | ||
[[文件:阻尼系数不同的范德珀尔振子机械能随时间变化图.png|缩略图|阻尼系数不同的范德珀尔振子机械能随时间变化图]] | [[文件:阻尼系数不同的范德珀尔振子机械能随时间变化图.png|缩略图|阻尼系数不同的范德珀尔振子机械能随时间变化图]] | ||
| 第49行： | 第69行： | ||
==模型的推导过程== | ==模型的推导过程== | ||
我们可通过如下步骤来探究<ref name="UConn_vanderpol_oscillator" />：考虑如图所示的RLC电路：它由电感 [math]\displaystyle{ L }[/math]、电容器 [math]\displaystyle{ C }[/math]、电阻 [math]\displaystyle{ R }[/math] 组成，我们假设电压源是电池（[math]\displaystyle{ E(t)=E }[/math]，常数）。当电池开关关闭时，根据[[基尔霍夫电压定律]]： | 我们可通过如下步骤来探究<ref name="UConn_vanderpol_oscillator" />：考虑如图所示的RLC电路：它由电感 [math]\displaystyle{ L }[/math]、电容器 [math]\displaystyle{ C }[/math]、电阻 [math]\displaystyle{ R }[/math] 组成，我们假设电压源是电池（[math]\displaystyle{ E(t)=E }[/math]，常数）。当电池开关关闭时，根据[[基尔霍夫电压定律]]： | ||
| + | |||
[math]\displaystyle{ L\dot{I}+RI+Q/C=E }[/math] | [math]\displaystyle{ L\dot{I}+RI+Q/C=E }[/math] | ||
| + | |||
电路中开始流动电流 [math]\displaystyle{ I(t) }[/math]，其中 [math]\displaystyle{ Q }[/math] 作为电容器中的电荷，其一阶导数等于 [math]\displaystyle{ I }[/math]。如果我们对等式两边求导，我们将得到式： | 电路中开始流动电流 [math]\displaystyle{ I(t) }[/math]，其中 [math]\displaystyle{ Q }[/math] 作为电容器中的电荷，其一阶导数等于 [math]\displaystyle{ I }[/math]。如果我们对等式两边求导，我们将得到式： | ||
| + | |||
[math]\displaystyle{ L\frac{d^2I}{dt^2}+R\frac{dI}{dt}+\frac{I}{LC}=0 }[/math] | [math]\displaystyle{ L\frac{d^2I}{dt^2}+R\frac{dI}{dt}+\frac{I}{LC}=0 }[/math] | ||
| + | |||
这是一个具有恒定系数的二阶线性常[[微分方程]]，它表示一个[[阻尼谐波振荡器]]。 | 这是一个具有恒定系数的二阶线性常[[微分方程]]，它表示一个[[阻尼谐波振荡器]]。 | ||
然而，范德珀尔考虑的电路显示的是有源元件（真空管阵列-半导体），而不是无源电阻。通过这种方式，当电流低时，半导体就像在系统中泵送能量一样，当电流过高时，半导体会抑制系统的能量。能量注入和能量吸收之间的相互作用导致电压和电流的周期性振荡。半导体的作用由函数 [math]\displaystyle{ I^2-\alpha }[/math] 所刻画，其中 [math]\displaystyle{ \alpha }[/math] 是电流的阈值水平，电流的方程变成下式： | 然而，范德珀尔考虑的电路显示的是有源元件（真空管阵列-半导体），而不是无源电阻。通过这种方式，当电流低时，半导体就像在系统中泵送能量一样，当电流过高时，半导体会抑制系统的能量。能量注入和能量吸收之间的相互作用导致电压和电流的周期性振荡。半导体的作用由函数 [math]\displaystyle{ I^2-\alpha }[/math] 所刻画，其中 [math]\displaystyle{ \alpha }[/math] 是电流的阈值水平，电流的方程变成下式： | ||
| + | |||
[math]\displaystyle{ L\dot{I}+(I^2-\alpha)I+Q/C=E }[/math] | [math]\displaystyle{ L\dot{I}+(I^2-\alpha)I+Q/C=E }[/math] | ||
| + | |||
为了简单起见，我们假设没有外部电压源。通过取方程两侧的导数，我们确实得到了范德珀尔振子的方程： | 为了简单起见，我们假设没有外部电压源。通过取方程两侧的导数，我们确实得到了范德珀尔振子的方程： | ||
| + | |||
[math]\displaystyle{ \ddot{I}-\frac{\alpha}{L}\left(1-\frac{3I^2}{\alpha}\right)\dot{I}+\frac{I}{LC}=0 }[/math] | [math]\displaystyle{ \ddot{I}-\frac{\alpha}{L}\left(1-\frac{3I^2}{\alpha}\right)\dot{I}+\frac{I}{LC}=0 }[/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] | ||
| + | |||
其中 [math]\displaystyle{ F }[/math] 是正弦驱动函数的振幅，[math]\displaystyle{ \omega }[/math] 是其角频率。此系统经历一系列分岔（Bifurcation）过程，这个系统从稳定的极限环走向混沌。 | 其中 [math]\displaystyle{ F }[/math] 是正弦驱动函数的振幅，[math]\displaystyle{ \omega }[/math] 是其角频率。此系统经历一系列分岔（Bifurcation）过程，这个系统从稳定的极限环走向混沌。 | ||
（注：受迫范德珀尔振子因引入时间维度而升维至三维及以上，打破了[[庞加莱-本迪克松定理]]的限制，从而触及[[混沌理论]]的层次。） | （注：受迫范德珀尔振子因引入时间维度而升维至三维及以上，打破了[[庞加莱-本迪克松定理]]的限制，从而触及[[混沌理论]]的层次。） | ||
| 第71行： | 第101行： | ||
第二，霍普夫分岔理论延续了范德珀尔对极限环稳定性的分析思路。 在范德珀尔的工作中，极限环的稳定性由非线性阻尼的能量调节机制保证：小振幅时注入能量、大振幅时耗散能量。在霍普夫分岔的一般理论中，这一物理图像被抽象为[[第一李雅普诺夫系数]]（first Lyapunov coefficient）的符号判定：当该系数为负时，分岔是超临界的，产生的极限环稳定；当该系数为正时，分岔是亚临界的，产生的极限环不稳定。对于范德珀尔振子，可以直接计算其[[第一李雅普诺夫系数]]为负，从而从一般理论的角度确认了范德珀尔通过物理分析得出的[[超临界分岔]]结论<ref name="Strogatz_Example_8.4.1" /><ref name="Perko_2013_Differential_Equations_and_Dynamical_Systems" />。 | 第二，霍普夫分岔理论延续了范德珀尔对极限环稳定性的分析思路。 在范德珀尔的工作中，极限环的稳定性由非线性阻尼的能量调节机制保证：小振幅时注入能量、大振幅时耗散能量。在霍普夫分岔的一般理论中，这一物理图像被抽象为[[第一李雅普诺夫系数]]（first Lyapunov coefficient）的符号判定：当该系数为负时，分岔是超临界的，产生的极限环稳定；当该系数为正时，分岔是亚临界的，产生的极限环不稳定。对于范德珀尔振子，可以直接计算其[[第一李雅普诺夫系数]]为负，从而从一般理论的角度确认了范德珀尔通过物理分析得出的[[超临界分岔]]结论<ref name="Strogatz_Example_8.4.1" /><ref name="Perko_2013_Differential_Equations_and_Dynamical_Systems" />。 | ||
第三，通过霍普夫分岔的[[正常型]]分析，可以清晰地看到范德珀尔振子的分岔本质。 具体而言，令 [math]\displaystyle{ \varepsilon=\mu }[/math]，并作变量变换 [math]\displaystyle{ u=\sqrt{\mu},x }[/math]，将范德珀尔方程转换为： | 第三，通过霍普夫分岔的[[正常型]]分析，可以清晰地看到范德珀尔振子的分岔本质。 具体而言，令 [math]\displaystyle{ \varepsilon=\mu }[/math]，并作变量变换 [math]\displaystyle{ u=\sqrt{\mu},x }[/math]，将范德珀尔方程转换为： | ||
| + | |||
[math]\displaystyle{ \ddot{u}+u+u^2\dot{u}-\varepsilon\dot{u}=0 }[/math] | [math]\displaystyle{ \ddot{u}+u+u^2\dot{u}-\varepsilon\dot{u}=0 }[/math] | ||
| + | |||
此式即为霍普夫分岔的标准[[正常型]]之一<ref name="Strogatz_Example_8.4.1" />。在这一变换后的变量 [math]\displaystyle{ u }[/math] 中，极限环的振幅正比于 [math]\displaystyle{ \sqrt{\varepsilon}=\sqrt{\mu} }[/math]，即当 [math]\displaystyle{ \mu }[/math] 从零增大时，极限环从原点”生长”出来，振幅按 [math]\displaystyle{ \sqrt{\mu} }[/math] 的规律增长——这正是超临界霍普夫分岔的标志性特征。值得注意的是，回到原始变量 [math]\displaystyle{ x=u/\sqrt{\mu} }[/math] 后，极限环的振幅为 [math]\displaystyle{ O(1) }[/math] 量级（小 [math]\displaystyle{ \mu }[/math] 时近似为 [math]\displaystyle{ 2 }[/math]），这说明范德珀尔振子的极限环在原始尺度上并非”小振幅”的，但其分岔结构（从平衡点”分叉”出极限环）在正常型尺度下完全符合霍普夫分岔的通用图像。 | 此式即为霍普夫分岔的标准[[正常型]]之一<ref name="Strogatz_Example_8.4.1" />。在这一变换后的变量 [math]\displaystyle{ u }[/math] 中，极限环的振幅正比于 [math]\displaystyle{ \sqrt{\varepsilon}=\sqrt{\mu} }[/math]，即当 [math]\displaystyle{ \mu }[/math] 从零增大时，极限环从原点”生长”出来，振幅按 [math]\displaystyle{ \sqrt{\mu} }[/math] 的规律增长——这正是超临界霍普夫分岔的标志性特征。值得注意的是，回到原始变量 [math]\displaystyle{ x=u/\sqrt{\mu} }[/math] 后，极限环的振幅为 [math]\displaystyle{ O(1) }[/math] 量级（小 [math]\displaystyle{ \mu }[/math] 时近似为 [math]\displaystyle{ 2 }[/math]），这说明范德珀尔振子的极限环在原始尺度上并非”小振幅”的，但其分岔结构（从平衡点”分叉”出极限环）在正常型尺度下完全符合霍普夫分岔的通用图像。 | ||
因此，霍普夫分岔理论不仅从一般性上解释了范德珀尔振子为何会在 [math]\displaystyle{ \mu=0 }[/math] 处产生极限环，还提供了定量工具（正常型变换、李雅普诺夫系数）来精确刻画分岔的性态。从霍普夫分岔的视角回望范德珀尔振子，我们可以看到：范德珀尔通过物理直觉和电路分析所发现的极限环涌现机制，本质上就是一个超临界霍普夫分岔过程；而霍普夫定理则将这一具体发现提升为适用于广泛动力系统的一般理论，使范德珀尔振子从”一个特殊的电路模型”升格为”理解分岔理论的教科书级原型”。 | 因此，霍普夫分岔理论不仅从一般性上解释了范德珀尔振子为何会在 [math]\displaystyle{ \mu=0 }[/math] 处产生极限环，还提供了定量工具（正常型变换、李雅普诺夫系数）来精确刻画分岔的性态。从霍普夫分岔的视角回望范德珀尔振子，我们可以看到：范德珀尔通过物理直觉和电路分析所发现的极限环涌现机制，本质上就是一个超临界霍普夫分岔过程；而霍普夫定理则将这一具体发现提升为适用于广泛动力系统的一般理论，使范德珀尔振子从”一个特殊的电路模型”升格为”理解分岔理论的教科书级原型”。 | ||
==量子物理领域——范德珀尔振子的[[哈密顿量]]== | ==量子物理领域——范德珀尔振子的[[哈密顿量]]== | ||
人们还可以为范德珀尔振荡器编写一个与时间无关的哈密顿形式，方法是使用辅助二阶非线性微分方程将其增强为四维自主[[动力系统]]，如下所示： | 人们还可以为范德珀尔振荡器编写一个与时间无关的哈密顿形式，方法是使用辅助二阶非线性微分方程将其增强为四维自主[[动力系统]]，如下所示： | ||
| + | |||
[math]\displaystyle{ \ddot{x}-\mu(1-x^2)\dot{x}+x=0 }[/math] | [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] | [math]\displaystyle{ \ddot{y}+\mu(1-x^2)\dot{y}+y=0 }[/math] | ||
| + | |||
请注意，由于 [math]\displaystyle{ x }[/math] 和 [math]\displaystyle{ y }[/math] 变量的时间演变之间的单向耦合，原始范德珀尔振荡器的动力学不受影响。可以通过这个方程组的[[哈密顿函数]] [math]\displaystyle{ H }[/math] 证明<ref name="Shah_Conservative_perturbation_theory_for_nonconservative_systems" />： | 请注意，由于 [math]\displaystyle{ x }[/math] 和 [math]\displaystyle{ y }[/math] 变量的时间演变之间的单向耦合，原始范德珀尔振荡器的动力学不受影响。可以通过这个方程组的[[哈密顿函数]] [math]\displaystyle{ H }[/math] 证明<ref name="Shah_Conservative_perturbation_theory_for_nonconservative_systems" />： | ||
| + | |||
[math]\displaystyle{ H(x,y,p_x,p_y)=p_xp_y+xy-\mu(1-x^2)y p_y }[/math] | [math]\displaystyle{ H(x,y,p_x,p_y)=p_xp_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] 分别是对应于 [math]\displaystyle{ x }[/math] 和 [math]\displaystyle{ y }[/math] 的共轭矩。原则上，这可能会导致范德珀尔振荡器的[[量子化]]，这样的[[哈密顿量]]也将具有时间依赖参数的极限周期系统的[[几何相位]]与相应的[[哈密顿量]]系统的[[汉奈角]]连接起来<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] 分别是对应于 [math]\displaystyle{ x }[/math] 和 [math]\displaystyle{ y }[/math] 的共轭矩。原则上，这可能会导致范德珀尔振荡器的[[量子化]]，这样的[[哈密顿量]]也将具有时间依赖参数的极限周期系统的[[几何相位]]与相应的[[哈密顿量]]系统的[[汉奈角]]连接起来<ref name="Chattopadhyay_Finding_the_Hannay_angle" />。 | ||
===关于量子振荡器=== | ===关于量子振荡器=== | ||
| 第97行： | 第134行： | ||
==催生混沌理论的诞生== | ==催生混沌理论的诞生== | ||
范德珀尔振子及受迫范德珀尔振子呈现的分岔性质，为最初被人类观测到的混沌现象之一，叩开了混沌世界的大门。1945年英国数学家卡特赖特（M.Cartwright）和李特伍德（J.Littlewood）在对受迫范德珀尔振子的研究中，在理论层面首次证明受迫范德珀尔方程存在[[奇异解]]，预示了混沌的可能性<ref name="Cartwright_Littlewood_On_non-linear_differential_equations_of_the_second_order" />。随后，日本学者上田睆亮（Ueda Yoshisuke）于1961年首次在[[模拟计算机]]上观察到混沌[[吸引子]]<ref name="Hayashi_Shibayama_Ueda_Quasi-periodic_Oscillations_in_Self-oscillatory_Systems_with_External_Force" />，随着后续研究的深入，混沌理论的面纱也逐渐被科学家们所揭晓。 | 范德珀尔振子及受迫范德珀尔振子呈现的分岔性质，为最初被人类观测到的混沌现象之一，叩开了混沌世界的大门。1945年英国数学家卡特赖特（M.Cartwright）和李特伍德（J.Littlewood）在对受迫范德珀尔振子的研究中，在理论层面首次证明受迫范德珀尔方程存在[[奇异解]]，预示了混沌的可能性<ref name="Cartwright_Littlewood_On_non-linear_differential_equations_of_the_second_order" />。随后，日本学者上田睆亮（Ueda Yoshisuke）于1961年首次在[[模拟计算机]]上观察到混沌[[吸引子]]<ref name="Hayashi_Shibayama_Ueda_Quasi-periodic_Oscillations_in_Self-oscillatory_Systems_with_External_Force" />，随着后续研究的深入，混沌理论的面纱也逐渐被科学家们所揭晓。 | ||
| − | + | =参考文献= | |
<ref name="Gerard-Lescuyer_Sur_un_paradoxe_electrodynamique">Gerard-Lescuyer, J. M. A. (1880). Sur un paradoxe électrodynamique.</ref> | <ref name="Gerard-Lescuyer_Sur_un_paradoxe_electrodynamique">Gerard-Lescuyer, J. M. A. (1880). Sur un paradoxe électrodynamique.</ref> | ||
<ref name="Nature_Musical_Arcs">Musical Arcs. (1901). Nature, 63(1640), 542–543.</ref> | <ref name="Nature_Musical_Arcs">Musical Arcs. (1901). Nature, 63(1640), 542–543.</ref> | ||
2026年10月22日 (四) 22:46的版本
范德珀尔振子是非线性动力学中的一个典型自激振荡模型，由荷兰物理学家在研究真空电子管振荡电路时提出，核心数学模型为一个二阶非线性阻尼振子模型。因阻尼系数不同，产生振动的不同动力学特性；同时，因振幅不同，阻尼的正负性、敛散性也出现差异，即出现分岔现象。 在动力学理论中，该模型可进一步延伸为非齐次的受迫范德珀尔振子及量子振荡器的模型。该模型作为典型的非线性动力学模型，也在机械振动、电子电器、生物医学等领域中广泛应用；此外，对受迫范德珀尔振子的深入研究也促进了混沌理论的诞生和发展。
概述及历史背景
范德珀尔振子（Van der Pol oscillator）是非线性动力学中最经典、最具代表性的自激振荡模型之一。它通过一个简明的二阶非线性微分方程，特别是其中关键的振幅相关非线性阻尼项，就能够展示线性振子所没有的稳定极限环、振幅不同而产生的分岔及其能量流向、自激振荡以及随非线性增强而出现的松弛震荡（亦称弛豫振荡）等核心非线性现象，因此常被视为理解极限环理论和非线性振荡的”原型系统”。 这一模型的诞生有着深刻的时代背景。19世纪末20世纪初，随着动力学理论与电子电气技术的持续发展，1880年Gérard-Lescuyer的串联电机[1]、1901年Duddell的音乐电弧[2]以及1907年de Forest发明的真空三极管[3]等装置，都展示了”自激振荡”现象。然而，面对这些自激振动现象，传统线性振动理论无法解释极限环及分岔现象，传统摄动分析方法也面临级数发散问题，始终缺失解释不同自激振动现象的统一理论体系。 在此背景下，时任飞利浦公司电气工程师的荷兰物理学家巴尔塔萨·范德珀尔（Balthasar van der Pol）于1920年采用非线性动力学分析方法，引入非线性阻尼概念，提出了这一二阶非线性模型[4]。该模型最初用于描述早期无线电设备中三极管电路的自持振荡现象，是首个被深入研究的”自激振荡”非线性模型。范德珀尔在研究该系统时开创性地提出了松弛振荡（张弛振荡，Relaxation Oscillation）这一概念，并系统阐述了此现象，同时开创了”非线性振荡”的系统性研究方法。他与Van der Mark于1927年在《Nature》上发表了关于”不规则噪声”的报告[5]，后被公认为混沌现象最早的实验观察之一，间接催生了后续混沌理论的诞生。 后来人们发现，这种机制并不局限于电子电路，而广泛出现在生物节律、神经活动、心脏起搏、机械振动等系统中，因此范德珀尔振子也成为从具体物理问题走向一般非线性动力学理论的一个标志性模型。 作为一个非保守非线性振荡器，范德珀尔振子的基本模型为一个二阶非线性微分方程：
[math]\displaystyle{ \ddot{x}-\mu(1-x^2)\dot{x}+x=0,\quad \mu>0 }[/math]
其中 [math]\displaystyle{ \mu }[/math] 为非线性阻尼系数，中间的非线性阻尼项 [math]\displaystyle{ -\mu(1-x^2)\dot{x} }[/math] 决定了系统的全部动力学特性，该模型将线性理论、分岔理论、极限环、快慢动力学乃至混沌的边缘整合了起来，实现了由传统线性振荡理论向非线性动力学、乃至于混沌理论的飞跃。
主要内容
数学模型
范德珀尔振子（Van der Pol oscillator）是研究极限环非线性系统的经典实例。作为一个耗散性非线性振荡系统，其数学模型为一个二阶非线性齐次方程：
[math]\displaystyle{ \ddot{x}-\mu(1-x^2)\dot{x}+x=0,\quad \mu>0 }[/math]
其中 [math]\displaystyle{ x }[/math] 为位置坐标（作为时间 [math]\displaystyle{ t }[/math] 的函数），[math]\displaystyle{ \mu }[/math] 是标量参数，二者决定阻尼的非线性和强度。 也可写为二维状态空间方程组的形式：
[math]\displaystyle{ \dot{x}=y,\quad \dot{y}=\mu(1-x^2)y-x }[/math]
其中：
（1）[math]\displaystyle{ x }[/math] 是位移（或电压等物理量）。
（2）[math]\displaystyle{ y=\dot{x} }[/math] 是速度，[math]\displaystyle{ \dot{y}=\ddot{x} }[/math] 是加速度。
（3）[math]\displaystyle{ \mu }[/math] 是非线性阻尼参数（[math]\displaystyle{ \mu\ge 0 }[/math]），控制着非线性效应的强弱。
当方程左侧为0的时候，模型为自由范德珀尔振子，不受外界噪声干扰；当方程右侧不为0时，模型变为受迫范德珀尔振子，可进一步展示分岔、混沌等更复杂的非线性现象。
物理意义与动力学行为
范德珀尔振子的动力学性质由三个相互关联的因素共同决定：阻尼系数的符号（[math]\displaystyle{ \mu }[/math] 的正负）、阻尼系数的大小（[math]\displaystyle{ |\mu| }[/math] 的强弱）以及振子位移的幅度（[math]\displaystyle{ |x| }[/math] 与 [math]\displaystyle{ 1 }[/math] 的相对大小）。阻尼系数的符号决定了系统是否发生分岔、是否存在极限环——当 [math]\displaystyle{ \mu }[/math] 从负值穿过零点变为正值时，系统经历霍普夫分岔，平凡平衡点由稳定变为不稳定，并产生稳定的极限环；阻尼系数的大小决定了极限环的形状和振动的非线性程度——[math]\displaystyle{ \mu }[/math] 较小时系统接近简谐振动，[math]\displaystyle{ \mu }[/math] 较大时呈现弛豫振荡；而位移幅度则决定了在一个振动周期内能量的实时分配——[math]\displaystyle{ |x|<1 }[/math] 时系统注入能量，[math]\displaystyle{ |x|>1 }[/math] 时系统耗散能量。三者相互耦合，共同塑造了范德珀尔振子从线性振荡到弛豫振荡、从稳定极限环到混沌的丰富动力学行为。下文将分别从这三个维度展开论述。
阻尼参数正负不同时的分岔性质
范德珀尔振子作为一个非线性振荡系统，其动力学特性与阻尼系数 [math]\displaystyle{ \mu }[/math] 的符号密切相关。将方程线性化于原点附近（即忽略 [math]\displaystyle{ x^2 }[/math] 项），可得线性近似方程 [math]\displaystyle{ \ddot{x}-\mu\dot{x}+x=0 }[/math] 其特征方程为 [math]\displaystyle{ \lambda^2-\mu\lambda+1=0 }[/math]，特征根为 [math]\displaystyle{ \lambda=(\mu\pm\sqrt{\mu^2-4})/2 }[/math]。由此可以看出，原点处的稳定性完全由 [math]\displaystyle{ \mu }[/math] 的符号决定：
当 [math]\displaystyle{ \mu<0 }[/math] 时，特征根的实部为负，原点是一个稳定的焦点（螺旋汇），相空间中所有轨迹均螺旋收敛于原点，系统最终静止，不存在极限环。此时若存在极限环，则必为不稳定的（排斥子），任何微小的扰动都会使系统远离该极限环。
当 [math]\displaystyle{ \mu=0 }[/math] 时，特征根为纯虚数 [math]\displaystyle{ \pm i }[/math]，原点退化为线性中心，系统退化为保守的简谐振荡器 [math]\displaystyle{ \ddot{x}+x=0 }[/math]，相轨迹为一族同心圆，既不发散也不收敛。[math]\displaystyle{ \mu=0 }[/math] 因此成为系统动力学性质发生突变的临界点——霍普夫分岔点。
当 [math]\displaystyle{ \mu>0 }[/math] 时，特征根的实部为正，原点变为不稳定的焦点（螺旋源），任何从原点附近出发的轨迹都会被”推开”，螺旋向外运动。然而，由于非线性阻尼项的存在，轨迹并不会无限发散，而是最终被”捕获”到一个稳定的闭合轨道——即稳定的极限环上。此时系统呈现自激振荡，振幅由非线性机制自动锁定。
上述 [math]\displaystyle{ \mu }[/math] 穿过零点时原点稳定性翻转、极限环随之产生或消失的现象，正是超临界霍普夫分岔的典型图景。值得注意的是，[math]\displaystyle{ \mu>0 }[/math] 与 [math]\displaystyle{ \mu<0 }[/math] 两种情形在时间反演下互为镜像：将 [math]\displaystyle{ t\to -t }[/math] 代入方程，[math]\displaystyle{ \mu }[/math] 的符号恰好翻转，因此一者为另一者的逆过程，体现了系统在时空意义上的对称性。从相图上看，当 [math]\displaystyle{ \mu>0 }[/math] 时，轨迹随机从平面内各点螺旋逼近极限环的位置（类似逆时针书写”6”字形，从“尾巴”出发，最后在下面的小圈处反复“画圈”），最终被极限环”捕获”并形成稳定循环；当 [math]\displaystyle{ \mu<0 }[/math] 时则相反，轨迹从极限环螺旋向内或向外随机收敛于平面内各点（类似顺时针书写“6”字形，从下面“小圈”左上角和“尾巴”衔接的地方出发，回到这里，再画出一条长长的“尾巴”）。分岔在此扮演了”入口”的角色——没有分岔，原点便是稳定的，所有轨迹终将被吸入口径为零的平衡点，极限环与后续的自持振荡便无从谈起。
阻尼参数大小不同时的振动非线性性质
当阻尼系数 [math]\displaystyle{ \mu\ge 0 }[/math] 时，系统呈现稳定或趋于稳定的振动。系统的动力学特性随 [math]\displaystyle{ \mu }[/math] 的大小不同而呈现三种典型情形，分别对应模型非线性的不同强度，且随着 [math]\displaystyle{ \mu }[/math] 的增加，极限环的非线性动力学特征也在不断增强：
（1）当 [math]\displaystyle{ \mu=0 }[/math] 时：方程退化为简谐振荡器 [math]\displaystyle{ \ddot{x}+x=0 }[/math]，系统做等幅正弦振荡。作为一个保守系统，没有能量耗散也没有能量注入。
（2）当 [math]\displaystyle{ 0<\mu\ll 1 }[/math] 时：非线性较弱，系统表现出接近正弦波的振荡，但振幅会缓慢调整，最终稳定在一个极限环上。此时极限环的形状接近圆形，振荡频率接近线性固有频率 [math]\displaystyle{ \omega\approx 1 }[/math]，收敛速度正比于 [math]\displaystyle{ \mu }[/math]。这一 regime 可用弱非线性摄动方法（如平均法、多尺度法）进行解析近似。
（3）当 [math]\displaystyle{ \mu\gg 1 }[/math] 时：非线性很强，系统表现出松弛振荡（Relaxation Oscillation）。波形不再是正弦波，而是呈现”缓慢充电、快速放电”的锯齿状或脉冲状特征：系统在 [math]\displaystyle{ |x|<1 }[/math] 区域内缓慢积累能量（准静态充电），一旦越过 [math]\displaystyle{ |x|=1 }[/math] 的阈值，非线性阻尼迅速转为正阻尼，系统在极短时间内释放能量（快速放电），随后又进入下一轮充电-放电循环。最终系统也将收敛于一个极限环，但该极限环在相空间中呈现明显的”圆角矩形”形状，左右两侧近似垂直的跳跃段对应于速度的瞬时剧变。这在生物神经元的动作电位模型中非常常见。
位移不同时的能量传递
在阻尼系数 [math]\displaystyle{ \mu }[/math] 为正且固定的条件下，系统的能量调节机制完全取决于位移 [math]\displaystyle{ x }[/math] 的大小，这一机制在一个振动周期内呈现出精细的周期性波动，具体表现为非线性阻尼项 [math]\displaystyle{ -\mu(1-x^2)\dot{x} }[/math] 的符号随位移幅值切换： 当位移较小（[math]\displaystyle{ |x|<1 }[/math]）时，[math]\displaystyle{ (1-x^2)>0 }[/math]，非线性阻尼项表现为负阻尼，系统向振荡注入能量，使小振幅运动增长。在相空间中，这对应于轨迹在极限环内侧的区域运动时，阻尼力与速度方向相同，对振子做正功，不断向系统输入能量，推动轨迹向外螺旋扩张。 当位移较大（[math]\displaystyle{ |x|>1 }[/math]）时，[math]\displaystyle{ (1-x^2)<0 }[/math]，非线性阻尼项表现为正阻尼，系统耗散能量，使大振幅运动衰减。在相空间中，这对应于轨迹在极限环外侧的区域运动时，阻尼力与速度方向相反，对振子做负功，从系统中抽取能量，迫使轨迹向内螺旋收缩。 在一个完整的振动周期中，振子位移大小随时间周期性变化，因此能量注入与耗散的机制也在周期性地交替：轨迹从极限环内侧出发时，[math]\displaystyle{ |x|<1 }[/math]，负阻尼注入能量，振幅被推大；随着振幅增大越过 [math]\displaystyle{ |x|=1 }[/math] 的阈值，正阻尼开始耗散能量，振幅被压小；当轨迹回到极限环内侧时，新一轮的注入-耗散循环又开始。正是这种”小振幅吸能、大振幅放能”的自调节机制，使得一个周期内注入的能量恰好等于耗散的能量，振幅被锁定在极限环上，系统实现稳定的自持振荡。 从能量守恒的角度看，极限环对应于相空间中一条闭合曲线，沿此曲线一周，非线性阻尼所做的总功为零——注入功与耗散功精确平衡。这一条件唯一确定了极限环的振幅：对于小 [math]\displaystyle{ \mu }[/math] 情形，极限环近似为半径为 [math]\displaystyle{ 2 }[/math] 的圆（即振幅 [math]\displaystyle{ A\approx 2 }[/math]）[6]；对于大 [math]\displaystyle{ \mu }[/math] 情形，极限环形状偏离圆形，但”一周期内净能量输入为零”的平衡条件依然成立，只是振幅的解析表达式需借助数值方法获得。
上述能量调节机制是范德珀尔振子能够维持稳定自激振荡的核心物理机制，也是解释霍普夫分岔产生极限环现象的关键——分岔决定了极限环的”诞生”（原点失稳、极限环从无到有），而能量调节机制则决定了极限环的”维持”（振幅被锁定、振荡持续不断）。二者共同构成了范德珀尔振子作为”自激振荡原型系统”的完整物理图像。
模型的推导过程
我们可通过如下步骤来探究[7]：考虑如图所示的RLC电路：它由电感 [math]\displaystyle{ L }[/math]、电容器 [math]\displaystyle{ C }[/math]、电阻 [math]\displaystyle{ R }[/math] 组成，我们假设电压源是电池（[math]\displaystyle{ E(t)=E }[/math]，常数）。当电池开关关闭时，根据基尔霍夫电压定律：
[math]\displaystyle{ L\dot{I}+RI+Q/C=E }[/math]
电路中开始流动电流 [math]\displaystyle{ I(t) }[/math]，其中 [math]\displaystyle{ Q }[/math] 作为电容器中的电荷，其一阶导数等于 [math]\displaystyle{ I }[/math]。如果我们对等式两边求导，我们将得到式：
[math]\displaystyle{ L\frac{d^2I}{dt^2}+R\frac{dI}{dt}+\frac{I}{LC}=0 }[/math]
这是一个具有恒定系数的二阶线性常微分方程，它表示一个阻尼谐波振荡器。 然而，范德珀尔考虑的电路显示的是有源元件（真空管阵列-半导体），而不是无源电阻。通过这种方式，当电流低时，半导体就像在系统中泵送能量一样，当电流过高时，半导体会抑制系统的能量。能量注入和能量吸收之间的相互作用导致电压和电流的周期性振荡。半导体的作用由函数 [math]\displaystyle{ I^2-\alpha }[/math] 所刻画，其中 [math]\displaystyle{ \alpha }[/math] 是电流的阈值水平，电流的方程变成下式：
[math]\displaystyle{ L\dot{I}+(I^2-\alpha)I+Q/C=E }[/math]
为了简单起见，我们假设没有外部电压源。通过取方程两侧的导数，我们确实得到了范德珀尔振子的方程：
[math]\displaystyle{ \ddot{I}-\frac{\alpha}{L}\left(1-\frac{3I^2}{\alpha}\right)\dot{I}+\frac{I}{LC}=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]
其中 [math]\displaystyle{ F }[/math] 是正弦驱动函数的振幅，[math]\displaystyle{ \omega }[/math] 是其角频率。此系统经历一系列分岔（Bifurcation）过程，这个系统从稳定的极限环走向混沌。 （注：受迫范德珀尔振子因引入时间维度而升维至三维及以上，打破了庞加莱-本迪克松定理的限制，从而触及混沌理论的层次。）
霍普夫分岔与范德珀尔振子的理论继承关系
范德珀尔在1920年代研究真空电子管振荡电路时，已经观察并系统描述了当阻尼参数改变符号时系统动力学行为发生的突变——稳定极限环的涌现、弛豫振荡的出现以及不同振荡模式之间的切换。这些工作实质上触及了现代分岔理论的核心概念，但范德珀尔的研究始终聚焦于具体的物理模型和电路实现，并未将其上升为一般的数学定理。 直到1942年，德国数学家赫尔曼·霍普夫（Hermann Hopf）在其论文中证明了如下一般性定理：对于一个依赖于参数的光滑动力系统，当平衡点的线性化矩阵有一对共轭复特征值的实部穿过虚轴（即从零变为正），且满足一定的非退化条件和横截条件时，必然在该平衡点附近产生一个极限环。1947年，安德罗诺夫（Andronov）与霍普夫合作发表了更完善的表述，后世将此结果称为安德罗诺夫-霍普夫分岔定理（Andronov-Hopf bifurcation theorem）[8][9]。 霍普夫定理与范德珀尔工作之间的继承关系体现在三个层面： 第一，范德珀尔提供了原型实例，霍普夫提供了通用框架。 范德珀尔振子是第一个被深入研究的、展示分岔产生极限环现象的非线性模型，它充当了霍普夫定理的”原型验证”——霍普夫定理将范德珀尔在特定电路中观察到的现象，推广到了任意满足线性化条件的有限维光滑动力系统。换言之，范德珀尔发现了”分岔产生极限环”这一现象，霍普夫则给出了”在什么条件下任何系统都会出现分岔产生极限环”的严格数学判据。 第二，霍普夫分岔理论延续了范德珀尔对极限环稳定性的分析思路。 在范德珀尔的工作中，极限环的稳定性由非线性阻尼的能量调节机制保证：小振幅时注入能量、大振幅时耗散能量。在霍普夫分岔的一般理论中，这一物理图像被抽象为第一李雅普诺夫系数（first Lyapunov coefficient）的符号判定：当该系数为负时，分岔是超临界的，产生的极限环稳定；当该系数为正时，分岔是亚临界的，产生的极限环不稳定。对于范德珀尔振子，可以直接计算其第一李雅普诺夫系数为负，从而从一般理论的角度确认了范德珀尔通过物理分析得出的超临界分岔结论[10][11]。 第三，通过霍普夫分岔的正常型分析，可以清晰地看到范德珀尔振子的分岔本质。 具体而言，令 [math]\displaystyle{ \varepsilon=\mu }[/math]，并作变量变换 [math]\displaystyle{ u=\sqrt{\mu},x }[/math]，将范德珀尔方程转换为：
[math]\displaystyle{ \ddot{u}+u+u^2\dot{u}-\varepsilon\dot{u}=0 }[/math]
此式即为霍普夫分岔的标准正常型之一[10]。在这一变换后的变量 [math]\displaystyle{ u }[/math] 中，极限环的振幅正比于 [math]\displaystyle{ \sqrt{\varepsilon}=\sqrt{\mu} }[/math]，即当 [math]\displaystyle{ \mu }[/math] 从零增大时，极限环从原点”生长”出来，振幅按 [math]\displaystyle{ \sqrt{\mu} }[/math] 的规律增长——这正是超临界霍普夫分岔的标志性特征。值得注意的是，回到原始变量 [math]\displaystyle{ x=u/\sqrt{\mu} }[/math] 后，极限环的振幅为 [math]\displaystyle{ O(1) }[/math] 量级（小 [math]\displaystyle{ \mu }[/math] 时近似为 [math]\displaystyle{ 2 }[/math]），这说明范德珀尔振子的极限环在原始尺度上并非”小振幅”的，但其分岔结构（从平衡点”分叉”出极限环）在正常型尺度下完全符合霍普夫分岔的通用图像。 因此，霍普夫分岔理论不仅从一般性上解释了范德珀尔振子为何会在 [math]\displaystyle{ \mu=0 }[/math] 处产生极限环，还提供了定量工具（正常型变换、李雅普诺夫系数）来精确刻画分岔的性态。从霍普夫分岔的视角回望范德珀尔振子，我们可以看到：范德珀尔通过物理直觉和电路分析所发现的极限环涌现机制，本质上就是一个超临界霍普夫分岔过程；而霍普夫定理则将这一具体发现提升为适用于广泛动力系统的一般理论，使范德珀尔振子从”一个特殊的电路模型”升格为”理解分岔理论的教科书级原型”。
量子物理领域——范德珀尔振子的哈密顿量
人们还可以为范德珀尔振荡器编写一个与时间无关的哈密顿形式，方法是使用辅助二阶非线性微分方程将其增强为四维自主动力系统，如下所示：
[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]
请注意，由于 [math]\displaystyle{ x }[/math] 和 [math]\displaystyle{ y }[/math] 变量的时间演变之间的单向耦合，原始范德珀尔振荡器的动力学不受影响。可以通过这个方程组的哈密顿函数 [math]\displaystyle{ H }[/math] 证明[12]：
[math]\displaystyle{ H(x,y,p_x,p_y)=p_xp_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] 分别是对应于 [math]\displaystyle{ x }[/math] 和 [math]\displaystyle{ y }[/math] 的共轭矩。原则上，这可能会导致范德珀尔振荡器的量子化，这样的哈密顿量也将具有时间依赖参数的极限周期系统的几何相位与相应的哈密顿量系统的汉奈角连接起来[13]。
关于量子振荡器
量子范德珀尔振荡器是经典范德珀尔振荡器的量子力学版本，使用林德布拉德方程来研究其量子动力学和量子同步[14]。请注意，上述具有辅助二阶方程的哈密顿方法产生无界相空间轨迹，因此无法构造物理上允许的束缚态，不能直接用于量子化范德珀尔振荡器。在弱非线性极限（即 [math]\displaystyle{ \mu\to 0 }[/math]）中，范德珀尔振荡器简化为Stuart-Landau方程。Stuart-Landau方程实际上描述了弱非线性极限中的一整类极限周期振荡器。经典Stuart-Landau方程的形式要简单得多，可以用Lindblad方程量化。 量子Stuart-Landau模型在量子同步研究中发挥了重要作用[15][16]。经典Stuart-Landau模型（[math]\displaystyle{ \mu\to 0 }[/math]）与更一般的极限周期振荡器（任意 [math]\displaystyle{ \mu }[/math]）之间的关系也在相应的量子模型中进行了数值证明。 然而，上述借助辅助变量构建哈密顿量的方法虽在数学上具有巧思，却止步于量子化的物理门槛——其产生的无界相空间轨迹无法构造物理上允许的束缚态。这恰恰揭示了经典非保守系统（如范德珀尔振荡器）量子化的本质困难：系统因非线性阻尼与外界环境存在持续的能量交换，必须将其视为一个开放量子系统。 要跨越这一鸿沟，理论物理学家转向了Lindblad主方程形式，通过引入增益和耗散两个独立的退相干通道，成功构建了量子范德珀尔振子的标准模型。近年来，随着离子阱、超导电路等量子模拟平台的飞速发展，该理论模型终于得以在实验室中精确”搭建”[14]。 2025年，中国科学技术大学林毅恒教授团队联合深圳大学及德国斯图加特大学等机构的研究人员，基于上述开放量子系统的理论框架，首次在离子阱实验中实现了量子范德珀尔振子。相关成果于2025年10月10日发表在国际学术期刊《Science Advances》上，演示了量子极限环的形成、量子同步及其增强等非经典行为[17]。
应用与实例
电子电气领域
范德珀尔振子在电子电气领域有着广泛而根基性的应用。其核心的极限环和非线性阻尼特性直接对应于含负阻元件的振荡电路行为。例如，在使用真空管的早期电路中，当系统在极限环附近受到外部周期信号驱动时，振荡频率会被”拖入”驱动信号的频率，两者形成同步——这一现象被称为频率牵引或注入锁定。 范德珀尔方程最初即源于对含非线性负阻元件（真空管、隧道二极管等）的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”的论文，建立了心脏的电子电路模型，将心跳视为松弛震荡[18]。 此后，改良的范德珀尔振荡器被广泛用于模拟心脏起搏器的动作电位时序。Grudziński与Żebrowski在2004年设计了修正的范德珀尔振荡器模型，通过在经典范德珀尔振荡器的相空间中增加一个鞍点和一个节点，并引入非对称阻尼项，使得起搏节点的放电频率可在不改变不应期长度的条件下独立调节[19]。 在神经科学中，Fitzhugh于1961年将范德珀尔方程与Hodgkin-Huxley模型结合，提出了FitzHugh-Nagumo模型（亦称Bonhoeffer-van der Pol模型）[20]。该模型的相平面中存在一个稳定的极限环，对应于神经元的周期性脉冲发放。研究表明，FitzHugh-Nagumo方程可以在很好的近似程度上被约化为一个受外部电流驱动的广义范德珀尔振荡器。
催生混沌理论的诞生
范德珀尔振子及受迫范德珀尔振子呈现的分岔性质，为最初被人类观测到的混沌现象之一，叩开了混沌世界的大门。1945年英国数学家卡特赖特（M.Cartwright）和李特伍德（J.Littlewood）在对受迫范德珀尔振子的研究中，在理论层面首次证明受迫范德珀尔方程存在奇异解，预示了混沌的可能性[21]。随后，日本学者上田睆亮（Ueda Yoshisuke）于1961年首次在模拟计算机上观察到混沌吸引子[22]，随着后续研究的深入，混沌理论的面纱也逐渐被科学家们所揭晓。
参考文献
[1] [2] [3] [4] [5] [8] [9] [6] [10] [11] [7] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22]
- ↑ 1.0 1.1 Gerard-Lescuyer, J. M. A. (1880). Sur un paradoxe électrodynamique.
- ↑ 2.0 2.1 Musical Arcs. (1901). Nature, 63(1640), 542–543.
- ↑ 3.0 3.1 De Forest, L. (1906). The Audion: A New Receiver for Wireless Telegraphy. Transactions of the American Institute of Electrical Engineers.
- ↑ 4.0 4.1 Van der Pol, B. (1920). A theory of the amplitude of free and forced triode vibrations. Radio Review, 1, 701-710.
- ↑ 5.0 5.1 Van der Pol, B., & Van der Mark, J. (1927). Frequency demultiplication. Nature, 120(3019), 363-364.
- ↑ 6.0 6.1 Strogatz, S. H. (2019). Nonlinear dynamics and chaos: with applications to physics, biology, chemistry, and engineering (2nd ed.). CRC Press. (Example 7.3.1)
- ↑ 7.0 7.1 Van der Pol Oscillator. Retrieved from https://www.phys.uconn.edu/~rozman/Courses/P2200_13F/downloads/vanderpol/vanderpol-oscillator-draft.pdf
- ↑ 8.0 8.1 Hopf, E. (1942). Abzweigung einer periodischen Lösung von einer stationären Lösung eines Differential systems. Berichte über die Verhandlungen der Sächsischen Akademie der Wissenschaften zu Leipzig, Mathematisch-Physikalische Klasse, 94, 1-22.
- ↑ 9.0 9.1 Andronov, A. A., & Hopf, E. (1947). Bifurcation of cycles in a certain class of differential equations. Recueil Mathématique [Mathématique] de la Nouvelle Série, 20(62), 189-203.
- ↑ 10.0 10.1 10.2 Strogatz, S. H. (2019). Nonlinear dynamics and chaos: with applications to physics, biology, chemistry, and engineering (2nd ed.). CRC Press. (Example 8.4.1)
- ↑ 11.0 11.1 Perko, L. (2013). Differential Equations and Dynamical Systems (3rd ed.). Springer.
- ↑ 12.0 12.1 Shah, T., Chattopadhyay, R., Vaidya, K., & Chakraborty, S. (2015). Conservative perturbation theory for nonconservative systems. Physical Review E, 92(6), 062927.
- ↑ 13.0 13.1 Chattopadhyay, R., Shah, T., & Chakraborty, S. (2018). Finding the Hannay angle in dissipative oscillatory systems via conservative perturbation theory. Physical Review E, 97(6), 062209.
- ↑ 14.0 14.1 14.2 Chia, A., Kwek, L. C., & Noh, C. (2020). Relaxation oscillations and frequency entrainment in quantum mechanics. Physical Review E, 102(4), 042213.
- ↑ 15.0 15.1 Walter, S., Nunnenkamp, A., & Bruder, C. (2014). Quantum Synchronization of a Driven Self-Sustained Oscillator. Physical Review Letters, 112(9), 094102.
- ↑ 16.0 16.1 Lee, T. E., & Sadeghpour, H. R. (2013). Quantum Synchronization of Quantum van der Pol Oscillators with Trapped Ions. Physical Review Letters, 111(23), 234101.
- ↑ 17.0 17.1 量子范德珀尔振子首次在离子阱实验中实现. Retrieved from https://www.cas.cn/syky/202510/t20251014_5085221.shtml
- ↑ 18.0 18.1 Van der Pol, B., & Van der Mark, J. (1928). LXXII. The heartbeat considered as a relaxation oscillation, and an electrical model of the heart. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 6(38), 763-775.
- ↑ 19.0 19.1 Grudziński, K., & Żebrowski, J. J. (2004). Modeling cardiac pacemakers with relaxation oscillators. Physica A: Statistical Mechanics and Its Applications, 336(1-2), 153-162.
- ↑ 20.0 20.1 Fitzhugh, R. (1961). Impulses and physiological states in theoretical models of nerve membrane. Biophysical Journal, 1(6), 445-466.
- ↑ 21.0 21.1 Cartwright, M. L., & Littlewood, J. E. (1945). 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. Journal of the London Mathematical Society, 20(3), 180-189.
- ↑ 22.0 22.1 Hayashi, C., Shibayama, H., & Ueda, Y. (1961). Quasi-periodic Oscillations in Self-oscillatory Systems with External Force. IECE Technical Report, Nonlinear Theory.