LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
%a,b,c 为系统的三个变量，用ODE45数值解算法
" v% x7 }) n5 u3 H! [. nclear
global a b c
" Z/ D' O7 G8 S: c* @9 ]a=10;b=8/3;c=22;
* g5 t; f( T! ~% Z8 n[T,X]=ode45('Lorenz',[0 100],[10;4;10]);
# F" Q, W/ d- L. r* Usubplot(2,2,1)
! c/ k/ |( q! D6 A$ Uplot(X(:,1),X(:,2))
7 s2 w, {" K+ U8 qsubplot(2,2,2)
plot(X(:,1),X(:,3))
subplot(2,2,3)
3 o2 `9 x* f- Z$ E3 Bplot(X(:,2),X(:,3))
1 ]0 m# }& k+ `' u4 b  `subplot(2,2,4)
plot3(X(:,1),X(:,2),X(:,3))
view(-20,60);
xlabel('x');ylabel('y');zlabel('z');

+ j( R1 q/ b8 c% X. M

+ v; W+ u% d" Vfunction dx=Lorenz(t,x)
global a b c
dx=[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x(2);x(1)*x(2)-b*x(3)];
1 Y& w6 ^5 E: _- ~; x! N) p" X! [

" M" H: S) _4 n, K; p画出来的图形附件！根据不同的参数画出来的图
  d' |! ~! S8 x4 w# w* K
# r- n6 v) s; h* u

天马绝尘

试试运气   

hngy

Lorentz混沌这么费体力

灼￡日

程序报错地方多

枫露之茗

鉴定完毕！

weixinmaths

guohuawin

声明一下：本人看贴和回贴的规则，好贴必看，精华贴必回。

ycy74_2005

我来了~~~~~~~~~ 闪人~~~~~~~~~~~~~~~~  

sigh

强烈支持。楼主万岁

阿里wudi

楼主，你写得实在是太好了。我惟一能做的，就只有把这个帖子顶上去这件事了

xiaoai

强人，佩服死了。呵呵，不错啊