QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 14010|回复: 20
打印 上一主题 下一主题

LORENZ吸引子

[复制链接]
字体大小: 正常 放大
yunfenger 实名认证       

1

主题

3

听众

42

积分

升级  38.95%

该用户从未签到

跳转到指定楼层
1#
发表于 2010-7-28 20:46 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
如下MATLAB程序。。。。MATLAB6.50 J# a0 k. f) X
%a,b,c 为系统的三个变量,用ODE45数值解算法. X) _# L5 U0 p8 |+ x
clear
" p# D5 J: a4 a* ?+ {# [2 [9 cglobal a b c
: v5 c& h3 o  U2 b1 na=10;b=8/3;c=22;
! Q$ \7 X/ q8 R) h' ^- d' a, x[T,X]=ode45('Lorenz',[0 100],[10;4;10]);1 q; J1 E+ P7 }
subplot(2,2,1)$ D: R  e7 x/ c/ Q5 i
plot(X(:,1),X(:,2))
3 g- ~4 v) f) j0 L+ Fsubplot(2,2,2)
" r& k& `9 f/ x. W. gplot(X(:,1),X(:,3))
! r. Q. b( ]  j& ]1 c7 I3 I; isubplot(2,2,3)- m! V  C' V+ P$ W3 a8 S1 w4 a
plot(X(:,2),X(:,3))! J0 i5 ~" y" B$ K
subplot(2,2,4)
% J  A# f6 V  a( c/ D. L" mplot3(X(:,1),X(:,2),X(:,3))
  D, _5 w0 b$ [% s4 p8 _* tview(-20,60);) o& T9 t5 n+ c3 H' x! r
xlabel('x');ylabel('y');zlabel('z');5 y4 q3 w& ~8 R
7 I: E$ T4 I# K1 @
- W+ E7 \% o) s+ j3 L( k
. J$ ^- _5 }3 r$ Y0 z0 G2 v6 |, Y+ J
function dx=Lorenz(t,x)
$ p1 J( j/ a6 Bglobal a b c
- O0 K2 t/ g$ r: ldx=[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x(2);x(1)*x(2)-b*x(3)];
  J% J: L, I5 Z: ~: q# n; ^: O. d! m7 p. O, ]+ g+ Y; f0 }; U, ~

' [1 N' S8 k+ U6 }; E5 a画出来的图形附件!根据不同的参数画出来的图" a$ d) f; r8 {9 G. L; U' N! s/ [0 n) n

. ]' _% u1 d) F$ N3 \

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 点

zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信

0

主题

2

听众

31

积分

升级  27.37%

该用户从未签到

新人进步奖

回复

使用道具 举报

0

主题

2

听众

30

积分

升级  26.32%

该用户从未签到

新人进步奖

回复

使用道具 举报

wy371tt1        

0

主题

2

听众

70

积分

升级  68.42%

该用户从未签到

新人进步奖

回复

使用道具 举报

0

主题

2

听众

70

积分

升级  68.42%

该用户从未签到

新人进步奖

回复

使用道具 举报

0

主题

2

听众

98

积分

升级  97.89%

该用户从未签到

新人进步奖

回复

使用道具 举报

547896321        

0

主题

2

听众

30

积分

升级  26.32%

该用户从未签到

新人进步奖

回复

使用道具 举报

mybaby        

0

主题

2

听众

30

积分

升级  26.32%

该用户从未签到

新人进步奖

回复

使用道具 举报

surfw        

0

主题

2

听众

97

积分

升级  96.84%

该用户从未签到

新人进步奖

回复

使用道具 举报

xiaoai        

0

主题

2

听众

98

积分

升级  97.89%

该用户从未签到

新人进步奖

回复

使用道具 举报

您需要登录后才可以回帖 登录 | 注册地址

qq
收缩
  • 电话咨询

  • 04714969085
fastpost

关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

手机版|Archiver| |繁體中文 手机客户端  

蒙公网安备 15010502000194号

Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

GMT+8, 2026-5-27 05:39 , Processed in 0.445630 second(s), 107 queries .

回顶部