matlab河南工业大学课件第6章.pptVIP

* 6.2 Simulink仿真技术与S函数 (2)M文件S-函数模板 Simulink为用户编写S-函数提供了各种模板文件，其中定义了S-函数完整的框架结构，用户可以根据自己的需要加以修改。 编写M文件S-函数时，推荐使用S-函数模板文件sfuntmpl.m。这个文件是一个完整的M文件S-函数，它包含1个主函数和6个子函数。 在主函数内，程序根据标志变量Flag，将执行流程转移到相应的子函数，即例程函数。Flag标志量作为主函数的参数由系统（Simulink引擎）调用时给出。 * 6.2 Simulink仿真技术与S函数 若要打开模板文件，可在MATLAB命令窗口输入“edit sfuntmpl”命令，模板如下： function [sys,x0,str,ts] = sfuntmpl(t,x,u,flag)switch flag, case 0, [sys,x0,str,ts]=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 error([Unhandled flag = ,num2str(flag)]);end * 6.2 Simulink仿真技术与S函数 function [sys,x0,str,ts]=mdlInitializeSizessizes = simsizes; % 生成sizes数据结构 sizes.NumContStates = 0; % 连续状态数，缺省为 0sizes.NumDiscStates = 0; % 离散状态数，缺省为 0sizes.NumOutputs = 0; % 输出量个数，缺省为 0sizes.NumInputs = 0; % 输入量个数，缺省为 0sizes.DirFeedthrough = 1; %是否存在直接馈通。1：存在；0：不存在，缺省为 1。sizes.NumSampleTimes = 1; % 采样时间个数，至少是一个sys = simsizes(sizes); % 返回sizes数据结构所包含的信息 * 6.2 Simulink仿真技术与S函数 % initialize the initial conditionsx0 = []; % 设置初始状态str = []; % 保留变量置空ts = [0 0]; % 采样时间：[采样周期 偏移量]，采样周期为0表示是连续系统计算导数例程子函数：给定 t，x，u， 计算连续状态的导数，用户应该在此给出系统的连续状态方程。该子函数可以不存在。 * 6.2 Simulink仿真技术与S函数 function sys=mdlDerivatives(t,x,u) sys = [];  % sys 表示状态导数，即dx。状态更新例程子函数：给定 t，x，u， 计算离散状态的更新。每个仿真步长必然调用该子函数，不论是否有意义。用户除了在此描述系统的离散状态方程外，还可以填入其它每个仿真步长都有必要执行的代码function sys=mdlUpdate(t,x,u) sys = []; % sys 表示下一个离散状态 即 x(k+1)。计算输出例程子函数：给定 t，x，u， 计算输出。该子函数必须存在，用户可以在此描述系统的输出方程 * 6.2 Simulink仿真技术与S函数 function sys=mdlOutputs(t,x,u) sys = []; % sys 表示下输出，即y。计算下一个采样时间，仅在系统是变采样时间系统时调用function sys=mdlGetTimeOfNextVarHit(t,x,u) sampleTime = 1; %设置下一次的采样时间是1s以后sys = t + sampleTime; %?sys 表示下一个采样时间点。?%仿真结束时要调用的例程函数，在仿真结束时调用，用户可以在此完成结束仿真所需的必要工作 * 6.2 S