新光滑多翅膀混沌系统研究

李文静，贾美美; Li Wenjing; Jia Meimei

网刊加载中。。。

使用Chrome浏览器效果最佳，继续浏览，你可能不会看到最佳的展示效果，

确定继续浏览么?

复制成功，请在其他浏览器进行阅读

新光滑多翅膀混沌系统研究

- ORCID：
李文静 ^1,2
- ORCID：
贾美美 ^1,2
✉

1. 内蒙古工业大学电力学院，呼和浩特 010080； 2. 内蒙古自治区机电控制重点实验室，呼和浩特 010051

最近更新：2021-08-30

DOI：10.6052/1672-6553-2020-072

目录contents

摘要

针对光滑多翅膀混沌系统不易于构造的问题，基于光滑代数函数提出两个新光滑多分段非线性函数，并用新光滑多分段非线性函数构造了新光滑单方向多翅膀S-M混沌系统和新光滑二方向网格多翅膀S-M混沌系统.然后，分析了两个新光滑多翅膀混沌系统的非线性动力学行为，主要包括相图、平衡点、不变性、耗散性、李亚普诺夫指数、分维数和庞加莱截面，并得到了多翅膀混沌吸引子的产生机理.其产生机理为指标2鞍焦平衡点用于产生翅膀.此外，通过单方向多翅膀S-M混沌吸引子和二方向网格多翅膀S-M混沌吸引子的电路实现，验证了理论推导和数值仿真的正确性.

关键词

非线性动力学行为; 多翅膀混沌系统; 代数函数; 多分段非线性函数; 电路实现

2020-04-19 收到第1稿，2020-06-30 收到修改稿.

引言

1963年，Lorenz提出了Lorenz混沌系统，其混沌吸引子的拓扑结构与蝴蝶的形状十分相似，是一种双翅膀混沌系统，该系统的提出掀起了混沌研究的热潮^［

1］.混沌在保密通信^{［参考文献 2-5}2-5］、图像加密^{［参考文献 6-8}6-8］等领域具有良好的应用前景.为了将混沌更好地应用于保密通信，人们长期致力于研究具有复杂拓扑结构和丰富动力学行为的混沌系统，例如单方向多翅膀混沌系统和多方向网格多翅膀混沌系统.由于单方向多翅膀混沌系统和多方向网格多翅膀混沌系统具有多个可调参数，即数量众多的翅膀能在相空间中呈现单个方向分布甚至多个方向分布的网格状图案，各个翅膀之间具有互相嵌套的拓扑结构，翅膀的数量、大小及网格状分布的图像可以由可调参数来进行调节，具有更复杂的动力学行为.这种复杂性体现在混沌吸引子的相轨迹或状态变量的取值能在多个翅膀之间随机的跳变，从而使得当翅膀的数量越多时，这种跳变的随机性就越大，更有利于混沌加密^{［参考文献 9

百度学术}9，10］.因此，基于双翅膀混沌系统构造多翅膀混沌系统具有重要意义.

2008年，Yu等^［

11］采用符号函数基于广义双翅膀Shimizu-Morioka（S-M）混沌系统构造了单方向多翅膀混沌系统和二方向网格多翅膀混沌系统.2010年，Yu等^{［参考文献 12

百度学术}12］采用符号函数基于广义双翅膀S-M混沌系统、双翅膀Lü混沌系统、双翅膀Chen混沌系统、双翅膀Sprott-1混沌系统和双翅膀Sprott-2混沌系统构造了相应的单方向多翅膀混沌系统.2012年，周欣等^{［参考文献 13

百度学术}13］采用符号函数基于广义双翅膀Lorenz混沌系统构造了二方向网格多翅膀混沌系统.2015年，Huang等^{［参考文献 14

百度学术}14］采用分段时滞函数基于双翅膀Lorenz混沌系统构造了二方向网格多翅膀混沌系统.2016年，孙克辉等^{［参考文献 15

百度学术}15］采用符号函数基于简化双翅膀Lorenz混沌系统构造了二方向网格多翅膀混沌系统.2018年，Zhang等^{［参考文献 16

百度学术}16］采用锯齿波函数和阶梯函数基于提出的双翅膀混沌系统构造了二方向网格多翅膀混沌系统.2018年，Yu等^{［参考文献 17

百度学术}17］采用阶梯切换函数基于双翅膀Lorenz混沌系统构造了三方向网格多翅膀混沌系统.

