一种具有冲激函数形式Lyapunov指数的多稳定性混沌系统的制作方法

本发明属于电子通信研究领域，具体涉及一种具有冲激函数形式lyapunov指数的多稳定性混沌系统。

背景技术：

混沌是指确定性系统产生的随机行为。从哲学上看，混沌是确定性与随机性的统一。混沌信号具有内在随机性、初值敏感性、遍历性和有界性等特点，能够产生类似白噪声的宽带信号，在信息加密、保密通信和混沌雷达等领域有着广泛的应用。混沌信号源是各类混沌应用系统的重要组成部分，研究开发新型混沌信号源对混沌理论的实用化至关重要。

1971年，华裔科学家蔡少棠教授根据电路变量的完备性，预测存在第四种基本电路元器件，并将其命名为忆阻器。直到2008年，忆阻器的实物才第一次被hp实验室所发现。由于忆阻器是非线性器件，且具有记忆功能，因此被广泛应用于许多领域，包括低功耗的闪存存储器、神经突触架构设计、神经网络构建、新型的混沌系统的构建等研究。将具有各种非线性特性的忆阻器引入经典混沌电路中,可以很容易地构造出新的混沌电路，且其动力学行为比原混沌系统更加复杂。

最近几年，混沌系统的多稳定性成为了人们的研究热点。它是指在相同的系统参数设置下，多种吸引子共存的现象。多稳定性是许多非线性系统中一种常见的现象。与一般的混沌系统相比，其具有更加复杂的动力学行为，在保密通信和图像视频加密等领域有着广泛的应用前景。

目前存在的多稳定性混沌系统，其lyapunov指数谱都是既有混沌状态也有非混沌状态，其状态过渡都是平缓的。但是本发明所提出的多稳定性混沌系统，系统参数设置为(x(0),1,0,0)时，其lyapunov指数谱除x(0)＝0外保持恒定，在x(0)＝0处发生突变，且其最大lyapunov指数类似于冲击函数的形式。值得注意的是，该系统的非混沌区域大约是十分逼近x(0)＝0，其突变程度十分接近冲激函数。这种剧烈的突变现象在混沌系统的研究领域有着十分重要的意义。新系统的动力学行为更加复杂，在信息加密、保密通信等领域有着巨大的应用前景。

技术实现要素：

本发明所要解决的技术问题是设计一种具有冲激函数形式lyapunov指数的多稳定性混沌系统，对其进行硬件电路实现。

为解决上述技术问题，本发明提供了一种具有冲激函数形式lyapunov指数的多稳定性混沌系统，设计了相应的硬件电路，具体如下：

一种具有冲激函数形式lyapunov指数的多稳定性混沌系统，由在一个三维混沌系统上引入一个忆阻器所构成的的四维忆阻混沌系统构成，所述四维忆阻混沌系统包括积分通道一、积分通道二、积分通道三和积分通道四，各积分通道中vx端口均连接在一起，各积分通道中vy端口均连接在一起，各积分通道中vz端口均连接在一起，各积分通道中-vx端口均连接在一起。

所述三维混沌系统所对应方程为：

其中，a、b、c和d都是常数，x、y和z是状态变量。

所述四维忆阻混沌系统所对应方程为：

其中，e和h是两个常数，w(w)是磁控忆阻模型的忆导函数，构建了电荷和磁通之间的关系，其表达式为：

w(w)＝f+3gw²

其中，f和g是两个正实数，w是状态变量。

所述积分通道一有3个输入端，分别为vx、vy、vz；

输入端vx串联一个电阻r7接于运算放大器u1的反相输入端；u1的反相输入端和输出端之间并联电阻r8，此时u1的输出端输出-vx；u1的输出端和运算放大器u2的反相输入端之间串联一个电阻r1；输入端vy和vz经乘法器m1相乘后串联一个电阻r2接于运算放大器u2的反相输入端；u2的反相输入端和输出端之间并联电容c1，此时u2的输出端输出vx；运算放大器u1和u2的同相输入端均接“地”。

所述积分通道二有4个输入端，分别为2个-vx、1个vy和1个vz；

