matlab 仿真求解方程根,基于MATLAB求解Rossler方程和模拟仿真.pdf

您所在位置：网站首页 > 海量文档

&nbsp>&nbsp计算机&nbsp>&nbspmatlab

基于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；

发表评论

请自觉遵守互联网相关的政策法规，严禁发布色情、暴力、反动的言论。

用户名:

验证码:

匿名?

发表评论