來源： http://blog.csdn.net/superdont

建立函數(shù)1：

function xdot=lorrenz(t,x)

for t=0:0.1:100;

if(cos(5.3*t)>=0)

p=1;

else

p=-1;

end

end

xdot=[-8/3*x(1)+p*x(2)*x(3);(-x(2)+x(3))*(25-10*cos(5.3*t));-p*x(1)*x(2)-(17.5*cos(5.3*t)+10.5)*x(2)+(13.3-14*cos(5.3*t))*x(3)];

建立文件2：

tic;

x0=[0;0;1e-10];tspan=[0,100];

[t,x]=ode45(@lorrenz,tspan,x0);

plot(t,x(:,1),'k');

r=length(t)-5000:length(t);

figure; plot(x(r,2),x(r,3),'k');

toc;

顯示相圖如下

輸出結(jié)果：

Elapsed time is 5.547000 seconds.

式中函數(shù)：

Tic

Toc

用于計(jì)算時(shí)間，例如：

for n = 1:100

A = rand(n,n);

b = rand(n,1);

tic

x = A/b;

t(n) = toc;

end

plot(t)

可輸出：

總結(jié)

以上是生活随笔為你收集整理的lorenz混沌系统matlab程序,Lorenz混沌系统相图的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。