LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
%a,b,c 为系统的三个变量，用ODE45数值解算法
5 {" k7 W8 I/ c$ p- y+ x$ A( qclear
' N- J! C" Z0 V( C# O3 D. Sglobal a b c
a=10;b=8/3;c=22;
6 Z0 @7 D9 [% y! _7 |8 P[T,X]=ode45('Lorenz',[0 100],[10;4;10]);
subplot(2,2,1)
7 J# S4 G& H+ Splot(X(:,1),X(:,2))
" ?4 Y+ w* ~. rsubplot(2,2,2)
plot(X(:,1),X(:,3))
subplot(2,2,3)
plot(X(:,2),X(:,3))
6 s1 d- p3 `( a8 Z; ?subplot(2,2,4)
; f# ^: S$ p4 t* Q# `: a' G9 Hplot3(X(:,1),X(:,2),X(:,3))
view(-20,60);
' [  x* p8 N* |2 e: m" O9 Cxlabel('x');ylabel('y');zlabel('z');



function 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)];


画出来的图形附件！根据不同的参数画出来的图

3 F# b; c/ u, B8 V& Z+ `

zhao6884269

试试运气啦~~~~~~~~~~~

lilijie_shangha

来报道!!!!!!!!!!!

wy371tt1

鉴定完毕！  

fengzhisheng106

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

想做就做

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

547896321

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

mybaby

哦～～

surfw

试试运气啦~~~~~~~~~~~

xiaoai

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