MATLAB、Simulink混沌理论仿真

1、毕业设计（论文）原创性声明本人郑重声明：所提交的毕业设计（论文），是本人在导师指导下，独立进行研究工作所取得的成果。除文中已注明引用的内容外，本毕业设计（论文）不包含任何其他个人或集体已经发表或撰写过的作品成果。对本研究做出过重要贡献的个人和集体，均已在文中以明确方式标明并表示了谢意。 论文作者签名： 日期： 年 月 日摘 要混沌在现代科学与工程学领域的应用十分广泛，混沌现象存在于自然界各个领域，包括通讯领域、气象学领域、生物学领域、医学诊断疾病等方面。学习混沌理论在未来的发展过程对我们是很有帮助的。在非线性的世界里，通过混沌理论洞察所有的非线性运动,对其进行控制和掌握。通过非线性电路对混沌系

2、统进行分析和理解，进而构造出符合二阶混沌系统的非线性电路和函数模型。Duffing方程就是典型的二阶非线性方程。运用MATLAB/Simulink对其混沌系统进行仿真实现，验证混沌系统的基本特性。关键词：混沌；非线性；Duffing方程; MATLAB/SimulinkABSTRACT Chaos widely used in modern science and engineering and chaos phenomenon exists in various fields of nature, including the communications field, the field o

3、f meteorology, biology, medical diagnosis of diseases. Learning Chaos Theory is very helpful to us in the development of this course in the future. In a nonlinear world, insight into the chaos theory, We can control and master non-linear movement. We analyze and understand the chaotic system via non

4、linear circuit, and then construct a second-order chaotic systems of nonlinear circuits and function model. Duffing equation is a typical second-order nonlinear equation. Using MATLAB/Simulink, we complete the chaotic system simulation and test the basic characteristics of chaotic systems.Key words：

5、Chaos；nonlinear；Duffing equation；MATLAB/Simulink目 录第一章 绪论11.1混沌理论11.2混沌的应用2第二章 二阶混沌系统的仿真实现52.1混沌系统52.1.1混沌产生的数学模型52.1.2 奇异吸引子与分形62.1.3 混沌系统的特征72.1.4 研究混沌的主要方法82.2 二阶混沌系统的实现9第三章 二阶非线性电路仿真实现153.1 Simulink仿真173.2 MATLAB语句命令演示模拟19第四章 结论22致 谢25参考文献26附录A27山东英才学院2012届本科生毕业设计(论文)第一章 绪论1.1混沌理论什么是混沌？现代科学意义上是很

6、难得出确切的定义，之所以这样是因为：到目前为止，还没有足够和统一数学定理可以将混沌理论完全表达出来，在数学理论的基础上通过混沌系统所表现出的普遍现象总结归纳出混沌的本质。对此，很多科学家给出很多观点费根包姆：“确定系统的内在随机运动。”洛仑兹：“确定性非周期流。”哈肯：“混沌性为来源于决定性方程的无规运动。”赫柏林：“没有周期性的有序。”1钱学森：“混沌是宏观无序、微观有序的现象。”等等。综上所述，我们可以对混沌的定义作出如下理解：混沌是指非线性系统在一定条件下所呈现的不可预测的随机现象；是将有序性与无序性融为一体的现象；其无序随机性不是来源于外部干扰，而是来源于内部的动力学方程中的非线性项，

7、正是由于非线性系统在一定的临界性条件下其对初值的敏感性表现出混沌现象，才导致内在的不稳定性的综合效果。通常我们把研究的对象称为系统，因此基于混沌的研究上我们把混沌称为非线性系统运动。正因为如此，我们所讨论的对象必然是非线性系统，或者确切地说是非线性动力系统。如直线函数就是一个最简单的线性函数，变量与自变量成一次方的函数关系就是在（x,y）平面中是一条直线。而函数y=f(x)对x的依赖关系高于一次，就像抛物线函数（其中y项是非线性项），那么这个函数所描述的系统就是“非线性系统”。可见，从函数构造的角度来说，非线性系统要比线性系统更复杂。线性系统与非线性系统的不同之处至少有两个方面。第一：线性系统

8、可以使用叠加原理，而非线性系统不能使用相关原理；第二：非线性系统对初值极敏感，而线性系统只对自变量有依赖关系。正如在运动形式上，线性现象一般表现为时空中的平滑运动，并可用简单的函数关系表示，而非线性现象则表现为从规则运动向不规则运动的转化和跃变。非线性系统是复杂不确定的，对它的进一步研究需要新的方法和思维方式。随着科学的理论研究，适时出现系统论、信息论、耗散结构、协同学等理论，成为研究非线性系统的主要工具。混沌理论成为非线性科学的主要研究对象。非线性系统在一定条件下，会表现出一些无规性，严格讲是貌似无规性；因为在这些貌似无规性中又会出现一定的规律性来，一般就称为系统出现了混沌状态，简而言之就说

9、出现了混沌，因此给混沌下一个定义的话，简而言之混沌就是系统的无规行为中的规律性。空间现象中简单与复杂、确定与随机的内在联系引出的混沌学使人们在观念和哲学方面发生了革命性的转变，因此，混沌理论、相对论和量子论被称为20世纪最伟大的科学革命。在自然科学发展史上，被称为“近代科学之父”的意大利物理学家与天文学家伽利略首先在科学实验的基础上融会贯通了数学、物理学和天文学三门知识，扩大、加深并改变了人们对物体运动和对宇宙的认识，开创了用实验来证实科学构想的方法。而在伽利略的天文学、力学的基础上，英国的伟大科学家牛顿完成了经典力学体系，发现并总结了运动定律，创立了万有引力理论，被我们现在用于解决物体之间的

