LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
- ?0 x2 G; n" U7 {3 h, L%a,b,c 为系统的三个变量，用ODE45数值解算法
% H7 ^" [" ?0 L, jclear
5 v# x6 `$ J, Y/ [1 G  Xglobal 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))
! r& ^+ V# T# _8 v# C( a+ z# y- dsubplot(2,2,2)
plot(X(:,1),X(:,3))
/ R- F3 Z3 R3 [$ O2 e. J! ^; msubplot(2,2,3)
plot(X(:,2),X(:,3))
2 D- ^3 n. w: j1 }subplot(2,2,4)
1 p& W* t1 U2 j8 i, @/ Y+ Fplot3(X(:,1),X(:,2),X(:,3))
, h! ^! E' T. V1 T7 l* w: Aview(-20,60);
4 D' }  X4 d/ _8 U6 x9 A! Jxlabel('x');ylabel('y');zlabel('z');

0 m! p! i9 Q9 |/ u

function dx=Lorenz(t,x)
# ?) Q2 [1 H0 _# O- @global a b c
' R2 k1 P; ?- Qdx=[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x(2);x(1)*x(2)-b*x(3)];
0 c  @- Z) d8 m1 a3 N8 k1 i
4 x% |3 `) D) B& F0 O* H, N
% ]" F, U  t0 Z" E6 Y! y" K# N9 [画出来的图形附件！根据不同的参数画出来的图
3 i9 A/ s8 F& H9 p' S  b% o0 F: }

zhao6884269

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

lilijie_shangha

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

wy371tt1

鉴定完毕！  

fengzhisheng106

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

想做就做

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

547896321

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

mybaby

哦～～

surfw

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

xiaoai

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