matlab解duffing方程

美好的日子

我的毕业设计题目就是混沌序列生成器，我是一个初学者，我们老师的要求是 用matlab解出每个混沌方程的解 然后动态的画出当参数确定时的混沌图形。哪位好心人帮助我一下，小弟先谢谢了。用matlab怎么解duffing方程并画出它的三维图形？然后怎么利用pause（n）来实现它的时间延迟？？我找了很多天的资料，用ODE45解了洛伦兹方程，也画出了图形，但是图形是一瞬间就形成了 怎么能求出每个 t 时刻的坐标 然后利用pause（n）来实现呢？代码大致如下：
先定义了一个LORENZED.M文件：function ydot=LORENZED(t,y)
                                                       global sigma r b
                                                       ydot=[b,0,y(2);0,-sigma,sigma;-y(2),r,-1]*y;

然后在另一个M文件中，按钮的回调函数代码： global sigma r b
                                                                          sigma=str2num(get(handles.edit1,'string'))         
                                                                           r=str2num(get(handles.edit2,'string'))
                                                                          b=str2num(get(handles.edit3,'string'))
        
                                                                           axis([10 40 -30 30 -30 30]);
                                                                             view(3)
                                                                                 hold on
                                                                            title('Lorenz ATTractor')
   
                                                                           [t,y]=ode45('LORENZED',[0,40],[0,0,0.00000000001]);
   
                                                                            figure;
                                                                               u=y(:,1);
                                                                                 v=y(:,2);
                                                                                w=y(:,3);
                                                                                   plot3(u,v,w)