10、运动问题和力学问题。牛顿的成就也带给了科学界哲学思想上的确定论，也就是认为一个系统的行为是有规的。一个系统，只要弄清楚它所处的环境，列出它的运动方程式，给定它的初始条件，则以后的运动状况都是可以计算出来的，这种系统就是确定性的系统，有规的系统；大量的科学实验可以证明了这一论断的正确性。而在我们分析系统的过程中，往往忽略排除了一些外界干扰和对方程过分的近似等，这就是我们相对的确定性的系统。而在后来庞加莱在研究三体问题中发现方程的解的状况非常复杂，以至于对于给定的初始条件，当时间趋于无穷时没有办法预测轨道，对于轨道的长时间行为时的不确定性，数学家和物理学家称之为混沌。庞加莱的发现可以说是混沌理论的

11、最早起源了。1972年12月29日，美国麻省理工学院教授、混沌学开创人之一E.N.洛伦兹在研究气候不能精确重演与长期天气预报者无能为力发现混沌系统对初始条件敏感性的一种依赖现象在美国科学发展学会第139次会议上发表了题为蝴蝶效应的论文，提出一个貌似荒谬的论断：在巴西一只蝴蝶翅膀的拍打能在美国得克萨斯州产生一个龙卷风，并由此提出了天气的不可准确预报性。随着对混沌理论更深层次的研究，人们对混沌的认识正在逐步深化。到目前为止，人们对混沌特征的理解主要体现在如下几个方面：混沌是对初始值的敏感；混沌是拓扑传递性以及周期点的稠密；混沌是内在的随机性；混沌是奇异吸引子；混沌是吸引子的杰作，是拉伸与折叠的变换

12、；混沌的签名是分形；混沌是周期3；混沌是正的李雅普诺夫指数；混沌是拓扑熵大于零；混沌是信息之源，是信息的膨胀；混沌是局部的不稳定和整体的稳定；混沌是简单和复杂的统一。混沌正在许多方面改变着人们的自然观：混沌是有序与无序的统一；混沌是确定性与概率性的统一；混沌是稳定性与不稳定性的统一；混沌是完全性与不完全性的统一；混沌是自相似性与非自相似性的统一；混沌是一个遵循辩证法规律的世界2。1.2 混沌的应用混沌在现代科技领域的应用十分广泛，在通讯领域内，通信在我们的生活中的作用越来越重要,随着电子商务的兴起，对保密通信提出了更高的要求。正是因为混沌信号最本质的特征是对初始条件极为敏感，并导致了混沌信号的

13、类随机特性。用它作为载波调制出来的信号当然也具有随机特性。因而，调制混沌信号即使被敌方截获，也很难被破译，这就为混沌应用于保密通信提供了有利条件。因此利用混沌进行保密通信是目前十分热门的研究课题。在混沌信号对初始条件极为敏感，并由此信号又具有整体稳定性,当我们用同一个混沌信号去驱动两个相同的系统时，两个系统的某些部分将产生同步化的行为，这就为混沌应用于保密通信提供了可行性。在气象学领域内，早在1904年，挪威气象学家Bjerknes就提出天气预报问题应提成大气运动方程组的初值问题。在近年的气象研究中，利用混沌进行中期预报的研究。由于气候系统是非线性系统，其初值问题的数值解是不确定的，研究气候状

14、态的特征就要研究混沌态的特征，研究气候系统的演变机制就要研究混沌态的变化。在这些研究中使用的数学工具主要是分形理论，如分数维、李雅普诺夫指数、标度指数和功率谱指数等。利用这些数学方法分别考察、分析气候状态特征量随控制变量的变化。在数学上把天气(气候)预报问题提成初值问题，即用动力学的方法进行预报，从认识论上讲就是把大气看成是确定论的系统，这在较短的时间尺度内是行得通的，而在时间较长的时候却是有问题的，主要是大气运动是非线性、强迫和耗散的。由这三大特点，可以得到一幅这样的图像：误差是随着时间呈指数增加的，初始场的作用随着时间是衰减的，必须考虑能量的补充和耗散。E.N.洛伦兹发现了“蝴蝶效应”，指

15、的就是初始场微小的不确定性的指数放大。这就提出了确定论预报的可预报性问题，中期数值天气预报逐日预报的可预报时限大约是两三周左右的时间。也就是说进行长期预报是不可能的。在生物学领域内的研究对我们自然界以及人类的发展具有深远的意义。混沌在生物学研究中的应用主要集中在生态学中的种群变化；医学诊断疾病等方面。这也许是最后的生命图景。从还原-综合还原高度整合。生命不仅是空间结构的范畴，更是时间的范畴，千千万万种物质的时空组合和演变，构成复杂而完美的动态生命体，将生命与非生命区分开来。如何解决生命的本质，很多科学家认为应从多角度入手，不仅仅是生物学的问题，而是一项系统工程，即未来的生命科学，必需包含有数学

