更改

其中 [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）是研究[[极限环]][[非线性系统]]的经典实例。作为一个[[耗散]]性非线性振荡系统，其数学模型为一个二阶非线性齐次方程：

当方程左侧为0的时候，模型为自由范德珀尔振子，不受外界噪声干扰；当方程右侧不为0时，模型变为[[受迫范德珀尔振子]]，可进一步展示[[分岔]]、[[混沌]]等更复杂的非线性现象。

当方程左侧为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{ |\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] 时系统[[耗散]]能量。三者相互耦合，共同塑造了范德珀尔振子从线性振荡到弛豫振荡、从稳定极限环到[[混沌]]的丰富动力学行为。下文将分别从这三个维度展开论述。

===阻尼参数正负不同时的分岔性质===

==阻尼参数正负不同时的分岔性质==

[[文件:阻尼系数不同的范德珀尔振子相图.png|缩略图|阻尼系数不同的范德珀尔振子振动相图]]

[[文件:阻尼系数不同的范德珀尔振子相图.png|缩略图|阻尼系数不同的范德珀尔振子振动相图]]

[[文件:不同起振点的范德珀尔振子.png|缩略图|同一阻尼系数，不同起振点的范德珀尔振子收敛于同一极限环]]

[[文件:不同起振点的范德珀尔振子.png|缩略图|同一阻尼系数，不同起振点的范德珀尔振子收敛于同一极限环]]

上述 [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] 的增加，极限环的非线性动力学特征也在不断增强：

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

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

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

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

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

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

在阻尼系数 [math]\displaystyle{ \mu }[/math] 为正且固定的条件下，系统的能量调节机制完全取决于位移 [math]\displaystyle{ x }[/math] 的大小，这一机制在一个振动周期内呈现出精细的周期性波动，具体表现为非线性阻尼项 [math]\displaystyle{ -\mu(1-x^2)\dot{x} }[/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{ (1-x^2)<0 }[/math]，非线性阻尼项表现为正阻尼，系统耗散能量，使大振幅运动衰减。在相空间中，这对应于轨迹在极限环外侧的区域运动时，阻尼力与速度方向相反，对振子做负功，从系统中抽取能量，迫使轨迹向内螺旋收缩。

当位移较大（[math]\displaystyle{ |x|>1 }[/math]）时，[math]\displaystyle{ (1-x^2)<0 }[/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]，常数）。当电池开关关闭时，根据[[基尔霍夫电压定律]]：

我们可通过如下步骤来探究<ref name="UConn_vanderpol_oscillator" />：考虑如图所示的RLC电路：它由电感 [math]\displaystyle{ L }[/math]、电容器 [math]\displaystyle{ C }[/math]、电阻 [math]\displaystyle{ R }[/math] 组成，我们假设电压源是电池（[math]\displaystyle{ E(t)=E }[/math]，常数）。当电池开关关闭时，根据[[基尔霍夫电压定律]]：