liu系统的混沌特性及其matlab仿真

1、Liu系统的混沌特性及其Matlab仿真 目 录 摘要 . 2 英文摘要 . 3 导言 . 4 第一章 混沌的定义及其相关理论 . 6 1.1 混沌学简史 . 7 1.2 混沌的定义 . 8 1.3 混沌的基本特征 . 8 1.4 混沌的主要研究方法 . 9 1.4.1 功率谱 . 9 1.4.2 李雅普诺夫（Lyapunov）指数 . 9 1.4.3 其它方法 . 10 1.5 混沌意义及应用 . 10 1.5.1 混沌的意义 . 10 1.5.2 混沌的应用及前景 . 10 第二章 Liu系统动力学行为的分析 . 11 2.1 基本动力学行为分析 . 11 2.1.1 对称性和不变性 .

2、12 2.1.2 耗散性和吸引子的存在性 . 12 2.1.3 平衡点及稳定性 . 13 2.2系统的参数影响 . 14 第三章 Liu系统的Matlab仿真 . 20 3.1 定义函数 . 20 3.2 ODE函数命令作图 . 21 3.2.1 吸引子图及其程序 . 21 3.2.2 时序图及其程序 . 21 3.2.3 相图及其程序 . 22 3.2.4 参数改变时的相图及其程序 . 22 第四章 Liu系统的功率谱分析 . 23 4.1.经典功率谱估计 . 23 4.1.1 Barlett 法 . 24 4.1.2 Welch 法 . 24 4.1.3 Nuttall法 . 25 4.2

3、 Liu 系统混沌的功率谱仿真 . 25 结论 . 27 参考文献 . 28 附录 . 29 附录（1） Liu系统的吸引子图，三维图程序 . 29 附录（2） Liu系统混沌时序图程序 . 29 附录（3） Liu系统混沌相图程序 . 30 附录（4） 经典功率谱分析方法比较 . 32 附录（5） Liu混沌系统功率谱MATLAB仿真程序 . 35 致谢 . 36 左自豪 指导老师：高心 西南民族大学电气信息工程学院电气041班 摘要：混沌现象几乎涉及到科学研究的每一个领域。物理力学的Lroenz模型、Rossler模型、Duffing系统，电子工程学的蔡氏双涡旋电路模型（连续动力学系统范例

4、），电力系统模型，生物学的Logistic映射，天体物理学的Henon映射（离散动力学系统范例）等等，这些范例表现出丰富的混沌行为。Liu混沌系统结构不同于以往的连续混沌系统，它是一类含有平方非线性项的三阶连续自治混沌系统1,2。本文采用相图，功率谱，和Lyapunove理论等研究混沌的一些方法,并借助MATLAB软件对之进行仿真研究,观察研究Liu混沌系统的基本动力学行为特性，良好的仿真结果验证了本文算法的有效性和快速性。 关键词：Liu混沌系统；系统动力学；吸引子；MATLAB The chaos of Liu system characteristic and its simulatio

5、n by Matlab Zuo Zihao Instructor: Gao Xin Southwest University for Nationalities Institute of Electrical and Information Engineering Electrical engineering and automation 041 classes Abstract: Chaotic phenomena related to scientific research in almost every area. Lroenz physical mechanics of the mod

6、el, Rossler model, Duffing system, electronic engineering Chuas double vortex circuit model (for example dynamic system), electric power system model, the Logistic map biology, astrophysics the Henon map (discrete dynamics Example system) and so on, these examples show the rich chaotic behavior. Liu

7、 chaotic system structure is different from the previous consecutive chaotic system, it is a kind of square with the third-order nonlinear of the chaotic system for autonomy. In this paper, phase diagram, the power spectrum, and Lyapunove Chaos Theory, and so on a number of methods and use of MATLAB

8、 simulation software for research, observation of Liu chaotic system dynamics of the basic characteristics of a good simulation results show the effective this algorithm And fast. Keywords: Liu chaotic systerm ; systerm dynamics; attractor; Matlab 导 言 自 20 世纪60 年代，洛伦兹发现混沌现象以来，混沌理论研究一直受到普遍的关注。随着混沌研究的

