yulei.m

LQR源程序

function a=yulei()
format compact
 Theta
   y0=[60 0 0 Theta*pi/180 0 0];%给函数赋初值
   tspan=[0 10];
     [t,y]=ode45('projectiledes',tspan,y0);
        subplot(2,3,1),plot(t,y(:,1))
        subplot(2,3,2),plot(t,y(:,2))
        subplot(2,3,3),plot(t,y(:,3))
        subplot(2,3,4),plot(t,y(:,4))
        subplot(2,3,5),plot(y(:,5),y(:,6))
      % plot(t,y(:,6))  
       


%--------------------------------------------------------------------------