整数阶和分数阶多涡卷混沌吸引子的电路实现

secret

[摘要]构造了一个由三角波函数驱动的新三维连续自治多涡卷整数阶和分数阶混沌系统，进行了相应的理论分析和数值仿真，设计并实现了相应的实验电路。新系统的涡卷平衡点均匀分布在x -y平面的网格点上，可以充分利用实验电路有源器件的有效工作区间，在保证有源器件运算精度下可产生更多涡卷数量的混沌吸引子。实验结果表明，新提出的整数阶和分数阶系统可产生多达81个涡卷的混沌吸引子。

    1.作品简介
多涡卷混沌吸引子生成的理论研究和电路实现已成为混沌研究领域的一个新方向。相比于传统的单涡卷和双涡卷混沌系统，多涡卷混沌系统呈现出更为复杂的结构和动力学行为，在保密通信和信息隐藏等领域具有较好的应用前景中。

为了能够实现在实际电路中产生更多涡卷数量的混沌信号，本文构造了一个新的混沌系统，使用三角波函数来驱动产生多涡卷混沌吸引子。与文献[2]所提出的三角波函数驱动的系统相比，本文所构造的系统产生的混沌吸引子均匀地分布在横坐标和纵坐标原点为中心的网格点上，这样可以充分利用运算放大器的有效工作区间。以此系统来设计电路，从而能产生涡卷数量更多的混沌信号输出。实验结果表明新构造的整数阶和分数阶系统都能产生81个涡卷的混沌吸引子。

指标要求:产生整数阶和分数阶49涡卷混沌吸引子和81涡卷混沌吸引子。


2.方案设计
为了获得多涡卷混沌吸引子，需要对系统方程中的非线性函数进行特殊设计，通常采用的函数构造方法有多转折点分段线性函数法、阶梯函数法、滞后函数法和饱和函数法等[2-5]。同样地，利用在整数阶多涡卷混沌系统中的构造方法，在分数阶混沌系统中也能获得相应的多涡卷混沌吸引子[6.7]。由于分数阶混沌系统的动力学特性要比整数阶混沌系统的更为复杂，因此，关于分数阶多涡卷混沌吸引子生成系统的理论研究和电路实现近几年成为了研究者们一个新的研究热点。

2.1整数阶三维多涡卷混沌系统
一个新的由三角波函数驱动实现的三维自治多涡卷混沌系统的状态方程可描述为:




其中μ是控制常数，fv(2)为三角波函数，其数学表达式为



式中A>0，B∈(0, A]称为三角波函数的变参数，N为正整数。



固定N= 1,A=1,可得系统(1)随参数B变化的分岔图如图1所示。由分岔图可知，系统(1)由Hopf分岔进入混沌状态，在B∈(0, 0.287]的区域内存在两维9涡卷混沌吸引子，如图2所示。同时从分岔图中可以观察到B值越小，涡卷分布越均匀。因此，除理论分析选择B值为0.1外，在以后的数值仿真和电路实验中，我们均选择B值为0.02。



当N= 1,A=1, B=0.1时，三角波函数fi (u)曲线如图3所示。当fi(u)= 0时有5个零点，对应u值分别是S:[- 2, 0, 2]和S:[-1, 1],这里s和S分别称为偶零点(或涡卷平衡点)和奇零点(或键波平衡点)。



    令dx_dy_dz=0.
    显然，对N=1,A=1,B=0.1,系统(1)有5x5个平衡点，其dtdt中3x3个平衡点位于z=0的xxy=[-2, 0,2]x[-2, 0, 2]平面网格_上，称之为偶平衡点Es; 其它2x2个平衡点位于z=0的xxy=[-1, 1]x[-1,1]平面网格上，称之为奇平衡点os。
    系统(1)在平衡点S*=(x*, y*, z*)处的Jacobian矩阵为:




这里



若u=2k(k=0, +1, +2...,+N),则有fN(u)=-1;或若u=(2k- lk1k)(k=+1, +2, ... +N),则有fN(u)=(A- B)/B=9。在平衡点S*附近线性化系统得到以下特征方程:



对偶平衡点S,有fj)(u)=-1,上述表达式变成为



注意到该三次多项式的系数在μ>0时均为正数，因此对所有h>0都有f(2)>0。数值计算表明当μ> 1时，方程式(6)有一个负实根和两个实部为正数的共轭复根。μ=1是Hopf分岔的临界点。
对奇平衡点0S，有f:(u)=9，方程式(5)变成为




