LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
) R+ h$ Z9 Y5 Z0 _%a,b,c 为系统的三个变量，用ODE45数值解算法
) s! T, e& ~+ ], Rclear
# \! [0 c4 {, u+ z0 \global a b c
$ s! ?1 g. W$ N  O! p2 ?% Ya=10;b=8/3;c=22;
" W6 ~7 W3 M% ~0 A8 a$ h- j, |[T,X]=ode45('Lorenz',[0 100],[10;4;10]);
subplot(2,2,1)
  U! z) Y% L4 f: Z8 O5 gplot(X(:,1),X(:,2))
subplot(2,2,2)
9 E# m9 z8 ?  P0 u: Y  \( f6 Vplot(X(:,1),X(:,3))
/ d' ~' u5 n# k! Usubplot(2,2,3)
- V6 t  ~0 p8 n% u7 T9 p  yplot(X(:,2),X(:,3))
subplot(2,2,4)
plot3(X(:,1),X(:,2),X(:,3))
' ]8 u2 l2 a) T& Y" lview(-20,60);
" y" \/ y, ?, t; ^; X2 oxlabel('x');ylabel('y');zlabel('z');

0 l8 I, {! @4 B9 {  C5 b( K
! h3 r3 m. z$ u+ w2 R" L! d: A) E: y
function dx=Lorenz(t,x)
3 O- W! j9 s9 X2 \global a b c
# k% i! v! D2 Q- F2 V" |0 Idx=[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x(2);x(1)*x(2)-b*x(3)];
" ^1 w" e, ~) x* X5 F2 g3 k
& E1 k% T4 p2 s: a
: F; j/ {: s1 h+ X# _+ P画出来的图形附件！根据不同的参数画出来的图
7 l. K# y8 \) z" A" w4 j

zhao6884269

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

lilijie_shangha

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

wy371tt1

鉴定完毕！  

fengzhisheng106

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

想做就做

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

547896321

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

mybaby

哦～～

surfw

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

xiaoai

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