基于simulink的离散系统

题目：编写 S－函数实现输出对输入的单位延迟，即 y( k+1) = u(k)。

解：单位延迟系统的状态方程可以表述为：x(k+1)=u(k),y(k+1)=x(k+1).

s-function的代码具体如下所示：

function [sys,x0,str,ts,simStateCompliance] = sfun_unitdelay(t,x,u,flag)

switch flag,
    
  case 0,
    [sys,x0,str,ts,simStateCompliance]=mdlInitializeSizes;

  case 1,
    sys=mdlDerivatives(t,x,u);

  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
    DAStudio.error('Simulink:blocks:unhandledFlag', num2str(flag));

end

function [sys,x0,str,ts,simStateCompliance]=mdlInitializeSizes

sizes = simsizes;
sizes.NumContStates  = 0;
sizes.NumDiscStates  = 1;
sizes.NumOutputs     = 1;
sizes.NumInputs      = 1;
sizes.DirFeedthrough = 0;
sizes.NumSampleTimes = 1;   % at least one sample time is needed
sys = simsizes(sizes);     %%将size结构传到sys中

x0  = 0;   %%初始状态向量，由传入的参数决定，没有为空  
str = [];
ts  = [0.1 0];  %%设置采样时间。第一位表示间隔的时间，连续系统为0。第二位表示偏移量。
                %%这里表示在S函数仿真开始后每隔0.1s运行一次
simStateCompliance = 'UnknownSimState';

function sys=mdlDerivatives(t,x,u)

sys = [];

function sys=mdlUpdate(t,x,u)

sys = u;

function sys=mdlOutputs(t,x,u)

sys = x;

function sys=mdlGetTimeOfNextVarHit(t,x,u)

sampleTime = 1;    %  Example, set the next hit to be one second later.
sys = t + sampleTime;

function sys=mdlTerminate(t,x,u)

sys = [];

simulink框图如下所示：

最后

以上就是落后灰狼为你收集整理的基于simulink的离散系统的全部内容，希望文章能够帮你解决基于simulink的离散系统所遇到的程序开发问题。

如果觉得靠谱客网站的内容还不错，欢迎将靠谱客网站推荐给程序员好友。