9、不断发展，人们开始把目光聚焦在控制混沌和利用混沌的研究上，控制和利用混沌的研究都是基于一些典型的混沌系统来进行的。 1999年Chen等采用线性反馈控制方法控制Lorenz混沌系统而发现了一种与Lorenz混沌系统类似但不拓扑等价的Chen混沌系统；2001年和2002年,吕金虎等人相继发现L混沌系统和连接上述三个混沌系统的统一混沌系统 沌系统的存在。 2004年刘崇新等人发现了新的Liu混沌系统。Liu系统是一个三维连续混沌系统，开展其动力学特性及应用研究具有重要的理论意义和实际价值。 通常对混沌系统一个完整的研究要使用到李雅普诺夫指数分析和李雅普诺夫指数谱，庞卡莱映射图，混沌吸引子图，混沌

10、相图，混沌时序图，混沌功率谱分析，混沌的电路实现等等。在对其动力学行为的分析中要考察研究其对称性，稳定性，不变性，平衡点和参数变化时的性质等等。 对于混沌系统的研究，我们要使用多种方法不同的角度来分析研究，同时在理论分析中结合仿真和电路设计实现，其中就常常用到MATLAB软件和EWB电子设计软件进行辅助验证分析结果。 混沌的产生，控制，同步等研究领域是目前最重要最前沿的研究体系。因为我们最重要的是实践也就是应用，目的是要应用转化为生产力。 自从美国海军实验室的Pecora和Carrol 提出了一种混沌同步方法2，以及随后在电子电路中首次观察到混沌同步现象以来，人们先后提出一系列有效的混沌同步控

11、制方法，例如，相互耦合同步,续变量反馈同步，自适应控制同步等，这些方法各有自己的使用范 围，其中由于线性反馈控制同步简单且容易在物理上实现而得到广泛的应用, 此方法主要是基于李雅谱诺夫稳定性理论，构造李雅谱诺夫函数，进行数学推导，得出同步控制参数的取值范围, 但是对于某些混沌系统，要构造其李雅谱诺夫函数并不是一件很容易的事情, 因此，研究采用其他方法来确定线性反馈控制同步的控制参数取值范围将具有重要的理论意义和实用价值。 近年来，随着人们对混沌现象的深入研究，对其动力学行为和基本特性的逐步了解，在图像处理、保密通信、电力电网动态分析和保护、机械振动分析与故障诊断、电子振荡发生器设计、信号检测与

12、信息处理等领域中已得到了有效的应用，随着混沌理论的不断发展和完善，混沌将会在更多的领域中得到广泛的应用在这些应用中，我们需要有目的地控制混沌，也需要有目的地生成混沌或者加强已存在的混沌行为。 122,72；2003年,Liu等发现了在三维连续自治混沌系统中能产生四螺旋混沌吸引子的混沌系统,并用实际的硬件电路证实了该混 近年来混沌系统已经从整数阶发展到分数阶了，其性质和应用将更广阔，更丰富。分数阶微积分是研究分数阶次的微积分算子特性及其应用的数学理论，它几乎与整数阶微积分理论具有同样长的发展历史，但由于分数阶微积分理论长期没有实际应用背景而发展缓慢),次指出了自然界及许多科学技术领域中存在大量的

13、分维数的事实，且在整数阶微积分与分数阶微积分理论描述的动力学系统之间存在着自相似现象，此后，作为分形几何和分维数的动力学基础，分数阶微积分才重新获得了新的发展而成为当前国际上的一个研究热点，并在电磁振荡、系统控制、材料力学等领域得到了有效的应用)近年来，分数阶混沌系统的电路实现及其应用已引起人们广泛的兴趣和深入的研究，整数阶微积分是分数阶微积分理论的特例，整数阶混沌系统都是对实际混沌系统的理想化处理分数阶微积分是整数阶微积分理论的推广，利用分数阶微积分算子能更准确地描述实际混沌系统的动力学特性) 特别是在最近，在Lorenz混沌系统、Chen 混沌系统、Chuas 混沌系统、Liu混沌系统以及

14、Rossler超混沌系统中，通过计算机数值仿真，发现当系统的阶数为分数时仍然出现混沌状态，且更能反映系统呈现的工程物理现象,促进了人们利用分数阶微积分理论更深入地研究混沌这一自然界普遍存在的物理现象，也促进了分数阶微积分理论的发展。 MATLAB是集数值运算、符号运算、数据可视化、数据图文字统一处理、系统动态仿真等功能于一体的数学软件具有很高的编程效率,在线性代数、矩阵分析、数值计算及优化、系统动力学、建模与仿真等领域中得到广泛应用3。MATLAB具有高效的数值计算和符号计算功能，能使用户从繁杂的数学运算分析中解脱出来，同时还拥有完备的图形处理功能，实现计算结果和编程的可视化，它的这些特点使得

