离散状态的S函数

 

 该函数实现的功能是将输入的波形延迟后输出在示波器中。

完整代码如下：

function [sys,x0,str,ts]=sfun_und(t,x,u,flag)
switch flag
    case 0,
        [sys,x0,str,ts]=mdlInitializeSizes;
    case 2,
        sys=mdlUpdate(t,x,u);
    case 3,
        sys=mdlOutputs(t,x,u);
    case 9,
        sys=[];
    otherwise 
        error(['Unhandled flag=',num2str(flag)]);    
end

function [sys,x0,str,ts]=mdlInitializeSizes(len)
sizes=simsizes;%返回一个变量，
sizes.NumContStates=0;%连续变量个数
sizes.NumDiscStates=1;%离散变量个数
sizes.NumOutputs=1;%输出个数
sizes.NumInputs=1;%输入个数
sizes.DirFeedthrough=0;%直接贯通，0或1，当时输出值直接依赖于同一时刻的输入时为1
sizes.NumSampleTimes=1;%采样时间
sys=simsizes(sizes);%返回值
x0=0;
str=[];
ts=[0.1 0];


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


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

 分析上述代码：

首先，主函数部分如下，程序在执行时会从这里开始：

function [sys,x0,str,ts]=sfun_und(t,x,u,flag)
switch flag
    case 0,
        [sys,x0,str,ts]=mdlInitializeSizes;
    case 2,
        sys=mdlUpdate(t,x,u);
    case 3,
        sys=mdlOutputs(t,x,u);
    case 9,
        sys=[];
    otherwise 
        error(['Unhandled flag=',num2str(flag)]);    
end

 

第一句，函数的输出部分为　[sys,x0,str,ts],其中sys是变量的返回值，当flag=3时表示输出结果值；x0是初始状态值，当flag=0时才有值；str是用来备用的，通常设置为空矩阵[]；ts是包含模块的采样时间和偏差值得两列矩阵。
ts=[0,0]时，表示连续采样；
ts=[-1,0]时，以驱动模块相同的速率运行；
ts=[0.2 0.1]时，系统仿真从0.1s开始，每隔0.2s运行一次。
接下来就是调用了

 

 

 接下来是初始化函数

function [sys,x0,str,ts]=mdlInitializeSizes(len)
sizes=simsizes;%开始调用simsizes函数
sizes.NumContStates=0;%连续变量个数
sizes.NumDiscStates=1;%离散变量个数
sizes.NumOutputs=1;%输出个数
sizes.NumInputs=1;%输入个数
sizes.DirFeedthrough=0;%直接贯通，0或1，当时输出值直接依赖于同一时刻的输入时为1
sizes.NumSampleTimes=1;%采样时间
sys=simsizes(sizes);%返回值
x0=0;
str=[];
ts=[0.1 0];

 

其中，sizes=simsizes;返回未初始化的sizes结构，然后下面开始定义sizes的结构属性

 

 初始化结束后，再将sizes结构中的信息传递给sys，即：sys=simsizes(sizes);

ts=[0.1 0]是指从0时刻起，每隔0.1s执行一次（10Hz）