吴锋

华南理工大学数学学院，广东 广州

收稿日期：2020年1月25日；录用日期：2020年2月14日；发布日期：2020年2月21日

本文基于segmented disc dynamo发电机，提出了一个新的四维混沌的广义发电机系统。作者研究了Pitchfork分叉、Hopf分叉、zero-zero Hopf分叉，给出了在不同条件下系统存在pitchfork分叉、Hopf分叉和周期解的证明，并给出了Pitchfork分叉和Hopf分叉的数值模拟。

关键词 :发电机模型，超混沌系统，Pitchfork分叉，Hopf分叉，Zero-Zero Hopf分叉

Copyright © 2020 by author(s) and beplay安卓登录

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

非线性科学自二十世纪下半叶以来，作为一门基础学科取得了迅猛的发展，尤其是其中的混沌现象更是成为当今科学前沿的热门话题，自1963年Lorenz引入第一个数学物理混沌模型 [ 1 ]，在那之后各个领域的科学家开始深入探索混沌的本质、路径、分叉和混沌系统的特性 [ 2 ] - [ 10 ]，在此基础上1979年，O. E. Rossler首次提出了超混沌的概念，并给出了具有一个以上正Lyapunov指数的四维超混沌系统，超混沌作为一种比混沌更复杂的动力学行为，在工程应用上又比混沌具有更大的应用潜力；众所周知，自治系统中产生超混沌的最小可能维数系统是4；因此，四维自治系统是超混沌研究和应用的主要方向，为了帮助理解天体物理间磁场的产生和逆转，过去广泛运用发电机模型 [ 11 ] - [ 15 ] 来进行研究，后来Moffa [ 16 ] 提出了一种新型发电机模型，该发电机模型可用以下常微分方程组来描述：

X ˙ = r ( Y − X ) (1)

Y ˙ = m X − ( 1 + m ) Y + X Z (2)

Z ˙ = g [ 1 + m X 2 − ( 1 + m ) X Y ] (3)

这里X(t)和Y(t)表示由于径向和方位角电流分布引起的磁通量，Z(t)表示圆盘的角速度，圆点表示相对于时间的微分，g表示施加的转矩，r和m为正常数，取决于电路的电性能。Moffat的发电机忽略了摩擦对转盘的影响，然而实际情况还要考虑轴承的摩擦，我们可以在(3)式后加上项“–vZ”来模拟这种摩擦，为方便，Moffat进行坐标变换：

X ¯ ′ = Y ¯ − X ¯ (5)

Y ¯ ′ = α ( X ¯ − Y ¯ ) + γ X ¯ + X ¯ Z ¯ (6)

Z ¯ ′ = − n Z ¯ + X ¯ ( X ¯ − β Y ¯ ) (7)

则方程(1)~(3)变为

其中， n = v / r , α = ( 1 + m ) / r , β = 1 + 1 / m , γ = ( g / v − 1 ) / r 。Moffat利用此方程组探究混沌的特性。

Pitchfork分叉和Hopf分叉是改变参数时从系统平衡点出现的两种分叉现象。局部分叉的研究有助于我们分析和阐明非线性系统如何从稳态演化到混沌和超混沌状态。因此，分叉点在研究平衡点的稳定性方面起到至关重要的作用。当平衡点的个数和稳定性发生变化时，且零点为非双曲点时，可能出现pitchfork分叉；当平衡点通过一对纯虚特征值来改变稳定性时，Hopf分叉将在动力系统平衡态的极限环中产生；当平衡点在特征值表现为一对纯虚特征值和零特征值时，将可能出现0-0 Hopf分叉。

本文主要研究的是基于Moffat提出的发电系统，对其广义化的四维新型超混沌系统及其基本性质进行研究，得到广义发电系统下pitchfork分叉、零Hopf分叉以及周期解在零零-Hopf平衡下的分叉问题。

令d为分叉参数，当 时，特征方程在奇点 处有特征值 ， ， ，

因此零点 是非双曲的，由中心流形定理和分叉理论，经过推导，可得Pitchfork分叉的以下定理。

定理1当 时，系统(8)在零平衡点 出现了一个Pitchfork分叉，而且当 时，仅有一个不稳定平衡点 在 的左侧；当 时，将出现三个平衡点， 将变成稳定的，另外两个不稳定的平衡点 在 的右侧。

证明：相关特征向量为 ， ， ，

根据 [ 17 ]，令

令 且 (13)

考虑变换

利用方程(13)、(14)将系统(8)转换为

其中

, ,

, ;

,

,

,

由中心流形定理可知，对于方程(15)存在一个中心流形，通过变量 和 [ 17 ] 可将其表示为

其中 和 充分小且 是 的全微分 ；要获得中心流形上的向量域，设

其中 表示为所有三阶项(例如 )和更高阶项(例如 )

又因为 ，则不变流形应满足

其中， ，

将方程(16)代入(17)，比较相同项的系数，我们可得到

因此方程(16)可改写为

将方程组(19)代入方程组(15)的 ，将向量场简化为中心流形，我们可得

其中

由平衡点分叉理论 [ 18 ]，容易证得以下条件是成立的，这说明系统(20)的平衡点 在 处发生了一个pitchfork分叉。

因此定理1得证

数值模拟为了验证定理1，设参数 根据定理1知道 ，系统(8)在 点产生了Pitchfork分叉，这说明在 附近平衡点个数和稳定性

