一个新三维混沌系统动力学分析与仿真研究.docVIP

一个新三维混沌系统动力学分析与仿真研究 　　摘 要：文章提出一个新的三维自治混沌系统，分析了系统的平衡点稳定性，计算出系统的Lyapunov指数，给出了系统变量的时域图、混沌吸引子相图和Lyapunov指数谱，理论分析证实新系统是混沌系统。并对此新系统进行了电路仿真实验，所采用的是电子工作平台electronic workbench（EWB）仿真软件，仿真结果再次表明新系统是混沌系统。 　　关键词：混沌吸引子；三维自治混沌系统；Lyapunov指数；EWB 　　引言 　　第一个混沌吸引子是1963年被Lorenz发现的，发现于一个三维自治混沌系统，此后，非线性科学研究的热点中便包含了对混沌理论的研究、新混沌系统的构造和混沌控制及其应用。许多新的混沌与超混沌系统被相继提出[1-7]，例如，chen系统的发现，此系统是一个在混沌系统反控制中与Lorenz系统并不拓扑等价的系统，是陈关荣等人发现的[1]；吕金虎等进一步发现了LV系统和链接Lorenz系统、Chen系统以及LV系统的统一混沌系统[2，3]；国内发布了有关新的离散与连续混沌系统的报道[4，5]；同时，还报道了两个不同的四维超混沌系统。新混沌与超混沌系统的提出与实现[6，7]，为混沌理论的深入研究和混沌的应用，特别是在混沌保密通信系统和微弱信号检测以及电力系统谐波抑制等领域的应用方面，提供技术支持，并为此奠定了理论基础。 　　在文章提出的系统之上，进行基本动力学特性的理论分析与数值仿真，诸如平衡点稳定性、混沌吸引子、耗散性和维数与Lyapunov指数等；并把该系统转化为实际的物理电路模型，所采用的是模块化设计法，并采用电子工作平台EWB软件对新的混沌系统进行电路仿真实验来进行验证。 　　1 新混沌系统及其基本动力学特性分析 　　文章提出的新三维混沌系统其状态方程为： 　　（1） 　　其中x，y，z为状态变量，a，b，c，d为系统参数，该系统存在一个非线性项y2，是一个三维二次自治系统。 　　系统中有一个混沌吸引子的条件是当a=2，b=0.18，c=2，d=2时， 对该系统采用MATLAB工具进行数值仿真实验，即可得到三维空间xyz相图、y变量的时域波形和y-z平面相图，例三维空间相图如图1（a）所示。 　　（a）混沌吸引子三维空间的相图 （b） Lyapunov指数谱 　　图1 　　为求系统的平衡点，令： 　　（2） 　　解方程组可求得系统的两个平衡点为P1（0，0，0）和P2（-1.91，-1.91，-3.82），在平衡点处，对应的Jacobian矩阵分别为： 　　（3） 　　由特征方程det[？姿I-J]=0，得J（P1），J（P2）的特征值分别为： 　　可见？姿1，？姿1为负实根，？姿2，？姿3和？姿2，？姿3是共轭复根，并且这两对复根均有正实部，故P1，P2平衡点是不稳定的鞍焦点且指标均为2。 　　系统（1）之散度： 　　（4） 　　当b=0.18，c=2时，？荦V0，因此系统（1）为耗散系统，即当t→∞时，小体积元-包含系统轨线以指数（-b-c+1）的速率收缩至零，所有系统轨线的渐近运动将会被固定在一个吸引子之内，其被限制于一个体积为零的集合。 　　混沌吸引子的相邻轨线之间具有一定的作用，它们具有彼此相排斥的趋势，此趋势会以指数的速率分离，对这种作用的描述采用Lyapunov指数，对Lyapunov指数的计算我们采用Jacobian矩阵法，可得此系统的分别为：？姿L1=0.2041，？姿L2=0，？姿L3=-1.3858，Lyapunov指数谱如图1（b）所示。可得Lyapunov维数： 　　（5） 　　由上对该系统相图、平衡点的稳定性、耗散性和Lyapunov指数与维数的理论分析结果得出，新系统为一混沌系统。 　　2 EWB电路仿真实验 　　为了确定电路是否可行，若可行电路可以优化，现在对所提出新系统（1）来进行电路设计，模拟乘法器、线性电阻、运算放大器和电容器等是电路的主要组成的部分，将系统转化为物理模型，转换之后对电路进行仿真实验，利用的是EWB仿真软件。 　　电路的各组成部分中用以做积分运算和加减运算的是运算放大器，用以实现系统的非线性项部分y2的是模拟乘法器，连接三个x，y，z的状态变量为一个整体。 　　设计过程中采用模块化的设计方法，主要考虑以下几个问题（1）有源器件的动态变化范围是否小于状态变量的取值变化范围，若不小于就不用处理，否则需要采用比例压缩变换操作。根据图l，变量的动态范围并未超出运算放大器TL082的线性范围（-13.5V，+13.5V），故不需要做变量比例变换。（2）更适合用于模块化设计的是反相加法或其运算电路，因为参数独立可调、便于调试是该类电路的主要特点，这是与同相加法或减