15、很快成为应用学科计算机辅助分析、设计、仿真等领域不可缺少的基础软件。 混沌理论研究的是非线性问题,难以用解析式表达,只能采用数值解法,而MATLAB在这方面便可展示其强大的潜能3。本文通过采用相图,时序图和功率谱,只借助MATLAB软件对Liu系统通向混沌的途径进行仿真研究,观察状态变量在时域和频域中的变化来了解系统的非线性特性,并通过调整其控制参数观察Liu系统动力学行为的演变过程验证其特性。 第一章 混沌定义及其相关理论 1.1混沌学简史 混沌一词译自英文chaos，chaos一词来自希腊文 ，其原意是指先于一切事物而存在的广袤虚无的空间，后来罗马人把混沌解释为原始的混乱和不成形的物质，而

16、宇宙的创造者就用这种物质创造出了秩序井然的宇宙4。 4,6自从牛顿三大定律及万有引力定律问世以来，确定论的思想就在人们心中根深蒂固，这种观点是伟大的法国数学家和自然哲学家拉普拉斯强烈主张的。在牛顿力学中，这 种信念是正确的，并且避免了任何可能的混合和含糊。但是，在真实世界里，初始状态的精确知识是得不到的，一个量不管测量得多么精确，我们总能要求测量得更精确些。尽管一般说来，我们没有能力知道这种精确知识，但通常我们假定，如果两个分别进行的实验的初始条件几乎相同，则最后结果亦将几乎相同。对于大多数具有光滑特性的常规系统，这种假设是正确的，但对于某些非线性系统，它是错误的，并且结果是确定性的混沌。 早

17、在上世纪末，伟大的法国数学家庞加莱就已经深刻地了解了这种可能性以及定量的结果。在他的科学与方法一书中写道：“初始条件中的微小差别会在最后现象中产生非常大差别的情况也可能发生，前者的微小误差将在后者中产生巨大的误差。预言变为不可能，而我们就有了偶然现象”。尽管庞加莱有惊人的洞察力，但直到本世纪年代初期，确定性混沌实际上仍然没有被仔细考察过。 交互式计算机的诞生终于为细致研究混沌提供了有力的工具。 1963年，气象学家洛仑兹根据牛顿定律建立了温度压强、压强和风速之间的非线性方程，他将该方程组在计算机上进行模拟实验，因嫌那些参数小数点后面的位数太多，输入时很繁琐，便舍去了几位，尽管舍去部分看来微不足

18、道，可是结果却大大出乎意料：该气象模型竟与没有舍去几位小数所得的气象模型大相径庭，变得完全不同。因此，洛仑兹断言：“长时期”天气预报是不可能的。在澳大利亚的一只蝴蝶偶然扇动翅膀所带来的微小气流，几星期后可能变成席卷北美佛罗里达洲的一场龙卷风。这就是天气系统的“蝴蝶效应”。 80年代后，混沌理论的研究一下子成为了热点，不仅是数学家、物理学家，而且生物学家、化学家、医学家、经济学家都不约而同地寻找不同形式的无规则性之间的联系。混沌之所以有如此大的吸引力，是因为它提供了把复杂的行为理解为有目的和有结构的某种行为的方法，而不是理解为外来的和偶然的行为。混沌是一种关于过程的科学，而不是关于状态的科学；是

19、关于演化的科学，而不是关于存在的科学，它使人们看到了运动演化中的生机和动力。 54 1.2 混沌的定义 混沌是学术界对非线性系统研究领域非常活跃的前沿课题,混沌现象几乎涉及到科学研究的每一个领域。他是指确定性系统中出现的一种类似随机过程的行为。一个非线性动力学系统,在系统参数达到一定匹配时便会出现混沌现象4。 混沌是J.Hadamard在19 世纪末研究Hamilton系统时发现的。研究热潮始于1963年Lorenz的三阶自治系统（Lorenz模型）。在Lroenz系统族中，常见的有以下三中典型模型2： Lorenz系统： Chens系统： x?a(y?x)y?cx?xz?yz?xy?bz (

20、a=10, b=8/3, c=28) x?a(y?x)y?(c?a)x?xz?cyz?xy?bz (a=35, b=3, c=28) L系统： x?a(y?x) y?xz?cy z?xy?bz ? (a=36, b=3, c=20) 混沌现象属于确定性系统且对初值极其敏感，具有稠密轨道的拓扑特征，并呈现多种“混乱无序却又颇有规则”的图像 存在的必要条件。 对于混沌的概念, 这个是很难确切地定义出来的,一般地,一个系统会受到小的扰动，非线性会放大这些扰动，系统一方面对于初始条件敏感依赖，另一方面又在有限范围内运动，就使那些初始状态和速度充分接近的轨道会以指数速度分离开来。由于轨道自己不能相交，所

