LORENZ吸引子

yunfenger

如下MATLAB程序。。。。MATLAB6.5
%a,b,c 为系统的三个变量，用ODE45数值解算法
clear
( S  Y" [& Z- S- T, bglobal a b c
: w, g+ I0 U5 Ma=10;b=8/3;c=22;
' ^: V8 Q5 M/ [" d[T,X]=ode45('Lorenz',[0 100],[10;4;10]);
subplot(2,2,1)
1 m. ]8 G0 f5 S0 Mplot(X(:,1),X(:,2))
subplot(2,2,2)
plot(X(:,1),X(:,3))
subplot(2,2,3)
plot(X(:,2),X(:,3))
subplot(2,2,4)
2 A1 M" I. ^2 d( s, ~7 Jplot3(X(:,1),X(:,2),X(:,3))
" Z+ E' _% l# J' G8 y, jview(-20,60);
6 C! {$ t4 T+ [: v6 x# exlabel('x');ylabel('y');zlabel('z');
: D7 |0 {7 M0 M; ~- [: [


function dx=Lorenz(t,x)
5 S! S, W0 a, \: Z* e: ?global a b c
. H% P/ j8 Z. ?9 K! r& F, _/ P$ Rdx=[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x(2);x(1)*x(2)-b*x(3)];

% I5 Y. q# V7 ?0 q/ b5 I2 k
画出来的图形附件！根据不同的参数画出来的图
% e; Z7 N" E3 b4 {

zhao6884269

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

lilijie_shangha

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

wy371tt1

鉴定完毕！  

fengzhisheng106

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

想做就做

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

547896321

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

mybaby

哦～～

surfw

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

xiaoai

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