16、成分。研究人体这一复杂系统，混沌学、甚至弦论都不可少，只有将这些能研究复杂体系的数理方法融入医学生物学才能带来真正意义上的生命科学的突破。因为生命的最基本本质是具有时空性、混沌性。生命的整体性，包括：（1）时间上整体，即生命活动的高度有序性；（2）空间上整体，即生命结构和物质相互作用、相互影响形成网络状整体。如形态上，人体结构和功能混沌调节机制（如免疫网络调控、心脏、肺、肠的分形结构学原理，心电的混沌产生与心脏普肯野氏纤维分形分布联系等等。在物质信号相互作用的动力学研究上，其相互作用不仅仅是激活、失活或抑制、促进这一简单的关系，而应包括复杂的数学过程，这种复杂的数学过程应是非线性的，很可能符合

17、混沌原理。这是因为这一点，我们可以通过研究生物学的混沌理论从而彻底了解生命意义。 混沌现象存在于自然界的所有领域，在人类社会中的经济领域同样存在着混沌现象。南美洲热带雨林中的彩蝶轻展双翅，北美大草原竟掀起了一场风暴，这是极言世界复杂性的蝴蝶效应。对此，美国赛纳尔公司（Cerner，纽约证交所代码CERN)首席执行官尼尔帕特森对此有刻骨铭心的认识。就因为他向公司400名中层经理发出的一份电子邮件竟让公司市值在短短三天时间内猛烈下跌了两成，逾3亿美元蒸发殆尽。类似的事情在经济学领域中数不胜数。由此而应运而生了经济混沌和经济波动的非线性动力学理论。虽然混沌现象的理论和实验研究在物理学、化学、生物学、

18、天体物理、气象学以及神经生理学等广泛领域获得重要进展，但在经济学中遇到严重困难。经济活动是人的行为，动力学系统的时间尺度和观察者相近。所以经济主体和经济结构随时间的演变难以忽略。时间序列的非稳态性质使目前常用的稳态时间序列分析方法难以应用。这是为什么经济混沌的研究比自然科学更为困难。 目前混沌理论的研究主要为宏观经济运动，因为发达国家的市场经济周期的观察积累了大量数据(Zarnowitz, 1992)。在各种代表性的美国宏观经济指数中发现经济混沌的普遍证据(Chen 1996a,b)。由于目前中国的经济统计数据的收集和整理还不够充分，例证主要取之于美国的数据。但其可能的应用对中国问题同样有潜在

19、的可能。混沌理论的应用在其他领域内有诸如此下方面的现象，天体运行的长期行为不可预测；电路中的混沌现象；非线性系统的控制；利用分形研究物质结构及性能。因此对于混沌理论的研究还有待进一步探索。本次设计的主要目标是通过对二阶混沌系统的仿真实现完成对混沌系统特征理解。论文的主要内容包括：第一章：绪论。介绍一下混沌概念以及混沌理论的应用，说明课题来源。第二章：二阶混沌系统的仿真实现。主要是对混沌系统、二阶非自治铁磁谐振电路的电路结构、数学模型进行分析，为电路仿真实现提供理论基础和依据。第三章：二阶非自治铁磁谐振电路仿真实现。详细说明二阶非自治铁磁谐振电路过程及运用Matlab软件对混沌系统进行模拟仿真结

20、果。第四章：结论。在二阶非自治铁磁谐振电路的基础上，上述运用Matlab软件对混沌系统进行模拟仿真结果，讨论了仿真的结果，并介绍了混沌电路研究的意义。对在仿真实现和理论分析过程中遇到的问题和困难进行总结分析，陈述了主要几个难题及其解决方法。第二章 二阶混沌系统的仿真实现2.1混沌系统2.1.1混沌产生的数学模型随着科学的发展和进步，数学模型的应用到两个经典的例子：一个例子是二体问题，从经典力学中的开普勒问题、相对论力学中的水星近日点进动，到量子力学中的氢原子和量子场论中的兰姆谱线位移，贯穿了经典和近代物理学的全部发展史。另一个例子是在显微镜下观察悬浮在水中的藤黄粉、花粉微粒，或在无风情形观察空

21、气中的烟粒、尘埃时都会看到的布朗运动，从爱因斯坦的直观处理和朗之万方程、福克普朗克方程、到涨落耗散定理的现代表述和随机过程的连续积分表示，引出了整个物理学中的概率论描述体系。这两个例子，一个为确定论，另一个则为概率论。同样混沌系统科学研究是基于数学模型进行理论研究，对于确定论系统中的随机性，即混沌现象，在此现象下也存在着特别的模型，这就是一维迭代过程。迭代函数是重复的与自身复合的函数，这个过程叫做迭代。迭代是非线性方程演化过程的有力工具。为了研究一个物理系统，把系统的状态用一组变量x,y,z.描述，它们都是时间t的函数，同一个系统还受某些可以调节的“控制参量”的影响。最简单的情景是固定一组参量

22、，把时间变化限制成等间隔的t,t+1,t+2.，看下一个时刻的系统状态如何依赖于当前状态。在只有一个状态变量x时，这个演化过程可以由一个非线性函数描述。Logistic方程是一个很简单的非线性的抛物线函数，可以表示为： (2.1.1.1)这个函数以X的值代入得到Y后，可以用Y的值作为新的X值得出新的Y值，这样的迭代可以不断地进行下去，因此这个函数可以采用以下的迭代形式： (2.1.1.2)其中是n次迭代时的值， 是n+1次迭代时的值。这里迭代前后的值都是X值，n次迭代得到的所有的值都可以在x轴上表示出来，因此它是一维迭代。这个迭代过程是通过抛物线函数把变换（或映射）成x的过程，所以它常被Log