21、以这些轨道只能在有限的空间内缠绕往复而形成非常复杂的形状,这就是混沌。 1.3 混沌的基本特征 混沌的基本特征有如下几点： (1) 具有类随机信号的特性。 (2) 对初始值非常敏感。 (3) 其相关函数类似于随机信号的相关函数，具有类似冲击函数的特性。 (4) 混沌信号的频谱与随机信号的频谱类似，表现为连续频谱。 (5) 混沌信号在相空间的吸引子表现为几何结构非常复杂具有分数维的奇怪吸引子。 (6)混沌吸引自具有正的李雅普诺夫（Lyapunov）指数。 混沌运动的最基本的特点是运动对初始条件极为敏感。两个很靠近的初始值所产生的轨迹，随时间的推移按指数方式分离。李雅普诺夫（Lyapunov）指数

22、是定量描述这一现象的量，在李雅普诺夫（Lyapunov）指数谱中，最小的Lyapunov指数，决定轨道收缩的的快慢；最大的Lyapnov指数则决定轨道发散即覆盖整个吸引子的快慢 轨线的分离问题，也就是“蝴蝶效应”的强弱的量。 2,776,7。因此混沌是一个关于过程的科学而不是一种关于状态的科学，是关于演化的科学而不是关于存在的科学，决定论规律的非线性是混沌运动。所有的指数之和可以认为是大体上表征轨道总的平均发散快慢。而单个的则是定量的表征相空间两相邻 由以上的介绍，我们可以判断混沌，判断噪音。下面介绍几种关于混沌时间序列判别的基本方法： （1）功率谱分析法 （2）主分量分析 （3）Poinca

23、re庞卡莱截面法 （4）Lyapumov指数 （5）C-C方法 （6）局部可变神经网络法 （7）指数衰减法 （8）频闪法 （9）代替数据法 1.4 混沌的主要研究方法。 对于混沌现象的客观反映就需要更严谨的数学描述以及更直观的物理现象。 1.4.1 功率谱 在实验分析方面，对混沌系统施行功率谱分析已是研究混沌的最有效、最直观的工具。 周期运动在功率谱中对应尖锋，混沌的特征是谱中出现噪声背景和宽锋。它是研究系统从分岔走向混沌的重要方法。 在很多实际问题中（尤其是对非线性电路的研究）常常只给出观测到的离散的时间序列X1, X2, X3,.Xn,那么如何从这些时间序列中提取前述的四种吸引子（零维不动

24、点、一维极限环、二维环面、奇怪吸引子）的不同状态的信息呢？ 我们可以运用数学上已经严格证明的结论，即拟合。我们将个采样值加上周期条件 Xn+i=Xi，则自关联函数（即离散卷积）为 然后对Cj完成离散傅氏变换，计算傅氏系数。 k说明第k个频率分量对Xi的贡献，这就是功率谱的定义3,8,。当采用快速傅氏变换算法后，可直接由Xi作快速傅氏变换，得到系数然后计算，由许多组Xi得一批Pk，求平均后即趋近前面定义的功率谱Pk。考虑到实际计算中，数据只能取有限个，谱也总以有限分辨度表示出来，从物理实验和数值计算的角度看，一个周期十分长的解和一个混沌解是难于区分的，这也正是功率谱研究的主要弊端3,8,9,10

25、,11。 1.4.2 李雅普诺夫（Lyapunov）指数 确定混沌区后，需要进一步对吸引子进行刻画。功率谱分析仍然有用，但更重要的是计算李亚普诺夫指数。 对初始条件的敏感依赖性是确定性系统混沌的关键特性。这意味着在相空间中相互靠近的两条轨线，随着时间的推移，它们将指数性的运动开。 李雅普诺夫指数的数值计算方法一般有如下四种： (1).定义法 (2).Wolf 方法（基准轨道法） (3).切空间法 (4).小数据量法 （M. T. Rosenstein，J. J. Collins & G. J. De Luca 1993，H. Kantz 1994）此方法对数据要求较低，计算精度有明显的

