离散忆阻混沌系统的Simulink 建模及其动力学特性分析*

扶龙香 贺少波 王会海 孙克辉

(中南大学物理与电子学院,长沙 410083)

为拓广离散忆阻器的研究与应用,基于差分算子,构建了具有平方非线性的离散忆阻模型,并实现了Simulink仿真.仿真结果表明,设计的忆阻器满足广义忆阻定义.将得到的离散忆阻引入三维混沌映射中,设计了一种新型四维忆阻混沌映射,并建立了该混沌映射的Simulink 模型.通过平衡点、分岔图、Lyapunov 指数谱、复杂度、多稳态分析了系统复杂动力学特性.本文从系统建模角度出发,构建离散忆阻与离散忆阻混沌映射,进一步验证了离散忆阻的可实现性,为离散忆阻应用研究奠定了基础.

1 引言

1971 年,蔡少棠教授[1]在研究基本电学物理量即电压V、电流I、电荷q、磁通φ之间的相互关系时,基于电路完备性考虑,预测存在除电阻、电感和电容外的第4 种无源非线性基本电路元件,用以描述电荷与磁通之间的关系,称其为忆阻器.2008 年,惠普(Hewlett-Packard,HP)实验室[2]利用TiO2和金属铂首次实现了世界上纳米级忆阻器,从此掀起了忆阻器研究热潮.近年来,忆阻器在许多应用领域都有深入研究,如图像处理[3]、存储器[4-6]、神经网络[7-9]、神经动力学[10,11]等.但由于HP 忆阻器纳米级工艺以及严格实验条件的限制,迄今为止,忆阻器尚未得到商用.为了推进对忆阻器等记忆元件的研究,建立合适的仿真模型至关重要.近年来,很多基于SPICE (simulation pro gram with integrated circuit emphasis)的各类忆阻仿真模型相继被提出[12-14].同时,Simulink 作为Matlab 中一种可视化建模工具,能对大多数非线性系统进行分析,具有流程清晰、易观察、效率高等诸多优点,便于更好地研究与分析非线性系统.因此,采用Matlab/Simulink 仿真工具对忆阻器进行建模值得进一步研究[15-17].

忆阻器是一种非线性元件,其在混沌系统方面的应用也是忆阻器的研究热点之一.在王春华等[18]对忆阻器相关应用总结中提到,可将忆阻器应用到混沌电路分为4 种方式,即用忆阻器替换蔡氏电路中非线性元件[19-21]、新增忆阻器[22-24]、忆阻器替换非Chua 系统中线性/非线性项的方法[25,26]、利用不同数学模型的忆阻器进行电路设计[27,28].本文将利用新增忆阻的方式,将离散忆阻引入到三维混沌映射中,构建一种新型四维忆阻混沌映射,并研究离散忆阻混沌系统建模方法.

本文内容安排如下.第2 节,基于离散忆阻数学模型建立离散忆阻Simulink 仿真模型;第3 节,将离散忆阻模型加入到三维Lorenz 混沌映射中,设计一新型混沌映射,并构建其Simulink 模型;第4 节,对该离散忆阻器及混沌系统Simulink 模型的应用进行分析,最后,对本文进行总结.

2 离散忆阻建模与Simulink 仿真

2.1 离散忆阻数学模型

类似于连续忆阻器模型,离散忆阻器的定义为[29]

其中,n=0,1,2,···,Δq(tn)=q(tn+1)-q(tn)为前向差分.因为

等号两边相加,则忆阻“内部电荷”计算式为

显然,q(t0)为离散忆阻内部初始电荷.在实际计算过程中,该离散忆阻器为一离散累加数学模型.一般地,可以将该离散忆阻改写为迭代式,即:

式中,i(tn)为离散输入电流信号,V(tn)为离散输出电压信号,qn为忆阻内部电荷,M(qn)为忆阻值.其中M(·)函数将采用平方函数[30],因此,该离散忆阻器模型可描述为