23、istic 映射或抛物线映射。随着参数，迭代过程显示：最终有无数点无规地出现，确定性系统中出现的类似随机的过程，系统进入了混沌状态。迭代函数的行为十分敏感地依赖于非线性的程度，即 的大小。一维非线性映射只对应于耗散系统是不可逆的，而二维映象在一维和高维之间许多方面起着衔接作用。所以说二维系统的混沌现象，出现在耗散系统和保守系统中。在保守系统中，由于系统的哈密顿函数H=常数（HT2-TO+V 式中的T2和TO分别为系统动能表示式中广义速度的二次项和零次项），系统存在一个能量积分，所以一维保守系统不可能出现混沌。二维哈密顿系统中研究较多的是所谓“标准映象”它出现在许多自由度为2的非线性振子理论中，

24、是带电粒子在环行磁场中运动的一种模型。 二维耗散系统中研究最多的一例是所谓埃农（Hénon）映射，令系统参数为b，只要b0，变换就是可逆的，b1时，它保持相体积不变，是一种保守系统。b<1对应耗散系统，b0则回到一维映射。埃农等人研究发现对于某些控制参数和初值，迭代结果迅速收敛到（x,y）平面上接近一维的“吸引子”上。这个吸引子很像是平滑曲线，但它具有宽度。如果取来吸引子的一小段不断放大，可以看到越来越小的尺度上重复出现近似的自相似结构。这是可以实际观察到的具有非整数维的“奇异吸引子”。 二维映象中在某些方面表现出了与一维线段映象不同的性质，比如：对b=1的埃农映象的研究给出的

25、结构普适常数8.7210和标度因子4.018均与一维单锋映象有所不同。但我们看到一维非线性迭代的普适特点在高维映象中仍然保持下来。2.1.2 奇异吸引子与分形保守系统由于相体积永远不变不存在吸引子，而耗散系统的相体积在演化过程中不断收缩，各种各样的运动在演化中逐渐衰亡，最后只剩下少数自由度决定的长时间行为，即：耗散系统的运动最终趋向维数比原始相空间低的极限集合，这个极限集合称为吸引子。其中包括平庸吸引子和奇异吸引子。（1）平庸吸引子 我们来考虑常微分方程解的极限集合，即相空间某一区域的点都取作初值时，这些轨道时的极限行为。极限集合的一些平庸情况是熟知的：零维不动点、一维极限环和二维环

26、面等。如果t时，系统趋向一个与时间无关的定常态，即相空间中的一个特定的点，这就是不动点。不动点是零维的吸引子。一维以上的系统原则上就可能具有不动点。如果t时，系统中剩下一个周期振动，这就是一维的吸引子极限环。只有在二维以上的相空间中，才可能出现极限环。通常极限环是由不动点发展起来的。当某个不动点在参数变变化过程中由稳定而失稳，新的稳定状态往往是围绕着原有不动点的周期运动，这个过程称为霍普夫分岔(hopf)。 （2）奇异吸引子 奇异吸引子是耗散系统混沌现象的另一个重要的特征。简单地说奇异吸引子就是相空间（对连续的动力学系统，至少是三维；对离散的动力学系统，至少是二维）的一个有限的区域内

27、，由无穷多个不稳定点集组成的一个集合体。奇异吸引子有两个最重要的特征：（1）对初始条件有敏感的依赖性。在初始时刻从这个奇异吸引子上任何两个非常接近的点出发的两条运动轨道，最终必会以指数的形式互相分离。由于混沌对初值极为敏感，它表现为局部不稳定。但对耗散系统而言，则又具有相体积收缩的特性，因而造成轨道无穷多次折迭往返。混沌轨道在相空间中“添满”有限的区域，形成奇异吸引子。实际上，它有内外两种趋向，一切吸引子之外的运动都向它靠拢，这是稳定的方向；而一切到达吸引子内的轨道都又相互排斥（指数式分离），对应为不稳定方向。正是这种整体趋向稳定而局部又极为不稳定的矛盾，导致了奇异吸引子的另一个更奇异的性质；

28、（2）它具有非常奇特的拓扑结构和几何形式。 奇异吸引子是具有无穷多层次自相似结构的、几何维数为非整数的一个集合体。为了描述奇异吸引子的这种奇特结构，曼德尔布罗特(Mandelbrot)最早（1975年）引进了分形（既其维数是非整数的对象）的概念。2.1.3 混沌系统的特征混沌系统的主要特征表现为下列几个方面（1）对初始条件的敏感性经典学说认为：确定性的系统只要初始条件给定，方程的解也就随之确定了。一个随时间确定性变化或具有微弱随机性的变化系统，称为动力系统，它的状态可由一个或几个变量数值确定。在动力系统中，两个几乎完全一致的状态经过充分长时间后会变得毫无一致，恰如从长序列中随机选取的

