matlab中S函数编辑流程
- 格式:docx
- 大小:538.92 KB
- 文档页数:10
matlab中s-function用法创建一个有1输入(2维),2输出(1维),3个参数,还有全局变量的s-function。
1.新建sfunction的C语言文件打开simulink,点击User-Defined Functions里面的S-Function Examples。
这个里面有多个语言版本的模板,有C,C++,Ada,Fortran和M语言的版本,其实都大同小异,只要了解几个函数就很容易使用了。
选择C语言的版本:从S-function模块中选择C-file S-functions里面的Basic C-MEX template。
打开后,另存为自己的模块名字,如test.c。
下面我们来分析代码:#define S_FUNCTION_NAME test//这里把文件名sfuntmpl_basic修改为test #define S_FUNCTION_LEVEL 2#include "simstruc.h"//程序里面要用到的头文件在这里引用,如“math.h”等。
float global_var; //定义全局变量static void mdlInitializeSizes(SimStruct *S){//这个函数用来设置输入、输出和参数的。
ssSetNumSFcnParams(S, 3); /*设置参数个数,这里为3 */if (ssGetNumSFcnParams(S) != ssGetSFcnParamsCount(S)) {return;}; ssSetNumContStates(S, 0);//设置连续状态的个数,缺省为0ssSetNumDiscStates(S, 0);//设置离散状态的个数,缺省为0;if (!ssSetNumInputPorts(S, 1)) return;//设置输入变量的个数,这里为1 ssSetInputPortWidth(S, 0, 2); //设置输入变量0的维数为2ssSetInputPortRequiredContiguous(S, 0, true); //设置input0的访问方式,true就是临近访问,这样指针的增量后就可以直接访问下个input端口了。
matlab中S函数编写matlab中S函数编写S函数是system Function的简称,用它来写自己的simulink模块,可以用matlab、C、C++、Fortran、Ada等语言来写。
C或C++等语言写的s函数还可以实现对硬件端口的操作,还可以操作windows API等。
Simulink的仿真有两个阶段:一个为初始化,这个阶段主要是设置一些参数,像系统的输入输出个数、状态初值、采样时间等;第二个阶段就是运行阶段,这个阶段里要进行计算输出、更新离散状态、计算连续状态等等,这个阶段需要反复运行,直至结束。
在matlab的workspace里打edit sfuntmpl(这是matlab提供的s函数模板),具体分析S函数的结构。
它的第一行:function [sys,x0,str,ts,simStateCompliance]=sfuntmpl(t,x,u,flag) 输入与输出变量的含义:t是采样时间,x是状态变量,u是输入,flag是仿真过程中的状态标志(以它来判断当前是初始化还是运行等);sys输出根据flag的不同而不同,x0是状态变量的初始值,str是保留参数(一般在初始化中将它置空就可以了,str=[]),ts是一个1×2的向量,ts(1)是采样周期,ts(2)是偏移量。
下面结合sfuntmpl.m中的代码来讲具体的结构:switch flag, %判断flag,看当前处于哪个状态case 0,[sys,x0,str,ts]=mdlInitializeSizes;flag=0表示处于初始化状态,sys是一个结构体,用它来设置模块的一些参数,各个参数详细说明如下:size = simsizes; %用于设置模块参数的结构体用simsizes来生成sizes.NumContStates = 0; %模块连续状态变量的个数sizes.NumDiscStates = 0; %模块离散状态变量的个数sizes.NumOutputs = 0; %模块输出变量的个数sizes.NumInputs = 0; %模块输入变量的个数sizes.DirFeedthrough = 1; %模块是否存在直接贯通(输入能直接控制输出)sizes.NumSampleTimes = 1;%模块的采样时间个数,至少是一个sys = simsizes(sizes); %设置完后赋给sys输出举个例子,考虑如下模型:dx/dt=fc(t,x,u) 也可以用连续状态方程描述:dx/dt=A*x+B*u x(k+1)=fd(t,x,u) 也可以用离散状态方程描述:x(k+1)=H*x(k)+G*u(k)y=fo(t,x,u) 也可以用输出状态方程描述:y=C*x+D*u设上述模型连续状态变量、离散状态变量、输入变量、输出变量均为1个,我们就只需改上面那一段代码为:sizes.NumContStates=1;sizes.NumDiscStates=1;sizes.NumOutputs=1;sizes.NumInputs=1;其他的可以不变。
matlab中s函数模板-回复MATLAB中的s函数模板是一种用于编写自定义模块的工具,它允许用户以自己的方式定义系统的动态行为。
本文将详细介绍s函数模板的具体操作步骤和相关注意事项。
首先,为了能够使用s函数模板,我们需要在MATLAB工作区创建一个新的模型或打开一个已有的模型。
然后,我们可以按照以下步骤进行s函数模板的编写。
第一步,创建s函数模板文件。
在MATLAB的命令窗口中输入以下命令:sfunctionname这将会创建一个名为"sfunctionname.m"的文件,其中"sfunctionname"是我们给s函数模板命名的名称。
第二步,打开s函数模板文件,并编辑其内容。
在"function"关键字下方的行中,可以定义输入输出参数以及任何需要的局部变量。
例如,以下是一个基本的s函数模板:function [sys,x0,str,ts] = sfunctionname(t,x,u,flag)switch flagcase 0 初始化模型[sys,x0,str,ts] = mdlInitializeSizes; 调用模型初始化函数case 1 计算模型状态方程sys = mdlDerivatives(t,x,u);case 2 计算模型输出方程sys = mdlOutputs(t,x,u);case 3 处理模型状态信息sys = mdlUpdate(t,x,u);case 4 处理模型输出信息sys = mdlGetTimeOfNextVarHit(t,x,u);case 9 处理模型示意图sys = mdlTerminate(t,x,u);otherwiseerror(['Unhandled flag = ',num2str(flag)]);endend在上述代码中,我们定义了四个输入参数t、x、u和flag,以及四个输出参数sys、x0、str和ts。
matlab中s函数编写心得(转)Part I:所谓s函数是system Function的简称, 用它来写自己的simulink模块. s函数可以用matlab、C、C++、Fortran、Ada等语言来写,这儿我只介绍怎样用matlab语言来写吧(主要是它比较简单)< xmlnamespace prefix ="o" ns="urn:schemas-microsoft-com:office:office" />先讲讲为什么要用s函数,我觉得用s函数可以利用matlab的丰富资源,而不仅仅局限于simulink提供的模块,而用c或c++等语言写的s函数还可以实现对硬件端口的操作,还可以操作windows API等先介绍一下simulink的仿真过程(以便理解s函数),simulink的仿真有两个阶段:一个为初始化,这个阶段主要是设置一些参数,像系统的输入输出个数、状态初值、采样时间等;第二个阶段就是运行阶段,这个阶段里要进行计算输出、更新离散状态、计算连续状态等等,这个阶段需要反复运行,直至结束.在matlab的workspace里输入edit sfuntmpl(这是matlab自己提供的s函数模板),我们看它来具体分析s函数的结构.1. 函数的函数头函数的第一行:function [sys,x0,str,ts]=sfuntmpl(t,x,u,flag) , 先讲输入与输出变量的含义:t是采样时间, x是状态变量, u是输入(是做成simulink模块的输入), flag是仿真过程中的状态标志(以它来判断当前是初始化还是运行等) sys输出根据flag的不同而不同(下面将结合flag来讲sys的含义), x0是状态变量的初始值, str是保留参数(mathworks公司还没想好该怎么用它, 一般在初始化中将它置空就可以了, str=[]), ts是一个1×2的向量, ts(1)是采样周期, ts(2)是偏移量2. 函数分析下面结合sfuntmpl.m中的代码来讲具体的结构:switch flag, %判断flag,看当前处于哪个状态case 0,[sys,x0,str,ts]=mdlInitializeSizes;// 解释说明flag=0表示当前处于初始化状态,此时调用函数mdlInitializeSizes进行初始化,此函数在该文件的第149行定义. 其中的参数sys是一个结构体,它用来设置模块的一些参数,各个参数详细说明如下size = simsizes;%用于设置模块参数的结构体用simsizes来生成sizes.NumContStates = 0; %模块连续状态变量的个数sizes.NumDiscStates = 0; %模块离散状态变量的个数sizes.NumOutputs = 0; %模块输出变量的个数sizes.NumInputs = 0; %模块输入变量的个数sizes.DirFeedthrough = 1; %模块是否存在直接贯通sizes.NumSampleTimes = 1; %模块的采样时间个数, 至少是一个sys = simsizes(sizes); %设置完后赋给sys输出举个例子,考虑如下模型:dx/dt=fc(t,x,u) 也可以用连续状态方程描述:dx/dt=A*x+B*u x(k+1)=fd(t,x,u) 也可以用离散状态方程描述:x(k+1)=H*x(k)+G*u(k)y=fo(t,x,u) 也可以用输出状态方程描述:y=C*x+D*u设上述模型连续状态变量、离散状态变量、输入变量、输出变量均为1个,我们就只需改上面那一段代码为(一般连续状态与离散状态不会一块用, 我这儿是为了方便说明):sizes.NumContStates=1;sizes.NumDiscStates=1;sizes.Nu mOutputs=1;sizes.NumInputs=1;其他的可以不变, 继续在mdlInitializeSizes函数中往下看:x0 = []; %状态变量设置为空,表示没有状态变量,以我们上面的假设,可改为x0=[0,0](离散和连续的状态变量我们都设它初值为0)str = []; %保留参数, 置[]就可以了, 没什么用ts = [0 0]; %采样周期设为0表示是连续系统, 如果是离散系统在下面的mdlGetTimeOfNextVarHit函数中具体介绍case 1,sys=mdlDerivatives(t,x,u);//flag=1表示此时要计算连续状态的微分, 即上面提到的dx/dt=fc(t,x,u)中的dx/dt, 找到193行的函数mdlDerivatives, 如果设置连续状态变量个数为0, 此处只需sys=[]就可以了, 按我们上述讨论的那个模型, 此处改成sys=fc(t,x(1),u)或sys=A*x(1)+B*u, 我们这儿x(1)是连续状态变量, 而x(2)是离散的, 这儿只用到连续的, 此时的输出sys就是微分case 2,sys=mdlUpdate(t,x,u);//flag=2表示此时要计算下一个离散状态, 即上面提到的x(k+1)=fd(t,x,u), 找到mdlUpdate函数, 它这儿sys=[]表示没有离散状态, 我们这儿可以改成sys=fd(t,x(2),u)或sys=H*x(2)+G*u;%sys即为x(k+1)case 3,sys=mdlOutputs(t,x,u);//flag=3表示此时要计算输出, 即y=fo(t,x,u), 找到218行的mdlOutputs函数. 如果sys=[]表示没有输出, 我们改成sys=fo(t,x,u)或sys=C*x+D*u %sys此时为输出ycase 4,sys=mdlGetTimeOfNextVarHit(t,x,u);//flag=4表示此时要计算下一次采样的时间, 只在离散采样系统中有用(即上文的mdlInitializeSizes中提到的ts设置ts(1)不为0), 连续系统中只需在mdlGetTimeOfNextVarHit函数中写上sys=[]. 这个函数主要用于变步长的设置, 具体实现大家可以用edit vsfunc看vsfunc.m这个例子case 9,sys=mdlTerminate(t,x,u);//flag=9表示此时系统要结束,一般来说写上在mdlTerminate函数中写上sys=[]就可, 如果你在结束时还要设置什么,就在此函数中写完了.Part II此外, s函数还可以带用户参数, 下面给个例子, 它和simulink下的gain模块功能一样function [sys,x0,str,ts] = sfungain(t,x,u,flag,gain)switch flag,case 0,sizes = simsizes;sizes.NumContStates = 0;sizes.NumDiscStates = 0;sizes.NumOutputs = 1;sizes.NumInputs = 1;sizes.DirFeedthrough = 1;sizes.NumSampleTimes = 1;sys = simsizes(sizes);x0=[];str=[];ts=[0,0];case 3,sys=gain*u;case {1,2,4,9},sys = [];end做好了s函数后, simulink --> user-defined function下拖一个S-Function到你的模型, 就可以用了. 在simulink -->user-defined function还有个s-Function Builder, 他可以生成用c 语言写的s函数. 或者在matlab的workspace下打sfundemos, 可以看到很多演示s函数的程序Part IIISIMULINK s-function的设计Simulink为用户提供了许多内置的基本库模块, 通过这些模块进行连接而构成系统的模型. 对于那些经常使用的模块进行组合并封装可以构建出重复使用的新模块, 但它依然是基于Simulink原来提供的内置模块.而Simulink s-function是一种强大的对模块库进行扩展的新工具.(一) s-function的概念s-function是一个动态系统的计算机语言描述, 在MATLAB里, 用户可以选择用m文件编写, 也可以用c或mex文件编写, 在这里只给大家介绍如何用m文件编写s-function.S-function提供了扩展Simulink模块库的有力工具, 它采用一种特定的调用语法, 使函数和Simulink解法器进行交互.S-function最广泛的用途是定制用户自己的Simulink模块. 它的形式十分通用, 能够支持连续系统、离散系统和混合系统.(二) 建立m文件s-function1. 使用模板文件:sfuntmp1.m, 其格式为[sys,x0]=function(t,x,u,flag). 该该模板文件位于MATLAB根目录下toolbox/simulink/blocks目录下模板文件里s-function的结构十分简单, 它只为不同的flag的值指定要相应调用的m文件子函数. 比如当flag=3时, 即模块处于计算输出这个仿真阶段时, 相应调用的子函数为sys=mdloutputs(t,x,u). 模板文件使用switch语句来完成这种指定, 当然这种结构并不唯一, 用户也可以使用if语句来完成同样的功能. 而且在实际运用时, 可以根据实际需要来去掉某些值, 因为并不是每个模块都需要经过所有的子函数调用.模板文件只是Simulink为方便用户而提供的一种参考格式, 并不是编写s-function的语法要求, 用户完全可以改变子函数的名称, 或者直接把代码写在主函数里, 但使用模板文件的好处是, 比较方便, 而且条理清晰.使用模板编写s-function, 用户只需把s-函数名换成期望的函数名称, 如果需要额外的输入参量, 还需在输入参数列表的后面增加这些参数, 因为前面的4个参数是simulink调用s-function时自动传入的.对于输出参数, 最好不做修改. 接下去的工作就是根据所编s-function要完成的任务, 用相应的代码去替代模板里各个子函数的代码即可.Simulink在每个仿真阶段都会对s-function进行调用. 在调用时, Simulink会根据所处的仿真阶段为flag传入不同的值, 而且还会为sys这个返回参数指定不同的角色. 也就是说尽管是相同的sys变量, 但在不同的仿真阶段其意义却不相同, 这种变化由simulink自动完成.m文件s-function可用的子函数说明如下:mdlInitializeSizes(flag=0) -- 定义s-function模块的基本特性, 包括采样时间、连续或者离散状态的初始条件和sizes数组mdlDerivatives(flag=1) -- 计算连续状态变量的微分方程mdlUpdate(flag=2) -- 更新离散状态、采样时间和主时间步的要求mdlOutputs(flag=3) -- 计算s-function的输出mdlGetTimeOfNextVarHit(flag=4) -- 计算下一个采样点的绝对时间, 这个方法仅仅是在用户在mdlInitializeSizes 里说明了一个可变的离散采样时间概括说来, 建立s-function可以分成两个分离的任务: 第一, 初始化模块特性包括输入输出信号的宽度, 离散连续状态的初始条件和采样时间. 第二, 将算法放到合适的s-function子函数中去。
S函数的简介及编写摘自恒润科技S-function的编写1. S函数模板编辑环境进入:在MATLAB主界面中直接输入:edit sfuntmpl即可弹出S函数模板编辑的M文件环境,修改即可。
在MATLAB主界面中直接输入:sfundemos,即可调出S 函数的许多编程例子。
2. S函数模板的相关基础:1)M文件S函数的引导语句为:xtflagfuFunction[psysxstrtsp,,12,...),,,0,],(,S函数默认的四个输入参数:t ,x ,u ,flagS函数默认的四个输出函数:sys ,x0 ,str ,ts各个参数的含义如下:T :代表当前的仿真时间,该输入决定了下一个采样时间;X :表示状态向量,行向量,引用格式:X(1),X(2)U :表示输入向量;Flag :控制在每一个仿真阶段调用哪一个子函数的参数,由SIMULINK在调用时自动取值;Sys :通用的返回变量,返回的数值决定Flag值,mdlUpdates里:列向量,引用格式:Sys(1,1),Sys(2,1);mdlOutputs里:行向量,引用格式:Sys =x、X0 :初始的状态值;列向量,引用格式:X0=[ 0;0;0 ]Str :空矩阵,无具体含义;Ts :包含模块采样时间与偏差的矩阵。
[period, offset]当Ts为-1时,表示与输入信号同采样周期。
2)S函数工作方式:Flag = 0时,调用mdlInitializeSizes函数,定义S函数的基本特性,包括采样时间,连续或者离散状态的初始条件与Sizes数组;Flag = 1时,调用mdlDerivatives函数,计算连续状态变量的微分方程;求所给表达式的等号左边状态变量的积分值的过程。
Flag = 2时,调用mdlUpdate函数,用于更新离散状态,采样时间与主时间步的要求;Flag = 3时,调用mdlOutputs函数,计算S函数的输出;Flag = 4时,调用mdlGetTimeOfNextVarHit函数,计算下一个采样点的绝对时间,这个方法仅仅就是使用户在mdlInitializeSize 里说明一个可变的离散采样时间;Flag = 9时,调用mdlTerminate函数,实现仿真任务的结束。
sfunction使用方法在matlab命令行中输入 edit sfuntmpl ,可以打开S-function的模板,里面对S-function的函数和变量作用均作了介绍。
S-function的基础是现代控制理论的状态定义,使用S-function需具备现代控制理论的基本知识为什么要有S-function•因为matlab中的模块太过标准化,无法满足多变的应用和调试需求,在以下情况下,S-function可以很方便的对系统进行建模。
•控制或受控系统具有复杂时变特性•控制系统在控制过程中自我迭代,存在需要保留和记录的变量•控制系统的输出较为复杂,难以用function和m-function表示•系统在生成、运行或结束时需要执行特殊操作S-function属性窗口介绍•S-function name :S-function的名字,随便写,自己认识即可•S-function parameters:S-function的模块参数,默认为空•S-function modules :S-function的模块,无需修改,采用系统默认模块即可也就是说,这三个参数,只需要修改第一个参数为模块命名即可•点击Edit,可以进入S-function的代码编辑界面S-function内部函数介绍•[sys,x0,str,ts] = functionName(t,x,u,flag) 这是函数的总入口,收到信号后,首先进入这个函数.这个函数包含一个switch语句,根据情况进入不同的子函数•[sys,x0,str,ts,simStateCompliance]=mdlInitializeS izes; S-function进行基本的设置,具体参数会在后面进行介绍。
相当于构造函数•sys=mdlDerivatives(t,x,u); 该函数仅在连续系统中被调用,用于产生控制系统状态的导数•sys=mdlUpdate(t,x,u); 该函数仅在离散系统中被调用,用于产生控制系统的下一个状态•sys=mdlOutputs(t,x,u); 产生(传递)系统输出•sys=mdlGetTimeOfNextVarHit(t,x,u) 获得下一次系统执行(next hit)的时间,该时间为绝对时间此函数仅在采样时间数组中指定变量离散时间采样时间[-2 0]时会被调用。
matlab的sfunction执行流程MATLAB的S-Function(System Function)是一种用于扩展Simulink模型功能的强大工具。
S-Function允许用户以C、C++或MATLAB语言编写自定义的动态系统。
以下是S-Function的基本执行流程:1.初始化(Initialization):•当Simulink模型开始仿真时,S-Function的初始化函数(如mdlInitializeSizes)首先被调用。
在这个阶段,S-Function定义了它的输入和输出数量、离散状态的数量、采样时间等。
•如果S-Function使用工作向量(work vector),它也会在这个阶段被分配。
2.计算采样时间(Compute Sample Times):•在初始化之后,Simulink会调用S-Function的mdlGetSampleTimes 函数来确定其采样时间。
这对于确定何时调用S-Function的其他函数至关重要。
3.开始仿真(Start of Simulation):•在仿真开始时,可以调用一个可选的S-Function函数(如mdlStart),用于执行任何需要在仿真开始时执行的初始化或预计算。
4.仿真循环(Simulation Loop):•在仿真的每个采样点,S-Function的主要功能函数(如mdlOutputs,mdlUpdate,mdlDerivatives)被调用。
▪mdlOutputs:计算S-Function的输出。
▪mdlUpdate:更新S-Function的离散状态。
▪mdlDerivatives:计算S-Function的连续状态导数(对于连续时间系统)。
•根据S-Function的类型(连续、离散或混合)和配置,这些函数可能会被以不同的方式调用。
5.结束仿真(End of Simulation):•在仿真结束时,可以调用一个可选的S-Function函数(如mdlTerminate),用于执行任何需要在仿真结束时执行的清理操作。
S函数编辑流程及相关1. S函数模板编辑环境进入:在MATLAB主界面中直接输入:edit sfuntmpl即可弹出S函数模板编辑的M文件环境,修改即可。
在MATLAB主界面中直接输入:sfundemos,即可调出S函数的许多编程例子。
2. S函数模板的相关基础:1)M文件S函数的引导语句为:txfpusysFunction[pxstrtsflag,,12,...),,,0,],(,S函数默认的四个输入参数:t ,x ,u ,flagS函数默认的四个输出函数:sys ,x0 ,str ,ts各个参数的含义如下:T :代表当前的仿真时间,该输入决定了下一个采样时间;X :表示状态向量,行向量,引用格式:X(1),X(2)U :表示输入向量;Flag :控制在每一个仿真阶段调用哪一个子函数的参数,由SIMULINK在调用时自动取值;Sys :通用的返回变量,返回的数值决定Flag值,mdlUpdates 里:列向量,引用格式:Sys(1,1),Sys(2,1);mdlOutputs里:行向量,引用格式:Sys =x.X0 :初始的状态值;列向量,引用格式:X0=[ 0;0;0 ]Str :空矩阵,无具体含义;Ts :包含模块采样时间和偏差的矩阵。
[period, offset]当Ts为-1时,表示与输入信号同采样周期。
2)S函数工作方式:Flag = 0时,调用mdlInitializeSizes函数,定义S函数的基本特性,包括采样时间,连续或者离散状态的初始条件和Sizes 数组;Flag = 1时,调用mdlDerivatives函数,计算连续状态变量的微分方程;求所给表达式的等号左边状态变量的积分值的过程。
Flag = 2时,调用mdlUpdate函数,用于更新离散状态,采样时间和主时间步的要求;Flag = 3时,调用mdlOutputs函数,计算S函数的输出;Flag = 4时,调用mdlGetTimeOfNextVarHit函数,计算下一个采样点的绝对时间,这个方法仅仅是使用户在mdlInitializeSize 里说明一个可变的离散采样时间;Flag = 9时,调用mdlTerminate函数,实现仿真任务的结束。
3)S函数仿真过程:①初始化:mdlInitializeSizes,初始化S函数●初始化SimStruct,包含了S函数的所有信息;●设置输入、输出端口数;●设置采样时间;●分配存储空间。
②数值积分:mdlDerivatives●用于连续状态的求解和非采样过零点;●如果存在连续状态,调用mdlDerivatives和mdlOutput两个子函数;●如果存在非采样过零点,调用mdlOutput和mdlZeroCrossings子函数,以定位过零点。
③更新离散状态:mdlUpdate④计算输出:mdlOutputs,计算所有输出端口的输出值。
⑤计算下一个采样时间点:mdlGetTimeOfNextVarHit⑥仿真结束:mdlTerminate,在仿真结束时调用。
3. S函数的编写:1)参数初始设定:初始化sizes结构,再调用simsizes函数;Sizes结构体:NumContStates:连续状态的个数NumDiscStates:离散状态的个数NumOutputs:输出变量的个数NumInputs:输入变量的个数DirFeedthrough:有无直接馈入,值为1时表示输入直接传到输出口NumSampleTimes:采样时间的个数,值为1时表示只有一个采样周期Simsizes函数的调用:sys = simsizes (sizes) ,即将sizes结构体中的信息传递给sys。
2)状态的动态更新:连续模块的状态更新由mdlDerivatives函数来进行;离散模块的状态更新由mdlUpdate函数来进行;3)输出信号的计算:计算出模块的输出信号,系统的输出仍然由sys变量返回。
4. M文件S函数的模块化:在动态系统仿真设计,分析中,用户可以使用S-Function 模块来调用S-函数。
1)S-Function模块是一个单输入单输出的模块,如果有多个输入与输出信号,可以使用Mux模块与Demux模块对信号进行组合和分离操作;2)在S-Function模块的参数设置对话框中,包含了调用的S函数名和用户输入的参数列表,如下图所示:3)S-Function模块是以图形的方式提供给用户一个调用S函数的接口,S函数中的源文件必须由用户自行编写;4)S-Function模块中的S-函数名和参数值列表必须与用户填写的S函数源文件的名称和参数列表完全一致,包括参数的顺序。
5. S-Function模块建立流程:创建一个有1输入(2维),2输出(1维),3个参数,还有全局变量的s-function。
1.新建sfunction的C语言文件打开simulink,点击User-Defined Functions里面的S-Function Examples。
这个里面有多个语言版本的模板,有C,C++,Ada,Fortran和M语言的版本,其实都大同小异,只要了解几个函数就很容易使用了。
选择C语言的版本:从S-function模块中选择C-file S-functions里面的Basic C-MEX template。
打开后,另存为自己的模块名字,如test.c。
下面我们来分析代码:#define S_FUNCTION_NAME test//这里把文件名sfuntmpl_basic修改为test#define S_FUNCTION_LEVEL 2#include "simstruc.h"//程序里面要用到的头文件在这里引用,如“math.h”等。
float global_var; //定义全局变量static void mdlInitializeSizes(SimStruct *S){//这个函数用来设置输入、输出和参数的。
ssSetNumSFcnParams(S, 3); /*设置参数个数,这里为3 */if (ssGetNumSFcnParams(S) != ssGetSFcnParamsCount(S)) {return;}ssSetNumContStates(S, 0);//设置连续状态的个数,缺省为0;ssSetNumDiscStates(S, 0);//设置离散状态的个数,缺省为0;if (!ssSetNumInputPorts(S, 1)) return;//设置输入变量的个数,这里为1ssSetInputPortWidth(S, 0, 2); //设置输入变量0的维数为2 ssSetInputPortRequiredContiguous(S, 0, true); //设置input0的访问方式,true就是临近访问,这样指针的增量后就可以直接访问下个input端口了。
ssSetInputPortDirectFeedThrough(S, 0, 1);// 设置输入端口的信号是否mdlOutputs 函数中使用,这儿设置为true。
if (!ssSetNumOutputPorts(S, 2)) return;//设置输出变量的个数ssSetOutputPortWidth(S, 0, 1);//设置输出变量0的维数为1维ssSetOutputPortWidth(S, 1, 1);//设置输出变量1的维数为1维ssSetNumSampleTimes(S, 1); //设置采样时间,此处为1s。
ssSetNumRWork(S, 0);//不管ssSetNumIWork(S, 0);ssSetNumPWork(S, 0);ssSetNumModes(S, 0);ssSetNumNonsampledZCs(S, 0);ssSetOptions(S, 0);//下面可以写全局变量的初始化程序global_var=1;}static void mdlInitializeSampleTimes(SimStruct *S)//暂时不管{ssSetSampleTime(S, 0, CONTINUOUS_SAMPLE_TIME);ssSetOffsetTime(S, 0, 0.0);}#define MDL_INITIALIZE_CONDITIONS /* Change to #undef to remove function */#if defined(MDL_INITIALIZE_CONDITIONS)static void mdlInitializeConditions(SimStruct *S)//暂时不管{}#endif /* MDL_INITIALIZE_CONDITIONS */#define MDL_START /* Change to #undef to remove function */#if defined(MDL_START)static void mdlStart(SimStruct *S)//暂时不管{}#endif /* MDL_START */static void mdlOutputs(SimStruct *S, int_T tid)//这里填入相关的运算、算法等{real_T *para1 = mxGetPr(ssGetSFcnParam(S,0));real_T *para2 = mxGetPr(ssGetSFcnParam(S,1));real_T *para3 = mxGetPr(ssGetSFcnParam(S,2));const real_T *u = (const real_T*) ssGetInputPortSignal(S,0);real_T *y1 = ssGetOutputPortSignal(S,0);real_T *y2 = ssGetOutputPortSignal(S,1);y1[0]=u[0]*para1[0]+u[1]*para2[0];y2[0]=u[1]*para3[0]+u[0]*para1[0];}#define MDL_UPDATE /* Change to #undef to remove function */#if defined(MDL_UPDATE)static void mdlUpdate(SimStruct *S, int_T tid){}#endif /* MDL_UPDATE */#define MDL_DERIV ATIVES /* Change to #undef to remove function */#if defined(MDL_DERIV ATIVES)static void mdlDerivatives(SimStruct *S){}#endif /* MDL_DERIV ATIVES */static void mdlTerminate(SimStruct *S)//这里需要把global变量全部初始化,否则下次运行程序时,全局变量还是之前的值。