| 第4行: |
第4行: |
| | 对于二维线性动力系统<math>\dot{\mathbf{x}}=A\mathbf{x}</math>,状态向量可写为<math>\mathbf{x}=(x,y)^T</math>。在相平面中,每一个状态点<math>(x,y)</math>都通过矩阵 A 对应一个速度向量<math>(\dot{x},\dot{y})^T</math>,表示系统在该点处的瞬时变化方向和变化速率。所有状态点上的速度向量共同构成系统的<nowiki>'''</nowiki>向量场<nowiki>'''</nowiki>。向量场直观展示了系统在二维状态空间中的流动趋势,是分析相轨迹、平衡点和稳定性的基础。 | | 对于二维线性动力系统<math>\dot{\mathbf{x}}=A\mathbf{x}</math>,状态向量可写为<math>\mathbf{x}=(x,y)^T</math>。在相平面中,每一个状态点<math>(x,y)</math>都通过矩阵 A 对应一个速度向量<math>(\dot{x},\dot{y})^T</math>,表示系统在该点处的瞬时变化方向和变化速率。所有状态点上的速度向量共同构成系统的<nowiki>'''</nowiki>向量场<nowiki>'''</nowiki>。向量场直观展示了系统在二维状态空间中的流动趋势,是分析相轨迹、平衡点和稳定性的基础。 |
| | ===相轨迹的几何形态=== | | ===相轨迹的几何形态=== |
| − | 在二维线性动力系统中,若两个状态变量分别表示位置和速度,则系统状态可写为 \mathbf{x}(t)=(x(t),v(t))^T,对应的二维状态空间称为<nowiki>'''</nowiki>相平面<nowiki>'''</nowiki>。给定初始状态 \mathbf{x}_0=(x_0,v_0)^T 后,状态点会随时间 t 在相平面中运动,其描绘出的曲线称为<nowiki>'''</nowiki>相轨迹<nowiki>'''</nowiki>。相轨迹用于直观表示系统从初始状态出发后的运动路径,并可用于判断系统是否趋于平衡、发散或呈现振荡行为。 | + | 在二维线性动力系统中,若两个状态变量分别表示位置和速度,则系统状态可写为<math>\mathbf{x}(t)=(x(t),v(t))^T</math>,对应的二维状态空间称为<nowiki>'''</nowiki>相平面<nowiki>'''</nowiki>。给定初始状态<math>\mathbf{x}_0=(x_0,v_0)^T</math>后,状态点会随时间 t 在相平面中运动,其描绘出的曲线称为<nowiki>'''</nowiki>相轨迹<nowiki>'''</nowiki>。相轨迹用于直观表示系统从初始状态出发后的运动路径,并可用于判断系统是否趋于平衡、发散或呈现振荡行为。 |
| | ===相图的定性分析=== | | ===相图的定性分析=== |
| − | 在二维线性动力系统中,研究重点不仅在于求出解析解 \mathbf{x}(t),也在于通过相平面中的整体轨迹分布判断系统的定性行为。通过观察相轨迹的几何形状、流动方向以及其与平衡点的关系,可以在不显式求解微分方程的情况下,分析系统的稳定性、振荡性和长期演化趋势。例如,无阻尼弹簧振子在相平面中的轨迹通常表现为围绕平衡点的闭合曲线,反映了系统能量守恒和周期振荡的特征。 | + | 在二维线性动力系统中,研究重点不仅在于求出解析解<math>\mathbf{x}(t)</math>,也在于通过相平面中的整体轨迹分布判断系统的定性行为。通过观察相轨迹的几何形状、流动方向以及其与平衡点的关系,可以在不显式求解微分方程的情况下,分析系统的稳定性、振荡性和长期演化趋势。例如,无阻尼弹簧振子在相平面中的轨迹通常表现为围绕平衡点的闭合曲线,反映了系统能量守恒和周期振荡的特征。 |
| | ===平衡点附近的相平面结构=== | | ===平衡点附近的相平面结构=== |
| − | 在二维线性动力系统中,满足 A\mathbf{x}=\mathbf{0} 的状态称为系统的<nowiki>'''</nowiki>平衡点<nowiki>'''</nowiki>或<nowiki>'''</nowiki>不动点<nowiki>'''</nowiki>。对于连续时间系统 \dot{\mathbf{x}}=A\mathbf{x},若系统矩阵 A 非奇异,则原点 (0,0) 是唯一平衡点;若 A 奇异,则可能存在一条直线或更高维子空间上的平衡点。平衡点附近相轨迹的走向决定系统的局部稳定性,例如轨迹趋向平衡点表示稳定,远离平衡点表示不稳定,而围绕平衡点旋转则可能对应振荡行为。 | + | 在二维线性动力系统中,满足 A\mathbf{x}=\mathbf{0} 的状态称为系统的<nowiki>'''</nowiki>平衡点<nowiki>'''</nowiki>或<nowiki>'''</nowiki>不动点<nowiki>'''</nowiki>。对于连续时间系统<math>\dot{\mathbf{x}}=A\mathbf{x}</math>,若系统矩阵 A 非奇异,则原点 (0,0) 是唯一平衡点;若 A 奇异,则可能存在一条直线或更高维子空间上的平衡点。平衡点附近相轨迹的走向决定系统的局部稳定性,例如轨迹趋向平衡点表示稳定,远离平衡点表示不稳定,而围绕平衡点旋转则可能对应振荡行为。 |
| | ===基于迹—行列式平面的分类=== | | ===基于迹—行列式平面的分类=== |
| − | 对于二维线性系统,局部几何结构可以由矩阵 [math]\displaystyle{ A }[/math] 的特征值完全刻画。设矩阵 [math]\displaystyle{ A }[/math] 的两个特征值为 [math]\displaystyle{ \lambda_1,\lambda_2 }[/math],其迹(Trace)与行列式(Determinant)分别为: | + | 对于二维线性系统,局部几何结构可以由矩阵<math>\displaystyle{ A }</math>的特征值完全刻画。设矩阵<math>\displaystyle{ A }</math>的两个特征值为<math>\displaystyle{ \lambda_1,\lambda_2 }</math>,其迹(Trace)与行列式(Determinant)分别为: |
| | | | |
| − | [math]\displaystyle{ T=\lambda_1+\lambda_2,\quad \Delta=\lambda_1\lambda_2 }[/math]
| + | <math>\displaystyle{ T=\lambda_1+\lambda_2,\quad \Delta=\lambda_1\lambda_2 }</math> |
| | | | |
| − | 根据特征方程 [math]\displaystyle{ \lambda^2-T\lambda+\Delta=0 }[/math],由求根公式可得: | + | 根据特征方程<math>\displaystyle{ \lambda^2-T\lambda+\Delta=0 }</math>,由求根公式可得: |
| | | | |
| − | <nowiki>[math]\displaystyle{ \lambda_{1,2}=\frac{T\pm\sqrt{T^2-4\Delta}}{2} }[/math]</nowiki> | + | <math>\displaystyle{ \lambda_{1,2}=\frac{T\pm\sqrt{T^2-4\Delta}}{2} }</math> |
| | | | |
| − | 因此,二维线性系统不动点的类型可统一由 [math]\displaystyle{ T }[/math]、[math]\displaystyle{ \Delta }[/math] 以及判别式 [math]\displaystyle{ T^2-4\Delta }[/math] 来决定。 | + | 因此,二维线性系统不动点的类型可统一由<math>\displaystyle{ T }</math>、<math>\displaystyle{ \Delta }</math>以及判别式<math>\displaystyle{ T^2-4\Delta }</math>来决定。 |
| | | | |
| | ==经典实例:二维谐振子系统<ref name=":4" />== | | ==经典实例:二维谐振子系统<ref name=":4" />== |
| − | 为了更直观地理解线性系统的降维与重构,可以以经典的无阻尼弹簧振子为例<ref>J. R. Taylor, Classical Mechanics. Sausalito, CA, USA: University Science Books, 2005.</ref>。设一个质量为 [math]\displaystyle{ m }[/math] 的物体连接在劲度系数为 [math]\displaystyle{ k }[/math] 的弹簧上,其偏离平衡位置的位移记为 [math]\displaystyle{ x }[/math]。根据牛顿第二定律与胡克定律,该系统满足以下二阶常微分方程: | + | 为了更直观地理解线性系统的降维与重构,可以以经典的无阻尼弹簧振子为例<ref>J. R. Taylor, Classical Mechanics. Sausalito, CA, USA: University Science Books, 2005.</ref>。设一个质量为<math>\displaystyle{ m }</math>的物体连接在劲度系数为<math>\displaystyle{ k }</math>的弹簧上,其偏离平衡位置的位移记为<math>\displaystyle{ x }</math>。根据牛顿第二定律与胡克定律,该系统满足以下二阶常微分方程: |
| | | | |
| − | [math]\displaystyle{ m\ddot{x}+kx=0 }[/math]
| + | <math>\displaystyle{ m\ddot{x}+kx=0 }</math> |
| | | | |
| − | 为了将其转化为动力系统分析中标准的一阶方程组形式,引入速度变量 [math]\displaystyle{ v=\dot{x} }[/math]。此时,原方程可降阶并改写为: | + | 为了将其转化为动力系统分析中标准的一阶方程组形式,引入速度变量<math>\displaystyle{ v=\dot{x} }</math>。此时,原方程可降阶并改写为: |
| | | | |
| − | [math]\displaystyle{ \dot{v}=-\frac{k}{m}x }[/math]
| + | <math>\displaystyle{ \dot{v}=-\frac{k}{m}x }</math> |
| | | | |
| | 联立上述两式,即可得到该二维线性系统的矩阵表达形式: | | 联立上述两式,即可得到该二维线性系统的矩阵表达形式: |
| | | | |
| − | [math]\displaystyle{ \begin{pmatrix} \dot{x}\\ \dot{v} \end{pmatrix}=\begin{pmatrix} 0 & 1\\ -\frac{k}{m} & 0 \end{pmatrix}\begin{pmatrix} x\\ v \end{pmatrix} }[/math]
| + | <math>\displaystyle{ \begin{pmatrix} \dot{x}\\ \dot{v} \end{pmatrix}=\begin{pmatrix} 0 & 1\\ -\frac{k}{m} & 0 \end{pmatrix}\begin{pmatrix} x\\ v \end{pmatrix} }</math> |
| | | | |
| − | 从该演化方程可以看出,要完全确定系统未来的状态,仅仅知道位置 [math]\displaystyle{ x }[/math] 这一物理量是不够的,还必须掌握系统当下的速度 [math]\displaystyle{ v }[/math]。因此,描述该系统的自然状态空间不再是一维的直线,而是由位置与速度构成的二维<nowiki>'''</nowiki>相平面<nowiki>'''</nowiki>(Phase plane),即 [math]\displaystyle{ (x,v) }[/math] 空间。 | + | 从该演化方程可以看出,要完全确定系统未来的状态,仅仅知道位置<math>\displaystyle{ x }</math>这一物理量是不够的,还必须掌握系统当下的速度<math>\displaystyle{ v }</math>。因此,描述该系统的自然状态空间不再是一维的直线,而是由位置与速度构成的二维<nowiki>'''</nowiki>相平面<nowiki>'''</nowiki>(Phase plane),即<math>\displaystyle{ (x,v) }</math>空间。 |
| | ==跨学科案例:罗密欧与朱丽叶的恋爱动力学<ref name=":4" />== | | ==跨学科案例:罗密欧与朱丽叶的恋爱动力学<ref name=":4" />== |
| | | | |
| 第40行: |
第40行: |
| | 二维线性系统不仅可以描述机械振动、电路演化等物理过程,也可以作为抽象模型刻画个体之间的相互作用。为了说明线性系统分类在不同学科中的普适性,非线性动力学研究者 Steven Strogatz 引入了一个经典的“恋爱动力学”(Love Dynamics)模型,用以描述两个人情感状态之间的线性反馈关系。 | | 二维线性系统不仅可以描述机械振动、电路演化等物理过程,也可以作为抽象模型刻画个体之间的相互作用。为了说明线性系统分类在不同学科中的普适性,非线性动力学研究者 Steven Strogatz 引入了一个经典的“恋爱动力学”(Love Dynamics)模型,用以描述两个人情感状态之间的线性反馈关系。 |
| | | | |
| − | 设 [math]\displaystyle{ R(t) }[/math] 表示时刻 [math]\displaystyle{ t }[/math] 罗密欧对朱丽叶的情感强度,[math]\displaystyle{ J(t) }[/math] 表示朱丽叶对罗密欧的情感强度,其中正值表示爱慕,负值表示厌恶。则两人的情感互动可表示为如下二维线性系统: | + | 设<math>\displaystyle{ R(t) }</math>表示时刻<math>\displaystyle{ t }</math>罗密欧对朱丽叶的情感强度,<math>\displaystyle{ J(t) }</math>表示朱丽叶对罗密欧的情感强度,其中正值表示爱慕,负值表示厌恶。则两人的情感互动可表示为如下二维线性系统: |
| | | | |
| − | [math]\displaystyle{ \dot{R}=aR+bJ }[/math]
| + | <math>\displaystyle{ \dot{R}=aR+bJ }</math> |
| | | | |
| − | [math]\displaystyle{ \dot{J}=cR+dJ }[/math]
| + | <math>\displaystyle{ \dot{J}=cR+dJ }</math> |
| | | | |
| | 写成矩阵形式即: | | 写成矩阵形式即: |
| | | | |
| − | [math]\displaystyle{ \begin{pmatrix} \dot{R}\\ \dot{J} \end{pmatrix}=\begin{pmatrix} a & b\\ c & d \end{pmatrix}\begin{pmatrix} R\\ J \end{pmatrix} }[/math]
| + | <math>\displaystyle{ \begin{pmatrix} \dot{R}\\ \dot{J} \end{pmatrix}=\begin{pmatrix} a & b\\ c & d \end{pmatrix}\begin{pmatrix} R\\ J \end{pmatrix} }</math> |
| | | | |
| − | 其中,参数 [math]\displaystyle{ a }[/math] 和 [math]\displaystyle{ d }[/math] 描述个体对自身当前情感状态的内在反应;参数 [math]\displaystyle{ b }[/math] 和 [math]\displaystyle{ c }[/math] 描述双方对彼此情感变化的外部响应方式。不同参数符号和大小的组合,对应于不同的互动模式和长期演化结果。这一模型说明,即使系统变量来自社会或心理意义上的抽象量,其数学分析仍然可以采用与物理系统相同的线性动力学框架。 | + | 其中,参数 <math>\displaystyle{ a }</math>和<math>\displaystyle{ d }</math>描述个体对自身当前情感状态的内在反应;参数<math>\displaystyle{ b }</math>和<math>\displaystyle{ c }</math>描述双方对彼此情感变化的外部响应方式。不同参数符号和大小的组合,对应于不同的互动模式和长期演化结果。这一模型说明,即使系统变量来自社会或心理意义上的抽象量,其数学分析仍然可以采用与物理系统相同的线性动力学框架。 |
| | ===特例分析:相爱相杀的循环=== | | ===特例分析:相爱相杀的循环=== |
| | 作为一个典型情形,考虑如下系统: | | 作为一个典型情形,考虑如下系统: |
| | | | |
| − | [math]\displaystyle{ \dot{R}=bJ,\quad \dot{J}=cR }[/math]
| + | <math>\displaystyle{ \dot{R}=bJ,\quad \dot{J}=cR }</math> |
| | | | |
| − | 其中,[math]\displaystyle{ b>0 }[/math],[math]\displaystyle{ c<0 }[/math]。这表示罗密欧主要依据朱丽叶的态度来调整自己:朱丽叶越爱他,他越热情;而朱丽叶则对罗密欧的接近表现出回避倾向:罗密欧越热情,她越倾向于后退。此时系统矩阵为: | + | 其中,<math>\displaystyle{ b>0 }[/math],[math]\displaystyle{ c<0 }</math>。这表示罗密欧主要依据朱丽叶的态度来调整自己:朱丽叶越爱他,他越热情;而朱丽叶则对罗密欧的接近表现出回避倾向:罗密欧越热情,她越倾向于后退。此时系统矩阵为: |
| | | | |
| − | [math]\displaystyle{ A=\begin{pmatrix} 0 & b\\ c & 0 \end{pmatrix} }[/math]
| + | <math>\displaystyle{ A=\begin{pmatrix} 0 & b\\ c & 0 \end{pmatrix} }</math> |
| | | | |
| | 计算其迹与行列式可得: | | 计算其迹与行列式可得: |
| | | | |
| − | [math]\displaystyle{ T=0,\quad \Delta=-bc>0 }[/math]
| + | <math>\displaystyle{ T=0,\quad \Delta=-bc>0 }</math> |
| | | | |
| | 根据二维线性系统的不动点分类,该系统的特征值为一对纯虚数,原点对应于一个<nowiki>'''中心'''</nowiki>(Center)。在相平面上,系统轨迹表现为围绕原点的闭合曲线。这意味着两人的情感状态既不会收敛到静止平衡,也不会无限发散,而是呈现出周期性的循环变化。 | | 根据二维线性系统的不动点分类,该系统的特征值为一对纯虚数,原点对应于一个<nowiki>'''中心'''</nowiki>(Center)。在相平面上,系统轨迹表现为围绕原点的闭合曲线。这意味着两人的情感状态既不会收敛到静止平衡,也不会无限发散,而是呈现出周期性的循环变化。 |