2016年，Zhou等^［

18］采用符号函数与光滑磁控忆阻器基于改进的双翅膀Lü混沌系统构造了单方向多翅膀混沌及超混沌系统.文献［19-21］也报道了如何采用忆阻器来构造多翅膀混沌及超混沌系统.文献［11-17］采用的符号函数、分段时滞函数、锯齿波函数、阶梯函数、阶梯切换函数都是非光滑函数，其构造的多翅膀混沌或超混沌系统都是非光滑系统.然而，构造光滑多翅膀混沌系统和超混沌系统更加具有挑战性^{［参考文献 19

百度学术}19］.那么，是否可以采用光滑非线性函数基于光滑双翅膀混沌系统构造光滑多翅膀混沌系统？本文将给出明确的答案.

自1980年，Shimizu和Morioka^［

22］提出双翅膀S-M混沌系统以来，目前鲜有文献报道基于双翅膀S-M混沌系统构造光滑多翅膀S-M混沌系统问题.需要注意的是：双翅膀S-M混沌系统是光滑系统；本文采用的代数函数是光滑非线性函数；本文构造的新单方向多翅膀S-M混沌系统和新二方向网格多翅膀S-M混沌系统都是光滑系统.

本文首先以双翅膀S-M混沌系统为基础，基于代数函数提出了两个新多分段非线性函数 $f_{N} (x)$ 和 $f_{M} (z)$ .然后，通过采用非线性函数 $f_{N} (x)$ 代替双翅膀S-M混沌系统中的平方项，构造出新单方向多翅膀S-M混沌系统.通过采用多分段非线性函数 $f_{M} (z)$ 代替新单方向多翅膀S-M混沌系统中的状态变量z，构造出新二方向网格多翅膀S-M混沌系统.随后，分析了新单方向多翅膀S-M混沌系统和新二方向网格多翅膀S-M混沌系统的非线性动力学行为，并得到了新单方向多翅膀S-M混沌吸引子和新二方向网格多翅膀S-M混沌吸引子的产生机理.其产生机理为指标2鞍焦平衡点用于产生翅膀.同时，对上述系统进行了电路实现.

1 单方向多翅膀S-M混沌系统的构造

首先引入代数函数，如式（1）所示

g (x) = \frac{x}{\sqrt[]{d + x^{2}}}

（1）

随着参数值 $d$ 的变化，g（x）∈［-1，1］.当 $d = 1$ 时，代数函数 $g (x)$ 变为具体的代数函数^［

23］.图1表示参数

d

取不同值时代数函数

g (x)

的曲线.基于代数函数

g (x)

，构造多分段非线性函数

f_{N} (x)

，可表示为

图1 代数函数 $g (x)$ 的曲线

Fig.1 Algebraic function $g (x)$

f_{N} (x) = F_{0} \cdot x^{2} - 3.55 - \sum_{i = 1}^{N} F_{i} \cdot [g (x - E_{i}) - g (x + E_{i}) + 2]

（2）

其中， $F_{0} = 2$ ， $F_{i} = 2 \cdot E_{i}$ ， $E_{i} = i + 1.2$ ， $i = 1,2, . . ., N$ （ $N < 4$ ）， $N$ 为正整数. $g (\cdot)$ 由式（1）定义，其中 $d = 0.001$ .经典双翅膀S-M混沌系统^［

11］可表示为

