MATLABS函数编写示范.doc

MATLAB S函数编写示范（MATLAB2010a）?? 2010-05-05 21:22:47|??分类： MATLAB之控制系统 |??标签：matlab??s-function??编写??教程??s??builder?? |字号大中小?订阅 1.M语言S函数 ????? 在命令行输入命令“edit sfuntmpl”即可打开M语言s函数模板。M语言S函数是以flag为标志进行运算的。flag==0时，初始化参数；flag==1时，计算微分；flag==2时，状态更新；flag==3时，计算输出；flag==4时适用于变步长结算；其他参数不常用。 ??????（1）flag==0时，调用mdlInitializeSizes函数，sizes结构及解释如下： ?????????????????????? sizes.NumContStates? = 0;???//连续变量个数?????????????????????? sizes.NumDiscStates? = 0;?? //离散变量个数?????????????????????? sizes.NumOutputs???? = 0;????//输出个数?????????????????????? sizes.NumInputs????? = 0;??? //输入个数?????????????????????? sizes.DirFeedthrough = 1;? //是否直通，当在输出里包含了输入时，为直通?????????????????????? sizes.NumSampleTimes = 1;? //采样时间 ?????????????????????? 对状态变量进行初始化,如果初始条件是0初始条件，则 ?????????????????????? x0=[0 0] ?????? （2）flag==1时，调用mdlDerivatives(t,x,u)函数计算微分，如有一系统为状态方程为x=2x+u，y=5x+u时，编写如下：sys=2*x+u; ?????? （3）flag==3时，计算输出，如上例系统，则编写如下：sys=5*x+u; ???????? ok,下面只需要将编写好的s函数名称写入s-function模块。 ???????? 例子：对如下系统（零初始条件） ????????????????x1=-3x1-x2+u; ??????????????? x2=2x1; ??????????????? y=x2;编写s函数如下： ?function [sys,x0,str,ts,simStateCompliance] = sfun_m_hyj(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));endfunction [sys,x0,str,ts,simStateCompliance]=mdlInitializeSizessizes = simsizes; sizes.NumContStates? = 2;sizes.NumDiscStates? = 0;sizes.NumOutputs???? = 1;sizes.NumInputs????? = 1;sizes.DirFeedthrough = 1;sizes.NumSampleTimes = 1;?? % at least one sample time is neededsys = simsizes(sizes);x0? = [0 0];str = [];ts? = [0 0];simStateCompliance = UnknownSimState;function sys=mdlDerivatives(t,x,u)sys = [-3*x(1)-x(2)+u;2*x(1)]; //注意此处传递的是x向量，包含了x1和x2，引用时使用x（1），x（2），下标从1开始funct