输入端vy串联一个电阻r3接于运算放大器u3的反相输入端；输入端-vx和vz经乘法器m2相乘后串联一个电阻r4接于运算放大器u3的反相输入端；输入端-vx串联一个忆阻器接于运算放大器u3的反相输入端；u3的反相输入端和输出端之间并联电容c2，此时u3的输出端输出vy；运算放大器u3的同相输入端接“地”。

所述积分通道三有3个输入端，分别为-vx、vy和vz；

输入端vz串联一个电阻r5接于运算放大器u4的反相输入端；输入端-vx和vy经乘法器m3相乘后串联一个电阻r6接于运算放大器u4的反相输入端；u4的反相输入端和输出端之间并联电容c3，此时u4的输出端输出vz；运算放大器u4的同相输入端接“地”。

所述积分通道四只有1个输入端-vx；

输入端-vx串联一个电阻ra接于运算放大器u5的反相输入端；u5的反相输入端和输出端之间并联电容c4，此时u5的输出端输出vw；u5的输出端vw经乘法器ma相乘作平方运算后输出vw²；vw²和-vx经乘法器mb相乘后串联电阻rb接于运算放大器u3的反相输入端；输入端-vx串联电阻rc也接于运算放大器u3的反相输入端。

本发明的有益效果在于：

本发明相比于一般的混沌系统，新系统结构简单，易于电路实现，且具有更加复杂的动力学特性，对研究共存多吸引子及其硬件电路实现有着重要的理论物理意义和工程应用价值。

附图说明

图1为积分通道一的电路图；

图2为积分通道二的电路图；

图3为积分通道三的电路图；

图4为系统参数设置为(x(0),1,0,0)，x(0)变化时，状态变量w的分岔图；

图5为系统参数设置为(x(0),1,0,0)，x(0)变化时的lyapunov指数谱；

图6为x-y的多稳定性混沌系统的数值仿真相轨图；

图7为x-z的多稳定性混沌系统的数值仿真相轨图；

图8为y-z的多稳定性混沌系统的数值仿真相轨图；

图9为x-w的多稳定性混沌系统的数值仿真相轨图；

图10为y-w的多稳定性混沌系统的数值仿真相轨图；

图11为z-w的多稳定性混沌系统的数值仿真相轨图；

图12为vx-vy的多稳定性混沌系统的pspice电路仿真结果图；

图13为vx-vz的多稳定性混沌系统的pspice电路仿真结果图；

图14为vy-vz的多稳定性混沌系统的pspice电路仿真结果图；

图15为vx-vw的多稳定性混沌系统的pspice电路仿真结果图；

图16为vy-vw的多稳定性混沌系统的pspice电路仿真结果图；

图17为vz-vw的多稳定性混沌系统的pspice电路仿真结果图。

具体实施方式

下面结合附图对本发明做进一步描述。

新提出的多稳定性混沌系统具有冲激函数形式的lyapunov指数，此外还具有离散的分岔图，可以产生两种结构不同且沿第四轴平行分布的无限多共存吸引子。

一种具有冲激函数形式lyapunov指数的多稳定性混沌系统，在一个三维混沌系统的基础上引入一个忆阻器，将其改造为一个新的四维忆阻混沌系统。通过相轨图、lyapunov指数谱和分岔图分析，研究了这个系统的混沌行为，分析表明新系统可以产生两种结构不同且沿第四轴平行分布的无限多共存吸引子，且具有冲激函数形式的lyapunov指数和离散的分岔图。

所述主电路如图1所示，包括：积分通道一、积分通道二、积分通道三和积分通道四。积分通道一有3个输入端，分别为1个“vx”、1个“vy”和1个“vz”，分别通过反相器和乘法器，再通过积分器最后输出“vx”；积分通道二有4个输入端，分别为2个“-vx”、1个“vy”和1个“vz”，先通过乘法器和忆阻器，再通过积分器最后输出“vy”；积分通道三有3个输入端，分别为1个“-vx”、1个“vy”和1个“vz”，通过乘法器和积分器最后输出“vz”；积分通道四为图2方框中的忆阻器部分，其只有1个输入端“-vx”，通过积分器输出“vw”，通过两个乘法器后再通过积分通道2中的积分器，最后输出“vy”；运算放大器u1、u2、u3、u4和u5的同相输入端均接“地”。