当μ>1时，方程式(7)有一个负实根和两个实部为正数的共轭复根。
理论分析和数值仿真表明由三角波函数驱动的系统(1)在μ>1的条件下会产生混沌行为。这里系统(1)的所有平衡点S*都是指数2的鞍点。
选择μ=2, A=1, B=0.02, N分别为3和4时，新提出系统是混沌的，能分别产生49涡卷和81涡卷的混沌吸引子，如图4所示。当N=3时，有7x7个涡卷平衡点和5x5个键波平衡点;而当N=4时，有9x9个涡卷平衡点和7x7个键波平衡点。从图6中可以看到，新构造系统的涡卷平衡点均匀分布在x -y平面的网格点上，每个涡卷平衡点都能产生一个涡卷混沌吸引子，这样吸引子运行轨道在较小的尺度范围内可以较充分地填满整个相空间。




2.2.分数阶三维多涡卷混沌系统
与常微分方程的数值求解不同，分数阶微分方程的数值仿真仍成问题。在分数阶混沌领域的相关文献中，提出了两种数值解析分数阶微分方程近似方法。
第一种方法是改进的Adams- Bashforth-Moulton算法[8]，是基于预估校正原理提出的。第二种方法，称之为频域近似9，在频域中近似分数阶系统的特性。这里采用频域近似法由电路仿真给出分数阶多涡卷混沌系统的动力学解析结果。


3.硬件电路设计
可以采用纯模拟电路方法，设计一个实现整数阶和分数阶多涡卷混沌系统的实验电路。选择TI公司的TL082CP运算放大器10完成整数阶和分数阶多涡卷混沌系统方    程中的加减运算、积分运算和三角波函数产生等功能。TL082CP运算放大器的特点有:低功耗，广泛通用模式和电压变化范围;输入偏置低和低补偿失调电流;有输出短路保护:高输入阻抗;内置频率补偿;速度快;输入电源电压通用等。其电源电压E=+15V,实验测得输出电压的饱和值为Vsur=+13.5 V。为了便于电路实验，所有电阻均采用普通色环电阻与精密可调电位器组合实现。



图6的实验电路为N = 4的三角波函数发生器。如果将图6中最上部的两个运算放大器与其它电路断开，那么所得到的电路为N = 3的三角波函数发生器。图7为基本的三维系统电路。

图7中的1/(R*C)为积分器的积分常数，同时也是时间尺度变换因子,固定R = 1002不变，改变C的大小，可改变时间尺度变换因子，从而可改变混沌信号的频谱范围。在    电路实验中,取R =1002, C= 1μF。图7中的电阻R7和Rs决定了μ的大小，其数学关系可以表示为μ=Rs/R7。当μ=2时，取R7=5kS, Rs= 10 k2。
图7实现的整数阶和分数阶系统电路是通过一-组 拨动开关进行转换的。当开关组向左拨动，与分数阶单元电路相连，就可把整数阶系统电路转变成了一个分数阶系统电路。



4.作品性能测试与分析
用实际电路进行混沌系统的实验验证具有重要的物理意义和应用价值1251。但在实际应用中，电路实现更多涡卷数量的混沌吸引子，受到实际元器件的许多限制。例如，有源器件如运算放大器的高精度运算的动态范围一般是较为有限的，在输入信号的动态范围变化较大的情况下，器件输出信号与输入信号之间的运算精度会降低，误差将增大，再加之各个运算放大器以及电路中其它元件(如电阻、电容等)参数的离散性，使得实际运算放大器难以实现较为精确的数学运算关系。这可能是目前很少有文献报道通过电路实验产生更多涡卷混沌吸引子的一一个主要原因。

按照上面设计的图6和图7电路原理图，我们制作了一个实际的实现整数阶和分    数阶多涡卷混沌系统的实验电路。

在整数阶(1, 1, 1)即a=β=y=1时，图7电路中三个通路的开关接通电容器C;在分数阶(0.9, 0.9, 0.9)即a=β=γ= 0.9时，图7电路中三个通路的开关接通图6(a)单元电路。在TL082运算放大器输出端x和-y上分别连接示波器两个通路。

从数字存储示波器上观察到的硬件电路实验结果分别如图8和图9所示。图8为整数阶和分数阶49涡卷混沌吸引子的电路实验输出，图9为整数阶和分数阶81涡卷混沌吸引子的电路实验输出。



5.结论
一般来说，一个混沌系统产生的混沌吸引子涡卷个数越多，其混沌特性就越好，用于信号调制和加密的保密性能也就越好，但由于运算放大器的动态范围是有限的，因此电路输出的涡卷个数也就受到了限制。本文通过构造一-个新的由三角波函数驱动的三维连续自治系统，实现了49涡卷和81涡卷的整数阶和分数阶混沌系统，进行了相应的理论分析和数值仿真，并用一个实际电路实现了整数阶及其分数阶的多涡卷混沌吸引子产生系统。硬件电路实验结果表明，本文提出的实现多涡卷混沌系统这一方案是可行和有效的。图10为电路实物照片。

满座衣冠胜雪

谢谢分享