regulation.m

氧乐果控制过程

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