数学建模社区-数学中国

标题: LORENZ吸引子 [打印本页]
作者: yunfenger    时间: 2010-7-28 20:46
标题: LORENZ吸引子
如下MATLAB程序。。。。MATLAB6.5  y( V; B6 H- Z
%a,b,c 为系统的三个变量，用ODE45数值解算法
: o6 \' q- @, [' nclear
: ^' X9 L. S9 u& c9 [& _global a b c
/ V* W( |1 _' R, z9 Wa=10;b=8/3;c=22;
5 f# A" f6 x2 `0 p[T,X]=ode45('Lorenz',[0 100],[10;4;10]);; Q$ g# ]1 i0 G9 q, g, r" ?
subplot(2,2,1)' i" h: ?8 K1 O( H9 ^* @. r
plot(X(:,1),X(:,2))
! E0 |3 w* a6 k' lsubplot(2,2,2)
+ m& `+ g8 e$ u! `5 Fplot(X(:,1),X(:,3))
" ~* n1 X) K) R2 B7 p4 wsubplot(2,2,3)
2 W; f# ?, f" u( C, S: uplot(X(:,2),X(:,3))
$ R  Q( g7 R1 w( Y1 \' N7 Z: Wsubplot(2,2,4)' f0 @5 l3 g6 s8 P/ m% D
plot3(X(:,1),X(:,2),X(:,3))1 f' E( }; L2 |0 I0 c* w! i
view(-20,60);  W5 ]* ^. V/ c# w5 k* r
xlabel('x');ylabel('y');zlabel('z');
- R4 Y* d; q( U7 ~% c' ?2 t0 s- O/ r) l' F# h

4 H1 b7 R7 C- M- E
$ H/ R' w4 ?* S9 q' xfunction dx=Lorenz(t,x)) w) @& M2 s9 W( [! v7 Z1 H4 }
global a b c
% x/ Y8 L% }1 C2 E, Z% {dx=[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x(2);x(1)*x(2)-b*x(3)];: f% K* s# s$ G# ?* |1 e. a7 O  o
; ?" U4 w; o# e1 O

7 V  b7 c6 v" @2 M. @% ~画出来的图形附件！根据不同的参数画出来的图0 }2 y* |+ Y- x

& r2 c3 i) @$ m6 @) ~

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
鉴定完毕！
: Q5 R" T% x5 ?* K4 U
作者: 灼￡日    时间: 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