数学建模社区-数学中国

标题: LORENZ吸引子 [打印本页]
作者: yunfenger    时间: 2010-7-28 20:46
标题: LORENZ吸引子
如下MATLAB程序。。。。MATLAB6.5
2 C9 h; i, y% _, ?4 k%a,b,c 为系统的三个变量，用ODE45数值解算法
, E2 R, m3 B$ n4 mclear5 y9 O0 M; b1 M1 g8 v. }) L& Y. `
global a b c
  Y1 x' d# ^3 _- l; W1 Wa=10;b=8/3;c=22;
8 j% t. [/ k8 c/ d[T,X]=ode45('Lorenz',[0 100],[10;4;10]);" D2 S6 s& S% D% ]* J: c
subplot(2,2,1)
3 U- u7 {( M7 W0 P. H6 X( ^+ V' W$ Qplot(X(:,1),X(:,2))# P$ Z0 x6 ^8 O+ g) A( Q
subplot(2,2,2)
8 W4 {" i1 d5 \  z+ R2 Q" iplot(X(:,1),X(:,3))
, D& D  Z( N5 y1 z# @' Zsubplot(2,2,3)
( ~; w' @, k2 t5 u2 X# x! ]plot(X(:,2),X(:,3))" P" _4 P; R3 q9 a9 ?) y/ h  K
subplot(2,2,4)
( s- C1 K0 b3 W; A8 F9 A2 L! eplot3(X(:,1),X(:,2),X(:,3)): ?: ~8 Z% R0 Y8 y  s% y
view(-20,60);
* M+ o0 t0 D/ y+ j/ Pxlabel('x');ylabel('y');zlabel('z');" |3 f, q" ?% h5 f! c1 }: g
! r" w6 @- r6 Q; {
! i" E. z3 C  c' T
5 }4 p9 E5 G6 i8 W: D
function dx=Lorenz(t,x)
  }- I! ~8 s5 W7 @$ hglobal a b c
- ?- d" Y) m# s" S/ Ndx=[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x(2);x(1)*x(2)-b*x(3)];$ _, ?0 Z/ _3 B4 z
) J8 z: f) \' G6 K' X

1 K) ^; M: f/ \: q画出来的图形附件！根据不同的参数画出来的图
7 }3 B, B2 T% J) q3 n& y9 y
) Q4 N5 [( Z( p. y

r602.fig

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

r1.fig

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

r8.fig

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

r15.fig

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

r24.fig

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

r30.fig

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

r40.fig

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

r50.fig

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

r60.fig

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

r80.fig

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

r146.fig

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

r150.fig

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

r168.fig

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

r180.fig

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

r200.fig

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

r220.fig

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

r240.fig

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

r242.fig

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

r244.fig

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

r302.fig

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

r343.fig

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

r402.fig

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

r502.fig

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
鉴定完毕！
# b7 U0 }0 B, J0 |# `9 {' g
作者: 灼￡日    时间: 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