数学建模社区-数学中国

标题: LORENZ吸引子 [打印本页]

---

作者: yunfenger    时间: 2010-7-28 20:46
标题: LORENZ吸引子
如下MATLAB程序。。。。MATLAB6.5
%a,b,c 为系统的三个变量，用ODE45数值解算法
- D6 T0 S& B' Z3 D, h$ G' z* Eclear
0 j& n. S. c( M# E) W4 a, _global a b c
a=10;b=8/3;c=22;
, K# Q8 ^- U, \$ X[T,X]=ode45('Lorenz',[0 100],[10;4;10]);
subplot(2,2,1)
+ h# ]% B3 t7 _# @7 `  nplot(X(:,1),X(:,2))
0 e  q5 q- \- w' {subplot(2,2,2)
6 m7 G6 b( O4 q% Aplot(X(:,1),X(:,3))
2 y% b6 K2 {, Y: r5 R& ]5 M" r, H2 N0 P  rsubplot(2,2,3)
plot(X(:,2),X(:,3))
9 X* D9 K1 E0 G' x. W8 U1 msubplot(2,2,4)
& i4 F) U7 S4 ^+ |" bplot3(X(:,1),X(:,2),X(:,3))
( F; `1 _' V2 c/ h! Lview(-20,60);
xlabel('x');ylabel('y');zlabel('z');
' S  m' x4 l+ J/ p$ [

! ^9 s% A7 \9 L# D
function dx=Lorenz(t,x)
$ N6 ]8 b9 B) z. [global a b c
dx=[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x(2);x(1)*x(2)-b*x(3)];
. I) W5 R3 G- v' p4 B1 ]8 I

1 }. g" Z6 E3 S, S5 I画出来的图形附件！根据不同的参数画出来的图
4 X$ I7 q8 z1 E6 \
( l4 H4 J' ~% G! t9 E. D

r602.fig

2010-7-28 20:45 上传

点击文件名下载附件

下载积分: 体力 -2 点

125.32 KB, 下载次数: 6, 下载积分: 体力 -2 点

r1.fig

2010-7-28 20:45 上传

点击文件名下载附件

下载积分: 体力 -2 点

108.77 KB, 下载次数: 2, 下载积分: 体力 -2 点

r8.fig

2010-7-28 20:45 上传

点击文件名下载附件

下载积分: 体力 -2 点

119.01 KB, 下载次数: 0, 下载积分: 体力 -2 点

r15.fig

2010-7-28 20:45 上传

点击文件名下载附件

下载积分: 体力 -2 点

126.88 KB, 下载次数: 0, 下载积分: 体力 -2 点

r24.fig

2010-7-28 20:45 上传

点击文件名下载附件

下载积分: 体力 -2 点

394.63 KB, 下载次数: 0, 下载积分: 体力 -2 点

r30.fig

2010-7-28 20:45 上传

点击文件名下载附件

下载积分: 体力 -2 点

37.27 KB, 下载次数: 0, 下载积分: 体力 -2 点

r40.fig

2010-7-28 20:45 上传

点击文件名下载附件

下载积分: 体力 -2 点

465.23 KB, 下载次数: 1, 下载积分: 体力 -2 点

r50.fig

2010-7-28 20:45 上传

点击文件名下载附件

下载积分: 体力 -2 点

58.53 KB, 下载次数: 0, 下载积分: 体力 -2 点

r60.fig

2010-7-28 20:45 上传

点击文件名下载附件

下载积分: 体力 -2 点

65.47 KB, 下载次数: 0, 下载积分: 体力 -2 点

r80.fig

2010-7-28 20:45 上传

点击文件名下载附件

下载积分: 体力 -2 点

663.95 KB, 下载次数: 0, 下载积分: 体力 -2 点

r146.fig

2010-7-28 20:45 上传

点击文件名下载附件

下载积分: 体力 -2 点

874.89 KB, 下载次数: 0, 下载积分: 体力 -2 点

r150.fig

2010-7-28 20:45 上传

点击文件名下载附件

下载积分: 体力 -2 点

888.51 KB, 下载次数: 0, 下载积分: 体力 -2 点

r168.fig

2010-7-28 20:45 上传

点击文件名下载附件

下载积分: 体力 -2 点

960.79 KB, 下载次数: 0, 下载积分: 体力 -2 点

r180.fig

2010-7-28 20:45 上传

点击文件名下载附件

下载积分: 体力 -2 点

988.35 KB, 下载次数: 0, 下载积分: 体力 -2 点

r200.fig

2010-7-28 20:45 上传

点击文件名下载附件

下载积分: 体力 -2 点

1.01 MB, 下载次数: 0, 下载积分: 体力 -2 点

r220.fig

2010-7-28 20:45 上传

点击文件名下载附件

下载积分: 体力 -2 点

1.05 MB, 下载次数: 0, 下载积分: 体力 -2 点

r240.fig

2010-7-28 20:45 上传

点击文件名下载附件

下载积分: 体力 -2 点

1.13 MB, 下载次数: 1, 下载积分: 体力 -2 点

r242.fig

2010-7-28 20:45 上传

点击文件名下载附件

下载积分: 体力 -2 点

58.76 KB, 下载次数: 0, 下载积分: 体力 -2 点

r244.fig

2010-7-28 20:45 上传

点击文件名下载附件

下载积分: 体力 -2 点

414.01 KB, 下载次数: 0, 下载积分: 体力 -2 点

r302.fig

2010-7-28 20:45 上传

点击文件名下载附件

下载积分: 体力 -2 点

86.22 KB, 下载次数: 0, 下载积分: 体力 -2 点

r343.fig

2010-7-28 20:45 上传

点击文件名下载附件

下载积分: 体力 -2 点

96.07 KB, 下载次数: 0, 下载积分: 体力 -2 点

r402.fig

2010-7-28 20:45 上传

点击文件名下载附件

下载积分: 体力 -2 点

100.38 KB, 下载次数: 0, 下载积分: 体力 -2 点

r502.fig

2010-7-28 20:45 上传

点击文件名下载附件

下载积分: 体力 -2 点

111.45 KB, 下载次数: 0, 下载积分: 体力 -2 点

---

作者: zhao6884269    时间: 2010-8-26 11:33
试试运气啦~~~~~~~~~~~

---

作者: lilijie_shangha    时间: 2010-8-26 11:34
来报道!!!!!!!!!!!

---

作者: wy371tt1    时间: 2010-8-26 11:35
鉴定完毕！  

---

作者: fengzhisheng106    时间: 2010-8-26 11:36
强人，佩服死了。呵呵，不错啊

---

作者: 想做就做    时间: 2010-8-26 11:38
强人，佩服死了。呵呵，不错啊

---

作者: 547896321    时间: 2010-8-26 11:46
我来了~~~~~~~~~ 闪人~~~~~~~~~~~~~~~~  

---

作者: mybaby    时间: 2010-8-26 13:52
哦～～

---

作者: surfw    时间: 2010-8-26 18:59
试试运气啦~~~~~~~~~~~

---

作者: xiaoai    时间: 2010-8-26 19:17
强人，佩服死了。呵呵，不错啊

---

作者: 阿里wudi    时间: 2010-8-27 00:00
楼主，你写得实在是太好了。我惟一能做的，就只有把这个帖子顶上去这件事了

---

作者: sigh    时间: 2010-8-27 12:00
强烈支持。楼主万岁

---

作者: ycy74_2005    时间: 2010-8-27 15:00
我来了~~~~~~~~~ 闪人~~~~~~~~~~~~~~~~  

---

作者: guohuawin    时间: 2010-8-27 20:00
声明一下：本人看贴和回贴的规则，好贴必看，精华贴必回。

---

作者: weixinmaths    时间: 2011-6-23 09:01

---

作者: 枫露之茗    时间: 2012-7-16 11:41
鉴定完毕！

---

作者: 灼￡日    时间: 2013-10-28 19:35
程序报错地方多

---

作者: hngy    时间: 2013-11-6 11:26
Lorentz混沌这么费体力

---

作者: 天马绝尘    时间: 2014-2-12 10:35
试试运气   

---

欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)

Powered by Discuz! X2.5