蔡氏电路MATLAB混沌仿真

1、蔡氏电路的 Matlab 混沌仿真研究班级：姓名：学号：摘要本文首先介绍非线性系统中的混沌现象， 并从理论分析与仿真计算两个方面细致研究了非线 性电路中典型混沌电路， 即蔡氏电路反映出的非线性性质。 通过改变蔡氏电路中元件的参数， 进 而产生多种类型混沌现象。 最后利用软件对蔡氏电路的非线性微分方程组进行编程仿真， 实现了 双涡旋和单涡旋状态下的同步，并准确地观察到混沌吸引子的行为特征。关键词 ：混沌；蔡氏电路； MATLAB 仿真AbstractThis paper introduces the chaos phenomenon in nonlinear circuits. Chua'

2、; scircuit 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 ihnuCa' s circuit On the platform of Matlab, mathematical model of Chua' scircuit was prog

3、rammed and simulated to acquire the synchronization of dual and single cochlear volume. Meanwhile, behavioral characteristics of chaos attractor were observed.Key words： chaos phenomenon；Chua' s circu；it Simulation1、引言混沌理论的基本思想起源于 20世纪初，完善于 20世纪 60年代后，发展壮大于 20世纪 80 年代，被认为是继相对论、 量子力学之后， 人类认识世界和改造

4、世界的最富有创造性的科学领域 第三次大革命。 混沌理论揭示了有序与无序的统一， 确定性与随机性的统一， 简单性与复杂性的 统一，稳定性与不稳定性的统一，完全性与不完全性的统一，自相似性与非相似性的统一，并成 为正确的宇宙观和自然哲学的里程碑。 今天，混沌理论又与计算机科学等相结合， 使人们对一些 久悬未解的难题的研究取得了突破性进展， 在探索、描述及研究客观世界的复杂性方面发挥了巨 大作用。混沌现象是非线性系统在特定条件下产生的特殊行为， 作为一种普遍存在的非线性现象， 混 沌的发现对科学的发展具有深远的影响。 混沌行为是确定性因素导致的类似随机运动的行为， 即 : 一个可由确定性方程描述的非

5、线性系统， 其长期行为表现为明显的随机性和不可预测性， 我们就 认为该系统存在混沌现象。 混沌具有三个特点 1、随机性，即混沌具有类似随机变量的杂乱表现； 2、遍历性，即能够不重复地历经系统的所有状态点； 3、规律性，即混沌是由确定性的迭代式产 生的。混沌还有一个很重要的性质：系统行为对初始条件非常敏感。物理、化学、生物学， 以及社会科学等等各个学科领域中都有混沌现象。近年来许多学者通 过非线性电路对混沌行为进行了广泛地研究，其中最典型的是由美国 Berkeley 大学的 Leon. O.Chua提出的蔡氏电路，它是能产生混沌行为的最小、最简单的三阶自治电路。自从 1983 年电路发现以来，它

6、一直是人们研究混沌现象的重要模型，在许多文献中都以蔡 氏电路为基础研究混沌现象。它的优点在于电路非常简单，但其非线性动力学行为却极其丰富。 本文对蔡氏电路的混沌特性进行了理论分析， 并通过仿真观察三阶自治动力系统的混沌双涡卷吸 引子和稳定周期轨道。2、蔡氏电路结构模型1983年，美籍华裔科学家蔡少棠教授首次提出了著名的蔡氏电路(Chua's circuit)。它是历史上第一例用电子电路来证实混沌现象的电路， 也是迄今为止在非线性电路中产生复杂动力学行为 的最为有效和较为简单的电路。通过改变蔡氏电路的拓扑结构或电路参数，可以产生倍周期分叉、 单涡卷、双涡卷吸引子、多涡卷吸引子等十分丰富的

7、混沌现象。因此可以说蔡氏电路开启了混沌 电子学的大门，人们已围绕它开展了混沌机理的探索、混沌在保密通信中的应用研究，并取得了 一系列丰硕的成果。自治动力系统产生混沌现象需要以下条件： 系统至少有三个状态变量，并且存在一定的非线 性环节。蔡氏电路使用三个储能元件和一个分段线性电阻， 电路如图1所示。可以把电路分为线 性部分和非线性部分。其中线性部分包括：电阻 R、电感L和两个电容Ci与C2；非线性部分只 有一个分段线性电阻Rm,其伏安特性如图2所示。LiLiC2+ 卄U2G 十 U1 RmkI图1蔡氏电路图2非线性电阻的伏安特性根据图1可以列写所示电路的三阶微分方程组为：C1d u1dtr(U2

8、 ui) f(ujd u2R(U1 U2)iL(1)L当dtU2其中,ir = f（U1） = f （Ur）,它是一个三段线性的分段线性函数：m°U1E(m)1 m°)uu1 Erf（5）m1u1u1EU1E(mh m°)uu1 E也可以写成:f(Ui)i丄(m12m0U1mo) U1u1I，U2E如果定义:miR ,m°R ,则原微分方程组（1）变为（y x f（x）bx其中:irf (Ui)axbx蔡氏电路中的非线性电阻又称为蔡氏二极管,u1 EC2C1，c2r2可米用多种方式实现。一种较简单的实现，它相当于两个非线性电阻 RN1和RN2的并联。图3

9、给出RN1和RN2的电路及其伏安特性。适当选取电阻参数值使 E2远大E1,也远大于蔡氏电路工作时uc1的变化范围，则在电路的工作范围内，RN2是一个线性负电阻，RN1和RN2并联后可实现图I中非线性电阻RN的伏安 特性。电路及其佻安伶性时仏电略凤算假安转牲图3两个非线性电阻的及其伏安特性3、蔡氏电路Matlab仿真通过电路参数的调整，我们可以从蔡氏电路中观察到丰富的非线性动态特性，以下我们详细地给出各种特性的Matlab仿真结果。在图1所示的电路中，选择电路参数为：Ci= 10nF, C2= 100nF, L = 18mH，改变电路各参数均可改变电路的非线性特性。以改变电阻参数为例，通过 Ma

10、tlab仿真电阻由小到大改变，初初值取x = 0.01, y = 0, z = 0;步长为0.01，共计算20000个点。当电阻为1.86 kQ时系统开始出现双 涡卷吸引子如图4所示。(a)吸引子在相平面iL-u(1)上的投影(c)吸引子在相平面u(1)-u(2)上的投影(d)吸引子在相平面t- u(1)上的投影图4电阻R = 1.86 kQ时的吸引子给出几种典型的吸引子在状态空间的投影如图4所示。当电阻为1934Q时系统开始出现螺旋吸引子，皿 g £4-202 亦 iI I10-6 -B-2 u 2! E Q 10R=1932 Q时的双涡卷吸引子(b) R=1933 Q时的双涡卷吸引子(c) R=193