248
个编辑
更改
跳到导航
跳到搜索
第1行：
第1行： − 鞍结点分岔（Saddle-Node Bifurcation），又称折叠分岔（Fold Bifurcation）、切线分岔（Tangent Bifurcation）或蓝天分岔（Blue Sky Bifurcation），是非线性系统中最常见的不稳定性机制之一。鞍结点分岔常作为一种局部分岔（Local Bifurcation）出现。+ − 在数学上，我们可以通过标准型来研究分岔的本质特征。对于一维动力系统，鞍结点分岔的最简微分方程形式（标准型）如下+ 第8行：
第8行： − + − + − 鞍结点分岔的关键性质在于，随着系统参数<math>r</math>的改变，相空间中原本存在的两个不动点（一个稳定，一个不稳定）相互靠近、碰撞并最终相互抵消而消失，导致系统原本的稳定状态不复存在。也可反过来看待，原本状态不稳定的系统因系统参数<math>r</math>的改变，使其先出现了一个半稳定不动点，最终出现了两个不动点（一个稳定，一个不稳定）。+ − + − + − + − + − + − + 第29行：
第29行： − 相图如下，展示了系统参数<math>r</math>在不同取值下的鞍结点分岔的三种不同状态，横坐标为<math>x</math>，纵坐标为<math>\dot{x}</math>。实心点为稳定不动点，空心点为不稳定不动点，半空心点为临界半稳定不动点，箭头为演化方向。+ − 分岔图如下，展示了系统参数<math>r</math>在不同取值下的不动点的分布情况。+ 第40行：
第40行： − 相图如下，展示了系统参数<math>r</math>在不同取值下的鞍结点分岔的三种不同状态，横坐标为<math>x</math>，纵坐标为<math>\dot{x}</math>。实心点为稳定不动点，空心点为不稳定不动点，半空心点为临界半稳定不动点，箭头为演化方向。+ − 分岔图如下，展示了系统参数<math>r</math>在不同取值下的不动点的分布情况。+ 第56行：
第56行： − + 第62行：
第62行： − + 第89行：
第89行： − + 第104行：
第104行： − + − + − + − 此时这块下降段区域的变化就是鞍结点分岔。+ − 相图如下+ − 与一维动力系统中的鞍结点分岔类似，二维鞍结点分岔同样描述了一个稳定不动点和一个不稳定不动点随参数变化相互靠近、碰撞并最终同时消失的过程。在环、无限周期分岔等其他非线性系统的参数讨论中也有出现鞍结点分岔的身影。+ 第125行：
第125行： − + − + − + − + − 三维以上的标准型与二维类似，仅在某个维度有分岔行为（类似<math>\dot x = \mu - x^2</math>），系统在其余维度上稳定（均为类似<math>\dot y=ay</math>的形式，其中参数小于零）。+ 第136行：
第136行： − x方向的不动点分析与一维情况完全类似。随着参数 <math>\mu</math> 的变化，二维相平面上的拓扑结构经历以下三个阶段：+ − + − 系统存在两个平衡点：+ − + − + − 动力学特征：相空间被不稳定不动点的分割，大部分轨线最终收敛于稳定不动点。+ − + − 两个不动点在原点 <math>(0,0)</math> 处发生碰撞并融合。此时雅可比矩阵具有一个零特征值（<math>\lambda_1=0</math>）和一个负特征值（<math>\lambda_2=-1</math>）。+ − + − 方程在实数域内无解，所有平衡点消失。+ − +
无编辑摘要
鞍结点[[分岔]]（Saddle-Node [[Bifurcation]]），又称折叠[[分岔]]（Fold [[Bifurcation]]）、切线[[分岔]]（Tangent [[Bifurcation]]）或蓝天[[分岔]]（Blue Sky [[Bifurcation]]），是[[非线性系统]]中最常见的不稳定性机制之一。鞍结点[[分岔]]常作为一种局部[[分岔]]（Local [[Bifurcation]]）出现。
=标准型和性质=
=标准型和性质=
在[[数学]]上，我们可以通过标准型来研究[[分岔]]的本质特征。对于一维[[动力系统]]，鞍结点[[分岔]]的最简[[微分方程]]形式（标准型）如下
{{NumBlk|:|
{{NumBlk|:|
<math display="block">
<math display="block">
</math>
</math>
|{{EquationRef|1}}}}
|{{EquationRef|1}}}}
* 其中<math>x \in \mathcal{R}</math>是状态变量，<math>r \in \mathcal{R}</math>为诱导分岔产生的外部参数。
* 其中<math>x \in \mathcal{R}</math>是[[状态变量]]，<math>r \in \mathcal{R}</math>为诱导[[分岔]]产生的外部参数。
* 对于实际研究中的更复杂的系统，如<math>\dot{x} = r - x^2 - e^x</math>，<math>x</math>的取值范围为应为不动点(<math>x=x^*</math>)的邻域，即<math>-\varepsilon <x<\varepsilon </math>，<math>r</math>的取值范围应为临界点<math>r=r_c</math>的邻域，即<math>-\delta <r<\delta </math>，否则无法忽略高阶项。
* 对于实际研究中的更复杂的系统，如<math>\dot{x} = r - x^2 - e^x</math>，<math>x</math>的取值范围为应为[[不动点]](<math>x=x^*</math>)的邻域，即<math>-\varepsilon <x<\varepsilon </math>，<math>r</math>的取值范围应为[[临界点]]<math>r=r_c</math>的邻域，即<math>-\delta <r<\delta </math>，否则无法忽略高阶项。
鞍结点[[分岔]]的关键性质在于，随着系统参数<math>r</math>的改变，[[相空间]]中原本存在的两个[[不动点]]（一个稳定，一个不稳定）相互靠近、碰撞并最终相互抵消而消失，导致系统原本的[[稳定状态]]不复存在。也可反过来看待，原本状态不稳定的系统因系统参数<math>r</math>的改变，使其先出现了一个半[[稳定不动点]]，最终出现了两个[[不动点]]（一个稳定，一个不稳定）。
==不动点和稳定性分析==
==[[不动点]]和稳定性分析==
以<math display="block">
以<math display="block">
\dot{x} = r - x^2
\dot{x} = r - x^2
</math>
</math>
为例
为例
===不动点判断===
===[[不动点]]判断===
令<math>f(x)=r-x^2</math>，令<math>f(x)=0</math>，此时<math>r</math>的范围在<math>\left [ 0,+ \infty \right ) </math>时，解<math>x^*=\pm \sqrt{r}</math>有实际意义。当<math>r<0</math>时，没有不动点；当<math>r=0</math>时，有一个不动点；当<math>r>0</math>时，有两个不动点。
令<math>f(x)=r-x^2</math>，令<math>f(x)=0</math>，此时<math>r</math>的范围在<math>\left [ 0,+ \infty \right ) </math>时，解<math>x^*=\pm \sqrt{r}</math>有实际意义。当<math>r<0</math>时，没有[[不动点]]；当<math>r=0</math>时，有一个[[不动点]]；当<math>r>0</math>时，有两个[[不动点]]。
===不动点稳定性分析===
===[[不动点]]稳定性分析===
<math>f'(x^*)=-2x^*</math>。当<math>r>0</math>时，<math>x^*=\sqrt{r}</math>，有<math>f'(x^*)<0</math>，是稳定不动点；<math>x^*=-\sqrt{r}</math>，有<math>f'(x^*)>0</math>，是不稳定不动点。<math>r=0</math>为分岔点，有<math>f'(x^*)=0</math>
<math>f'(x^*)=-2x^*</math>。当<math>r>0</math>时，<math>x^*=\sqrt{r}</math>，有<math>f'(x^*)<0</math>，是[[稳定不动点]]；<math>x^*=-\sqrt{r}</math>，有<math>f'(x^*)>0</math>，是不[[稳定不动点]]。<math>r=0</math>为[[分岔]]点，有<math>f'(x^*)=0</math>
==相图与分岔图==
==[[相图]]与[[分岔]]图==
对于<math display="block">
对于<math display="block">
</math>
</math>
[[相图]]如下，展示了系统参数<math>r</math>在不同取值下的鞍结点[[分岔]]的三种不同状态，横坐标为<math>x</math>，纵坐标为<math>\dot{x}</math>。实心点为[[稳定不动点]]，空心点为不[[稳定不动点]]，半空心点为[[临界]]半[[稳定不动点]]，箭头为演化[[方向]]。
[[File:Phase Portrait for Saddle-Node Bifurcation 1.png|800px]]
[[File:Phase Portrait for Saddle-Node Bifurcation 1.png|800px]]
[[分岔]]图如下，展示了系统参数<math>r</math>在不同取值下的[[不动点]]的分布情况。
[[File:Saddle-Node Bifurcation Diagram 01.png|600px]]
[[File:Saddle-Node Bifurcation Diagram 01.png|600px]]
\dot{x} = r - x^2
\dot{x} = r - x^2
</math>
</math>
[[相图]]如下，展示了系统参数<math>r</math>在不同取值下的鞍结点[[分岔]]的三种不同状态，横坐标为<math>x</math>，纵坐标为<math>\dot{x}</math>。实心点为[[稳定不动点]]，空心点为不[[稳定不动点]]，半空心点为[[临界]]半[[稳定不动点]]，箭头为演化[[方向]]。
[[File:Phase Portrait for Saddle-Node Bifurcation 2.png|800px|]]
[[File:Phase Portrait for Saddle-Node Bifurcation 2.png|800px|]]
[[分岔]]图如下，展示了系统参数<math>r</math>在不同取值下的[[不动点]]的分布情况。
[[File:Saddle-Node Bifurcation Diagram 02.png|600px]]
[[File:Saddle-Node Bifurcation Diagram 02.png|600px]]
可以看作是一阶参数与二阶自变量的线性组合。
可以看作是一阶参数与二阶自变量的线性组合。
===推导过程===
===推导过程===
将<math>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>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">
</math>
</math>
|{{EquationRef|1}}}}
|{{EquationRef|1}}}}
其中<math>f\left(x^{*}, r_{c}\right)=0</math>，且因为<math>x^*</math>是不动点，且由相切情况可得<math>\left. \frac {\partial f}{\partial x}\right|_{\left(x^{*} ; r_{c}\right)}=0</math>，再将<math>\left. \frac {\partial f}{\partial r}\right|_{\left(x ; r_{c}\right)}</math>视为常数<math>a</math>，将<math>\left. \frac {1}{2} \frac{\partial^{2} f}{\partial x^{2}}\right|_{\left(x ; r_{c}\right)}</math>视为常数<math>b</math>。忽略后面的高阶项，即可得到<math>\dot{x} = a(r-r_c)+b(x-x^*)^2</math>
其中<math>f\left(x^{*}, r_{c}\right)=0</math>，且因为<math>x^*</math>是[[不动点]]，且由相切情况可得<math>\left. \frac {\partial f}{\partial x}\right|_{\left(x^{*} ; r_{c}\right)}=0</math>，再将<math>\left. \frac {\partial f}{\partial r}\right|_{\left(x ; r_{c}\right)}</math>视为常数<math>a</math>，将<math>\left. \frac {1}{2} \frac{\partial^{2} f}{\partial x^{2}}\right|_{\left(x ; r_{c}\right)}</math>视为常数<math>b</math>。忽略后面的高阶项，即可得到<math>\dot{x} = a(r-r_c)+b(x-x^*)^2</math>
=案例：云杉蚜虫爆发模型=
=案例：云杉蚜虫爆发模型=
令<math>\frac{\mathrm{d} x}{\mathrm{d} \tau } =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>r>0</math>。此时<math>x^*=0</math>通常为不[[稳定不动点]]，因此分析实际问题时一般关注<math>x^*>0</math>的[[不动点]]的情况，即
<math display="block">
<math display="block">
因<math>k</math>是一个慢变量，所以将其视为一个常数（上图中<math>k=10</math>），将<math>r</math>（上图纵坐标）视为主要控制参数。
因<math>k</math>是一个慢变量，所以将其视为一个常数（上图中<math>k=10</math>），将<math>r</math>（上图纵坐标）视为主要控制参数。
* 当<math>r</math>很小时，可看到两条曲线仅有一个交点，且为一个稳定不动点。
* 当<math>r</math>很小时，可看到两条曲线仅有一个交点，且为一个[[稳定不动点]]。
* 令<math>r</math>逐渐增大到某个临界值<math>r_1</math>时，<math>y_1(x)</math>与<math>y_2(x)</math>的下降段相切，出现了一个半稳定不动点。
* 令<math>r</math>逐渐增大到某个[[临界]]值<math>r_1</math>时，<math>y_1(x)</math>与<math>y_2(x)</math>的下降段相切，出现了一个半[[稳定不动点]]。
* 继续增大<math>r</math>，<math>y_1(x)</math>与<math>y_2(x)</math>的下降段出现了两个不动点（一个不稳定，一个稳定）。
* 继续增大<math>r</math>，<math>y_1(x)</math>与<math>y_2(x)</math>的下降段出现了两个[[不动点]]（一个不稳定，一个稳定）。
此时这块下降段区域的变化就是鞍结点[[分岔]]。
[[相图]]如下
[[File:Saddle Node case phase portrait.png|600px]]
[[File:Saddle Node case phase portrait.png|600px]]
=二维系统（以上）=
=二维系统（以上）=
与一维[[动力系统]]中的鞍结点[[分岔]]类似，二维鞍结点[[分岔]]同样描述了一个[[稳定不动点]]和一个不[[稳定不动点]]随参数变化相互靠近、碰撞并最终同时消失的过程。在[[环]]、[[无限周期分岔]]等其他[[非线性系统]]的参数讨论中也有出现鞍结点[[分岔]]的身影。
==标准型==
==标准型==
</math>
</math>
* <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>的形式，其中参数小于零）。
==稳定性分析==
==稳定性分析==
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> [[方向]]排斥。
[[动力学]]特征：[[相空间]]被不[[稳定不动点]]的分割，大部分轨线最终收敛于[[稳定不动点]]。
2. 临界分岔点 (<math>\mu = 0</math>)：
2. [[临界]][[分岔]]点 (<math>\mu = 0</math>)：
两个[[不动点]]在原点 <math>(0,0)</math> 处发生碰撞并融合。此时[[雅可比矩阵]]具有一个零[[特征值]]（<math>\lambda_1=0</math>）和一个负[[特征值]]（<math>\lambda_2=-1</math>）。
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>）。[[系统状态]]在穿过该区域时会发生显著的[[临界]][[迟滞]]现象，形成所谓的“瓶颈效应”。
=参考文献=
=参考文献=
<references />
<references />