积分通道一中，输入端“vx”串联一个电阻r7接于运算放大器u1的反相输入端；u1的反相输入端和输出端之间并联电阻r8，此时输出端输出“-vx”；u1的输出端和运算放大器u2的反相输入端之间串联一个电阻r1；输入端“vy”和“vz”经乘法器m1相乘后串联一个电阻r2接于运算放大器u2的反相输入端；u2的反相输入端和输出端之间并联电容c1，此时u2的输出端输出“vx”；运算放大器u1和u2的同相输入端均接“地”。

积分通道二中，输入端“vy”串联一个电阻r3接于运算放大器u3的反相输入端；输入端“-vx”和“vz”经乘法器m2相乘后串联一个电阻r4接于运算放大器u2的反相输入端；输入端“-vx”串联一个忆阻器(即方框内电路)接于运算放大器u3的反相输入端；u3的反相输入端和输出端之间并联电容c2，此时u3的输出端输出“vy”；运算放大器u3的同相输入端接“地”。

积分通道三中，输入端“vz”串联一个电阻r5接于运算放大器u4的反相输入端；输入端“-vx”和“vy”经乘法器m3相乘后串联一个电阻r6接于运算放大器u4的反相输入端；u4的反相输入端和输出端之间并联电容c3，此时u3的输出端输出“vz”；运算放大器u4的同相输入端接“地”。

积分通道四(即方框中的忆阻器电路)中，输入端“-vx”串联一个电阻ra接于运算放大器u5的反相输入端；u5的反相输入端和输出端之间并联电容c4，此时u5的输出端输出“vw”；u5的输出端“vw”经乘法器ma相乘作平方运算后输出“vw²”；“vw²”和“-vx”经乘法器mb相乘后串联电阻rb接于运算放大器u3的反相输入端；输入端“-vx”串联电阻rc也接于运算放大器u3的反相输入端；u3的反相输入端和输出端之间并联电容c2，此时u3的输出端输出“vy”；运算放大器u3的同相输入端接“地”。

本发明提供一种具有冲激函数形式lyapunov指数的多稳定性混沌系统。主要思想为，在一个三维混沌系统的基础上引入一个忆阻器，将其改造为一个新的四维忆阻混沌系统，通过相轨图、lyapunov指数谱和分岔图分析，研究了这个系统的混沌行为。分析表明新系统可以产生两种结构不同且沿第四轴平行分布的无限多共存吸引子，且具有冲激函数形式的lyapunov指数和离散的分岔图。

所述三维混沌系统所对应方程为：

其中a、b、c和d都是常数，x、y和z是状态变量。

在上述三维混沌系统中引入一个忆阻器，将其改造为一个四维忆阻混沌系统，改造后的四维混沌系统所对应方程为：

其中，e和h是两个正常数。w(w)是磁控忆阻模型的忆导函数，构建了电荷和磁通之间的关系，其表达式如下式所示。

w(w)＝f+3gw²

其中，f和g是两个正实数，w是状态变量。

新系统的基本动力学分析

(1)对称性

对称性广泛存在于具有偶数个吸引子的混沌系统中。忆阻混沌系统(2)与原三维混沌系统(1)的对称性一致，即系统(2)在变换下保持不变性。这意味着状态空间中的吸引子必须相对于z轴对称。

(2)平衡点和稳定性

可以看出提出的多稳定性混沌系统的平衡点仅取决于x，y和z，与w不相关。该系统具有一个线性平衡点：

o＝{(x,y,z,w)x＝y＝z＝0,w＝k}

其中k是任意常数。

将系统(2)在原点o线性化，可以得到其雅可比矩阵：

系统的特征方程如下式所示。

λ(λ-a)(λ+c)(λ+d)＝0

可以求出系统的特征根：

λ1＝0,λ2＝a,λ3＝-c,λ4＝-d

