sfun08832_1.m

来自「MATLAB 6.5 从基础讲解其功能」· M 代码 · 共 52 行

function [sys,x0,str,ts] = sfun08832_1(t,x,u,aa,lb,ub,xi)
switch aa,
case 0,
    [sys,x0,str,ts]=mdlInitializeSizes(xi);
case 1,
    sys=mdlDerivatives(t,x,u,lb,ub);
case 2,
    sys=mdlUpdate(t,x,u);
case 3,
    sys=mdlOutputs(t,x,u);
case 4,
    sys=mdlGetTimeOfNextVarHit(t,x,u);
case 9,
    sys=mdlTerminate(t,x,u);
otherwise
    error(['Unhandled aa = ',num2str(aa)]);
end

function [sys,x0,str,ts]=mdlInitializeSizes(xi)
sizes = simsizes;
sizes.NumContStates  = 1;							%一个连续状态变量
sizes.NumDiscStates  = 0;
sizes.NumOutputs     = 1;							%一个输出
sizes.NumInputs      = 1;							%一个输入
sizes.DirFeedthrough = 0;
sizes.NumSampleTimes = 1; 

sys = simsizes(sizes);
x0  = xi;
str = [];
ts  = [0 0];

function sys=mdlDerivatives(t,x,u,lb,ub)
if (x<=lb & u<0)|(x>=ub & u>0)
    sys = 0;
else
    sys=u;
end

function sys=mdlUpdate(t,x,u)
sys = [];

function sys=mdlOutputs(t,x,u)
sys =x;

function sys=mdlGetTimeOfNextVarHit(t,x,u)
sampleTime = 1;    
sys = t + sampleTime;

function sys=mdlTerminate(t,x,u)
sys = [];