蔡氏电路MATLAB混沌仿真.docx

蔡氏电路的 Matlab 混沌仿真研究班级姓名学号摘要本文首先介绍非线性系统中的混沌现象，并从理论分析与仿真计算两个方面细致研究了非线 性电路中典型混沌电路，即蔡氏电路反映出的非线性性质通过改变蔡氏电路中元件的参数，进 而产生多种类型混沌现象最后利用软件对蔡氏电路的非线性微分方程组进行编程仿真，实现了 双涡旋和单涡旋状态下的同步，并准确地观察到混沌吸引子的行为特征关键词：混沌；蔡氏电路；MATLAB仿真AbstractThis paper introduces the chaos phenomenon in nonlinear circuits. Chu'a s circuit was a typical chaos circuit, thus theoretical analysis and simulation was made to research it. Many kinds of chaos phenomenon on would generate as long as one component parameter was altered in Chua's circuit．On the platform of Matlab, mathematical model of Chua's circuit was programmed and simulated to acquire the synchronizationof dual and single cochlear volume. Meanwhile, behavioral characteristics of chaos attractor were observed.Key words： chaos phenomenon； Chua's circuit； Simulation1、引言混沌理论的基本思想起源于20世纪初，完善于 20世纪60年代后，发展壮大于 20世纪80 年代，被认为是继相对论、量子力学之后，人类认识世界和改造世界的最富有创造性的科学领域 第三次大革命。

混沌理论揭示了有序与无序的统一，确定性与随机性的统一，简单性与复杂性的 统一，稳定性与不稳定性的统一，完全性与不完全性的统一，自相似性与非相似性的统一，并成 为正确的宇宙观和自然哲学的里程碑今天，混沌理论又与计算机科学等相结合，使人们对一些 久悬未解的难题的研究取得了突破性进展，在探索、描述及研究客观世界的复杂性方面发挥了巨 大作用混沌现象是非线性系统在特定条件下产生的特殊行为，作为一种普遍存在的非线性现象，混 沌的发现对科学的发展具有深远的影响混沌行为是确定性因素导致的类似随机运动的行为，即: 一个可由确定性方程描述的非线性系统，其长期行为表现为明显的随机性和不可预测性，我们就 认为该系统存在混沌现象混沌具有三个特点 1、随机性，即混沌具有类似随机变量的杂乱表现； 2、遍历性，即能够不重复地历经系统的所有状态点；3、规律性，即混沌是由确定性的迭代式产 生的混沌还有一个很重要的性质：系统行为对初始条件非常敏感物理、化学、生物学，以及社会科学等等各个学科领域中都有混沌现象近年来许多学者通 过非线性电路对混沌行为进行了广泛地研究，其中最典型的是由美国 Berkeley 大学的 Leon. O.Chua 提出的蔡氏电路，它是能产生混沌行为的最小、最简单的三阶自治电路。

自从 1983 年电路发现以来，它一直是人们研究混沌现象的重要模型，在许多文献中都以蔡 氏电路为基础研究混沌现象它的优点在于电路非常简单，但其非线性动力学行为却极其丰富 本文对蔡氏电路的混沌特性进行了理论分析，并通过仿真观察三阶自治动力系统的混沌双涡卷吸 引子和稳定周期轨道2、蔡氏电路结构模型1983年，美籍华裔科学家蔡少棠教授首次提出了著名的蔡氏电路(Chua' s circuit)它是 历史上第一例用电子电路来证实混沌现象的电路，也是迄今为止在非线性电路中产生复杂动力学 行为的最为有效和较为简单的电路通过改变蔡氏电路的拓扑结构或电路参数，可以产生倍周期 分叉、单涡卷、双涡卷吸引子、多涡卷吸引子等十分丰富的混沌现象因此可以说蔡氏电路开启 了混沌电子学的大门，人们已围绕它开展了混沌机理的探索、混沌在保密通信中的应用研究，并 取得了一系列丰硕的成果自治动力系统产生混沌现象需要以下条件：系统至少有三个状态变量，并且存在一定的非线 性环节蔡氏电路使用三个储能元件和一个分段线性电阻，电路如图 1 所示可以把电路分为线性部分和非线性部分其中线性部分包括：电阻R、电感L和两个电容C 与 C2；非线性部分只有 一个分段线性电阻R，其伏安特性如图2所示。

mR图 1 蔡氏电路图 2 非线性电阻的伏安特性根据图 1 可以列写所示电路的三阶微分方程组为：壮；1 = R(u 2-u1) 一 / (u1)C du2 =丄(u — u ) — i2 d t R 1 2 L(1)diL—i = —udt 2其中，i 二 f (u) = f (u),r 1 r它是一个三段线性的分段线性函数:ir = f ("1) = 1mu + E(m — m )0 1 1 0mu11mu + E(m — m )0 1 1 0u = u > Eilu1 < Eu = u < E1也可以写成： f (u ) =1m u + (m0 1 2 1—m ) + E — |u + E0 1L-E如果定义：T = t,a = m R ,C R 12则原微分方程组(1) 变为b=m0RC— cC1(2)dx在=a (y - x - /(x))(3)dy = x — y + zd t芽=-0ybx + a 一 b其中：I=f(u1)=1axbx 一 a + bx > 1|x| < 1x = 一1蔡氏电路中的非线性电阻又称为蔡氏二极管，可采用多种方式实现一种较简单的实现，它 相当于两个非线性电阻R和R的并联。

图3给出R和R的电路及其伏安特性N1 N2 N1 N2适当选取电阻参数值使E2远大E1,也远大于蔡氏电路工作时uc1的变化范围，则在电路的 工作范围内,RN2是一个线性负电阻,RN1和RN2并联后可实现图l中非线性电阻RN的伏安特性打£;电路及其伏安轉桂图 3 两个非线性电阻的及其伏安特性3、蔡氏电路Mat lab仿真通过电路参数的调整，我们可以从蔡氏电路中观察到丰富的非线性动态特性，以下我们详细地给出各种特性的Matlab仿真结果在图1所示的电路中，选择电路参数为：C二10nF, C= 100nF, L = 18mH，改变电路各参数12均可改变电路的非线性特性以改变电阻参数为例，通过Matlab仿真电阻由小到大改变，初初值取x二0.01, y = 0，z = 0；步长为0.01，共计算20000个点当电阻为1.86 kQ时系统开始出现双涡卷吸引子如图4所示a)吸引子在相平面iL- u(1)上的投影(b)吸引子在相平面iL- u(2)上的投影0 1 2 3 5(c)吸引子在相平面u(1)- u(2)上的投影(d)吸引子在相平面t- u(1)上的投影图4电阻R二1.86 kQ时的吸引子给出几种典型的吸引子在状态空间的投影如图4所示。

当电阻为1934Q时系统开始出现螺旋吸引子a) R=1932Q时的双涡卷吸引子(b) R=1933Q时的双涡卷吸引子(c) R=1934Q时的螺旋吸引子(d) R=1965Q时的螺旋吸引子图5吸引子在相平面iL- u(1)上的投影4、总结由上述分析及仿真结果可知，虽然蔡氏电路非常简单，但其非线性动力学行为却极其丰富，当选择适当的电路参数时，其动态特性出现混沌现象，此时，奇怪吸引子具有双涡卷结构，出现 连续的功率谱，分形特征可用分数维来描述。