29、两个状态那样，这种系统被称为敏感地依赖于初始条件，这就是系统对初值的敏感，还有混沌的敏感表现在一些控制参数的变化。（2）整体稳定局部不稳定稳定性是有关扰动现象的。如果一个动力系统中发生轻微的变化，这个系统还会保持它的运动状态，保持它的能力和属性。混沌的整体稳定性指一个微小的扰动也不会改变系统原有的性能。在混沌系统中，局部不稳定性表现在混沌对初值的敏感依赖性，一个微小的初值变化就会引起系统局部的不稳定。（3）奇异吸引子及其分形奇异吸引子将混沌运动的特征初始条件的敏感性和确定性的随机直观地反映出来。在耗散系统当中，当连续流在收缩体积时，一边沿这些地方压缩，另一边又沿其他地方延伸。不过连续流是固定在

30、一个有界的区域内，这种伸缩和折叠过程会使运动轨道在奇异吸引子上产生混沌运动。可见，奇异吸引子是轨道不稳定和耗散系统相体积收缩两种因素的内在性质同时发生的现象。它的几何特性由分形来刻画，具有大尺度与小尺度之间的相似性，具有无穷无尽自相似的精细图案，具有分数维数。分形的形状是一些难以用传统的几何学来描述的极度不规则的图形；分形存在着很小的比例精密的细节结构；分形的维数大于等于它的拓扑维；分形产生于迭代过程具有自相似性，这种自相似性可以是严格的，也可以是近似的或统计意义上的。分形和混沌是同一种规律的不同表现，这种统一的规律反映在空间分布上表现为分形，出现在时间分布上表现为混沌。（4）分岔当系统的一些

31、控制参数发生变化时，新的定常状态解、周期解、拟周期解或者是混沌解就会分叉出来，其中相轨迹图发生拓扑结构的突变，分岔理论是非线性解定性行为数学理论，失稳是发生分岔的物理前提，分岔后，系统的不同状态便会有了突变，经过不断的分岔，最终达到的状态就是混沌理论的研究对象。（5）遍历性及有界性混沌运动的轨迹经历混沌吸引子内每一个状态点的地方，不重复，不紊乱。混沌的有界性最好的证明是奇异吸引子，混沌的运动轨迹虽说有点乱，但它始终在一个确定的区域里，有一定的规律性。（6）普适性若将第n倍周期分岔（或混沌带合并）时对应的参数m记为mn，则相继两次分岔（或合并）的间隔之比趋于同一个常数：=4.6692016091

32、0299067，它是一个普适常数：一类具有相同的单峰映射性质的函数中的任何一个，在沿倍周期分岔的道路进入混沌时，都会出现同一个；在沿倍周期分岔的道路进入混沌的过程中，不仅在周期区内分岔序列按速率收敛，在混沌区中的倒分岔序列也以同样的速率收敛。并证明了此种结构所具有的定量特征有着普适性，既出现于不同的非线性系统之中，又反映于同一系统的不同层次。普适性有结构普适性和测度普适性两种3。结构普适性指出无论是指数函数或是三角函数，只要是单峰映射，那么函数表现出来的结构与有着某种共同的数学性质的非线性动力系统的逻辑斯蒂方程所表现出来的结构相同，都为复杂的分岔结构。同样都是经倍周期分岔进入混沌状态。测度普适

33、性指在沿倍周期分岔进入混沌的过程中隐含着一种深刻的规律，它以常数的形式表现出来。倍周期分岔序列具有一个确定的收敛速率。费根鲍姆普适常数d 的数值只与系统的某种非线性性质有关，而与各个系统的其他具体细节无关，反映出混沌演化过程中所存在的一种普适性，说明混沌内部存在着一定的统一规律，是混沌内在规律性的另一个侧面反映，为认识和研究混沌提供了坚实的基础。2.1.4 研究混沌的主要方法我们在研究混沌的主要方法：混沌现象的客观反映就需要更严谨的数学描述以及更直观的物理现象。在实验分析方面，对混沌系统施行功率谱分析是研究混沌的最有效、最直观的工具。对采集的信号进行傅里叶变换，得到相关的频谱。周期运动在功率谱

34、中对应尖锋，混沌的特征是谱中出现“噪声背景”和宽锋。同时我们确定混沌区后，需要进一步对吸引子进行刻画。功率谱分析仍然有用，但更重要的是计算李雅普诺夫指数。对初始条件的敏感依赖性是确定性系统混沌的关键特性。这意味着在相空间中相互靠近的两条轨线，随着时间的推移，它们将指数性的运动开。就拿差分方程或离散映射所描述的系统Xn+1=f(Xn)来讲，若初始值X01和X02相差d0=X01-X02，那么n次迭代后其中称为李雅普诺夫指数，即代表相邻点之间距离的平均辐射率。李雅普诺夫指数的大小刻划了吸引子的动力学，反映了系统产生或消除不确定因素的速率。初始不确定性经多少时间将覆盖整个吸引子由最大的指数决定，而对

35、吸引子摄动的渐近消失则由最小的指数控制。在研究混沌系统中确定奇异吸引子的各种维数、确定混沌系统的所谓柯尔莫哥洛夫熵、对周期驱动系统较易实现的分频采样法以及前面介绍的取庞加莱截面的方法等等。理论分析中则大量使用泛函分析，借用相变理论中的重正化群方法。比如费根鲍姆推出的普适常数及标度因子就是解一个反映自相似结构的泛函方程形式的重正化群方程。另外，分析系统的梅尔尼科夫(Melnikov)函数也是确定混沌的一个重要方法。通向混沌的道路有三条：倍周期分岔，阵发混沌和准周期进入混沌。与之对应的是非线性方程中三种不同类型的分岔倍周期分岔、切分岔和霍普夫分岔。2.2 二阶混沌系统的实现本次毕业设计混沌系统的实

