“非线性中心”的版本间的差异
| 第1行: | 第1行: | ||
{{DISPLAYTITLE:非线性中心}} | {{DISPLAYTITLE:非线性中心}} | ||
| − | + | 在非线性动力学中,'''非线性中心(Nonlinear Center)'''是指相平面上某一孤立不动点附近的所有轨迹均为闭合曲线(周期轨道)的动力学现象。在线性系统中,若雅可比矩阵的特征值为纯虚数,则该不动点为线性中心;然而,对于非线性系统,特征值为纯虚数并不能保证该点一定是中心——微小的非线性项就可能将中心转变为稳定或不稳定的'''焦点(Focus)'''。 | |
| − | + | 为了判定非线性系统的中心稳定性,需要借助系统的附加结构性质。两类重要的系统——'''保守系统'''与'''可逆系统'''——提供了判定非线性中心的有力工具。这两种方法分别从能量守恒的几何性质和时间反演的对称性出发,确立了非线性中心存在的充分条件。 | |
| − | |||
| − | + | == '''保守系统中的非线性中心''' == | |
| + | === 定义 === | ||
| + | 在经典力学与非线性动力学中,若一个二阶系统 <math>\dot{\mathbf{x}} = \mathbf{f}(\mathbf{x})</math>(其中 <math>\mathbf{x} \in \mathbb{R}^2</math>)存在一个光滑的实值函数 <math>E(\mathbf{x})</math>,使得沿任意轨迹 <math>\mathbf{x}(t)</math> 均有 <math>\frac{dE}{dt} = 0</math>,且 <math>E(\mathbf{x})</math> 在任何开集上不为常数,则称该系统为'''保守系统''',<math>E(\mathbf{x})</math> 称为'''守恒量'''或'''第一积分'''<ref>STROGATZ S H. 非线性动力学与混沌[M]. 孙梅, 汪小帆, 译. 2版. 北京: 机械工业出版社, 2017</ref>。 | ||
| − | + | === 定理(保守系统的非线性中心)=== | |
| + | 考虑系统 <math>\dot{\mathbf{x}} = \mathbf{f}(\mathbf{x})</math>,其中 <math>\mathbf{x} = (x, y) \in \mathbb{R}^2</math>,<math>\mathbf{f}</math> 是连续可微的。假设存在一个守恒量 <math>E(\mathbf{x})</math>,且 <math>\mathbf{x}^*</math> 是一个'''孤立不动点'''(即在其某个邻域内没有其他不动点)。若 <math>\mathbf{x}^*</math> 是 <math>E(\mathbf{x})</math> 的一个'''局部极小值'''(或局部极大值),则所有充分接近 <math>\mathbf{x}^*</math> 的轨迹都是闭合的,即 <math>\mathbf{x}^*</math> 是一个'''非线性中心'''。 | ||
| − | + | '''证明思路:''' 由于能量沿轨迹守恒,每条轨迹都位于 <math>E(\mathbf{x})</math> 的某一等高线上。在局部极值点附近,等高线是闭合的。又因假设 <math>\mathbf{x}^*</math> 是孤立不动点,故在该点的小邻域内,等高线上不存在其他不动点,因此轨迹必定沿着闭合等高线运动,形成周期轨道。 | |
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| + | [[文件:保守系统的非线性中心示例图.png|缩略图|图1:双阱系统能量等高线图及其相平面投影|right|300px]] | ||
| − | == | + | === 示例:双阱势系统 === |
| − | + | 考虑一个单位质量粒子在双阱势 <math>V(x) = -\frac{1}{2}x^2 + \frac{1}{4}x^4</math> 中的无阻尼运动。运动方程为 <math>\ddot{x} = x - x^3</math>,等价系统为: | |
| + | :<math>\dot{x} = y, \quad \dot{y} = x - x^3.</math> | ||
| + | 系统具有守恒量(总能量): | ||
| + | :<math>E(x, y) = \frac{1}{2}y^2 - \frac{1}{2}x^2 + \frac{1}{4}x^4.</math> | ||
| + | 不动点为 <math>(0,0)</math> 和 <math>(\pm 1, 0)</math>。在 <math>(0,0)</math> 处,能量曲面的等高线呈现鞍点特征;而在 <math>(\pm 1, 0)</math> 处,<math>E(x, y)</math> 达到局部极小值。根据定理,这两个不动点均为非线性中心,附近存在一族闭轨(周期轨道)。 | ||
| + | == '''可逆系统中的非线性中心''' == | ||
| + | ===定义=== | ||
| + | 如果存在一个相空间到自身的映射 <math>R(\mathbf{x})</math>,满足 <math>R^2(\mathbf{x}) = \mathbf{x}</math>(即 <math>R</math> 是一个对合变换),使得系统 <math>\dot{\mathbf{x}} = \mathbf{f}(\mathbf{x})</math> 在变换 <math>t \to -t</math>,<math>\mathbf{x} \to R(\mathbf{x})</math> 下保持不变,则称该系统为'''可逆系统'''。 | ||
| + | 最常见的二维情形是系统具有关于 <math>x</math> 轴的反射对称性:即变换 <math>t \to -t</math>,<math>y \to -y</math> 保持系统不变。这等价于要求方程 <math>\dot{x} = f(x,y)</math>,<math>\dot{y} = g(x,y)</math> 中的 <math>f</math> 关于 <math>y</math> 是奇函数(<math>f(x,-y) = -f(x,y)</math>),且 <math>g</math> 关于 <math>y</math> 是偶函数(<math>g(x,-y) = g(x,y)</math>)。 | ||
| + | === 定理(可逆系统的非线性中心)=== | ||
| + | 假设原点 <math>\mathbf{x}^* = \mathbf{0}</math> 是一个连续可微可逆系统的'''线性中心'''。若系统是可逆的,则原点也是一个'''非线性中心''',即所有充分接近原点的轨迹都是闭合的。 | ||
| + | '''证明思路:''' 考虑一条从正 <math>x</math> 轴出发靠近原点的轨迹。由于线性部分是中心,该轨迹将绕原点旋转并最终与负 <math>x</math> 轴相交。利用可逆性,将轨迹关于 <math>x</math> 轴反射并反转时间方向,将得到一条具有相同端点但箭头相反的轨迹。这两条轨迹共同构成一个闭合曲线。因此原点被一族闭轨包围。 | ||
| − | == | + | === 示例:无阻尼钟摆 === |
| − | + | 无阻尼钟摆的非线性本质通过小角度的近似值<math>sin𝜃 ≈ 𝜃</math>被忽略。利用相图的分析钟摆旋转到顶部这一大角度的区域,在缺少阻力和外部吸引力的情况下,钟摆由下式决定: | |
| − | + | :<math>\frac{d^2\theta}{dt^2} + \frac{g}{L}\sin\theta = 0.</math> | |
| − | = | + | 令 <math>\omega = \sqrt{g/L}</math>,<math>\tau = \omega t</math> 进行无量纲化,得: |
| − | + | :<math>\ddot{\theta} + \sin\theta = 0.</math> | |
| − | + | 写成相平面系统(<math>v = \dot{\theta}</math>): | |
| − | + | :<math>\dot{\theta} = v, \quad \dot{v} = -\sin\theta.</math> | |
| − | + | 该系统在变换 <math>\tau \to -\tau</math>,<math>v \to -v</math> 下保持不变,因此是可逆的。不动点 <math>(\theta^*, v^*) = (0,0)</math> 处的雅可比矩阵为 <math>\begin{pmatrix}0 & 1 \\ -1 & 0\end{pmatrix}</math>,特征值为 <math>\pm i</math>,故为线性中心。根据可逆系统的非线性中心定理,<math>(0,0)</math> 也是一个非线性中心,对应钟摆在最低点的微小摆动。此外,系统也具有守恒量 <math>E = \frac{1}{2}v^2 - \cos\theta</math>,且 <math>(0,0)</math> 是 <math>E</math> 的局部极小值点,这从保守系统的角度也印证了该结论。 | |
| − | |||
| − | |||
| − | |||
| − | |||
[[文件:钟摆.png|缩略图|右|图2:无阻尼钟摆|right|250px]] | [[文件:钟摆.png|缩略图|右|图2:无阻尼钟摆|right|250px]] | ||
| − | + | ||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
2026年4月17日 (五) 17:22的版本
在非线性动力学中,非线性中心(Nonlinear Center)是指相平面上某一孤立不动点附近的所有轨迹均为闭合曲线(周期轨道)的动力学现象。在线性系统中,若雅可比矩阵的特征值为纯虚数,则该不动点为线性中心;然而,对于非线性系统,特征值为纯虚数并不能保证该点一定是中心——微小的非线性项就可能将中心转变为稳定或不稳定的焦点(Focus)。
为了判定非线性系统的中心稳定性,需要借助系统的附加结构性质。两类重要的系统——保守系统与可逆系统——提供了判定非线性中心的有力工具。这两种方法分别从能量守恒的几何性质和时间反演的对称性出发,确立了非线性中心存在的充分条件。
保守系统中的非线性中心
定义
在经典力学与非线性动力学中,若一个二阶系统 [math]\displaystyle{ \dot{\mathbf{x}} = \mathbf{f}(\mathbf{x}) }[/math](其中 [math]\displaystyle{ \mathbf{x} \in \mathbb{R}^2 }[/math])存在一个光滑的实值函数 [math]\displaystyle{ E(\mathbf{x}) }[/math],使得沿任意轨迹 [math]\displaystyle{ \mathbf{x}(t) }[/math] 均有 [math]\displaystyle{ \frac{dE}{dt} = 0 }[/math],且 [math]\displaystyle{ E(\mathbf{x}) }[/math] 在任何开集上不为常数,则称该系统为保守系统,[math]\displaystyle{ E(\mathbf{x}) }[/math] 称为守恒量或第一积分[1]。
定理(保守系统的非线性中心)
考虑系统 [math]\displaystyle{ \dot{\mathbf{x}} = \mathbf{f}(\mathbf{x}) }[/math],其中 [math]\displaystyle{ \mathbf{x} = (x, y) \in \mathbb{R}^2 }[/math],[math]\displaystyle{ \mathbf{f} }[/math] 是连续可微的。假设存在一个守恒量 [math]\displaystyle{ E(\mathbf{x}) }[/math],且 [math]\displaystyle{ \mathbf{x}^* }[/math] 是一个孤立不动点(即在其某个邻域内没有其他不动点)。若 [math]\displaystyle{ \mathbf{x}^* }[/math] 是 [math]\displaystyle{ E(\mathbf{x}) }[/math] 的一个局部极小值(或局部极大值),则所有充分接近 [math]\displaystyle{ \mathbf{x}^* }[/math] 的轨迹都是闭合的,即 [math]\displaystyle{ \mathbf{x}^* }[/math] 是一个非线性中心。
证明思路: 由于能量沿轨迹守恒,每条轨迹都位于 [math]\displaystyle{ E(\mathbf{x}) }[/math] 的某一等高线上。在局部极值点附近,等高线是闭合的。又因假设 [math]\displaystyle{ \mathbf{x}^* }[/math] 是孤立不动点,故在该点的小邻域内,等高线上不存在其他不动点,因此轨迹必定沿着闭合等高线运动,形成周期轨道。
示例:双阱势系统
考虑一个单位质量粒子在双阱势 [math]\displaystyle{ V(x) = -\frac{1}{2}x^2 + \frac{1}{4}x^4 }[/math] 中的无阻尼运动。运动方程为 [math]\displaystyle{ \ddot{x} = x - x^3 }[/math],等价系统为:
- [math]\displaystyle{ \dot{x} = y, \quad \dot{y} = x - x^3. }[/math]
系统具有守恒量(总能量):
- [math]\displaystyle{ E(x, y) = \frac{1}{2}y^2 - \frac{1}{2}x^2 + \frac{1}{4}x^4. }[/math]
不动点为 [math]\displaystyle{ (0,0) }[/math] 和 [math]\displaystyle{ (\pm 1, 0) }[/math]。在 [math]\displaystyle{ (0,0) }[/math] 处,能量曲面的等高线呈现鞍点特征;而在 [math]\displaystyle{ (\pm 1, 0) }[/math] 处,[math]\displaystyle{ E(x, y) }[/math] 达到局部极小值。根据定理,这两个不动点均为非线性中心,附近存在一族闭轨(周期轨道)。
可逆系统中的非线性中心
定义
如果存在一个相空间到自身的映射 [math]\displaystyle{ R(\mathbf{x}) }[/math],满足 [math]\displaystyle{ R^2(\mathbf{x}) = \mathbf{x} }[/math](即 [math]\displaystyle{ R }[/math] 是一个对合变换),使得系统 [math]\displaystyle{ \dot{\mathbf{x}} = \mathbf{f}(\mathbf{x}) }[/math] 在变换 [math]\displaystyle{ t \to -t }[/math],[math]\displaystyle{ \mathbf{x} \to R(\mathbf{x}) }[/math] 下保持不变,则称该系统为可逆系统。 最常见的二维情形是系统具有关于 [math]\displaystyle{ x }[/math] 轴的反射对称性:即变换 [math]\displaystyle{ t \to -t }[/math],[math]\displaystyle{ y \to -y }[/math] 保持系统不变。这等价于要求方程 [math]\displaystyle{ \dot{x} = f(x,y) }[/math],[math]\displaystyle{ \dot{y} = g(x,y) }[/math] 中的 [math]\displaystyle{ f }[/math] 关于 [math]\displaystyle{ y }[/math] 是奇函数([math]\displaystyle{ f(x,-y) = -f(x,y) }[/math]),且 [math]\displaystyle{ g }[/math] 关于 [math]\displaystyle{ y }[/math] 是偶函数([math]\displaystyle{ g(x,-y) = g(x,y) }[/math])。
定理(可逆系统的非线性中心)
假设原点 [math]\displaystyle{ \mathbf{x}^* = \mathbf{0} }[/math] 是一个连续可微可逆系统的线性中心。若系统是可逆的,则原点也是一个非线性中心,即所有充分接近原点的轨迹都是闭合的。 证明思路: 考虑一条从正 [math]\displaystyle{ x }[/math] 轴出发靠近原点的轨迹。由于线性部分是中心,该轨迹将绕原点旋转并最终与负 [math]\displaystyle{ x }[/math] 轴相交。利用可逆性,将轨迹关于 [math]\displaystyle{ x }[/math] 轴反射并反转时间方向,将得到一条具有相同端点但箭头相反的轨迹。这两条轨迹共同构成一个闭合曲线。因此原点被一族闭轨包围。
示例:无阻尼钟摆
无阻尼钟摆的非线性本质通过小角度的近似值[math]\displaystyle{ sin𝜃 ≈ 𝜃 }[/math]被忽略。利用相图的分析钟摆旋转到顶部这一大角度的区域,在缺少阻力和外部吸引力的情况下,钟摆由下式决定:
- [math]\displaystyle{ \frac{d^2\theta}{dt^2} + \frac{g}{L}\sin\theta = 0. }[/math]
令 [math]\displaystyle{ \omega = \sqrt{g/L} }[/math],[math]\displaystyle{ \tau = \omega t }[/math] 进行无量纲化,得:
- [math]\displaystyle{ \ddot{\theta} + \sin\theta = 0. }[/math]
写成相平面系统([math]\displaystyle{ v = \dot{\theta} }[/math]):
- [math]\displaystyle{ \dot{\theta} = v, \quad \dot{v} = -\sin\theta. }[/math]
该系统在变换 [math]\displaystyle{ \tau \to -\tau }[/math],[math]\displaystyle{ v \to -v }[/math] 下保持不变,因此是可逆的。不动点 [math]\displaystyle{ (\theta^*, v^*) = (0,0) }[/math] 处的雅可比矩阵为 [math]\displaystyle{ \begin{pmatrix}0 & 1 \\ -1 & 0\end{pmatrix} }[/math],特征值为 [math]\displaystyle{ \pm i }[/math],故为线性中心。根据可逆系统的非线性中心定理,[math]\displaystyle{ (0,0) }[/math] 也是一个非线性中心,对应钟摆在最低点的微小摆动。此外,系统也具有守恒量 [math]\displaystyle{ E = \frac{1}{2}v^2 - \cos\theta }[/math],且 [math]\displaystyle{ (0,0) }[/math] 是 [math]\displaystyle{ E }[/math] 的局部极小值点,这从保守系统的角度也印证了该结论。
参考资料
- ↑ STROGATZ S H. 非线性动力学与混沌[M]. 孙梅, 汪小帆, 译. 2版. 北京: 机械工业出版社, 2017