式中,α和β为忆阻器参数,令k=1.

2.2 Simulink 仿真实现

根据离散忆阻器的数学模型,构建离散忆阻器的Simulink 模型,设计的模型如图1 所示.该模型将离散忆阻数学模型中的加减乘除分类整合,且由离散时间求和模块结合前向差分算子实现电荷q的迭代,其中可通过改变Constant 与Gain1 的取值来改变α与β的值,Scope1 显示的是忆阻器输入电压vn与流经其电流in之间的关系.首先,对该离散忆阻Simulink 模型进行验证,这里取α=1,β=0.000002.发现在正弦电流的作用下,忆阻器电流-电压特性曲线见图2(a),具有迟滞环特性.其次,分别改变输入正弦信号I(tn)=A0sin(2πwtn)的幅值和频率以讨论输入信号幅值和频率对该忆阻模型的影响.取不同的幅值A0=1,1.5,2 A,其中ω=0.005 rad/s,结果如图2(b)所示.可见,在输入信号频率保持不变的情况下,该忆阻模型的滞回曲线随着输入信号幅值的增大而增大.保持输入信号幅值A0=1.5 A 不变,改变输入信号频率,即ω分别取0.1,0.01,0.003 rad/s,结果如图2(c)所示,可见在输入信号幅值一定的情况下,该忆阻模型的滞回曲线随着输入信号频率增大而减小,最终退化为线性器件.经以上分析,根据忆阻器特性判定原则可知,所构建的离散忆阻Simulink 模型满足忆阻定义[31],可用于进一步的应用研究.

图1 离散忆阻Simulink 模型Fig.1.Discrete memristor Simulink model.

图2 离散忆阻器的电流-电压特性曲线 (a)电流-电压特性曲线;(b)幅值变化;(c)频率变化Fig.2.Current-voltage characteristic curves of discrete memristor:(a)Current-voltage characteristic curve;(b)amplitude change;(c)frequency change.

3 离散忆阻混沌系统及其Simulink仿真

3.1 离散忆阻混沌系统

Sprott[32]报道了一种三维混沌映射,并称之为Lorenz 混沌映射,其动力学方程为

采用结构化设计方法,该映射可以被设计为图3 所示的结构图.可见,该系统存在3 个延时模块,1 个加法器和1 个乘法器.

图3 三维Lorenz 混沌映射(6)系统框图Fig.3.Structure diagram of the three-dimensional Lorenz chaotic map (6).

引入离散忆阻,得到一新型三维忆阻混沌映射,其定义为

其中,a为系统参数,用于控制忆阻的输出值对系统动力学行为的影响.根据上述系统方程(7),绘制引入忆阻器后的混沌系统框图,如图4 所示.可见,离散忆阻增加在xn与yn的乘积项上.

图4 加入忆阻器后的三维混沌映射(7)系统框图Fig.4.Structure diagram of the three-dimensional chaotic map with memristor (7).

结合离散忆阻数学模型,令:

进而可以得到系统的四维表达式为

根据(9)式绘制设计的忆阻混沌系统框图,如图5 所示,其中阴影部分描述的是离散忆阻的结构.

图5 四维忆阻混沌映射(9)系统框图Fig.5.Structure diagram of the four-dimensional memristor chaotic map (9).

3.2 离散忆阻混沌系统Simulink 仿真

根据四维忆阻混沌系统,建立该系统Simulink仿真模型,如图6 所示,其中忆阻模块采用图1 所示的模型.模型中增益Gain2 表示系统参数a的值,Scope2 显示的是系统xn,yn,zn时间序列,并可改变3 个延时器以及忆阻模块中离散时间求和模块中的初值,从而分别改变系统初始状态x0,y0,z0,w0值.

图6 四维忆阻混沌映射(9)Simulink 模型Fig.6.Simulink model of thefour-dimensional memristor chaotic map (9).

