基于Matlab非线性电路的混沌仿真计算机课设

1、 计算机课程设计报告基于 Matlab 的非线性电路混沌实验仿真姓 名:学 院: 班 级: 指导老师:摘 要混沌是指发生在确定系统中的貌似随机的不规则运动。 而混沌 对于非线性动力学的研究有着非常重要的作用,本文结合非线性 电路的混沌的课堂教学,设计了 Matlab/Simulink仿真实验 , 研究蔡 氏电路的模拟仿真过程。本文首先通过对非线性电路(蔡氏电路 与非线性动力学进行了阐述和分析,建立非线性动力学方程,然 后利用 Matlab/Simulink软件进行仿真,研究系统波形图、单吸引 子、双吸引子、相面图以及在不同的非线性电阻的导纳下的不同 形状。达到了预期的实验效果,基于 Matla

2、b/Simulink对非线性电 路混沌的仿真对学生对非线性电路实验混沌适应实验的理解有着 较大的参考价值。关键词 非线性电路 混沌现象 Matlab/Simulink仿真目 录摘要 .3 1 混沌的概述 .4 1.1 混沌现象的概述 .41.2 混沌电路综述 .52 混沌理论基础 .5 2.1 混沌的基本定义 .5 2.2 混沌的基本特征 .52.3 混沌理论的基本概念 .73 蔡氏电路的分析与仿真 .8 3.1 蔡氏电路的分析 .93.2计算机仿真 .104 结论 .15 致谢 .16 参考书目 .161混沌的概述混沌是指发生在确定系统中的貌似随机的不规则运动, 长期以 来, 人们在认识和描

3、述运动时, 大多只局限于线性动力学描述方法, 即确定 的运动有一个完美确定的解析解。 但是自然界在相当多情况下, 非线性现象 却起着很大的作用。 1963年美国气象学家 Lorenz 在分析天气预报模型时, 首先发现空气动力学中混沌现象,该现象只能用非线性动力学来解。于是, 1975年混沌作为一个新的科学名词首先出现在科学文献中。从此,非线性 动力学迅速发展, 并成为有丰富内容的研究领域。 该学科涉及非常广泛的科 学范围从电子学到物理学,从气象学到生态学,从数学到经济学等。与我们通常研究的线性科学不同, 混沌学研究的是一种非线性科学, 而 非线性科学研究似乎总是把人们对 “ 正常 ” 事物 “

4、 正常 ” 现象的认识转向对 “ 反 常 ” 事物 “ 反常 ” 现象的探索。例如,孤波不是周期性振荡的规则传播; “ 多媒 体 ” 技术对信息贮存、压缩、传播、转换和控制过程中遇到大量的 “ 非常规 ” 现象产生所采用的 “ 非常规 ” 的新方法; 混沌打破了确定性方程由初始条件严 格确定系统未来运动的 “ 常规 ” ,出现所谓各种 “ 奇异吸引子 ” 现象等。混沌来自于非线性动力系统, 而动力系统又描述的是任意随时间发展变 化的过程, 并且这样的系统产生于生活的各个方面。 举个例子, 生态学家对 某物种的长期性态感兴趣, 给定一些观察到的或实验得到的变量 (如捕食者 个数、气候的恶劣性、食

5、物的可获性等等 ,建立数学模型来描述群体的增 减。如果用 Pn 表示 n 代后该物种极限数目的百分比,则著名的 “ 罗杰斯蒂 映射 ” :Pn+1=kP(1-Pn (其中 k 是依赖于生态条件的常数, “ n+1” 是脚标 可以用于在给定 Po , k 条件下,预报群体数的长期性态。如果将常数 k 处 理成可变的参数 k ,则当 k 值增大到一定值后, “ 罗杰斯蒂映射 ” 所构成的动 力系统就进入混沌状态。混沌(Chaos 也作混沌,指确定性系统产生的一种对初始条件具有 敏感依赖性的回复性非周期运动。浑沌与分形 (fractal和孤子 (soliton是非线性科学中最重要的三个概念。 浑沌理

6、论隶属于非线性科学, 只有非线性系 统才能产生浑沌运动。据 1991年出版的浑沌文献总目统计,已收集到 与浑沌研究有直接关系的书 269部、 论文 7157篇。 到 1996年底, 还不断有新 的浑沌研究成果发表。科学史上只有量子力学的攻坚热情可与之媲美。1.2 混沌电路综述混沌理论是现代非线性科学的一个重要分支,混沌应用是一个 全新的非线性利于,从混沌理论的研究到实际应用,混沌电路起 着重要的纽带作用,随着非线性力学的飞速发展,是的混沌理论 应用与实践成为了可能,混沌逐渐应用与电子电路,混沌通信等 诸多领域,因此,无论是在理论研究还是工程应用上,对混沌电 路域混沌系统的研究都是很有必要的。由

7、于混沌的初值的敏感性,使其在电路设计中对于案件的精度 要求比较高,因此从硬件出发研究混沌现象比较困难,而且某些 元器件仅在理论上成立,因此目前大都从理论上进行研究。 不过,经过多年的发展,人们在混沌电路的实现已经取得了可 喜的成绩,诸如蔡氏电路,可以较为方便的观察到单吸引子,双 吸引子等现象。2、混沌理论基础由于混沌系统的奇异性和复杂性至今尚未被人们彻底了解,因 此至今混沌还没有一个统一的定义,不过经过多年的发展,从不 同的理论观点出发,揭示了混沌现象的一些本质。基于混沌的“蝴蝶效益” ,洛伦兹(E.N.Lorenz 教授于 1963年大 气科学 杂志上发表了 “决定性的非周期流”一文,阐述了

8、在气候不能精确 重演与长期天气预报者无能为力之间必然存在着一种联系, 这就是非周期性 与不可预见性之间的关系。 洛伦兹在计算机上用他所建立的微分方程模拟气 候变化的时候, 偶然发现输入的初始条件的极细微的差别, 可以引起模拟结果的巨大变化。 洛伦兹打了个比喻, 即我们在文首提到的关于在南半球巴西 某地一只蝴蝶的翅膀的偶然扇动所引起的微小气流, 几星期后可能变成席卷 北半球美国得克萨斯州的一场龙卷风,这就是天气的“蝴蝶效应” 。从现象上看,混沌貌似是随机的不可预测的,但是混沌与随 机有着本质的区别,混沌运动是有确定的物理规律引起的是源于 内在特性的外在表现,因此又称为确定性混沌,下面就混沌的特

9、性加以介绍:混沌理论是近代非线性动力学中重要的组成部分,虽然混沌的定义 多繁复杂,但混沌还是有自己的一些与其他非线性系统所没有的基本特征, 具体表现为如下:(1对初始条件的敏感性经典学说认为:确定性的系统只要初始条件给定, 方程的解也就随之确 定了。 一个随时间确定性变化或具有微弱随机性的变化系统, 称为动力系统, 它的状态可由一个或几个变量数值确定。 在动力系统中, 两个几乎完全一致 的状态经过充分长时间后会变得毫无一致, 恰如从长序列中随机选取的两个 状态那样, 这种系统被称为敏感地依赖于初始条件, 这就是系统对初值的敏 感,还有混沌的敏感表现在一些控制参数的变化。1972年洛伦兹在华盛顿

10、科学进步协会上的报告上指出:“在巴西的一只 蝴蝶拍打翅膀会引发得克萨斯州的一场龙卷风” 。 这就是著名的 “蝴蝶效应” 。 这句话的意思是说任意一个微小的扰动可能会引起世界另一边天气的变化, 这种微小的扰动如同蝴蝶扇一下翅膀, 都有可能发生巨大的改变。 这一现象 的指出就是对混沌初值敏感性的最好的诠释。(2整体稳定局部不稳定稳定性是有关扰动现象的。 如果一个动力系统中发生轻微的变化, 这个 系统还会保持它的运动状态, 保持它的能力和属性。 混沌的整体稳定性指一 个微小的扰动也不会改变系统原有的性能。一个系统并不能只是绝对的稳 定, 还要有局部的稳定, 这样这个系统才能进化。 局部不稳定性表现在

11、混沌 对初值的敏感依赖性,一个微小的初值变化就会引起系统局部的不稳定。(3奇怪吸引子及其分形奇怪吸引子将混沌运动的特征初始条件的敏感性和确定性的随机直观 地反映出来。 在耗散系统当中, 当连续流在收缩体积时, 一边沿这些地方压 缩, 另一边又沿其他地方延伸。 不过连续流是固定在一个有界的区域内, 这 种伸缩和折叠过程会使运动轨道在奇怪吸引子上产生混沌运动。 可见, 奇怪吸引子是轨道不稳定和耗散系统相体积收缩两种因素的内在性质同时发生 的现象。它的几何特性由分形来刻画, 具有大尺度与小尺度之间的相似性, 具有 无穷无尽自相似的精细图案, 具有分数维数。 分形的形状是一些难以用传统 的几何学来描述

12、的极度不规则的图形; 分形存在着很小的比例精密的细节结 构; 分形的维数大于等于它的拓扑维; 分形具有自相似性, 这种自相似性可 以是严格的, 也可以是近似的或统计意义上的; 分形一般都产生于迭代过程 这些规则。 分形和混沌是同一种规律的不同表现, 这种统一的规律反映在空 间分布上表现为分形,出现在时间分布上表现为混沌。(4分岔(Bifurcation 当系统的一些控制参数发生变化时, 新的定常状态解、 周期解、 拟周期 解或者是混沌解就会分叉出来, 其中相轨迹图发生拓扑结构的突变, 分岔理 论是非线性解定性行为数学理论, 失稳是发生分岔的物理前提, 分岔后, 系 统的不同状态便会有了突变,

13、经过不断的分岔, 最终达到的状态就是混沌理 论的研究对象。(5遍历性及有界性混沌运动的轨迹经历混沌吸引子内每一个状态点的地方, 不重复, 不紊 乱。混沌的有界性最好的证明是奇怪吸引子,混沌的运动轨迹虽说有点乱, 但它始终在一个确定的区域里,有一定的规律性。(6普适性结构普适性指出无论是指数函数或是三角函数, 只要是单峰映射, 那么 函数表现出来的结构与有着某种共同的数学性质的非线性动力系统的逻辑 斯蒂方程所表现出来的结构相同, 为复杂的分岔结构。 同样都是经倍周期分 岔进入混沌状态。 测度普适性指在沿倍周期分岔进入混沌的过程中隐含着一 种深刻的规律, 它以常数的形式表现出来。 倍周期分岔序列具

14、有一个确定的 收敛速率。费根鲍姆普适常数 的数值只与系统的某种非线性性质有关,而与各 个系统的其他具体细节无关,反映出混沌演化过程中所存在的一种普适性, 说明混沌内部存在着一定的统一规律,是混沌内在规律性的另一个侧面反 映,为认识和研究混沌提供了坚实的基础。2.3 混沌理论的基本概念1. 混沌运动确定性系统中局限于有限相空间的高度不稳定的运动, 混沌电路分析及其 在保密通信中的应用研究所谓轨道高度不稳定, 是指近邻的轨道随时间的发 展会指数的分离。 由于这种不稳定性, 系统的长时间行为会显示出来某种混 乱性。2. 相空间在连续动力系统中 , 用一组一阶微分方程描述运动 , 以状态变量 (或 状

15、态向量 为坐标轴的空间构成系统的相空间。系统的一个状态用相空间的 一个点表示 , 通过该点有唯一的一条积分曲线。3. 混沌运动是确定性系统中局限于有限相空间的高度不稳定的运动。所谓轨道高 度不稳定 , 是指近邻的轨道随时间的发展会指数地分离。 由于这种不稳定性 , 系统的长时间行为会显示出某种混乱性。4. 分形和分维 :分形是 n 维空间一个点集的一种几何性质 , 该点集具有无限精细的 结构 , 在任何尺度下都有自相似部分和整体相似性质 , 具有小于所在空间 维数 n 的非整数维数。分维就是用非整数维分数维来定量地描述分形 的基本性质。5. 不动点 :又称平衡点、定态。不动点是系统状态变量所取

16、的一组值 , 对于这些 值系统不随时间变化。在连续动力学系统中 , 相空间中有一个点 0x , 若满 足当 t 时 , 轨迹 0( x t x , 则称 0x 为不动点。6. 吸引子 :指相空间的这样的一个点集 s (或一个子空间 , 对 s 邻域的几乎任 意一点 , 当 t 时所有轨迹线均趋于 s, 吸引子是稳定的不动点。7. 奇异吸引子 :又称混沌吸引子 , 指相空间中具有分数维的吸引子的集合。 该吸引集 由永不重复自身的一系列点组成 , 并且无论如何也不表现出任何周期性。 混 沌轨道就运行在其吸引子集中。8. 分叉和分叉点 :又称分岔或分支。指在某个或者某组参数发生变化时 , 长时间动力

17、学 运动的类型也发生变化。 这个参数值 (或这组参数值 称为分叉点 , 在分叉点 处参数的微小变化会产生不同性质的动力学特性 , 故系统在分叉点处是结 构不稳定的。9. 周期解 :对于系统 1( n n x f x += , 当 n 时 , 若存在 n i n x x += , 则称 该系统有周期 i 解 。不动点可以看作是周期为 1的解 , 因为它满足 1n n x x +=。10. 初值敏感性:对初始条件的敏感依赖是混沌的基本特征,也有人用它来定义混沌:混沌系统是其终极状态极端敏感地依赖于系统的初始状态的系统。 敏感依赖 性的一个严重后果就在于,使得系统的长期行为变得不可预见。3 蔡氏电路

18、的分析与仿真蔡氏电路是著名的非线性混沌电路,结构简单,但却出现 了双涡卷奇怪吸引子和极其丰富的混沌力学行为。蔡氏混沌电路(Chua s chaotic circuit是美籍华裔科学家蔡 绍棠与 1984年提出的一种典型的混沌电路,蔡氏电路是一种物理 结构和数学模型简单的混沌系统。实验电路如下图 3.1所示, 图中只有一个非线性元件 NR , 它是一 个有源非线性负阻元件, 电感 L 和电容 2c 组成一个损耗可以忽略的振 荡回路;可变电阻 V R 和电容器 1C 串联将振荡器产生的正弦信号移相输出。 较理想的非线性元件 NR 是一个三段分段的线性元件, 3.2所示的是 该电阻的伏安特性曲线,从

19、特性曲线上显示加在此非线性元件上 的电压与通过它的电流极性是相反的。由于加在此元件上的电压 增加时,通过它的电流却减小,因而将此元件称为非线性负阻元 件。 3-1 蔡氏电路图 3-1中非线性电阻 NR 式压控型非线性电阻, 它具有三段分段的 伏安特性;是它的非线性部分 ,其中 E E I E V E m m V m v g C C C C C c -+=-+-+=111101101 (21 ( 3-2蔡氏电路负阻元件的三段曲线根据图 3-1得出蔡氏电路的非线性动力学方程-=+=2212211211 (C L L C C C C C C C V dt di L i V V G dt dV C V

20、 g V V G dt dV C根据上述方程和 (1c V g 的表达式得出下面方程:则 , , , , 选择参数:带入方程式,的表达式,令 根据上述方程式和 , 0. 1, 8. 0, 5. 0m 7. 071191, , x (101322111=-=-=E m G L C C x i X V V V g L C C c-=+-=+-+-=233212112117. 07. 07. 00. 10. 1(35. 13. 68. 1x dt dx x x x dtdx x x x x dt dx此时电路有混沌解,对上述方程选取非线性项的系数为一个变量 a ,即如下 式表示:-=+-=+-+-=

21、233212112117. 07. 07. 00. 10. 1(3. 68. 1d x dt dx x x x dt dx x x a x x dt x在不同的参数 a 值对方程应用 Matlab 仿真。(1在 Matlab 编辑器中建立方程式的微分方程文件;(2通过编写程序文件解微分方程;(3由解出的方程三个变量的数值绘制波形图和吸引子;取 该 混 沌 电 路 的 典 型 状 态 进 行 仿 真 , 选 择 在 参 数 a=1.05, a=1.15,a=1.25,a=1.35下, 及电路处于周期震荡、 倍频震荡、 单涡旋混沌和双 涡旋混沌状态下进行仿真,得到如下结果。当 a=1.05时,波形

22、图,混沌吸引子及投影图形 图 3-4 a=1.05时混沌吸引子极其在各个平面的投影图形 当 a=1.15时,波形图,混沌吸引子及投影图形 时,波形图,混沌吸引子及投影图形 时,波形图,混沌吸引子及投影图形 4、结论通过分析和仿真,得出如下结论:1. 蔡氏电路的结构和蔡氏混沌电路结构完全行对称。2. 蔡氏对偶混沌电路通过计算机仿真, 在混沌振荡过程中出现了与 蔡氏电路同样的双涡卷混沌奇怪吸引子,改变其参数和初始值, 电路呈现出丰富的混沌动力学行为。次换的电路也证实了蔡氏混 沌电路所描述的非线性动态方程,正确的揭示了自然界的一种工 程物理现象的非线性激励。3. 具有负电容的三阶自制混沌电路与蔡氏电

23、路的非线性动态方程 是一致的,其电路结构是众多蔡氏混沌电路中的一种新结构。 4. 具有负电容的三阶自制混沌电路在混沌振荡过程中, 出现了与蔡氏电路同样的双涡卷混沌奇怪吸引子,改变其参数和初始值,电 路呈现出丰富的混沌动力学行为。南昌航空大学计算机课程设计报告 参考文献 1 周建兴 岂兴明.MATLAB 从入门到精通. 人民邮电出版社 2008.11 2 朱全水 龚勇清.大学物理实验. 科学出版社 2011.8 3 The Essence of Chaos混沌的本质 E.N. 洛仑兹， 气 ） 象出版社，北京 1997 4 王正林 刘 明.精通 MATLAB 7.电子工业出版社，北京 2007

24、5 于万波. 混沌的计算实验与分析.科学出版社 2011.8 6 刘秉正. 非线性动力学与混沌基础. 东北师范大学出版 社.1994.35-38 致谢 在本次设计中感谢陈常婷老师在设计中对我们的指导和帮 助！在整个创新实践过程中，陈老师以严谨求实的态度，兢兢业 业的工作作风，开拓进取的科研精神影响着我、激励着我，特此 感谢。 附录 试验程序： function yprime=bxchual(t,y %电路非线性动力学方程 yprime=-1.8*y(1+6.3*y(2-1.35*(abs(y(1-1-abs(y(1+ 1; 0.7*y(1-0.7*y(2+y(3; -7*y(2; %方程中非线性系数g=1.05,1.15,1.25,1.35 tspan=0,200; 16 %二元时间向量 南昌航空大学计算机课程设计报告 y0=-0.5;0.4;0.1; %微分方程的初始值 t,y=ode45('bxchual