yulei.asv

LQR源程序

global DeltaCav DeltaFinE Theta
  Theta=0;DeltaCav=0;DeltaFinE=0;
   y0=[60 0 0 Theta*pi/180 0 0];%给函数赋初值
   tspan=[0 10];
     [t,y]=ode45('projectiledes',tspan,y0,[],DeltaCav,DeltaFinE,Theta);
        
        subplot(2,3,5),plot(y(:,5),y(:,6))
      % plot(t,y(:,6))  
       


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