当参数a，c和d都是正实数时，特征根λ3和λ4都为负，而特征根λ2总为正。因此，该系统存在一个正实数根、一个零根和两个负实数根，即系统(2)具有不稳定的鞍点。

(3)耗散性

系统(2)的耗散性可以由下式所示。

当参数a、c和d满足a-c-d＜0时，系统(2)是耗散的。

系统的数值仿真

根据图1所示一种具有冲激函数形式lyapunov指数的多稳定性混沌电路，利用仿真软件平台，可以对所描述的系统进行数值仿真分析。当系统参数a＝4，b＝6，c＝20，d＝5，e＝0.01，f＝1，3g＝0.1，h＝0.1时，设置初值为(x(0),1,0,0)，给出x(0)在[-10⁴,10⁴]区间内的状态变量w的分叉图和lyapunov指数谱如图4、图5所示。

初值分别设置为(1,1,0,0)，(±10,1,0,0)，(±50,1,0,0)，(±100,1,0,0)系统(2)可以产生七个不同的共存吸引子，如图6-图11所示。其中，绿色对应初值(1,1,0,0)，蓝色和红色分别对应初值(10,1,0,0)，(-10,1,0,0)，粉色和青色分别对应初值(50,1,0,0)，(-50,1,0,0)，黄色和黑色分别对应初值(100,1,0,0)，(-100,1,0,0)。可以发现，系统的相轨图与上面的分叉图和lyapunov指数谱具有良好的一致性。

系统的电路实现

(1)对多稳定性混沌系统进行变量比例压缩变换。电源的供电电压是±15v，运算放大器的饱和电压是±13.5v，乘法器的电压范围在±10v之间。混沌吸引子各变量的动态范围可能超出了元件的饱和电压。因此，将系统(2)状态变量压缩10倍，可以得到：

(2)对多稳定性混沌系统进行时间尺度变换，系统的无量纲方程可表示为：

(3)根据变换后的系统状态方程搭建电路，可以得到：

其中，vx，vy，vz和vw分别是各电容上的电压。可求出相应的电阻和电容表达式。c1＝c2＝c3＝c4＝c，r1＝r/a，r2＝r/10b，r3＝r/c，r4＝r/10，r5＝r/d，r6＝r/10，ra＝r/h，rb＝r/(e*3g*100)，rc＝r/ef。

令r＝100kω和c＝10000nf。四维忆阻混沌系统的系统参数为a＝4，b＝6，c＝20，d＝5，e＝0.01，f＝1，3g＝0.1，h＝0.1，系统初值设置为(1,1,0,0)。因此，相应的电阻取值如下所示。r1＝25kω，r2＝1.67kω，r3＝5kω，r4＝10kω，r5＝20kω，r6＝10kω，ra＝1000kω，rb＝1000kω，rc＝10000kω。电容c1和c2的初始电压设置为0.1v，其它电容初始电压保持为0v。值得注意的是，混沌系统状态变量压缩了10倍，系统的初始值也应压缩10倍。电路仿真结果图如图12-图17所示。

本发明公开了一种具有冲激函数形式lyapunov指数的多稳定性混沌电路。主要思想为，在一个三维混沌电路的基础上引入一个忆阻器，将其改造为一个新的四维忆阻混沌电路。电路由基于运算放大器与电阻或电容所构成的反向、加法和积分等运算模块，以及模拟乘法器与运算放大器所构成的非线性运算模块组成。该电路在一定的系统参数下可以产生两种结构不同且沿第四轴平行分布的无限多共存吸引子，且具有冲激函数形式的lyapunov指数和离散的分岔图。与一般混沌系统相比，其具有更加复杂的动力学特性。同时对研究共存多吸引子及其硬件电路实现有着重要的理论物理意义和工程应用价值。

本发明实现的一种具有冲激函数形式lyapunov指数的多稳定性混沌系统，其结构简单，可以产生两种结构不同且沿第四轴平行分布的无限多共存吸引子，并且具有冲激函数形式的lyapunov指数，此外还具有离散的分岔图。