36、现是基于二阶非线性电路，即非线性连续系统的混沌控制方法，自控制反馈连续控制方法4。非线性混沌系统的输出信号与输入信号的自反馈耦合，或者从系统外部强迫注入某一周期信号，或者直接将系统自身的输出信号取出一部分经过一定的时间延迟后再反馈到原混沌系统中去。作为控制信号，通过调节控制因子及控制信号的大小实现稳定控制。一个非线性电路产生周期、拟周期或混沌振荡，必须满足一定的电路参数条件。同一个非线性电路不同的参数，其解也不会一样。当非线性电路的参数发生变化，引起电路解的性质发生质的变化，例如由平衡点解变为周期振荡解,这种解的质的变化就称为分岔5。而这是由于非线性电路中的分岔现象与非线性电路中产生拟周期解与

37、混沌有密切的关系。通过Duffing方程5可以实现共振现象、调和振动、次调和振动、拟周期振动、概周期振动、奇异吸引子和混沌现象(或随机过程)的简单数学模型。因此，在研究非线性振动理论中Duffing方程具有重要的意义。它的标准形式为: 其中，为阻尼系数。g(x)是含有三次方项的非线性函数，f(x,t)为一个周期函数。Duffing方程通常作如下分类：（1）假设g(x)满足超线性条件:,则称Duffing方程是超线性的;（2）假设g(x)满足次线性条件:,则称Duffing方程是次线性的;（3）假设g(x)满足半线性条件:,则称Duffing方程式半线性的。Duffing方程的形式虽然简单但解法

38、很复杂，对其有个近似解法是小参数法，即二阶常微分方程：式中，是一个无量纲的正数，其值很小；是=0时无阻尼线性电路系统的谐振频率，是其变量的非线性函数，可以显示时间t，这说明方程式非自治的，此后假定与t有关的项是，正弦波形电源引入，电路会有周期稳态解，其周期和电源周期相同，这就说明当电路系统的非线性程度不严重即小，周期解会接近正弦变化。当接近于电源频率时，电路接近谐振。即便和电源频率两者相差很大，也可能出现某高次谐波谐振的情况。这样一来就忽略电感器的非线性特性，成为线性电路。而Duffing方程系统作为一个典型的非线性振动系统，尽管是从简单物理模型中得出来的非线性振动模型。在实际的工程中，在数学

39、理论基础上通过转化为该方程来解决许多非线性振动问题，如船的横摇运动、化学键的破坏等，一些张力扰动模型和动力学方程也与Duffing系统基本相似，因此Duffing系统被广泛地应用到尖锐碰摩转子的故障检测、微弱周期信号检测、电力系统周期振荡分析、周期电路系统的模拟与控制等实际工程中。关于Duffing系统还有许多问题尚未彻底研究清楚，如Duffing方程的分数谐波振动、超谐波振动、组合振动等等，而且研究结果中规律性的成果可以推广到其他类似系统。因此从某种角度来说，对非线性Duffing系统的研究是研究许多复杂动力学系统的基础。典型的Duffing方程的具体形式为：其中和为状态变量，rcos(t)

40、为周期驱动，M为阻尼系数，-x(t)+x3 (t)为非线性激励。一般为了便于对系统的分析，同时达到降阶的目的，将系统描述如下：在时域中，用或关系曲线描述状态变量x(t)或x(t)随时间t的变化规律称为时序图。以为横轴和为纵轴所构成的平面称为相平面。在相平面上做关系曲线，表征系统状态的相点随时间变化的轨迹称为相迹图。为了实现二阶混沌系统的仿真，采用铁磁谐振电路，铁芯线圈是典型的非线性电感器、含有铁芯线圈、电容器和交流电源的电路。电源电动势是时间t的函数，因此这是非自治电路。本文主要讨论Duffing方程形式：用一+R种并联LC铁磁混沌振荡电路实现：图2.1 并联LC铁磁混沌振荡电路设图中电阻R和

41、电容C是正常数，非线性电感器的磁链和电流是非线性关系，电源的电动势是正弦电动势是正弦电动势，则电路方程为其中是磁链，。对方程进行归一化处理，令，则。电路图的状态方程可化为上式的等号右方仍是频率为的正弦量。假定非线性电感器的的关系是简单的函数关系来简化计算，设曲线关于原点对称，做变量代换，令，将状态方程改写成 该式被称为Duffing方程。 混沌振动是一种由确定性系统产生对于初始条件极为敏感而具有随机性和长期预测不可能性的往复非周期运动。考虑一个由非线性磁链和线性电感组成的电流和磁链系统在简谐激励作用下的受迫振动。电感电流与磁通链数的关系满足：则系统的状态方程为 设上述系统的参数为：取两组差别很

42、小的初值： 图2.2 两个初值下的时间历程从图中可以看到，随着时间的变化，开始很接近的两个信号越来越分开，这就是系统对初值的敏感性，这也是混沌振动的一个基本特征，这也使得混沌振动具有长期不可预测性仅有初值敏感性还不能称为混沌振动，混沌振动还必须是往复的非周期性运动，这是非线性系统的又一个特征。混沌振动的往复非周期性可以用相平面图的几何方法表示出来。图2.3 相平面图当周期运动的周期很长时，仅根据相平面图难以区分周期振动和混沌振动，而Poincare映射的几何方法能更好地刻划混沌振动的往复非周期性。Poincare映射就是按系统激励的周期采样，将相轨线离散化为相点，用较少的数据获得较多信息。不同

