【S-Function】离散状态变量变量

%**************************************************************************
% x1(k+1)=x1(k)+2*x2(k)
% x2(k+1)=3x1(k)+2*u(k)
% y(k+1)=2x1(k)+u(k)
%**************************************************************************

function [sys,x0,str,ts]=ls2ztfc(~,x,u,flag)
switch flag,
case 0                         % 调用初始化函数
   [sys,x0,str,ts] = mdlInitializeSizes();
case 2                         % 调用离散状态的更新函数
   sys = mdlUpdates(x,u);
case 3                         % 调用输出量的计算函数
   sys = mdlOutputs(x,u);
case {1,4,9}                 % 未使用的 flag 值
   sys = [];
otherwise                      % 处理错误
   error(['Unhandled flag = ',num2str(flag)]);
end;
% 当 flag 为 0 时进行整个系统的初始化
function [sys,x0,str,ts] = mdlInitializeSizes()
% 首先调用 simsizes 函数得出系统规模参数 sizes, 并根据离散系统的
% 实际情况设置 sizes 变量
sizes = simsizes;              % 读入初始化参数模板
sizes.NumContStates = 0;       % 连续状态变量个数：0
sizes.NumDiscStates = 2;       % 离散状态变量个数：2
sizes.NumOutputs = 2;          % 输出两个量：跟踪信号和微分信号
sizes.NumInputs = 1;           % 系统输入信号一路
sizes.DirFeedthrough = 1;      % 输入直接传到输出口
sizes.NumSampleTimes = 1;      % 单个采样周期