设置该模型系统参数a=0.9,初始状态为x0=0.5,y0=0.5,z0=0.5,w0=0.4,忆阻参数α=1,β取值分别为—0.1,—0.02,—0.000002.将示波器中的图像保存到WorkSpace,再利用Matlab 绘制出图像,得到x-y平面的吸引子相图,如图7(a)—(c)所示.可见,随着β的取值的不同,吸引子大小形状产生较大变化.设置系统初始状态为x0=0.5,y0=0.5,z0=0.5,w0=0.4,忆阻参数α=1,β=—0.001,改变系统参数a,分别取值为0.25,0.5,0.9,得到的xn-yn平面吸引子相图,如图7(d)—(f)所示,表明参数a也能够改变吸引子大小与形状.绘制图7 中各个状态下对应的Lyapunov 指数,如图8所示.可见,系统参数以及忆阻参数的变动能够改变系统的动力学行为.

图7 四维忆阻混沌映射(9)的Simulink 仿真吸引子图 (a)β=—0.1,超混沌状态;(b)β=—0.02,超混沌状态;(c)β=—0.000002,混沌状态;(d)a=0.25,混沌状态;(e)a=0.5,混沌状态;(f)a=0.9,混沌状态Fig.7.Simulink simulation results of the four-dimensional memristor chaotic map (9):(a)β=—0.1,hyperchaotic;(b)β=—0.02,hyperchaotic;(c)β=—0.000002,chaos;(d)a=0.25,chaos;(e)a=0.5,chaos;(f)a=0.9,chaos.

图8 四维忆阻混沌映射(9)的Simulink 仿真吸引子图对应Lyapunov 指数谱 (a)β=—0.1,超混沌状态;(b)β=—0.02,超混沌状态;(c)β=—0.000002,混沌状态;(d)a=0.25,混沌状态;(e)a=0.5,混沌状态;(f)a= 0.9,混沌状态Fig.8.Simulink simulation attractor diagram of the four-dimensional memristor chaotic map (9)corresponds to Lyapunov Exponent spectra:(a)β=—0.1,hyperchaotic;(b)β=—0.02,hyperchaotic;(c) β=—0.000002,chaos;(d)a=0.25,chaos;(e)a=0.5,chaos;(f)a=0.9,chaos.

4 动力学特性分析

4.1 平衡点及其稳定性分析

令(9)式左右两边相等,即可求出系统平衡点为

其中,C可以为任意值,也就是说w(e)可以取任意值,即系统具有无限多的平衡点.

对系统进行线性化分析,得其Jacobian 矩阵为

代入平衡点后矩阵为

经计算得其特征值为λi=[1,—1,0.5+0.866i,0.5—0.866i].可见,系统存在不稳定平衡点,可产生混沌,即验证了前述Simulink 仿真结果.

4.2 分岔图与Lyapunov 指数谱分析

系统(9)具有无限多不稳定的平衡点,表明此离散忆阻混沌系统具有丰富的动力学行为,为了进一步研究系统的复杂动力学特性,分析了当a在[0,1]范围变化时,Lyapunov 指数谱和系统的分岔图.

Lyapunov 指数是相空间中邻近轨迹发散或收敛的平均指数率,出现正的Lyapunov 指数时,表明系统是混沌的.初始状态设置为x0=0.5,y0=0.5,z0=0.5,w0=0.4,忆阻参数α=1,β=—0.001,迭代次数为N=10000,通过Matlab 程序仿真绘制xn随参数a变化的分岔图和Lyapunov 指数谱,如图9 所示.可见,Lyapunov 指数谱与分岔图吻合,两图都表明了系统能产生复杂的动力学行为.当a∈[0,0.26)时,Lyapunov 指数有两个大于零的值,但其值较小,且分岔图形状与a∈[0.26,1]时有明显差别.因此,进一步分析a的取值与该区间xn序列值及对应的Lyapunov 指数,即a=0.1,忆阻参数和初值保持不变,绘制xn序列值及对应Lyapunov 指数如图10 所示.通过xn序列发现此时系统处于非周期状态,存在混沌成分,且Lyapunov 指数值很小,处于弱混沌状态.当a=0.9 时,Lyapunov 指数中有两个正数,且值相对较大,结合分岔图可以发现该系统处于超混沌状态.通过分岔图还可以发现在超混沌区间中出现了周期窗口,表明系统的动力学行为对参数a的取值比较敏感.