43、的Poincare映射图形对应着不同的运动形态。当m=1.0,ek=0.5,k=1.0,ef=7.5,w=1.0,n=0具有分形结构的Poincare映射图 图2.4 Poincare映射图第三章 二阶非线性电路仿真实现MATLAB 产品家族是美国 MathWorks公司开发的用于概念设计，算法开发，建模仿真，实时实现的理想的集成环境。由于其完整的专业体系和先进的设计开发思路，使得 MATLAB 在多种领域都有广阔的应用空间，特别是在 MATLAB 的主要应用方向科学计算、建模仿真以及信息工程系统的设计开发上已经成为行业内的首选设计工具，MATLAB用户广泛的分布在航空航天，金融财务，机械化工

44、，电信，教育等各个行业。在MATLAB产品家族中，MATLAB工具箱是整个体系的基座，它是一个语言编程型（M语言）开发平台，提供了体系中其他工具所需要的集成环境（比如M语言的解释器）。同时由于MATLAB对矩阵和线性代数的支持使得工具箱本身也具有强大的数学计算能力。 MATLAB产品体系的演化历程中最重要的一个体系变更是引入了Simulink，用来对动态系统建模仿真。其框图化的设计方式和良好的交互性，对工程人员本身计算机操作与编程的熟练程度的要求降到了最低，工程人员可以把更多的精力放到理论和技术的创新上去。在MATLAB/Simulink基本环境之上，MathWorks公司为用户提供了丰富的扩

45、展资源，这就是大量的Toolbox和Blockset。从1985年推出第一个版本以后的近二十年发展过程中，MATLAB已经从单纯的Fortran数学函数库演变为多学科，多领域的函数包，模块库的提供者。用户在这样的平台上进行系统设计开发就相当于已经站在了巨人的肩膀上，众多行业中的专家、精英 们的智慧结晶可以信手拈来。同时，MATLAB开放的体系结构允许用户在平台上进行自由扩展，目前在全世界范围内已经有大量的商业的或者免费的MATLAB二次开发产品发布。用户购买一套MATLAB，获得的是世界范围的专家支持。而对于用户自己开发的算法包，MATLAB也提供了包括Compiler应用发布和Web网络发布

46、在内的众多方式的发布途径，使得用户一方面能够充分地利用MATLAB的算法资源形成技术成果，同时又可以有效的保护自己的知识产权。在这样一个产品体系中，我们可以看到，由于MATLAB及其丰富的Toolbox资源的支持，使得用户可以方便的进行具有开创性的建模与算法开发工作，并通过MATLAB强大的图形和可视化能力反映算法的性能和指标。所得到的算法则可以在Simulink环境中以模块化的方式实现，通过全系统建模，进行全系统的动态仿真以得到算法在系统中的动态验证。 但是这样一个开发流程总是欠缺和工程实现的有效连接，系统级的设计产物无法和硬件产品直接挂钩。工程师无法直接应用 MATLAB/Simulink

47、 的宝贵资源。为了改善设计流程中的这一缺陷， MATLAB 产品体系中加入了连接工程实现的桥梁 实时代码生成工具 Real-Time Workshop （ RTW ）。 RTW 使用户可以直接将 Simulink 框图模型转化为实时标准 C 代码，进而为快速原型系统、半物理仿真系统或者产品提供设计输入。 MATLAB是国际上最流行的科学与工程计算的软件，MATLAB是以复数矩阵作为基本编程单元的一种高级程序设计语言，它可以用来各种矩阵的运算与操作，并且有较强的绘图能力。MATLAB是一个高度集成的软件系统，它集科学与工程计算、计算机仿真、图形可视化、图像处理和多媒体处理于一身，提供了实用的Wi

48、ndows图形界面设计方法，使用户能设计出友好的图形界面。说到MATLAB，Simulink是该语言的重要组成部分，基于框图的可视化建模与仿真模拟方法，为复杂的工程系统建模与仿真提供了有利的工具和思路。Simulink和外界硬件环境的直接接口搭建实现了纯数字仿真与半实物仿真及实时控制，MATLAB/Simulink为工程在信息科技上的应用提供了一个新的平台6。现代科学与工程的进程都离不开数学。科学家往往关心的是数学问题的解析解和解的存在性证明，而在此不同的是工程技术人员关心的是如何求得解。在数学函数问题上有的没有解析解，有的求取解析解很麻烦，有的只有数值解。MATLAB为数学问题的解法提供了极

49、大的方便，因而被各个领域应用，MATLAB为解决数学问题提供了一组配套齐全、结构严谨的命令，在力学领域，常用有限元法求解偏微分方程，在航空、航天与自动控制领域，常用数值线性代数与常微分方程的数值解法解决问题，在工程与非工程系统的计算机仿真中，常用到各种差分方程、常微分方程的数值解法等等，MATLAB语言可以求解几乎所有的数值问题。在 MATLAB 上面执行一些数列或矩阵的运算非常方便，而它的方程式结构有点类似 BASIC或 C 程式，写起来十分平易近人，而函数调用很简单，再加以目前各种视窗上的版本都已推出，编辑程式，执行、观看结果和列印，都可轻松的透过视窗的切换及下拉式功能来完成，更值得特别一

