数学建模社区-数学中国

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

作者: yunfenger    时间: 2010-7-28 20:46
标题: LORENZ吸引子
如下MATLAB程序。。。。MATLAB6.5
2 v* h8 C. e. T%a,b,c 为系统的三个变量,用ODE45数值解算法. O' H/ z+ e5 F, s* j0 Q
clear
0 w" c  _7 M; r" v  Hglobal a b c" n5 o. n$ E# V
a=10;b=8/3;c=22;  g8 `+ ^6 Y/ A+ h( [0 z) H4 `
[T,X]=ode45('Lorenz',[0 100],[10;4;10]);6 ?  _& O9 O& l  @( Z
subplot(2,2,1)! o, y; a- [2 E0 {. F% ~4 P0 w
plot(X(:,1),X(:,2))( `" i+ L3 P; `+ [2 m5 `1 i
subplot(2,2,2)0 f7 P/ C; @3 [0 b" D( N) n! w8 P
plot(X(:,1),X(:,3))# o; T' {# b' ~1 \& }; F. a. ]
subplot(2,2,3)
2 j6 _9 Q- @) B  H0 q1 G5 eplot(X(:,2),X(:,3))
4 k$ E0 B  C9 p/ x, c$ X. ksubplot(2,2,4)1 {4 c% ]- ~1 T& W$ C( J5 v( f8 x
plot3(X(:,1),X(:,2),X(:,3))
! f; l0 [8 k4 o0 y* b) v- I; {view(-20,60);
( {3 h) @0 J; U& `& axlabel('x');ylabel('y');zlabel('z');5 Z0 N. C4 p+ L; C9 I  i

! G: a3 p8 N7 O& [6 y2 Y& h7 O0 A, t( _
, p! K9 b! S) g
function dx=Lorenz(t,x)' i. b( B+ q! u" `; \# H$ I
global a b c
# o# s7 I0 g8 L) m" \dx=[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x(2);x(1)*x(2)-b*x(3)];
5 {# q( b0 Z% |0 ^
' s8 G. h* ~, c+ {. ]0 X6 \2 m2 a7 E; t6 L: C! E
画出来的图形附件!根据不同的参数画出来的图8 \4 C  m9 u# V' s- p4 G
* u$ \0 v" ?0 n) X8 u) M# ^

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
鉴定完毕!
2 y2 _$ n/ s6 z) [) H
作者: 灼£日    时间: 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