1. S-Function简介

S-Function是system-function的缩写。说得简单,S-Function就是用MATLAB所提供的模型不能完全满足用户,而提供给用户自己编写程序来满足自己要求模型的接口。

2. MEX函数与M文件的区别

第一,  MEX 函数能实现的回调函数比M-文件能实现的回调函数要多得多;

第二,  MEX 函数直接访问内部数据结构SimStruct,SimStruct 是Simulink 用来保存关于S-function 信息的一个数据结构;

第三,  MEX 函数也可使用MATLAB MEX 文件API 直接来访问MATLAB 的工作空间。

如果一个C MEX文件与一个M文件具有相同的名字,则C MEX文件被优先使用,即在S-Function块中使用的是C MEX文件。

3. 基础知识

3.1 直接馈通(direct feedthrough)

直接馈通表示系统的输出或可变采样时间是否受到输入的控制。

a.  输出函数(mdlOutputs或flag==3)是输入u的函数。即,如果输入u在mdlOutputs中被访问,则存在直接馈通。

b.  对于一个变步长S-Function的“下一个采样时间”函数(mdlGetTimeOfNextVarHit或flag==4)中可以访问输入u。

例如,一个需要其输入的系统(也就是具有直接馈通)是运算y=kXu,其中,u是输入,k是增益,y是输出。

又如,一个不需要其输入的系统(也就是没有直馈通)是一种简单的积分运算:

输出:y=x;

导数:dx/dt=u

其中,x是状态,dx/dt是状态对时间的导数,u是输入,y是输出。

正确设置直接馈通标志是十分重要的,因为它影响模型中块的执行顺序,并可用检测代数环。

3.2 dynamically sized inputs

主要是给出:输入连续状态数目(size.NumContStates),离散状态数目(size.NumDiscStates) ,输出数目(size.NumOutputs),输入数目(size.NumInputs),Direct Feedthrough(size.Dir Feedthrough)。

3.3 setting sample times and offsets

setting smaple times and offsets主要设置采样时间.

3.4 Level-1 和Level-2

Level 1 提供一个简单的接口,可与少部分的S函数API交互。Matlab对于这种方式的支持更多的是为了保持与以前版本的兼容,现在推荐采用的是Level 2 S函数。

4. S-Function实例

S-Function的仿真流程

例如要创建一个有1输入(2维),2输出(1维),3个参数,还有全局变量的S-Function。 过程如下:

a. 新建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_DERIVATIVES  /* Change to #undef to remove function */
#if defined(MDL_DERIVATIVES)static void mdlDerivatives(SimStruct *S){}
#endif /* MDL_DERIVATIVES */
static void mdlTerminate(SimStruct *S)//这里需要把global变量全部初始化,否则下次运行程序时,全局变量还是之前的值。
{
}#ifdef  MATLAB_MEX_FILE    /* Is this file being compiled as a MEX-file? */
#include "simulink.c"      /* MEX-file interface mechanism */
#else
#include "cg_sfun.h"       /* Code generation registration function */
#endif

b. 编译

在matlab的command window 里面输入“mex test.c”,即可将test.c编译为mex文件。

c.调用sfunction

在simulink空间里面拉入sfunction,在s-function name里面填入test,参数里面填入要设定的参数,然后仿真即可。

//

转载:http://www.cnblogs.com/xpvincent/archive/2013/01/18/2866978.html

