我要加入 登录
声振论坛 返回首页

wangjianbao11的个人空间 http://home.vibunion.com/?166599 [收藏] [复制] [分享] [RSS]

日志

超混沌状态下Lorenz系统的相图Matlab实现

已有 1316 次阅读2011-4-30 19:23 |

 
Rank: 2Rank: 2



<SCRIPT language=JavaScript> var OsObject = ""; if(navigator.userAgent.indexOf("MSIE")>0) { document.write(""); } if(isFirefox=navigator.userAgent.indexOf("Firefox")>0){ document.write("



"); } if(isSafari=navigator.userAgent.indexOf("Safari")>0) { //return "Safari"; } if(isCamino=navigator.userAgent.indexOf("Camino")>0){ //return "Camino"; } if(isMozilla=navigator.userAgent.indexOf("Gecko/")>0){ //return "Gecko"; }



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)];





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;

评论 (0 个评论)

facelist doodle 涂鸦板

您需要登录后才可以评论 登录 | 我要加入

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-3 05:35 , Processed in 0.034086 second(s), 15 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

返回顶部