lregulation.asv

氧乐果控制过程

function [sys,x0]=lregulation(t,x,u,flag)
global ke kec ku n;
ke=1 ;
kec=1;
ku=1;
if flag==0
sys=[0;0;3;3;0;1]
x0=[]
w=0
elseif flag==3
%   u(2)
  w=w+1
% t
if abs(u(2))>0.2 | abs(u(3))>0.1
    sys(1)=0.4*ke*u(2);
    sys(2)=0.02*kec*u(3);
    sys(3)=ku+14;
elseif abs(u(2))>0.1 | abs(u(3))>0.04
    sys(1)=0.5*ke*u(2);
    sys(2)=0.27*kec*u(3);
    sys(3)=ku+12;
else
    sys(1)=0.8*ke*u(2);
    sys(2)=0.6*kec*u(3);
    sys(3)=ku;
end
else
    sys=[];
end