| 第1行: |
第1行: |
| − | 鞍结点[[分岔]](Saddle-Node [[Bifurcation]]),又称折叠[[分岔]](Fold [[Bifurcation]])、切线[[分岔]](Tangent [[Bifurcation]])或蓝天[[分岔]](Blue Sky [[Bifurcation]]),是[[非线性系统]]中相对简单且常见的不稳定性机制之一。 | + | 鞍结点[[分岔]](Saddle-Node Bifurcation),又称折叠[[分岔]](Fold Bifurcation)、切线[[分岔]](Tangent Bifurcation)或蓝天[[分岔]](Blue Sky Bifurcation),是[[非线性系统]]中相对简单且常见的不稳定性机制之一。 |
| | | | |
| | =历史= | | =历史= |
| − | 鞍结点分岔的早期研究可追溯到庞加莱(Jules Henri Poincaré)时期。庞加莱最早提出了“分岔”的概念。他重点研究了这样一种现象:当外部条件(即参数)发生改变时,系统原本可能存在的几种不同状态(即方程不同的解的分支)会逐渐靠近、合并,并最终完全消失。他当时描绘的这一过程,与鞍结点分岔的情况非常相似。<ref name="Poincare1885">H. Poincare (1885). sur l'i~ouilibre d'ui~e masse fluide animi~e d'un mouvemeht de rotation</ref> | + | 鞍结点分岔的早期研究可追溯到庞加莱(Jules Henri Poincaré)时期。庞加莱最早提出了“[[分岔]]”的概念。他重点研究了这样一种现象:当外部条件(即参数)发生改变时,系统原本可能存在的几种不同状态(即方程不同的解的分支)会逐渐靠近、合并,并最终完全消失。他当时描绘的这一过程,与鞍结点分岔的情况非常相似。<ref name="Poincare1885">H. Poincare (1885). sur l'i~ouilibre d'ui~e masse fluide animi~e d'un mouvemeht de rotation</ref> |
| | | | |
| | 鞍结点分岔的名称最早由苏联学者安德罗诺夫(A. A. Andronov)与列昂托维奇(E. A. Leontovich)给出,并且在他们的研究中,系统地讨论了参数变化对系统不动点的影响。<ref name="Andronov1957">A. A. Andronov, E. A. Leontovich (1957). On the generation of limit cycles from a loop of a separatrix and from the separatrix of the state of equilibrium of saddle-knot type</ref> | | 鞍结点分岔的名称最早由苏联学者安德罗诺夫(A. A. Andronov)与列昂托维奇(E. A. Leontovich)给出,并且在他们的研究中,系统地讨论了参数变化对系统不动点的影响。<ref name="Andronov1957">A. A. Andronov, E. A. Leontovich (1957). On the generation of limit cycles from a loop of a separatrix and from the separatrix of the state of equilibrium of saddle-knot type</ref> |
| | | | |
| | =标准型和性质= | | =标准型和性质= |
| − | 在[[数学]]上,我们可以通过标准型来研究[[分岔]]的本质特征。对于一维[[动力系统]](参考[[一维动力学]]),鞍结点[[分岔]]的最简[[微分方程]]形式(标准型)如下<ref name="Strogatz1994">Strogatz, S. H. (1994). ''Nonlinear Dynamics and Chaos.'' Reading, MA: Perseus Books.</ref>
| + | 在数学上,我们可以通过标准型来研究[[分岔]]的本质特征。对于一维[[动力系统]](参考[[一维动力学]]),鞍结点[[分岔]]的最简微分方程形式(标准型)如下<ref name="Strogatz1994">Strogatz, S. H. (1994). ''Nonlinear Dynamics and Chaos.'' Reading, MA: Perseus Books.</ref> |
| | {{NumBlk|:| | | {{NumBlk|:| |
| | <math display="block"> | | <math display="block"> |
| 第18行: |
第18行: |
| | ** <math>x</math> 的取值范围应为[[不动点]] <math>x^*</math> 的 <math>\varepsilon</math>-邻域:<math>x \in (x^* - \varepsilon, x^* + \varepsilon)</math>,其中 <math>\varepsilon > 0</math> 且充分小。 | | ** <math>x</math> 的取值范围应为[[不动点]] <math>x^*</math> 的 <math>\varepsilon</math>-邻域:<math>x \in (x^* - \varepsilon, x^* + \varepsilon)</math>,其中 <math>\varepsilon > 0</math> 且充分小。 |
| | ** <math>r</math> 的取值范围应为[[临界]]值 <math>r_c</math> 的 <math>\delta</math>-邻域:<math>r \in (r_c - \delta, r_c + \delta)</math>,其中 <math>\delta > 0</math> 且充分小。 | | ** <math>r</math> 的取值范围应为[[临界]]值 <math>r_c</math> 的 <math>\delta</math>-邻域:<math>r \in (r_c - \delta, r_c + \delta)</math>,其中 <math>\delta > 0</math> 且充分小。 |
| − | ** 只有在此微小邻域内,[[泰勒展开]]后的低阶项才能主导系统的定性行为。具体可参见下文中的推广形式。 | + | ** 只有在此微小邻域内,泰勒展开后的低阶项才能主导系统的定性行为。具体可参见下文中的推广形式。 |
| | | | |
| − | 鞍结点[[分岔]]的关键性质在于,随着系统参数<math>r</math>的改变,[[相空间]]中原本存在的两个[[不动点]](一个稳定,一个不稳定)相互靠近、碰撞并最终相互抵消而消失,导致系统原本的[[稳定状态]]不复存在。也可反过来看待,原本状态不稳定的系统因系统参数<math>r</math>的改变,使其先出现了一个半[[稳定不动点]],最终出现了两个[[不动点]](一个稳定,一个不稳定)。
| + | 鞍结点分岔的关键性质在于,随着系统参数<math>r</math>的改变,[[相空间]]中原本存在的两个[[不动点]](一个稳定,一个不稳定)相互靠近、碰撞并最终相互抵消而消失,导致系统原本的[[稳定状态]]不复存在。也可反过来看待,原本状态不稳定的系统因系统参数<math>r</math>的改变,使其先出现了一个半稳定不动点,最终出现了两个不动点(一个稳定,一个不稳定)。 |
| | | | |
| | ==[[不动点]]和稳定性分析== | | ==[[不动点]]和稳定性分析== |
| 第38行: |
第38行: |
| | | | |
| | 当<math>r>0</math>时, | | 当<math>r>0</math>时, |
| − | * 分支<math>x^*=\sqrt{r}</math>,有特征值<math>f'(x^*)=-2r<0</math>,微小扰动会不断衰减,是稳定[[不动点]]; | + | * 分支<math>x^*=\sqrt{r}</math>,有特征值<math>f'(x^*)=-2r<0</math>,微小扰动会不断衰减,是稳定不动点; |
| − | * 分支<math>x^*=-\sqrt{r}</math>,有<math>f'(x^*)=2r>0</math>,微小扰动会被放大,是不稳定[[不动点]]。 | + | * 分支<math>x^*=-\sqrt{r}</math>,有<math>f'(x^*)=2r>0</math>,微小扰动会被放大,是不稳定不动点。 |
| | 当 <math>r=0</math>,有<math>f'(x^*)=0</math>,此处为[[分岔]]点。 | | 当 <math>r=0</math>,有<math>f'(x^*)=0</math>,此处为[[分岔]]点。 |
| | | | |
| 第60行: |
第60行: |
| | \dot{x} = r + x^2 | | \dot{x} = r + x^2 |
| | </math> | | </math> |
| − | 情况类似,只是[[方向]]发生了调转。
| + | 情况类似,只是方向发生了调转。 |
| | | | |
| | [[File:相图:鞍结点分岔标准型2.png|800px]] | | [[File:相图:鞍结点分岔标准型2.png|800px]] |
| 第69行: |
第69行: |
| | | | |
| | ==推广形式== | | ==推广形式== |
| − | 为了更好地发现在复杂的[[动力系统]]中,常以局部[[分岔]]的方式出现的鞍结点[[分岔]],我们需要了解它的推广形式。<ref name="Strogatz1994" /> | + | 为了更好地发现在复杂的[[动力系统]]中,常以局部分岔的方式出现的鞍结点分岔,我们需要了解它的推广形式。<ref name="Strogatz1994" /> |
| | {{NumBlk|:| | | {{NumBlk|:| |
| | <math display="block"> | | <math display="block"> |
| 第77行: |
第77行: |
| | 可以看作是一阶参数与二阶自变量的线性组合。 | | 可以看作是一阶参数与二阶自变量的线性组合。 |
| | ===推导过程=== | | ===推导过程=== |
| − | 将<math>\dot x=f</math>视为关于<math>x</math>与<math>r</math>的函数<math>f(x,r)</math>,对<math>x</math>,<math>r</math>在<math>x^*</math>与<math>r_c</math>处分别进行[[泰勒展开]]得到(其中<math>r_c</math>是[[临界状态]]时的<math>r</math>值) | + | 将<math>\dot x=f</math>视为关于<math>x</math>与<math>r</math>的函数<math>f(x,r)</math>,对<math>x</math>,<math>r</math>在<math>x^*</math>与<math>r_c</math>处分别进行泰勒展开得到(其中<math>r_c</math>是[[临界状态]]时的<math>r</math>值) |
| | {{NumBlk|:| | | {{NumBlk|:| |
| | <math display="block"> | | <math display="block"> |
| 第93行: |
第93行: |
| | | | |
| | <math display="block"> | | <math display="block"> |
| − | \begin{aligned} | + | \begin{aligned} |
| | \dot{x} & =r-x-\mathrm{e}^{-x} \\ | | \dot{x} & =r-x-\mathrm{e}^{-x} \\ |
| | & =r-x-\left[1-x+\frac{x^{2}}{2!}+\cdots\right] \\ | | & =r-x-\left[1-x+\frac{x^{2}}{2!}+\cdots\right] \\ |
| | & =(r-1)-\frac{x^{2}}{2}+\cdots | | & =(r-1)-\frac{x^{2}}{2}+\cdots |
| − | \end{aligned} | + | \end{aligned} |
| | </math> | | </math> |
| | | | |
| 第128行: |
第128行: |
| | </math> | | </math> |
| | | | |
| − | * 令<math>\frac{\mathrm{d} x}{\mathrm{d} \tau } =0</math>。生态学背景下,一般假设<math>r>0</math>。此时<math>x^*=0</math>通常为不稳定[[不动点]],因此分析实际问题时一般关注<math>x^*>0</math>的[[不动点]]的情况,即 | + | * 令<math>\frac{\mathrm{d} x}{\mathrm{d} \tau } =0</math>。生态学背景下,一般假设<math>r>0</math>。此时<math>x^*=0</math>通常为不稳定不动点,因此分析实际问题时一般关注<math>x^*>0</math>的不动点的情况,即 |
| | | | |
| | <math display="block"> | | <math display="block"> |
| 第143行: |
第143行: |
| | 因<math>k</math>是一个慢变量,所以将其视为一个常数(上图中<math>k=10</math>),将<math>r</math>(上图纵坐标)视为主要控制参数。 | | 因<math>k</math>是一个慢变量,所以将其视为一个常数(上图中<math>k=10</math>),将<math>r</math>(上图纵坐标)视为主要控制参数。 |
| | | | |
| − | * 当<math>r</math>很小时,可看到两条曲线仅有一个交点,且为一个稳定[[不动点]],记为不动点1。 | + | * 当<math>r</math>很小时,可看到两条曲线仅有一个交点,且为一个稳定不动点,记为不动点1。 |
| | | | |
| − | * 令<math>r</math>逐渐增大到第一个[[临界]]值<math>r_1</math>时(约为<math>0.384</math>),<math>y_1(x)</math>与<math>y_2(x)</math>的下降段相切,出现了半[[稳定不动点]]。 | + | * 令<math>r</math>逐渐增大到第一个临界值<math>r_1</math>时(约为<math>0.384</math>),<math>y_1(x)</math>与<math>y_2(x)</math>的下降段相切,出现了半稳定不动点。 |
| | | | |
| − | * 继续增大<math>r</math>,<math>y_1(x)</math>与<math>y_2(x)</math>的下降段出现了两个[[不动点]](一个不稳定,一个稳定),记为不动点2、不动点3。以上是该系统第一个鞍结点分岔。 | + | * 继续增大<math>r</math>,<math>y_1(x)</math>与<math>y_2(x)</math>的下降段出现了两个不动点(一个不稳定,一个稳定),记为不动点2、不动点3。以上是该系统第一个鞍结点分岔。 |
| | | | |
| − | * 再接着增大<math>r</math>到第二个[[临界]]值<math>r_2</math>(约为<math>0.56</math>),<math>y_1(x)</math>与<math>y_2(x)</math>在单峰附近相切,原本不动点1与不动点2相遇出现了一个半[[稳定不动点]]。 | + | * 再接着增大<math>r</math>到第二个[[临界]]值<math>r_2</math>(约为<math>0.56</math>),<math>y_1(x)</math>与<math>y_2(x)</math>在单峰附近相切,原本不动点1与不动点2相遇出现了一个半稳定不动点。 |
| | | | |
| | * 继续增大<math>r</math>,<math>y_1(x)</math>与<math>y_2(x)</math>只剩下不动点3。以上是该系统的第二个鞍结点分岔。 | | * 继续增大<math>r</math>,<math>y_1(x)</math>与<math>y_2(x)</math>只剩下不动点3。以上是该系统的第二个鞍结点分岔。 |
| 第162行: |
第162行: |
| | | | |
| | ==物理意义== | | ==物理意义== |
| − | 上面的相图和分岔图描述了当外部环境条件(参数<math>r</math>)发生缓慢且[[连续]]的变化时,系统的最终归宿([[稳定不动点]]的位置和数量)会发生怎样的改变。 | + | 上面的相图和分岔图描述了当外部环境条件(参数<math>r</math>)发生缓慢且连续的变化时,系统的最终归宿([[稳定不动点]]的位置和数量)会发生怎样的改变。 |
| | | | |
| | * 当参数<math>r</math>在较小数值范围内变动(<math>0<r<0.384</math>):系统只有一个[[不动点]],环境的改变只引起[[系统状态]]的微小调整。比如气温稍微升高,蚜虫数量只是略微增加,依然会被捕食者维持在较低数量的稳定态。 | | * 当参数<math>r</math>在较小数值范围内变动(<math>0<r<0.384</math>):系统只有一个[[不动点]],环境的改变只引起[[系统状态]]的微小调整。比如气温稍微升高,蚜虫数量只是略微增加,依然会被捕食者维持在较低数量的稳定态。 |
| 第174行: |
第174行: |
| | | | |
| | =二维系统(以上)= | | =二维系统(以上)= |
| − | 与一维[[动力系统]]中的鞍结点[[分岔]]类似,二维及以上的鞍结点[[分岔]]同样描述了一个稳定[[不动点]]和一个不稳定[[不动点]]随参数变化相互靠近、碰撞并最终同时消失的过程。 | + | 与一维[[动力系统]]中的鞍结点[[分岔]]类似,二维及以上的鞍结点[[分岔]]同样描述了一个稳定不动点和一个不稳定不动点随参数变化相互靠近、碰撞并最终同时消失的过程。 |
| | | | |
| | ==二维标准型== | | ==二维标准型== |
| 第182行: |
第182行: |
| | </math><ref name="Strogatz1994" /> | | </math><ref name="Strogatz1994" /> |
| | | | |
| − | * <math>x</math>与<math>y</math>为[[状态变量]]。 | + | * <math>x</math>与<math>y</math>为状态变量。 |
| − | * <math>\mu</math>为[[分岔参数]]。 | + | * <math>\mu</math>为分岔参数。 |
| | * 一式描述了系统在<math>x</math>维度上的[[分岔]]行为。 | | * 一式描述了系统在<math>x</math>维度上的[[分岔]]行为。 |
| − | * 二式描述了系统在<math>y</math>维度上的收缩行为(<math>\lambda_y = -1</math>),表明系统在非[[分岔]][[方向]]上是稳定的。 | + | * 二式描述了系统在<math>y</math>维度上的收缩行为(<math>\lambda_y = -1</math>),表明系统在非[[分岔]]方向上是稳定的。 |
| | | | |
| | 三维以上的标准型与二维类似,仅在某个维度有[[分岔]]行为(类似<math>\dot x = \mu - x^2</math>),系统在其余维度上稳定(均为类似<math>\dot y=ay</math>的形式,其中参数<math>a<0</math>)。 | | 三维以上的标准型与二维类似,仅在某个维度有[[分岔]]行为(类似<math>\dot x = \mu - x^2</math>),系统在其余维度上稳定(均为类似<math>\dot y=ay</math>的形式,其中参数<math>a<0</math>)。 |
| 第193行: |
第193行: |
| | | | |
| | | | |
| − | x[[方向]]的[[不动点]]分析与一维情况完全类似。随着参数 <math>\mu</math> 的变化,二维相平面上的[[拓扑结构]]经历以下三个阶段:
| + | x方向的[[不动点]]分析与一维情况完全类似。随着参数 <math>\mu</math> 的变化,二维相平面上的[[拓扑结构]]经历以下三个阶段: |
| | | | |
| | 1. 双[[不动点]]共存阶段 (<math>\mu > 0</math>): | | 1. 双[[不动点]]共存阶段 (<math>\mu > 0</math>): |
| − | * 稳定[[不动点]] <math>(\sqrt{\mu}, 0)</math>:在两个特征[[方向]]上均表现为吸引,对应系统的一个[[稳态]]。 | + | * 稳定[[不动点]] <math>(\sqrt{\mu}, 0)</math>:在两个特征方向上均表现为吸引,对应系统的一个[[稳态]]。 |
| − | * 不稳定[[不动点]] <math>(-\sqrt{\mu}, 0)</math>:在 <math>y</math> [[方向]]吸引,在 <math>x</math> [[方向]]排斥。 | + | * 不稳定[[不动点]] <math>(-\sqrt{\mu}, 0)</math>:在 <math>y</math> 方向吸引,在 <math>x</math> 方向排斥。 |
| | [[动力学]]特征:[[相空间]]被不稳定[[不动点]]分割,大部分轨线最终收敛于稳定[[不动点]]。 | | [[动力学]]特征:[[相空间]]被不稳定[[不动点]]分割,大部分轨线最终收敛于稳定[[不动点]]。 |
| | | | |
| 第204行: |
第204行: |
| | | | |
| | 3. [[不动点]]消失与[[幽灵吸引子]] (<math>\mu < 0</math>): | | 3. [[不动点]]消失与[[幽灵吸引子]] (<math>\mu < 0</math>): |
| − | 方程在[[实数]]域内无解,所有[[不动点]]消失。
| + | 方程在实数域内无解,所有[[不动点]]消失。 |
| | | | |
| − | * [[幽灵吸引子]] (Ghost [[Attractor]]):虽然[[不动点]]不再存在,但在原点附近区域,[[向量场]]流速极慢(<math>\dot{x} \approx 0</math>)。[[系统状态]]在穿过该区域时会发生显著的[[临界迟滞]]现象,形成所谓的“[[瓶颈效应]]”。 | + | * [[幽灵吸引子]] (Ghost Attractor):虽然[[不动点]]不再存在,但在原点附近区域,[[向量场]]流速极慢(<math>\dot{x} \approx 0</math>)。[[系统状态]]在穿过该区域时会发生显著的[[临界迟滞]]现象,形成所谓的“[[瓶颈效应]]”。 |
| | | | |
| | =相关概念= | | =相关概念= |
| | | | |
| − | 鞍结点[[分岔]]的形式在所有[[分岔]]现象里是最为普遍的一种。其在系统中常作为一种[[局部分岔]](Local [[Bifurcation]])出现。高维系统会有[[环上的鞍结点分岔]]、[[无限周期分岔]]等其他更复杂的[[分岔]]类型,而在某些局部上可以将其看成是一维的鞍结点[[分岔]]来进行分析,所以掌握一维鞍结点[[分岔]]的性质是非常基本和重要的。 | + | 鞍结点[[分岔]]的形式在所有[[分岔]]现象里是最为普遍的一种。其在系统中常作为一种[[局部分岔]](Local Bifurcation)出现。高维系统会有[[环上的鞍结点分岔]]、[[无限周期分岔]]等其他更复杂的[[分岔]]类型,而在某些局部上可以将其看成是一维的鞍结点[[分岔]]来进行分析,所以掌握一维鞍结点[[分岔]]的性质是非常基本和重要的。 |
| | | | |
| | 从[[突变论]](catastrophe theory)的角度看<ref>Arnold, V. I., Afrajmovich, V. S., & Il’yashenko, Y. S. (2013). Dynamical Systems V: Bifurcation Theory and Catastrophe Theory. Springer. https://books.google.com/books?id=8-LtCAAAQBAJ</ref>,鞍结点[[分岔]]实际上是[[折叠突变]](fold catastrophe)在[[动力系统]]中的具体表现形式。[[突变论]]由 [[René Thom]] 建立,其核心思想是用[[奇点理论]]研究[[系统状态]]随控制参数[[连续]]变化而发生的不[[连续]]跃迁。在 Thom 的七种基本突变类型中,折叠突变是最简单的一类,通常由一个[[状态变量]]和一个控制参数描述。 | | 从[[突变论]](catastrophe theory)的角度看<ref>Arnold, V. I., Afrajmovich, V. S., & Il’yashenko, Y. S. (2013). Dynamical Systems V: Bifurcation Theory and Catastrophe Theory. Springer. https://books.google.com/books?id=8-LtCAAAQBAJ</ref>,鞍结点[[分岔]]实际上是[[折叠突变]](fold catastrophe)在[[动力系统]]中的具体表现形式。[[突变论]]由 [[René Thom]] 建立,其核心思想是用[[奇点理论]]研究[[系统状态]]随控制参数[[连续]]变化而发生的不[[连续]]跃迁。在 Thom 的七种基本突变类型中,折叠突变是最简单的一类,通常由一个[[状态变量]]和一个控制参数描述。 |
| | | | |
| − | 因此,在[[奇点理论]]框架下,折叠突变描述的是平衡解[[集合]]几何结构的折叠[[奇点]],而鞍结点[[分岔]]则是该几何结构在[[动力系统]]中的[[动力学]]实现。二者在[[数学]]结构上由相同的[[奇点]]类型控制,只是描述视角分别来自[[突变论]]的[[势函数]]几何和[[分岔理论]]的[[相空间]][[动力学]]。这种对应关系揭示了[[分岔理论]]与[[突变论]]之间的深层联系:许多基本[[分岔]]都可以理解为 Thom 基本突变在[[动力系统]]中的具体体现。 | + | 因此,在[[奇点理论]]框架下,折叠突变描述的是平衡解[[集合]]几何结构的折叠[[奇点]],而鞍结点[[分岔]]则是该几何结构在[[动力系统]]中的[[动力学]]实现。二者在数学结构上由相同的[[奇点]]类型控制,只是描述视角分别来自[[突变论]]的[[势函数]]几何和[[分岔理论]]的[[相空间]][[动力学]]。这种对应关系揭示了[[分岔理论]]与[[突变论]]之间的深层联系:许多基本[[分岔]]都可以理解为 Thom 基本突变在[[动力系统]]中的具体体现。 |
| | | | |
| | =参考文献= | | =参考文献= |
| | <references /> | | <references /> |