\{\begin{array}{l} \dot{x} = y \\ \dot{y} = (1 - z) \cdot x - a \cdot y \\ \dot{z} = x^{2} - b \cdot z \end{array}

（3）

其中，a=0.75，b=0.45是参数. x，y和z是状态变量.

为了构造单方向多翅膀S-M混沌系统，用多分段非线性函数 $f_{N} (x)$ 代替双翅膀S-M混沌系统（3）中的 $x^{2}$ 项，可得单方向多翅膀S-M混沌系统（4）.

\{\begin{array}{l} \dot{x} = y \\ \dot{y} = (1 - z) \cdot x - a \cdot y \\ \dot{z} = f_{N} (x) - b \cdot z \end{array}

（4）

单方向多翅膀S-M混沌系统（4）可以产生4翅膀S-M混沌吸引子，6翅膀S-M混沌吸引子和8翅膀S-M混沌吸引子.当初值为 $(1,1, 1)$ ，仿真时间为1000 s，4翅膀S-M混沌吸引子在x-z平面的相图如图2所示.4翅膀S-M混沌吸引子的Multisim电路实现结果如图3所示.对比图2和图3可知，电路实现结果与数值仿真结果一致.其他单方向多翅膀S-M混沌吸引子的电路实现结果也与对应的数值仿真结果一致，这里不再给出.

图2 4翅膀S-M混沌吸引子在x-z平面的相图

Fig.2 Phase diagrams of the 4-wing S-M chaotic attractor in the x-z plane

图3 4翅膀S-M混沌吸引子的电路实现

Fig.3 Circuit implementation of the 4-wing S-M chaotic attractor

2 单方向多翅膀S-M混沌系统的非线性动力学行为分析

2.1　单方向多翅膀S-M混沌吸引子的平衡点

为了研究单方向多翅膀S-M混沌吸引子的产生机理，对平衡点 $Q_{\pm n} (x_{\pm n}, y_{\pm n}, z_{\pm n})$ 进行分析.令此混沌系统（4）中的 $\dot{x} = \dot{y} = \dot{z} = 0$ ，可得

\{\begin{array}{l} \dot{x} = y = 0 \\ \dot{y} = (1 - z) \cdot x - a \cdot y = 0 \\ \dot{z} = f_{N} (x) - b \cdot z = 0 \end{array}

（5）

由式（5）可得

\{\begin{array}{l} y_{\pm n} = 0 \\ z_{\pm n} = 1 \\ f_{N} (x_{\pm n}) = b \cdot z_{\pm n} \end{array}

（6）

根据式（6）可得单方向多翅膀S-M混沌系统（4）的平衡点 $Q_{\pm n} (x_{\pm n}, y_{\pm n}, z_{\pm n})$ ， $n = 0,1, 2, . . ., N$ .图4表示当 $N = 2$ 时的多分段非线性函数 $f_{N} (x)$ .可以看到在x轴方向上有6个交点（ $x_{\pm 0}$ ， $x_{\pm 1}$ 和 $x_{\pm 2}$ ），在图中用‘·’表示.这6个交点可产生6个翅膀.

图4 当 $N = 2$ 时，多分段非线性函数 $f_{_{N}} (x)$

Fig.4 Multisegment nonlinear function $f_{_{N}} (x)$ with $N = 2$

6翅膀S-M混沌吸引子的平衡点Q_±_n（x_±n，y_±n，z_±_n）和特征值 $λ$ ，如表1所示.由表1可知，特征值 $λ$ 有一个负实根和一对实部为正的共轭复根.这表明平衡点（ $Q_{\pm 0}$ ， $Q_{\pm 1}$ 和 $Q_{\pm 2}$ ）为指标2^［

9］鞍焦平衡点.每个指标2鞍焦平衡点可产生一个相对应的翅膀，因此6个指标2鞍焦平衡点可产生6个翅膀.

表1 6翅膀S-M混沌吸引子的平衡点

Q_{\pm n} (x_{\pm n}, y_{\pm n}, z_{\pm n})

和特征值

λ

Table 1 Equilibrium

Q_{\pm n} (x_{\pm n}, y_{\pm n}, z_{\pm n})

and eigenvalues

λ

of the 6-wing S-M chaotic attractor

$Q_{\pm n} (x_{\pm n}, y_{\pm n}, z_{\pm n})$	$λ$
$Q_{_{\pm 0}} (\pm 1.4151,0, 1)$	$- 2.4228, 0.6114 \pm i 1.7107$
$Q_{_{\pm 1}} (\pm 2.4698,0, 1)$	$- 3.3481, 1.0741 \pm i 2.5251$
$Q_{_{\pm 2}} (\pm 4.5044,0, 1)$	$- 4.1143, 1.4572 \pm i 3.1945$

2.2　单方向多翅膀S-M混沌吸引子的不变性和耗散性

单方向多翅膀S-M混沌系统（4）在 $(x, y, z) \to (- x, - y, z)$ 的变换下保持不变，即单方向多翅膀S-M混沌系统（4）关于z轴对称.通过体积 $V = [\dot{x}, \dot{y}, \dot{z}]$ 的散度来计算单方向多翅膀S-M混沌系统（4）的耗散性，即

\nabla V = \frac{\partial \dot{x}}{\partial x} + \frac{\partial \dot{y}}{\partial y} + \frac{\partial \dot{z}}{\partial z} = - a - b = - 1.2 < 0

（7）

因此，单方向多翅膀S-M混沌系统（4）是耗散的.

2.3　单方向多翅膀S-M混沌吸引子的李亚普诺夫指数、分维数和庞加莱截面

本节以6翅膀S-M混沌吸引子为例来研究李亚普诺夫指数、分维数和庞加莱截面.当初值为 $(1,1, 1)$ ，仿真时间为3000 s，6翅膀S-M混沌吸引子的李亚普诺夫指数如图5所示.该图中有一个正的指数 $L E_{1}$ ，这表明单方向多翅膀S-M混沌系统（4）是混沌的.由 $L E_{1} = 0.348291$ ， $L E_{2} = 0.002137$ 和 $L E_{3} = - 1.550427$ ，可计算分维数^［

9］

D_{L} = 2.226020

.可知分维数是介于2和3之间的一个分数值，这表明单方向多翅膀S-M混沌系统（4）是混沌的.

图5 6翅膀S-M混沌吸引子的李亚普诺夫指数

Fig.5 Lyapunov exponents of the 6-wing S-M chaotic attractor

由图6可知，6翅膀S-M混沌吸引子的庞加莱截面由成片密集点组成，这表明单方向多翅膀S-M混沌系统（4）是混沌的.

图6 6翅膀S-M混沌吸引子的庞加莱截面

Fig.6 Poincaré section of the 6-wing S-M chaotic attractor

3 二方向网格多翅膀S-M混沌系统的构造

本节利用代数函数 $g (x)$ 构造多分段非线性函数 $f_{M} (z)$ ，并利用多分段非线性函数 $f_{M} (z)$ 构造二方向网格多翅膀S-M混沌系统.此系统能产生二方向网格多翅膀S-M混沌吸引子.利用代数函数 $g (x)$ 构造多分段非线性函数 $f_{M} (z)$ ，其可以表示为

\begin{array}{l} f_{M} (z) = z + 1 - G_{0} \cdot {g (z) + [g (z - U_{i}) \\ + g (z + U_{i})]} \end{array}

（8）

其中， $G_{0} = 1.5$ ， $U_{i} = 2 i G_{0}$ ， $i = 1,2, . . ., M$ （ $M < 2$ ）， $M$ 为整数. $g (\cdot)$ 由式（1）定义，其中 $d = 0.001$ .用多分段非线性函数 $f_{M} (z)$ 代替单方向多翅膀S-M混沌系统（4）中的状态变量z，可得二方向网格多翅膀S-M混沌系统（9）.

\{\begin{array}{l} \dot{x} = y \\ \dot{y} = [1 - f_{M} (z)] \cdot x - a \cdot y \\ \dot{z} = f_{N} (x) - b \cdot f_{M} (z) \end{array}

（9）

其中， $a = 0.75$ ， $b = 0.45$ 是系统参数. $f_{N} (x)$ 和 $f_{M} (z)$ 分别由式（2）和式（8）定义. $(2 N + 2) \times (2 M + 2)$ 网格多翅膀S-M混沌系统在x-z平面的相图如图7所示.图7（a），7（b）和7（c）分别表示当 $N = 1,2, 3$ ， $M = 0$ 时， $4 \times 2$ ， $6 \times 2$ 和 $8 \times 2$ 网格多翅膀S-M混沌吸引子的相图.图7（d），7（e）和7（f）表示当 $N = 1,2, 3$ ， $M = 1$ 时， $4 \times 4$ ， $6 \times 4$ 和 $8 \times 4$ 网格多翅膀S-M混沌吸引子的相图.

(a) $4 \times 2$ 网格多翅膀S-M混沌吸引子( $N = 1$ ， $M = 0$ )

(a) $4 \times 2$ -grid multi-wing S-M chaotic attractor ( $N = 1$ ， $M = 0$ )

(b) $6 \times 2$ 网格多翅膀S-M混沌吸引子(N=2，M=0)

(b) $6 \times 2$ -grid multi-wing S-M chaotic attractor(N=2，M=0)

(d) $4 \times 4$ 网格多翅膀S-M混沌吸引子( $N = 1$ ， $M = 1$ )

(d) $4 \times 4$ -grid multi-wing S-M chaotic attractor( $N = 1$ ， $M = 1$ )

(e) $6 \times 4$ 网格多翅膀S-M混沌吸引子( $N = 2$ ， $M = 1$ )

(e) $6 \times 4$ -grid multi-wing S-M chaotic attractor( $N = 2$ ， $M = 1$ )

(f) $8 \times 4$ 网格多翅膀S-M混沌吸引子( $N = 3$ ， $M = 1$ )

(f) $8 \times 4$ -grid multi-wing S-M chaotic attractor( $N = 3$ ， $M = 1$ )

图7 $(2 N + 2) \times (2 M + 2)$ 网格多翅膀S-M混沌吸引子在x-z平面的相图

Fig.7 Phase diagrams of the $(2 N + 2) \times (2 M + 2)$ -grid multi-wing S-M chaotic attractors in the x-z plane

$4 \times 2$ 网格多翅膀S-M混沌吸引子的Multisim电路实现结果如图8所示.对比图7（a）和图8可知，电路实现结果与数值仿真结果一致.其他二方向网格多翅膀S-M混沌吸引子的电路实现结果也与对应的数值仿真结果一致，这里不再给出.

图8 $4 \times 2$ 网格多翅膀S-M混沌吸引子的电路实现

Fig.8 Circuit implementation of the $4 \times 2$ -grid multi-wing S-M chaotic attractor

4 二方向网格多翅膀S-M混沌吸引子的非线性动力学行为分析

4.1　二方向网格多翅膀S-M混沌吸引子的平衡点、不变性和耗散性

本节通过分析平衡点 $Q_{\pm n, \pm m} (x_{\pm n}, y_{\pm n}, z_{\pm m})$ 来研究二方向网格多翅膀S-M混沌系统的产生机理.令二方向网格多翅膀S-M混沌系统（9）中 $\dot{x} = \dot{y} = \dot{z} = 0$ ，可得

\{\begin{array}{l} \dot{x} = y = 0 \\ \dot{y} = [1 - f_{M} (z)] \cdot x - a \cdot y \\ \dot{z} = f_{N} (x) - b \cdot f_{M} (z) = 0 \end{array} = 0

（10）

由式（10）可得

\{\begin{array}{l} y_{\pm n} = 0 \\ f_{M} (z_{\pm m}) = 1 \\ f_{N} (y_{\pm n}) = b \cdot f_{M} (z_{\pm m}) \end{array}

（11）

由式（11）得到二方向网格多翅膀S-M混沌系统（9）的平衡点 $Q_{\pm n, \pm m} (x_{\pm n}, y_{\pm n}, z_{\pm m})$ ， $n = 0,1, 2, . . ., N$ ， $m = 0,1, 2, \dots M$ .图9表示当 $M = 1$ 时的多分段非线性函数 $f_{M} (z)$ .可以看到在z轴方向上有4个交点（ $z_{\pm 0}$ 和 $z_{\pm 1}$ ），在图中用‘·’表示.这4个交点可产生4个翅膀.

图9 当 $M = 1$ 时，多分段非线性函数 $f_{M} (z)$

Fig.9 Multisegment nonlinear function $f_{M} (z)$ with $M = 1$

$6 \times 4$ 网格多翅膀S-M混沌吸引子的平衡点 $Q_{\pm n, \pm m} (x_{\pm n}, y_{\pm n}, z_{\pm m})$ 和特征值 $λ$ ，如表2所示.由表2可知，特征值 $λ$ 有一个负实根和一对实部为正的共轭复根.这表明平衡点（ $Q_{\pm 0, \pm 0}$ ， $Q_{\pm 1, \pm 0}$ ， $Q_{\pm 2, \pm 0}$ ， $Q_{\pm 0, \pm 1}$ ， $Q_{\pm 1, \pm 1}$ ， $Q_{\pm 2, \pm 1}$ ）为指标2^［

9］鞍焦平衡点.每个指标2鞍焦平衡点可产生一个相对应的翅膀，因此24个指标2鞍焦平衡点可产生24个翅膀.

表2

6 \times 4

网格多翅膀S-M混沌吸引子的平衡点

Q_{\pm n, \pm m} (x_{\pm n}, y_{\pm n}, z_{\pm m})

和特征值

λ

Table 2 Equilibrium

Q_{\pm n, \pm m} (x_{\pm n}, y_{\pm n}, z_{\pm m})

and eigenvalues

λ

of the

6 \times 4

-grid multi-wing S-M chaotic attractor

$Q_{\pm n, \pm m} (x_{\pm n}, y_{\pm n}, z_{\pm m})$	$λ$
$Q_{\pm 0, \pm 0} (\pm 1.4151,0, \pm 1.5)$	$- 2.4221, 0.6113 \pm i 1.7101$
$Q_{\pm 1, \pm 0} (\pm 2.5285,0, \pm 1.5)$	$- 3.3471, 1.0738 \pm i 2.5243$
$Q_{\pm 2, \pm 0} (\pm 3.5761,0, \pm 1.5)$	$- 4.1131, 1.4568 \pm i 3.1935$
$Q_{\pm 0, \pm 1} (\pm 1.4151,0, \pm 4.4996)$	$- 2.4224, 0.6113 \pm i 1.7104$
$Q_{\pm 1, \pm 1} (\pm 2.5285,0, \pm 4.4996)$	$- 3.3476, 1.0739 \pm i 2.5247$
$Q_{\pm 2, \pm 1} (\pm 3.5761,0, \pm 4.4996)$	$- 4.1137, 1.4569 \pm i 3.1940$

二方向网格多翅膀S-M混沌系统（9）的不变性和耗散性计算方法与单方向多翅膀S-M混沌系统（4）相同，这里不再详细讨论.

4.2　二方向网格多翅膀S-M混沌吸引子的李亚普诺夫指数、分维数和庞加莱截面

本节以 $6 \times 4$ 网格多翅膀S-M混沌吸引子为例来分析李亚普诺夫指数、分维数和庞加莱截面.当初值为 $(1,1, 1)$ ，仿真时间为3000 s， $6 \times 4$ 网格多翅膀S-M混沌系统的李亚普诺夫指数如图10所示.该图中显然存在一个正的李亚普诺夫指数 $L E_{1}$ ，这表明二方向网格多翅膀S-M混沌系统（9）是混沌的.由 $L E_{1} = 0.453331$ ， $L E_{2} = - 0.020405$ 和 $L E_{3} = - 1.567040$ ，可计算分维数^［

9］

D_{L} = 2.276269

.可知分维数是介于2和3之间的一个分数值，这表明二方向网格多翅膀S-M混沌系统（9）是混沌的.

图10 $6 \times 4$ 网格多翅膀S-M混沌吸引子的李亚普诺夫指数

Fig.10 Lyapunov exponents of the 6×4 -grid multi-wing S-M chaotic attractor

由图11可知， $6 \times 4$ 网格多翅膀S-M混沌吸引子的庞加莱截面是由成片的密集点构成，这表明二方向网格多翅膀S-M混沌系统（9）是混沌的.

图11 $6 \times 4$ 网格多翅膀S-M混沌吸引子的庞加莱截面

Fig.11 Poincaré section of the $6 \times 4$ -grid multi-wing S-M chaotic attractor

5 动力学行为演变方向及复杂度分析

以6翅膀S-M混沌系统（ $a = 0.75$ ， $b \in [0,2]$ ）为例，分析李亚普诺夫指数如图12所示及分岔图如图13所示.从图12-图13可知，随着参数 $b$ 的变化，6翅膀S-M混沌系统由混沌态（ $0 \leq b < 1$ ）演变为周期态（ $1 \leq b \leq 2$ ）.

图12 随参数 $b$ 变化的李亚普诺夫指数

Fig.12 Lyapunov exponents with varying parameter $b$

图13 分岔图

Fig.13 Bifurcation diagram

对经典双翅膀S-M混沌系统（3）、单方向多翅膀S-M混沌系统（4）以及二方向网格多翅膀S-M混沌系统（9），分别从SE复杂度^［

24］的最大值和平均值进行比较分析，结果如表3所示.由表3可知，对于二方向网格多翅膀S-M混沌系统（9），SE复杂度的平均值比经典双翅膀S-M混沌系统（3）、单方向多翅膀S-M混沌系统（4）的大.对于单方向多翅膀S-M混沌系统（4），SE复杂度的最大值和平均值比经典双翅膀S-M混沌系统（3）的大.由此可知对经典双翅膀S-M混沌系统（3）多翅膀化后，所得单方向多翅膀S-M混沌系统（4）以及二方向网格多翅膀S-M混沌系统（9）复杂度增加.需要注意的是，复杂度特性对比研究中，单方向多翅膀S-M混沌系统（4）以6翅膀S-M混沌系统（

a = 0.75

，

b \in [0,2]

）为例，二方向网格多翅膀S-M混沌系统（9）以6×4网格多翅膀S-M混沌系统（

a = 0.75

，

b \in [0,2]

）为例.

表3 S-M混沌系统(3)，(4)，(9)复杂度特性对比

Table 3 Comparison of complexity characteristics of S-M chaotic systems (3), (4), (9)

S-M chaotic system	SE_max	SE
Double-wing S-M chaotic system (3)	0.6167	0.1957
One-direction multi-wing S-M chaotic system (4)	0.8634	0.3178
Two-direction grid multi-wing S-M chaotic system (9)	0.8593	0.4031

6 结论

本文基于光滑代数函数构造了两个新光滑多分段非线性函数 $f_{N} (x)$ 和 $f_{M} (z)$ .然后，采用这两个新光滑多分段非线性函数构造了新光滑单方向多翅膀S-M混沌系统和新光滑二方向网格多翅膀S-M混沌系统.翅膀的个数由指标2鞍焦平衡点的个数来确定.需要说明的是，该光滑代数函数也可以被用于构造其他新光滑多分段非线性函数.此外，本文构造的新光滑单方向多翅膀S-M混沌系统和新光滑二方向网格多翅膀S-M混沌系统可应用于保密通信，图像加密等领域.

参考文献

Lorenz E N. Deterministic nonperiodic flow. Journal of the Atmospheric Sciences，1963，20：130~141 [百度学术]

李雄杰，周东华.一种基于强跟踪滤波的混沌保密通信方法.物理学报，2015，64（14）：140501（Li X J，Zhou D H. A method of chaotic secure communication based on strong tracking filter. Acta Physica Sinica，2015，64（14）：140501（in Chinese）） [百度学术]

Changchien S K，Huang C K，Nien H H，et al. Synchronization of the chaotic secure communication system with output state delay. Chaos， Solitons and Fractals，2009，39：1578~1587 [百度学术]

颜森林.激光混沌并行串联同步及其在中继器保密通信系统中的应用.物理学报，2019，68（17）：170502 [百度学术]

Yan S L. Chaotic laser parallel series synchronization and its repeater applications in secure communication. Acta Physica Sinica，2019，68（17）：170502（in Chinese） [百度学术]

Liu S Q，Jiang N，Zhao A K，et al. Secure optical communication based on cluster chaos synchronization in semiconductor lasers network. IEEE Access，2020，8：11872~11879 [百度学术]

庄志本，李军，刘静漪，等.基于新的五维多环多翼超混沌系统的图像加密算法.物理学报，2020，69（4）：040502 [百度学术]

Zhuang Z B，Li J，Liu J Y，et al. Image encryption algorithm based on new five-dimensional multi-ring multi-wing hyperchaotic system. Acta Physica Sinica，2020，69（4）：040502（in Chinese） [百度学术]

石航，王丽丹.一种基于压缩感知和多维混沌系统的多过程图像加密方案.物理学报，2019，68（20）：200501（Shi H，Wang L D. Multi-process image encryption scheme based on compressed sensing and multi-dimensional chaotic system. Acta Physica Sinica，2019，68（20）：200501（in Chinese）） [百度学术]

李桂珍，任晓芳.基于DNA合成图像和混沌映射的图像加密算法.控制工程，2018，25（7）：1278~1284 [百度学术]

Li G Z，Ren X F.Research of image encryption algorithm based on DNA image synthesis and chaotic mapping. Control Engineering of China，2018，25（7）：1278~1284（in Chinese） [百度学术]

禹思敏.混沌系统与混沌电路：原理，设计及其在通信中的应用.西安：西安电子科技大学出版社，2011（Yu S M. Chaotic systems and chaotic circuits：Principle，design and its application in communications. Xi'an：Xi'an Electronic Science＆Technology University Press，2011（in Chinese）） [百度学术]

Zhang C X，Yu S M. A novel methodology for constructing a multi-wing chaotic and hyperchaotic system with a unified step function switching control. Chinese Physics B，2016，25（5）：050503-1~050503-14 [百度学术]

Yu S M，Tang W K S，Lv J H，et al. Generation of n×m-wing Lorenz-like attractors from a modified Shimizu-Morioka model. IEEE Transactions on Circuits and Systems-II：Express Briefs，2008，55（11）：1168~1172 [百度学术]

Yu S M，Tang W K S，Lv J H，et al. Generating 2n-wing attractors from Lorenz-like systems. International Journal of Circuit Theory and Applications，2010，38：243~258 [百度学术]

周欣，王春华，郭小蓉.一个新的网格多翅膀混沌系统及其电路实现.物理学报，2012，61（20）：200506（Zhou X，Wang C H，Guo X R. A new grid multi-wing chaotic system and its circuit implementation. Acta Physica Sinica，2012，61（20）：200506（in Chinese）） [百度学术]

Huang Y，Zhang P，Zhao W F. Novel grid multiwing butterfly chaotic attractors and their circuit design. IEEE Transactions on Circuits and Systems-II：Express Briefs，2015，62（5）：496~500 [百度学术]

孙克辉，傅元理.简化Lorenz系统多翅膀混沌吸引子的设计与电路实现.动力学与控制学报，2016，14（5）：395~400 [百度学术]

Sun K H，Fu Y L. Design and circuit implementation of the simplified Lorenz multi-wing chaotic attractor. Journal of Dynamics and Control，2016，14（5）：395~400（in Chinese） [百度学术]

Zhang S，Zeng Y C，Li Z J，et al. A novel grid multiwing chaotic system with only non-hyperbolic equilibria. Pramana：Journal of Physics，2018，90（5）：63~71 [百度学术]

Yu N，Wang Y W，Liu X K，et al. 3D grid multi-wing chaotic attractors. International Journal of Bifurcation and Chaos，2018，28（4）：1850045 [百度学术]

Zhou L，Wang C H，Zhou L L. Generating hyperchaotic multi-wing attractor in a 4D memristive circuit. Nonlinear Dynamics，2016，85：2653~2663 [百度学术]

Zhou L，Wang C H，Zhou L L. Generating four-wing hyperchaotic attractor and two-wing，three-wing， and four-wing chaotic attractors in 4D memristive system. International Journal of Bifurcation and Chaos，2017，27（2）：1750027 [百度学术]

Zhou L，Wang C H，Zhou L L. A novel no-equilibrium hyperchaotic multi-wing system via introducing memristor. International Journal of Circuit Theory and Applications，2018，46：84~98 [百度学术]

Hu X Y，Liu C X，Liu L. Multi-scroll hidden attractors and multi-wing hidden attractors in a 5-dimensional memristive system. Chinese Physics B，2017，26（11）：110502 [百度学术]

Shimizu T， Morioka N. On the bifurcation of an asymmetric one in a simple model. Physics Letters A，1980，76（3-4）：201~204 [百度学术]

Nattagit J，Tachibana M，Wimol S U. Robustification of a one-dimensional generic sigmoidal chaotic map with application of true random bit generation. Entropy，2018，20（136）：1~15 [百度学术]

孙克辉，贺少波，何毅，等.混沌伪随机序列的谱熵复杂性分析.物理学报，2013，62（1）：010501 [百度学术]

Sun K H， He S B， He Y， et al. Complexity analysis of chaotic pseudo-random sequences based on spectral entropy algorithm. Acta Physica Sinica，2013，62（1）：010501（in Chinese） [百度学术]

新光滑多翅膀混沌系统研究

摘要

关键词

引言

1 单方向多翅膀S-M混沌系统的构造