您所在位置:网站首页 > 海量文档
 > 计算机 > matlab
基于MATLAB求解Rossler方程和模拟仿真.pdf2页
本文档一共被下载:次,您可全文免费在线阅读后下载本文档。
下载提示
1.本站不保证该用户上传的文档完整性,不预览、不比对内容而直接下载产生的反悔问题本站不予受理。
2.该文档所得收入(下载+内容+预览三)归上传者、原创者。
3.登录后可充值,立即自动返金币,充值渠道很便利
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;
发表评论
请自觉遵守互联网相关的政策法规,严禁发布色情、暴力、反动的言论。
用户名:
验证码:
匿名?
发表评论