LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
%a,b,c 为系统的三个变量，用ODE45数值解算法
clear
- B, P0 U) X, Cglobal a b c
a=10;b=8/3;c=22;
[T,X]=ode45('Lorenz',[0 100],[10;4;10]);
0 c3 ^2 D7 U: M! Z; Msubplot(2,2,1)
# t* N" P6 {  k3 n: M5 P- tplot(X(:,1),X(:,2))
2 ~& W! i" \, J) R$ E$ z, Jsubplot(2,2,2)
plot(X(:,1),X(:,3))
subplot(2,2,3)
  R. r" Y1 d2 `" @) J; cplot(X(:,2),X(:,3))
, z8 m* _2 v, e" }( t% P8 l, A: m  Qsubplot(2,2,4)
* @) q0 v  v$ D+ X& Q' ~plot3(X(:,1),X(:,2),X(:,3))
view(-20,60);
; J) q" g. P, w# H! u# Fxlabel('x');ylabel('y');zlabel('z');


# v( E! `1 L  d# F
function dx=Lorenz(t,x)
global a b c
$ |* p5 m9 f8 G/ m5 N7 Ndx=[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x(2);x(1)*x(2)-b*x(3)];
* \) K. q" D5 {# D% U
0 ?9 C5 B! b9 `. x, m' }% k+ V: U
画出来的图形附件！根据不同的参数画出来的图
1 H( a$ e5 O- B8 [; k+ n" l: z

zhao6884269

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

lilijie_shangha

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

wy371tt1

鉴定完毕！  

fengzhisheng106

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

想做就做

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

547896321

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

mybaby

哦～～

surfw

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

xiaoai

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