更改

范德珀尔振子（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" />，后被公认为[[混沌]]现象最早的实验观察之一，间接催生了后续[[混沌理论]]的诞生。

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

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

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

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

[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时，模型变为[[受迫范德珀尔振子]]，可进一步展示[[分岔]]、[[混沌]]等更复杂的非线性现象。

==物理意义与动力学行为==

==物理意义与动力学行为==

[[文件:阻尼系数不同的范德珀尔振子相图.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{ \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 }[/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|缩略图|阻尼系数不同的范德珀尔振子机械能随时间变化图]]

==模型的推导过程==

==模型的推导过程==

我们可通过如下步骤来探究<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）过程，这个系统从稳定的极限环走向混沌。

（注：受迫范德珀尔振子因引入时间维度而升维至三维及以上，打破了[[庞加莱-本迪克松定理]]的限制，从而触及[[混沌理论]]的层次。）

（注：受迫范德珀尔振子因引入时间维度而升维至三维及以上，打破了[[庞加莱-本迪克松定理]]的限制，从而触及[[混沌理论]]的层次。）

第二，霍普夫分岔理论延续了范德珀尔对极限环稳定性的分析思路。 在范德珀尔的工作中，极限环的稳定性由非线性阻尼的能量调节机制保证：小振幅时注入能量、大振幅时耗散能量。在霍普夫分岔的一般理论中，这一物理图像被抽象为[[第一李雅普诺夫系数]]（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" />。

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

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

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

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

范德珀尔振子及受迫范德珀尔振子呈现的分岔性质，为最初被人类观测到的混沌现象之一，叩开了混沌世界的大门。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>