图9 四维忆阻混沌映射(9)分岔图与Lyapunov 指数谱 (a)分岔图;(b)Lyapunov 指数谱Fig.9.Bifurcation and Lyapunov exponent spectra of the four-dimensional memristor chaotic map (9):(a)Bifurcation diagram;(b)Lyapunov exponent (LE)spectra.

图10 弱混沌系统动力学分析 (a)xn 序列图;(b)Lyapunov 指数谱Fig.10.Dynamic analysis of weakly chaotic system:(a)xn sequence diagram;(b)Lyapunov exponent spectra.

初始状态设置为x0=0.5,y0=0.5,z0=0.5,w0=0.4,N=10000,忆阻参数α=1,β=—0.001,系统参数a分别取0.25,0.5,绘制四维忆阻混沌映射(9)随初值z0变化的Lyapunov 指数谱,结果如图11(a)和图11(b)所示,原三维Lorenz 混沌映射(6)的Lyapunov 指数谱,如图11(c)所示,对比可见,加入忆阻器后系统动力学行为变得更加复杂了,由混沌状态进入了超混沌状态.

图11 随初值z0 变化的Lyapunov 指数谱 (a)四维忆阻混沌映射(9),a=0.25;(b)四维忆阻混沌映射(9),a=0.5;(c)三维Lorenz 混沌映射(6)Fig.11.Lyapunov exponent spectra with initial value z0:(a)Four-dimensional memristor chaotic map (9),a=0.5;(b)four-dimensional memristor chaotic map (9),a=0.5;(c)three-dimensional Lorenz chaotic map (6).

4.3 复杂度分析

复杂度是衡量一个混沌序列与随机序列接近程度的指标,复杂度越大意味着混沌序列越随机[33].绘制当初始状态设置为x0=0.5,y0=0.5,z0=0.5,w0=0.4,忆阻参数α=1,β=—0.001,迭代次数N=10000,a在[0,1]之间的系统样本熵复杂度[34],如图12 所示.可见系统样本熵复杂度呈上升趋势,且在a=0.26 后,系统进入高复杂度状态,系统中的混沌序列越来越接近随机序列,再次表明a值约为0.26,系统处于混沌状态.同时,可以发现该混沌系统高复杂度区域较宽,且与对应的分岔图和Lyapunov 指数谱匹配.

图12 四维忆阻混沌映射(9)样本熵复杂度Fig.12.Sample entropy (SampEn)complexity of the fourdimensional memristor chaotic map (9).

4.4 多稳态分析

设系统初始状态为x0=0.5,z0=0.5,w0=0.4,N=10000,系统参数a=0.9,忆阻参数α=1,β=—0.001,忆阻初值y0分别取值为1,0.5,0.25,绘制xn-yn平面的吸引子共存图,如图13(a)所示.取系统参数a=0.9,初始状态设置为x0=0.5,y0=0.5,w0=0.4,N=10000,忆阻参数α=1,β=—0.001.忆阻初值z0分别取值为1,0.5,0.25,绘制xn-yn平面的吸引子共存图,如图13(b)所示.设置系统初始状态为x0=0.5,y0=0.5,z0=0.5,N=10000,系统参数a=0.9,忆阻参数α=1,β=—0.001,忆阻初值w0分别取值为1,0.5,0.25,绘制x-y平面的吸引子共存图,如图13(c)所示.可见设计的离散忆阻系统存在吸引子共存现象.设系统参数a=0.9,图14 分别表示随初值y0,z0,w0变化的分岔图.可见,随着系统初值的变化,各分岔图都出现了较多的周期窗口且吸引子大小也在不断变化,进一步表明系统存在无限多共存吸引子.

