mjjqqq 发表于 2009-7-8 15:39

超混沌Lorenz系统的相图怎么实现

看论坛上一个代码，但是这个代码实现不了啊，请教

mjjqqq 发表于 2009-7-8 17:17

系统如下

比如dx/dt=a(y-x)+yz;dy/dt=bx-xz+w;dz/dt=-cz+x^2;dw/dt=-dy;用ode45求解器怎样求解，我试了试都不行，请高手指教
论坛里的一个程序如下：
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=;tspan=;
    =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;

无水1324 发表于 2009-7-8 19:05

麻烦先看一下matlab基础！

mjjqqq 发表于 2009-7-8 22:09

:@( 自己搞明白了

[ 本帖最后由 mjjqqq 于 2009-7-8 22:37 编辑 ]
页: [1]
查看完整版本: 超混沌Lorenz系统的相图怎么实现