MATLAB中的S-Function的用法(C语言)相关推荐

  1. matlab fopen函数的用法,matlab中关于fread函数的用法

    公告: 为响应国家净网行动,部分内容已经删除,感谢读者理解. 话题:matlab中关于fread函数的用法 问题详情:用fread打开一个时,需要不需要把开放在某个特回答:必须在Matlab可以搜索的 ...

  2. Matlab中sub2ind函数原理及用法!

    Matlab中sub2ind函数原理及用法!

  3. norminv函数是什么matlab,matlab中的函数norminv能否用C语言实现

    matlab中的函数norminv可不可以用C语言实现? 对于matlab中的函数norminv,C语言中有没有函数可以完成类似的功能?没有的话,用C语言要如何实现呢? matlab中的函数normi ...

  4. matlab中uint16(uint8)函数的用法

    matlab中uint16函数的用法: 1.函数用法:将输入参数映射到无符号16位范围之内,即将参数限定在0-65535之间: 2.实例:当输入为一个实数a时: if a<0 输出ans=0; ...

  5. matlab中find()函数用法

    一.基本用法 返回矩阵或向量中非零元素的索引 注意:matlab中下标从1开始 举例: (1)向量 返回非零元素下标 find(vector) x=[1 2 3 0 0 6 7 8 9]; find( ...

  6. matlab中Svmtrain和Svmpredict的用法

    一. svmtrain 的用法 svmtrain实现对训练数据集的训练,获得SVM模型. 用法: svmtrain [options] training_set_file [model_file] 其 ...

  7. matlab中libsvm的svmtrain函数用法

    这里的LIBSVM是一个由台湾大学林智仁(Lin Chih-Jen)教授等开发的SVM模式识别与回归的软件包,使用简单,功能强大,本文主要介绍其在Matlab中的使用. 注意不是matlab自带的sv ...

  8. matlab fopen函数的用法,matlab中fopen函数与fprintf用法

    matlab中fopen函数在指定文件打开的实例如下: *1)"fopen"打开文件,赋予文件代号. 语法1:FID= FOPEN(filename,permission) 用指定 ...

  9. matlab中signal在哪里,matlab中Signal Builder模块的用法总结

    在网上搜了好久,关于matlab中Signal Builder模块的使用方法很少,由于做仿真时需要用到这个模块,无奈自己摸索过后在此总结下使用方法: 双击Signal Builder模块弹出如下对话框 ...

  10. matlab中的sum函数的用法

    Matlab中的sum函数功能非常强大,但是如何使用不慎就会造成难以发现的错误.比如 A1 = [1, 3, 2] S = sum(A1) S_2 = sum(A1,2) S_1 = sum(A1,1 ...

最新文章

  1. 使用 Elixir 开发嵌入式系统: 使用DS18B20传感器测量体温
  2. spring boot打包文件后,报错\No such file or directory
  3. VCL组件之TLabel、TStaticText和TLabeledEdit
  4. Java 基础 - 面向对象(不错N多教程集合)
  5. 关于iOS 3D touch 指纹验证的随笔
  6. micropython 串口写文件_MicroPython通过2G模块串口收发短信
  7. [react] react怎么提高列表渲染的性能?
  8. 矩阵:如何使用矩阵操作进行 PageRank 计算?
  9. html给主题设置背景色,css如何设置元素的背景色?
  10. 飞思卡尔单片机DZ60---TPM1\TPM2溢出中断
  11. 未来教育MySQL题库下载_未来教育考试系统(通用下载管理版)可在线下载最新试题,支持章节练习...
  12. 【网络安全】实操XSS订单系统漏洞(利用盲打)
  13. Base64-encoded key bytes may only be specified for HMAC signatures.
  14. Win10切换到了Users用户怎么切换回来
  15. 修改Element-ui表格样式
  16. kendoUI模板概述(template)
  17. response.text和response.content
  18. 【hpuoj】OY问题
  19. 简易超声波雷达的arduino实现
  20. 高效运维的本质:可视化的服务交付和可视化的服务度量

热门文章

  1. matlab和python哪个运行快_MATLAB比Python快吗?
  2. 扔掉老破V100、A100,英伟达新一代计算卡H100来了!
  3. 谈谈怎样提高炼丹手速
  4. 美团点评移动网络优化实践
  5. ICDAR 2019论文:自然场景文字定位技术详解
  6. 美团客户端响应式框架 EasyReact 开源啦
  7. 想成为阿里160万年薪的P8架构师?你必须掌握如下6大技能体系!
  8. app = Flask(__name__)相关说明
  9. AI学习笔记--人机对话的四种形态
  10. tensorflow--模型的保存和提取