LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
%a,b,c 为系统的三个变量，用ODE45数值解算法
clear
& ?8 e, e' K7 j- L% Gglobal a b c
a=10;b=8/3;c=22;
5 C  k8 z  d, {5 |. k[T,X]=ode45('Lorenz',[0 100],[10;4;10]);
subplot(2,2,1)
plot(X(:,1),X(:,2))
subplot(2,2,2)
% |/ S# B5 U: C: g, `plot(X(:,1),X(:,3))
+ E/ G/ w. H, Y: s2 C! Q; Hsubplot(2,2,3)
plot(X(:,2),X(:,3))
subplot(2,2,4)
' j& v% w1 A1 e6 J" ~- tplot3(X(:,1),X(:,2),X(:,3))
view(-20,60);
xlabel('x');ylabel('y');zlabel('z');
, F" B' T/ z* B3 u5 a. c
' y6 b: |5 A; h# \. c
8 p; H% o( h5 D3 R/ K% ]. K1 {
6 T  }  ?( V8 ]' j- Sfunction dx=Lorenz(t,x)
5 a, g1 Q1 [4 a  Z* vglobal a b c
* |0 L1 P1 T; N7 ]* H5 Kdx=[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x(2);x(1)*x(2)-b*x(3)];
8 f# _! g8 e% J) O" l0 ?$ T% b% u
" q: I' r3 J  i, C, R8 Z
/ L1 K' B2 [6 u( J7 w画出来的图形附件！根据不同的参数画出来的图
- }3 y4 |% q( }; n' p' L
1 s  |" ?; ^, }6 @; B, `- F; s

zhao6884269

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

lilijie_shangha

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

wy371tt1

鉴定完毕！  

fengzhisheng106

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

想做就做

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

547896321

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

mybaby

哦～～

surfw

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

xiaoai

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