LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
%a,b,c 为系统的三个变量，用ODE45数值解算法
& t- d9 n! q  j$ A$ l: w3 |% `& wclear
, f' {/ I9 N4 N* u1 W; Y! Oglobal a b c
a=10;b=8/3;c=22;
[T,X]=ode45('Lorenz',[0 100],[10;4;10]);
subplot(2,2,1)
plot(X(:,1),X(:,2))
8 i/ @- P1 m$ Q0 `' J. |7 `' Hsubplot(2,2,2)
" f. Q) f' I2 Eplot(X(:,1),X(:,3))
subplot(2,2,3)
plot(X(:,2),X(:,3))
/ z4 p% x5 _. m* a+ j8 x4 `subplot(2,2,4)
plot3(X(:,1),X(:,2),X(:,3))
view(-20,60);
xlabel('x');ylabel('y');zlabel('z');

: `1 N: s- v  n0 c1 S0 _( \1 w

# I: ?, G* z! W( g3 `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)];
' v, c, J* Q- X  s! ?

画出来的图形附件！根据不同的参数画出来的图
, F) x" @2 ~  n9 [/ R5 ?+ n
6 \4 q- Z; h6 O

zhao6884269

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

lilijie_shangha

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

wy371tt1

鉴定完毕！  

fengzhisheng106

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

想做就做

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

547896321

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

mybaby

哦～～

surfw

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

xiaoai

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