发生了变化。在 左侧仅有一个不稳定的平衡点 ，在它的右侧有两个不稳定的平衡点 和一个稳定的平衡点 。图2所示为 在35~48之间的分叉图。显然，由图2得知分叉方向在 处，且 时 是不稳定的，当 时， 变为稳定的且 出现为不稳定的；这与定理1结论相同。

图2. 系统(8)在d = d⁰= 42处的Pitchfork分叉

令d分叉参数，在研究系统(8)的平衡点 的Hopf分叉之前，先介绍Hopf分叉的基本理论：

考虑一个向量场

若向量场将发生Hopf分叉，则要同时满足以下条件 [ 19 ] - [ 26 ] ：

a) 满足 ，雅各比矩阵 有一对纯虚根特征值 和 ，其余 个特征值实部非零。

b) 横截性条件：特征值实部满足 。

接下来，对Hopf分叉的两个条件逐一验证。

首先，由方程(11)通过计算得知当

系统(8)在 点的雅各比矩阵线性化后有一对纯虚根特征值和两个非零实特征值

其中 (23)

因此，条件a)成立。

其次，验证条件b)，已知 为方程(12)的根，对 中的d求偏导，得

则

横截条件b)满足。因此系统(8)在 处发生Hopf分岔，而分岔点的稳定性由第一个Lyapunov系数 所决定，故可得以下定理。

定理2对于系统(8)，( ， ， )，该系统在 处的第一个Lyapunov系数 可表示为

1) 如果 ，系统(8)在 处发生Hopf分岔，分岔点是不稳定的，而且当 ，靠近 时，系统在稳定奇点 附近存在一个不稳定的极限环。

2) 如果 ，系统(8)在 处发生Hopf分岔，分岔点是稳定的，而且当 ，靠近 时，系统在不稳定奇点 附近存在一个稳定的极限环。

证明：由(10)式的雅各比矩阵可知，当 时，

它的特征值为 ，

已知 ，

,

, ,

,

则根据公式

即可计算得到第一个Lyapunov系数 如定理中所示。

数值模拟 取 ，则 ， ，如图3所示，此时系统(8)得到一个初始条件为 的不稳定极限环。

图3. 系统(8)的初始条件 在满足参数 下的不稳定极限环

对于(8)可以考虑变形为以下微分形式：

其中 且 是充分小的正数， 是二阶连续偏导函数且各自定义域值域为 ；第一个变量为周期T， 是 上的一个开集。根据平均理论，研究系统(26)中T周期解的分叉问题，其主要是假设无干扰系统

有一个子流形。假设 是系统(27)在 处的周期解，那么其线性周期解 可以记为

令 为线性微分系统(28)的基本矩阵。假设存在一个开集V满足 ，使得对于每一个 ，都有一个系统(27)的周期解 。其中集合 是系统(27)的周期解集合，也被称之为系统(27)的同步集。因此，关于T周期解的分叉问题周期解的 全包含在集合 中，这将由以下定理给出。

定理3 [ 26 ] 假设有一个有界开集V满足 使得对于每一个 的解 是T周期的，则以下函数

满足陈述：

1) 若存在 满足 且 ，则存在一个系统(26)的周期解 使得当 时有

2) 周期解 的稳定性类型由雅各比矩阵 的特征值给出。

该定理证明过于繁琐，具体可见 [ 26 ]。

定理4令 且 ，则坐标原点为广义发电系统(8)的0-0-Hopf平衡点。

证明：易计算得系统(8)在 的特征方程为

当参数满足 ， 时，该方程的根依次为 ，定理4得证。

在求出系统(8)在原点处有0-0-Hopf分叉后，接下来则需要证明在该平衡处是否有周期解的存在，定理5给出了周期解的存在性和稳定性的满足条件。

定理5 令 且

当 时，广义发电系统(8)在原点平衡处有一个0-0-Hopf分叉，在这个平衡处有两个周期解，且当 时它们是稳定解；其中 具体值可参见附录所示。

证明：经过变形 ，其中 是充分小的正数，

是非零实数，则广义发电系统(8)变为

此外，进一步缩放变量，令 ，然后再次用 表示 ，则系统(29)变形为

借助定理3中描述的平均方法，对系统(30)进行研究；

考虑无扰动系统，

研究(31)的初值问题，得出其解为 满足

其中 ，

,

,

,

,

,

,

,

,

值得注意的是，系统(31)所有满足 解 都是周期的，含有相同周期 ，对于系统(31)的解 ，可以将线性化系统 的基本矩阵 写成如下形式：

其中 ，计算定理3中所建立的函数

其中 数值可见附录。

解非线性方程 ，得到三个解：

解 对应于在原点的平衡 ，对于另外两个解 满足：

因此，根据定理3，系统(30)存在两个周期为 的解 和 使得当 时有 和 成立。因为系统(30)是系统(29)的变形，所以系统(30)的周期解 对应于系统(29)的周期解 。

最后为了确定这两个周期解的稳定性类型，需要计算雅各比矩阵 ，不难得出矩阵 的特征多项式为

其中 数值可见附录。

令

根据Routh-Hurwitz准则，在 和 情况下，该特中多项式(32)的所有根均为负实部，即两个周期解 和 在原点生成是稳定的，定理5得证。

吴 锋. 一个四维混沌的广义发电机系统的动力学分析Dynamic Analysis of a 4D Hyperchaotic Generalized Segmented Disc Dynamo System[J]. 理论数学, 2020, 10(02): 111-127. https://doi.org/10.12677/PM.2020.102018