非线性报告-姚威doc.docVIP

非线性电路理论及应用课程报告 蔡氏混沌电路的MATLAB仿真及分析 姓名：姚威 学号：3115312022 专业：电气工程 2015年11月一 引言 混沌是非线性系统所独有且广泛存在的一种非周期运动形式, 其覆盖面涉及到自然科学和社会科学的几乎每一个分支。1972年12月29日，美国麻省理工学院教授、混沌学开创人之一E.N.洛伦兹在美国科学发展学会第139次会议上发表了题为《蝴蝶效应》的论文，提出一个貌似荒谬的论断：在巴西一只蝴蝶翅膀的拍打能在美国得克萨斯州产生一个龙卷风，并由此提出了天气的不可准确预报性。为什么会出现这种情况呢？这是混沌在作怪！ “混沌”译自英语中“chaos”一词，原意是混乱、无序，在现代非线性理论中，混沌则是泛指在确定体系中出现的貌似无则的、类随机的运动。’s Circuit），它是混沌发展史上的一个重要里程碑。蔡氏电路时目前众多混沌电路中最具代表性的一种，它的主要特点是能够产生双涡卷混沌吸引子，其混沌动力学行为已分别被数学分析、数值模拟和硬件实验所证实，并得到了严格的数学证明，成为理论和实践研究混沌的一个范例。 二 MATLAB中的龙格—库塔（Runge-Kutta）实现 MATLAB（Matrix Laboratory）是MathWorks公司开发的，目前国际上最流行应用最广的科学与工程计算机软件之一。MATLAB 软件以矩阵运算为基础，把计算，可视化，程序设计等有机的融合在一起，具有出色的数值计算能力和强大的图形处理功能。 基于Runge—Kutta法，MATLAB提供了求解微分方程数值解的函数，一般调用格式是： 其中fname 是定义的函数文件名，该函数文件必须返回一个列向量。Tspan形式是[t0,tf]，表示求解区间，y0是初始状态向量。这两个函数分别采用“二阶，三阶Runge—Kutta法”和“四阶，五阶Runge—Kutta法”，并采用自适应的求解方法，即当解的变化较慢时采用较大的步长，从而使计算速度很快，当解的变化较快时步长会自动变小长，从而使计算精度很高。在MATLAB中，一般选取四阶的龙格库塔方法。 三 单涡旋的变形蔡氏电路混沌系统的仿真与分析 1983年美国电学专家蔡少棠 (L. O. Chua)首次提出了著名的蔡氏电路 (Chua ‘s circuit) ,它是混沌发展史上的一个重要里程碑. 它使人们从被动的研究混沌现象向主动的设计和控制混沌迈出了关键的一步.蔡氏电路是目前众多混沌电路中最具典型代表性的一种。对于单涡旋的变形蔡氏电路的微分方程组为： 当＝5.5，＝7.4， ＝0.25， ＝0.1初始值（x,y,z）=(10,-50,1);时出现混沌状态，如图所示： 3-1 单涡旋蔡氏电路时域响应 3-2 单涡旋蔡氏电路相图 其初值的敏感性，分别仿真初始值（x,y,z）=(10,-50,1)和（x,y,z）=(10,-50,1)在经过短时间t=10和较长时间t=100的蔡氏电路系统图，如下图所示， 3-3 单涡旋蔡氏电路当初始值不同时的相图 如图3-3可以看出，随着初始值的改变，混沌震荡逐渐减弱。 四 多涡旋的蔡氏电路混沌系统的仿真与分析 对于多涡旋的情况，例如可以用多项式 产生三涡卷蔡氏混沌吸引子，其无量纲归一化状态方程可以写为： 其中，参数＝12.8，＝19.1，实验中固定b＝－1.1，c＝0.45。 初始值（x,y,z）=（0.025,-0.022,0.8）; 改变参数a，观察该电路模型的进入混沌状态的过程 4-1 a=0.3时蔡氏电路的时域响应 4-2 a=0.3时蔡氏电路的相图 a=0.3时，为单吸引子混沌震荡。 4-3 a=0.76时蔡氏电路的时域响应 4-4 a=0.76时蔡氏电路的相图 a=0.76时，为三吸引子混沌震荡。 4-5 a=0.79时蔡氏电路的时域响应 4-6 a=0.79时蔡氏电路的相图 a=0.79时，为双吸引子混沌震荡。 4-7 a=0. 9时蔡氏电路的时域响应 4-8 a=0.9时蔡氏电路的相图 a=0.9时，为单吸引子混沌震荡。 4-9 a=2.1时蔡氏电路的时域响应 4-10 a=2.1时蔡氏电路的相图 a=2.1时，为单吸引子混沌震荡. 由以上几张图可以看出，随着蔡氏电路中，系统参数a的改变，混沌震荡改变很明显，可能出现单吸引子混沌震荡、双吸引子混沌震荡和三吸引子混沌震荡。 五 具有负电容的三阶自治蔡氏混沌电路的分析 电路如图5-1所示，此电路是一个三阶自治电路，其中含有一个电压控制的非线性电阻元