数学建模社区-数学中国

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

作者: yunfenger    时间: 2010-7-28 20:46
标题: LORENZ吸引子
如下MATLAB程序。。。。MATLAB6.5" Z, H# Q9 D9 ^6 ], [
%a,b,c 为系统的三个变量,用ODE45数值解算法" u( j& Z3 l2 k8 i1 Q! y
clear
% F  H" x4 B2 r+ iglobal a b c9 y% C4 i7 l1 c- P9 M5 A
a=10;b=8/3;c=22;, I0 U2 @7 G* x: @# S" }
[T,X]=ode45('Lorenz',[0 100],[10;4;10]);2 V2 y+ V  m  T5 e
subplot(2,2,1)
% C4 T" K$ o; ^- @8 \' Hplot(X(:,1),X(:,2)). R9 d& Y3 Z: L4 c7 `" d: e; K
subplot(2,2,2)* ~3 z# g* ]* z4 x# x
plot(X(:,1),X(:,3))
) d; w4 `7 G* y( M6 qsubplot(2,2,3)1 U: V2 y' J+ w0 x0 f5 {
plot(X(:,2),X(:,3))" ]/ h: g* L- P# N* n
subplot(2,2,4)
% w4 A/ {: K* M  N  W6 l* yplot3(X(:,1),X(:,2),X(:,3))2 n- q, M" I, q9 P3 _( b& X  h
view(-20,60);8 o8 Y& J! |+ s9 J5 B" Y" G
xlabel('x');ylabel('y');zlabel('z');2 x9 U- S2 j2 X6 K
. o/ H- ]( P7 `. C4 Q) m4 D. a
4 U* `( _' y: \3 k' K9 x
5 B4 `4 {8 B( ?. k/ U
function dx=Lorenz(t,x)0 ~. u0 B) l, u& L/ m. P! |
global a b c
2 n+ _* i$ c, \9 p* e" Tdx=[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x(2);x(1)*x(2)-b*x(3)];
$ }3 [5 l. e7 `1 r8 z+ i4 b& q$ @- Z
9 W" q, G* B9 J5 t5 T' J3 P1 H" P: z! m2 u+ t4 |$ e
画出来的图形附件!根据不同的参数画出来的图3 O4 V) J6 K) _  Q
( Z* C2 ^: U7 Z$ ^7 Z+ L/ R

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
鉴定完毕!) [' B- O! N3 Q. L/ s( ~

作者: 灼£日    时间: 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