26、改善，而且比较稳健。 该李雅普诺夫指数的大小刻划了吸引子的动力学，反映了系统产生或消除不确定因素的速率。初始不确定性经过多少时间将覆盖整个吸引子由最大的指数决定，而对吸引子摄动的渐近消失则由最小的指数控制。 1.4.3 其它方法 描述混沌程度的方法有很多种，而且各有利弊。其中较主要的还有如确定奇怪吸引子的各种维数、确定混沌系统的所谓柯尔莫哥洛夫熵、对周期驱动系统较易实现的分频采样法以及前面介绍的取庞加莱截面的方法等等 似结构的泛函方程形式的重正化群方程 也是确定混沌的一个重要方法。 1.5 混沌意义及应用 1.5.1 混沌的意义 混沌现象几乎涉及到科学研究的每一个领域。物理力学的Lroenz模

27、型、Rossler模型、Duffing系统，电子工程学的蔡氏双涡旋电路模型（连续动力学系统范例），电力系统模型，生物学的Logistic映射，天体物理学的Henon映射（离散动力学系统范例）。这些范例表现出丰富的混沌行为。 1.5.2 混沌的应用及前景 根据混沌系统提取的非线性时间序列对系统的未来进行预测。如在电力系统短期负荷预测中的应用，电机转子剩余寿命的预报，计算机软件失效的预测，在天气、水文预报中的应用以及证券市场股价波动的混沌分析。近年来，随着人类对混沌吸引子现象的不断深入探索，对其动力学行为和基本特征的逐步了解，混沌在工程领域得到了很大的发展，混沌控制和混沌同步研究的也有突破性进展，

28、因此大大推进了混沌的应用研究。 （1） 改善和提高激光器的功效（性能和效率）。利用混沌控制技术研制品质优异、单色性好、高功率的激光装置。 （2） 混沌信号在通信中的应用，包括混沌保密通信、混沌载波数字通信、混沌序列调频扩频和直接序列扩频通信、混沌参分多址通信等领域。 （3）在其它科技领域中的可能应用，利用时空混沌作为信息源研制超高容量的动态信息存储器；未来能源受控热核聚变中的等离子体的混沌、揣流的控制约束；混沌控制与混沌同步与生命科学如神经网络、脑科学、心脏等领域的研究；人们在众多领域发现混沌是有用的或者有着巨大的应用前景。 当混沌有害时，抑制混沌动力学行为（混沌控制）；当混沌有用时，刻意产生

29、或加强混沌动力学行为（混沌反控制）；混沌的研究主要包括理论研究、混沌模型的建立和混沌信号在工程学上的应用三方面的内容：混沌发生器（算法、硬件）；分数阶混沌的研究；混沌信号的处理、控制和应用。 42。理论分析中则大量使用泛函分析，借用相变理论中的重正化群方法。比如费根鲍姆推出的普适常数及标度因子就是解一个反映自相。另外，分析系统的梅尔尼科夫(Melnikov)函数 第二章 Liu系统动力学行为的分析 混沌系统有着复杂的动力学行为，目前人们已知的混沌吸引子并不多，Liu混沌系统是一类含有平方非线性项的混沌系统,其系统的动力学方程为: ?x?a(y?x) y?bx?kxz z?cz?hx2 其中a,

30、b,k,c,h为系统参数,当a=10,b=40,k=1,c=2.5,h=4时,系统(1)处于混沌状态，系统的混沌吸引子和功率谱如下图所示。由于Liu系统是一个新的混沌系统,开展其动力学特性及应用研究具有重要的理论意义和实际价值。利用Jacobia方法计算李雅普诺夫指数得： 2,7?(1)? ?L?1.64328,?L?0,?L?14.42 123 固定k=1,h=4不变，因此,从系统的李雅普诺夫指数,可知该系统为混沌系统。 Liu系统混沌吸引子图 Liu系统混沌功率谱图 2.1 基本动力学行为分析 Liu系统与Lorenz系统和Chen系统具有许多相似的基本性质。 2.1.1 对称性和不变性

