基于MATLAB求解Rossler方程和模拟仿真.pdfVIP

2009年第 1期 牡丹汪师范学院学报 (自然科学版) No．1，2009 (总第 66期) JournalofMudanjiangNormalUniversity TotalNO66 基于MATLAB求解R0ssler方程和模拟仿真 陈永胜 (吉林师范大学数学学院，吉林 四平 136000) 摘 要 ：给 出用 MATLAB求解 ROssler方程 的一种 方法，用立体 图形动 态显示 出吸 引子 ，应用 MATIAB 仿 真工具进行混沌模拟 ，并对仿真结果作 了说明和讨论． 关键词 ：R6ssler方程 ；吸引子 ；MATIAB；混沌 ；仿 真 [中图分类~,]0175 [文献标识码]A [文章编号]1003—6180(2009)01--0009--02 R6ssler方程是 由ROssler在 1976年建立的 plot(t，y(：，3))％显示 y(3)即z与时间的 一 个三维系统，一个常微分方程 (ODE)系统用解 关系图 析法精确求解是不可能的，只能用计算机数值计 subplot(2，2，4) 算 ，主要有欧拉法 、亚当法和龙格一库塔法等．但用 plot(t，y) 显示三分量的关系图，即吸引子 这些方法都需要编写一段很长的程序 ，若步长取 figure； 得不恰 当，带来 的误差也很大．利用数学工具 plot3(y(：，1)，y(：，2)，y(：，3)) MATIAB求解微分方程却是一种较为方便的方 (3)求解 结果．在 MATLAB 窗 VI中执行 法，不仅求解速度快，且具有通用性． Rossler1．m 文件 ，结果见图 1，图 2． 1 R6ssler方程求解 R6ssler方程如下 ： dx — d—t一 一 y-- ’ 鲁一+o．2， (1) 图 1 图 2 dz ． d 一0．2--5．7+ z． t (4)动态显示吸引子的绘制过程．用 edit命 这是一个 自治的方程组 ，求解过程： 令建立一个命令文件 Rossler2．m ，内容为 (1)建立 自定义函数．用 edit命令建立 自定 [t，y]一ode45(Rossler，[0，lOO]，[o，0，o])； 义函数名为Rossler．m ，内容为 clf functiondx=Rossler(t．x) axis([一3O，3O，一3O，3O，一20，20]) 设定 dx=l—X(2)一X(3)；x(1)+0．2*x(2)； 合适的坐标范围 0．2+ (x(1)一5．7)*x(3)．-1： view([3o，60])；％设定较好的观察角度 (2)用ode45命令求僻．用 edit命令建立一个 holdon 保持坐标轴不变 命令文件 Rossler1．m ，内容为 comet3(y(：，1)，y(：，2)，Y(：，3))；9／6显 clf 示吸引子绘制过程 xO=Eo；0；