LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
/ S* ?$ V2 W+ X$ g%a,b,c 为系统的三个变量，用ODE45数值解算法
clear
global a b c
a=10;b=8/3;c=22;
2 ]3 g/ z7 p" z3 u4 S! O0 {[T,X]=ode45('Lorenz',[0 100],[10;4;10]);
subplot(2,2,1)
plot(X(:,1),X(:,2))
subplot(2,2,2)
# l  h! r$ o- f" rplot(X(:,1),X(:,3))
% y; }. \5 ]! t/ Q  u- i: `$ Psubplot(2,2,3)
- P" o5 p% _3 s( G  eplot(X(:,2),X(:,3))
- _6 \- J" G6 K% D- Asubplot(2,2,4)
plot3(X(:,1),X(:,2),X(:,3))
# x) X  P! p. @6 o; ]view(-20,60);
# b* t, O" M7 E3 Vxlabel('x');ylabel('y');zlabel('z');
, d5 F- q) k8 x1 c) i' m3 y% V
; E6 W6 u9 {3 j& t" ?4 e' J
! o, X, [* C3 h0 G; _& r, K7 d
function dx=Lorenz(t,x)
global a b c
" z3 l$ p4 @: O7 c, M8 W' b1 Adx=[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x(2);x(1)*x(2)-b*x(3)];
7 j' ^) ~/ v4 r- V; L4 _" ~
; R( t/ W/ C& |" `
画出来的图形附件！根据不同的参数画出来的图
) p# x$ B, L7 H* J
0 B! O# b# T6 U- t

zhao6884269

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

lilijie_shangha

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

wy371tt1

鉴定完毕！  

fengzhisheng106

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

想做就做

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

547896321

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

mybaby

哦～～

surfw

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

xiaoai

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