31、由于系统(1)模型在变换( x , y , z ) (- x,- y ,z)下的不变性，即系统的图象关于z轴具有对称性，且这种对所有的系统参数保持不变。 2.1.2 耗散性和吸引子的存在性 由于 ?x?y?z?V? ?a?c(2) 当a+c>0时,则系统(1)是耗散的,且以指数形式收敛 dV?(a?c)?edt(3) 即体积元V0在时刻t时收缩为体积元V0e?(a?c)t,这意味着,当t时包含系统轨线的每个体积元以指数率-a-c收缩到零。因此,所有系统轨线最终会被限制在一个体积为零的集合上,且它渐进运动固定在一个吸引子上，说明了吸引子是存在的。 2.1.3 平衡点及稳定性 令系统(1)式

32、的右边等于零,即 a(y?x)?0bx?kxz?0?cz?hx2?0 (4) 得三个平衡点为： S0?(0,0,0)S?S? ,b)?(?,?,b) 22 由其平衡点可知,a,c不变且都为正数,当b>0时,系统有三个平衡点;b<0时,系统仅有一个平衡点,可见在b=0处,则出现叉式分岔。 将系统(1)在平衡点处线性化,得其雅可比（Jacobi）矩阵为： ?a? J?b?z0 ?8x 0? 其特征方程为： a00 0? ?x0?c? (5) f() =(+ c)( ?2 + a + az0 - ab) + 8ax02 = 0 ( 6 ) 把S0=(0,0,0)代入(6)式,得其特征多项

33、式为 f0(?)?(?c)(?2?a?ab)?0 其特征根为： （7） ?a?a; ?3?. ? 1 = -c; ?2? 22 若b>0,则a2?4ab?a2,从而?1<0, ?2>0, ?3<0，则S0=(0,0,0)平衡点是不稳定的且为三维空间中的一个鞍点。 22a?4ab?a若b<0,则,从而?1<0, ?2<0, ?3<0.则S0=(0,0,0)平衡点是（1）唯一的平衡点且是稳定的汇。 把S?代入(6)式,得其特征多项式为: f?(?)?(a?c)?ac?2abc?0 32 （8） 由于a,c为正数,若b>0,则(8)式的所有系数均

34、为正,因此,对任意>0,则都有 f?(?)>0。因此,要使平衡点不稳定,则(8)式必须有两个正实部的复共轭特征根。若b=0,则(8)式有三个特征根: ?1=0, ?2?a, ?3=-c.当b从正趋近于零时,则的高阶项变的很小,忽略不计,则?1趋近于-2b,因此,在极限状态的情况下,稳定性将失去。在b从零逐渐增加的过程中,仅当Re(?)?0时,系统不稳定性才会发生。此时，(8)式的两个特征根则为 ?1,2?i,其中w为实数。 从而可得: ?1+?2+?3=-(a+c),因此, ?3=-a-c,把它代入(8)式,有-ac(a+c)+2abc=0,即 bh=(a+c)/2,若b<b

