LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
%a,b,c 为系统的三个变量，用ODE45数值解算法
0 _4 i1 _& h# U- q( \7 `. {clear
global a b c
a=10;b=8/3;c=22;
[T,X]=ode45('Lorenz',[0 100],[10;4;10]);
7 M! ^0 c9 R( Tsubplot(2,2,1)
7 u, P2 s+ n4 w2 P' n- @8 i& Wplot(X(:,1),X(:,2))
, r6 C/ y/ @& H* U% k1 gsubplot(2,2,2)
plot(X(:,1),X(:,3))
. ?* @% I( U0 l) hsubplot(2,2,3)
plot(X(:,2),X(:,3))
subplot(2,2,4)
% C$ J( R# ?/ t7 `+ @plot3(X(:,1),X(:,2),X(:,3))
' u6 t! N$ B! ~; `% S# Z$ f$ f1 n1 Cview(-20,60);
6 C3 G1 H7 \0 m. B% Vxlabel('x');ylabel('y');zlabel('z');



function dx=Lorenz(t,x)
0 B0 x! N+ s0 \; Zglobal a b c
% {2 i7 j: o$ }" ?3 T8 }; W5 Z0 Vdx=[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x(2);x(1)*x(2)-b*x(3)];

. ]- Z0 _: J" \0 c/ }
画出来的图形附件！根据不同的参数画出来的图
  S( _$ A5 k( f* m3 b
% B0 X; e: i2 d* y& S

zhao6884269

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

lilijie_shangha

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

wy371tt1

鉴定完毕！  

fengzhisheng106

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

想做就做

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

547896321

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

mybaby

哦～～

surfw

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

xiaoai

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