图13 xn-yn 平面上共存吸引子 (a)初值y0=1,0.5,0.25;(b)初值z0=1,0.5,0.25;(c)初值w0=1,0.5,0.25Fig.13.Coexisting attractors in the xn-yn plane:(a)y0=1,0.5,0.25;(b)z0=1,0.5,0.25;(c)w0=1,0.5,0.25.

图14 随初值变化的分岔图 (a)初值y0 变化,x0=0.5,z0=0.5,w0=0.4;(b)初值z0 变化,x0=0.5,y0=0.5,w0=0.4;(c)初值w0 变化,x0=0.5,y0=0.5,z0=0.5Fig.14.Bifurcation diagrams with the initial values:(a)Initial value y0 changes,x0=0.5,z0=0.5,w0=0.4;(b)initial value z0 changes,x0=0.5,y0=0.5,w0=0.4;(c)initial value w0 changes,x0=0.5,y0=0.5,z0=0.5.

对于不同系统初值条件下,共存的无限多吸引子的吸引盆可以描绘吸引子与初始条件的关系[35].系统参数分别为a=0.25,0.5,0.9,初始值为x0=0.5,y0=0.5,N=10000,忆阻参数α=1,β=—0.001,z0-w0初值平面上的吸引盆如图15,图中不同颜色代表不同的吸引子类型.可见,总体上,a较大时吸引子类型变化更加复杂.根据图15 中不同a值吸引盆,分析对应的随初值z0变化的样本熵复杂度如图16 所示,发现a=0.9 时,复杂度值总体上明显高于a=0.25 和a=0.5 时的值,且对初值的变化更加敏感,进一步验证了图15 中的结果.

图15 初值平面上吸引盆 (a)a =0.25;(b)a =0.5;(c)a =0.9Fig.15.Attraction basins in the planes constructed by two different initial conditions:(a)a =0.25;(b)a =0.5;(c)a =0.9.

图16 在不同参数a 下随初值z0 变化的样本熵复杂度 (a)a =0.25;(b)a =0.5;(c)a =0.9Fig.16.Sample entropy complexity with initial value z0 under different system parameter a:(a)a =0.25;(b)a =0.5;(c)a =0.9.

原系统随初值z0变化的复杂度如图17 所示,对比图16,再次验证加入忆阻器后系统复杂度变大,且高复杂度的区域更宽,表明设计的忆阻混沌系统具有较原系统更高的复杂度.

图17 三维Lorenz 混沌映射(6)样本熵复杂度Fig.17.Sample entropy complexity of the three-dimensional Lorenz chaotic map (6).

5 结论

本文在离散忆阻数学模型基础上构建了离散忆阻Simulink 模型.在Simulink 模型中借助离散时间求和模块实现电荷的迭代,以体现离散忆阻器独特的记忆特性.仿真结果发现,设计的忆阻器满足忆阻定义,表明本文从仿真的角度实现了忆阻器.同时,将此离散忆阻器引入到三维Lorenz 混沌映射中,实现了一种新型的四维忆阻混沌映射,并建立了该新型四维忆阻混沌映的Simulink 模型.通过分析该离散混沌映射的平衡点及稳定性,发现其具有无穷多平衡点.通过分析该混沌映射分岔图、Lyapunov 指数谱、系统复杂度、吸引子共存、吸引盆以及对应初值平面复杂度,研究了其复杂的动力学特性.数值仿真分析结果表明新型四维忆阻混沌映具有较原系统更丰富的动力学行为以及更高的复杂度.本文的研究进一步验证了离散忆阻模型的可实现性及其潜在的应用价值,为离散忆阻在非线性系统中的应用研究奠定了基础.