35、h,则S?和S?均为稳定的汇.若b>bh,则S?和S?均为不稳定的鞍焦点。 若b=bh时,两个复共轭特征根为?1,2?aci。 由于 ?2ac ?(b)?2 3?2(a?c)?ac 因此 （9） ac ?(0)?Re?(bh)?0,2 ac?(a?c) (a?c?0.2 ac?(a?c) (10) ?(0)?Im?(bh)? (11) 所以,根据霍普夫（Hopf）分岔的定理,则当b=bh时,系统在平衡点S?和S?处将出现霍普夫(Hopf)分岔,bh表示霍普夫（Hopf）分岔点。 2.2系统的参数影响 Liu系统的动力学行为是由控制参数决定的。为观察系统的动力学特性,可采用相图、功率谱、有

36、选择地研究控制参数，当参数变化时研究Liu系统的演化过程。 从上面分析可知,随着系统参数的改变,系统平衡点的稳定性将会发生变化,从而该系统也将处于不同的状态。下面利用仿真,分析各个系统参数变化时,系统的变化情况。 （1） 当a、c不变,b改变时，依据下面的b取不同值的相图可知： 当b< bh?6.25时,系统为稳定状态如（b=5的xz图）所示； 1 当b> b?6.25 时,系统为混沌状态如（b=7和b=20的xz图）所示。 b取不同值时的相图(分别依次取b=5,b=7,b=20) 很显然,系统在b= bh?6.25时,发生了霍普夫（Hopf）分岔。同时由图（b=7的xz图和b=2

37、0的xz图）比较还可以看出,随着b的增大,系统的混沌性增强。 （2） 当a,b不变,c改变时，依据下面的c取不同值的相图可知： 当0<c8时,系统为混沌状态； 当8<c<10.3时,系统为周期状态； 当10.3c<10.7时,系统为拟周期状态； 当10.7c<30时,系统为周期状态。 C取不同的值时的相图（分别依次取c=6,c=9,c=10.5,c=20） （3）当b、c不变,a改变时, 依据下面的c取不同值的相图可知: 当0<a4时,系统为周期状态； 当4<a14时,系统为混沌状态； 当14<a18时,系统为周期状态； 当18<a42.5

38、时,系统为混沌状态； 当42.5<a50时,系统为周期状态。 a取不同的值时的相图(分别依次取a=2,a=8,a=17,a=20,a=45) 综上所述，参数变化对系统影响的分析可以得到以下的结论： (1)Liu混沌系统在b=0,出现叉式分岔。 (2)Liu混沌系统在b=(a+c)/2,出现霍普夫分岔。 (3)Liu混沌系统当a,c不变,b改变时,随着b的增大,系统从稳定状态变为混沌状态,且系统的混沌性增强。 (4)Liu混沌系统当b,c不变,a改变时,随着a的增大,周期状态和混沌状态交替出现。 (5)Liu混沌系统当a,b不变,c改变时,随着c的增大,系统则出现混沌、周期和拟周期三种状态

39、。 第三章 Liu系统的Matlab仿真分析 Liu系统的动力学行为是由控制参数a、b、c决定的。为观察系统的动力学特性,可采用相图、功率谱、关联维数等方法有选择地研究控制参数，当参数变化时研究Liu系统的演化过程。 相空间就是由研究的物理量本身作为坐标分量所构成的广义空间,系统的任意状态相当于相空间中的一个点,系统状态随时间变化的过程对应于点在相空间中的变化,所有点的集合便构成了相图。非线性系统随时间的演变将趋向于维数比原来相空间低的极限集合即吸引子。通常的简单吸引子有不动点、极限环和环面,随系统参数的改变简单吸引子可发展为奇怪吸引子。像这种当控制参数变化到某个临界值时使系统的动力学行为更易

40、控制。 Matlab是一种既可交互使用又能解释执行的计算机语言。在研究非线性方程的解中使用方便直观，功能强大。通过Matlab仿真可以帮助我们很好的分析，研究和理解混沌系统的动力学行为特性9。 下面使用仿真分析Liu系统的动力学行为特性： 用Matlab实现Liu混沌系统的吸引子（状态图），结果见第一章系统混沌吸引子图。 3.1定义函数 为了方便表达将x,y,z表示为x(1),x(2),x(3)也即是列向量x中的3个分量，编译函数程序文件lu19pride_.m,具体程序如下： %设x=x(1),y=x(2),z=x(3); function xdot=liu19pride(t,x); a=1

41、0; b=40; k=1; c=2.5; h=4; xdot=a*(x(2)-x(1);b*x(1)-k*x(1)*x(3);-c*x(3)+h*x(1)2; 其中，a,b,k,c,h是Liu混沌系统的参数，参数值可以在Matlab装载m程序文件夹中直接修改和设置，保存便能实现参数的改变，方便编程仿真分析，此段程序主要是对系统函数方程的描述，定义Liu混沌系统的方程组，方程中x,y,z是对时间t的导数，因此是用Matlab求解常微分方程的解图像曲线，便要使用到ODE命令，在此采用Rung-Kutta(4,5) 公式的ode45调用函数。 3 3.2 ODE函数命令作图 3.2.1 吸引子图及其程序 在命令窗口编译下面程序实现Liu混沌系统的吸引子图，用plot3函数创建系统三维图型，其具体的图形见第一章系统吸引子图，它的具体编译程序见附录（1）。 3.2.2 时序图及其程序 时序图是各方程在时间域的时间序列图，其具体的编译程序见附录（2）： 程序中使用 axis函数设置图形的坐标范围，便于得到体现系统特征的仿真图形，从而更好地方便分析和研究系统的特性。而subplot函数则是选择绘图位置，用于在图形窗口上绘制多个子图。如