50、提是MATLAB 的工具箱 (TOOLBOX) 有 SIGNAL 及 IMAGE 两项法宝，对于我们学习、数位信号处理的过理中，能够提供完备的的辅助。唯一较遗憾的是目前 MATLAB 上面建立的程式仍然较难直接与一些界面(如影像处理卡)相互沟通，因此如果我们想利用它来从事实际上的影像处理，必须先在其他工作环境下，将影像抓取进来，存成图档后再进入MATLAB 中，将图档调用出来作深入的分析。另一点美中不足的地方是，虽然 MATLAB 在从事一些本身的内部建函式运算时速度很快，但是如果是执行我们所建立的一些外部建函式运算时速度却相当慢，因此，如果从事语音及影像分析时，大量的计算工作将使得电脑花费不

51、少处理的时间，这一点也使得它变得较为不切实际，而无法直接运用于数位信号处理。对于N阶微分方程初值问题，由于函数及其直至(N-1)阶导数在某自变量点的值已知，所以由泰勒级数展开可以算出新的函数及导数值。在MATLAB具体利用其命令来解初值问题，步骤如下：1.根据在方程的规律，定理和公式列出微分方程和相应的初始始条件；2.运用变量替换，把一个高阶方程写成一阶微分方程组，初始条件也要做相应的替换；3.根据变换后的一阶微分方程组，编写计算导数的M-file；4.使编好的ode函数文件和变换后的初值微分方程解算命令调用，运行后即可得到在指定时间区间上的数值解。3.1 Simulink仿真MATILAB中

52、的Simulink是一个动态系统建模仿真和分析的软件包，它是一种基于MATLAB的框图设计环境，支持线性系统和非线性系统，可以用连续采样时间、离散采样时间或两种混合的采样时间进行建模，它也支持多速率系统，也就是系统中的不同部分具有不同的采样速率6。Simulink中包括许多实现不同功能的模块库，选择不同的模块建模就能模拟出不同的系统。Simulink可以建立飞仙新模型，演示非线性参数和输入幅值变化的反正结果。Duffing方程的仿真模型图如图所示。图3.1 Simulink仿真模型图取初值(0，0),且当ek=0.1,ef=92.825时，x,y的仿真波形如图所示。图3.2 Simulink仿

53、真50S的x,y的波形图图3.3 信号x,y的平面图因为混沌系统对信号参数敏感，下面探究信号对参数的敏感性，取初值(0，0)：当ek=0.1,ef=88时图3.4 信号x,y的平面图当ek=0.1,ef=80时图3.5 信号x,y的平面图当ek=0.1,ef=40时图3.6 信号x,y的平面图从这些图中我们可以看到当参数变化时，系统特性发生变化，但各信号并没有强的规律性，这反映了确定系统中的不确定性的行为特征。3.2 MATLAB语句命令演示模拟Duffing方程为了用matlab求解，将这里的x，y表示为y(1)，y(2)（1）建立自定义函数duffingfunction dy=fangch

54、eng(t,y)global d; %定义全局变量，实现传递大数据的参数，实现函数调用global w;global f;f=1;w=1; dy=zeros(2,1); %设定参数矩阵 dy(1)=y(2);dy(2)=-y(1)-y(1)3-d*y(2)+f*cos(w*t);（2）编辑出各变量之间的关系,运用MATLAB中的ode常微分求解函数指令。ode是Matlab专门用于解微分方程的功能函数。该求解器有变步长（variable-step）和定步长（fixed-step）两种类型。不同类型有着不同的求解器，其中ode45求解器属于变步长的一种，采用Runge-Kutta算法；和他采用相

55、同算法的变步长求解器还有ode23。ode45表示采用四阶，五阶Runge-Kutta单步算法,截断误差为(x)3。解决的是Nonstiff(非刚性)常微分方程。function duffingglobal d; %定义全局变量，实现传递大数据的参数，实现函数调用global f;global w;x0=0.1;v0=0.1; k=1.5,1.35,1.15;for j=1:3 d=k(j); t,u=ode45(fangcheng,0:2*pi/300:200*pi,x0,v0); u(20000:300:30000,1)figuresubplot(2,2,1)plot(t,u(:,1)ax

56、is(0 100 -1.6 1.6)title('Trajectory'); %绘制函数轨迹图xlabel('t');ylabel('x');subplot(2,2,2)plot(u(:,1),u(:,2),'MarkerSize',10)title('Phase-space Diagram'); %绘制相位空间图xlabel('x');ylabel('da/dt');subplot(2,2,3)plot(u(20000:300:30000,1),u(20000:300:30000,

57、2),'r.')title('Poincare Section'); %绘制庞加莱截面图xlabel('x');ylabel('dx/dt');subplot(2,2,4)s=plot(0,sin(x0),0,cos(x0),'LineWidth',4,'erasemode','xor');axis(-1.5 1.5 -1.5 1.5)title('频率周期图'); for i=25000:30000 set(s,'xData',0,sin(u(i,1),'yData',0,cos(u(i,1); drawnow endend运行得出的各种图形图3.7 1周期下的庞加莱截面图3.8 2周期下的庞加莱截面 图3.10 3周期下的庞加莱截面第四章 结 论对于二阶混沌系统的仿真实现，采用Duffing方程，通过