lregulation.m

氧乐果控制过程

function [sys,x0]=lregulation(t,x,u,flag)
global ke kec ku n s uu;

s=0;
if flag==0
sys=[0;0;3;3;0;1]
x0=[]
n=0
ke=1 
kec=1
ku=1
uu=[]
% u(2)
elseif flag==3
  u(2)
  n=n+1
  uu(n)=n
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