蔡氏电路matlab仿真报告

蔡氏电路仿真分析学院：电气工程学院班级：硕6036姓名：张东海学号：31163120531目录1.基本分析..............................................................................................22.MATLAB仿真........................................................................................52蔡氏电路蔡氏电路是著名的非线性混沌电路，结构简单，但却出现双涡卷奇怪吸引子和及其丰富的混沌动力学行为。1.基本分析蔡氏电路是一个典型的混沌电路，最早由著名华裔科学家、美国加州大学蔡少堂教授设计。他证明了在满足以下条件时能够产生混沌现象。(1)非线性元件不少于1个；(2)线性有效电阻不少于1个；(3)储能元件不少于3个。根据以上条件，在图1.1中给出蔡氏电路方框图。图中R为线性有效电阻，L、C1、C2为储能元件，RN为非线性元件。图2.2给出非线性电阻伏安特性曲线。图1.1蔡氏电路方框图图1.2非线性电阻伏安特性曲线对于图2.1提出的蔡氏电路，其状态方程推导如下12112122121()()1()(1)CCCCCCCLLCduCuugudtRduCuuidtRdiLudt其中函数1()Cgu是分段线性函数，其形式为：311111()()()2CbCabCCguGuGGuEuE作变量代换：1222221,,,,1CCLuuixyzEEEGCCtGCCLGGR式（1）可以写为如下形式()(2)dxyxfxddyxyzddzyd式（2）即是蔡氏电路的标准方程形式。其中()fx可表示为如下形式1010101(),1(),1(),1mxmmxfxmxxmxmmx其中01,abmGEmGE蔡氏电路的三个状态方程式在状态空间的三个子空间为101={(,,)|1}={(,,)|1}={(,,)|1}DxyzxDxyzxDxyzx在状态空间的三个子空间内分别具有唯一平衡点如下1011(,0,),(0,0,0),(,0,).PkkDQDPkkD其中，1011mmkm在P、1P和Q处的雅可比矩阵分别为：41(1)011100PPmJJ，0(1)011100QmJ取10，15，01.2m，10.6m，则在P、1P处的特征值为一个实数值和一对共轭复数值。其中的一个实数值为-5.3938P而一对共轭复数值为j0.1969j3.3294PP可见0P，而0P。在Q处的特征值也为一个实数值和一对共轭复数值。其中的一个实数值为3.1201P而一对共轭复数值为j-1.0601j2.9140PP可见0P，而0P。因此，所有的平衡点P、1P和Q均为鞍焦点。2.MATLAB仿真在对蔡氏电路进行建模过程中，要用到加法器、积分器、增益模块、常数模块等单元，这些可以方便地从Simulink标准库模块中直接提取。除此之外，还需用到显示x和y信号相图的平面图形显示器（XYGraph）模块和示波器模块；最后，把模型顺序连接起来，得到用基本模块元素建立的数学仿真模型，如下图所示。5取0.1,0.1,0.1xyz为初值，作为系统的初始值,在t=[0,300]的时间范围内求解系统的运动轨迹，其投影相图及时域波形如下:图2.1x-y-z立体相图6图2.2x时域波形图2.3y时域波形7图2.4z时域波形图图2.7x-y平面相图8图2.6y-z平面相图图2.5x-z平面相图从仿真结果图可以，蔡氏电路的正规化状态方程描述了一个连续时间系统，这个系统在所给参数和初值的条件下可以产生双涡卷吸引子的混沌现象。改变参数和初值，还可以产生其它很多有趣的混沌现象。利用系统平衡点处的线性化矩阵，可以定性分析系统的动力学行为，以便寻找能使系统产生混沌的参数。